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