[R-sig-ME] Loglikelihood

[bbonit at tin.it]

 m1 <- matrix(c(597.1903,-0.97,-0.97,60.05023),nrow=2,byrow=TRUE)
D1 <-as.matrix(bdiag(m1, m1, m1))
V1 <-Z%*% D1 %*% t(Z) + R 
D1

         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]
[1,] 597.1903 -0.97000   0.0000  0.00000   0.0000  0.00000
[2,]  -0.9700 60.05023   0.0000  0.00000   0.0000  0.00000
[3,]   0.0000  0.00000 597.1903 -0.97000   0.0000  0.00000
[4,]   0.0000  0.00000  -0.9700 60.05023   0.0000  0.00000
[5,]   0.0000  0.00000   0.0000  0.00000 597.1903 -0.97000
[6,]   0.0000  0.00000   0.0000  0.00000  -0.9700 60.05023

insted of 
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]
[1,] 597.1903 -0.97000  -0.9700 -0.97000  -0.9700 -0.97000
[2,]  -0.9700 60.05023  -0.9700 -0.97000  -0.9700 -0.97000
[3,]  -0.9700 -0.97000 597.1903 -0.97000  -0.9700 -0.97000
[4,]  -0.9700 -0.97000  -0.9700 60.05023  -0.9700 -0.97000
[5,]  -0.9700 -0.97000  -0.9700 -0.97000 597.1903 -0.97000
[6,]  -0.9700 -0.97000  -0.9700 -0.97000  -0.9700 60.05023


but this dont affect the decomposition into ll (loglike) and dont change my loglike 
	[[alternative HTML version deleted]]