GAUSS에서 R 코드로 변환

Question

GAUSS로 작성된 MLE 견적서를 사용하여 R 코드를 다시 작성하려고합니다. 사용하지도 않고 GAUSS 자체도 사용하지 않았습니다. 그것). 코드에는 약간 혼란스러운 선이 있습니다.

invsig: scalar or m-by-m matrix with inverse of sigma 

을 나는 조각으로 조각을 작동하도록 코드를 받고 일하고 있지만, 내 첫 번째 : 검증 된 가우스 코드에서

는 상태 (의견의 일부) "입력"아래 줄이 질문은 비교적 간단합니다. 여기

나를 혼란 가우스 조각입니다 :

... 
local m, k, tobs, invsig 
m = rows(y);k = rows(x); tobs = rows(dat) 
invsig= eye(m)*invsig 
... 

나는 이것이 "입력"invsig를 곱한 행렬은 이해하지만 시뮬레이션 예제의 코드에 부착 된 로그 파일 (작동 할 수), invsig의 스칼라 값을 사용하여 프로그램을 시작할 수 있습니다. IE : 초기 값을 invsig = 1로 설정

R에서는 이것이 작동하지 않습니다.

y.mat <- rep(rexp(3)) 
x.mat <- matrix(rexp(36), 12, 3) 
myfct <- function(x,invsig){ 
    m <- nrow(x) 
    invsig <- diag(m)%*%invsig 
    return(invsig) 
} 
t1 <- myfct(x.mat, 1)  ##Non-conformable error 
t2 <- myfct(x.mat, y.mat) ##Works 

나는 내가 가우스와 R 사이의 변환에 뭔가를 놓친 거지되는 질문 R.에서 점점 오전 비 정합 오류를 이해 : 여기에 간단한 "테스트"코드는이를 얻으려고하는 것입니다? GAUSS는 도움말 매뉴얼을 온라인으로 읽는 동안 개별 기호 (*/+ -)를 사용하여 행렬 연산을 수행하고 요소 방식으로 일을 수행 할 때 "." 각 작업 전에. 그래서 저에게 GAUSS 코드는 행렬 곱셈 (R에서 %*%)을 수행한다고 말하면서 단순 함수가 시도하는 것입니다.

모든 의견이나 제안을 보내 주시면 대단히 감사하겠습니다.

Answer 1

문제는 3x3 매트릭스에 1x1 매트릭스를 매트릭스 곱할 수 없다는 것입니다. 나는 첫 번째 경우에 diag(m)하여 곱셈에 남아했지만, invsig 이미 매트릭스의 경우 물론이 불필요 편집

myfct <- function(x,invsig){ 
    if (is.matrix(invsig)) return(invsig) 
    m <- nrow(x) 
    return(diag(invsig, nrow = m)) 
} 

같은 것을 권 해드립니다.