Y<-matrix(c(400, 380, 410, 350, 420, 360, 390, 390, 430, 370), nrow=10) print (Y) X<-matrix (0,nrow=10, ncol=4) X [1:3,1]<-1 X [4:6,2]<-1 X [7:10,3]<-1 X [1,4]<-1 X [3,4]<-1 X [5,4]<-1 X [7,4]<-1 X [9,4]<-1 print(X) XT <- t(X) print (XT) Z1 <- matrix (0,nrow = 10, ncol = 16) Z1 [1,7]<-1 Z1 [2,8]<-1 Z1 [3,9]<-1 Z1 [4,10]<-1 Z1 [5,11]<-1 Z1 [6,12]<-1 Z1 [7,13]<-1 Z1 [8,14]<-1 Z1 [9,15]<-1 Z1 [10,16]<-1 print (Z1) Z1T <- t(Z1) print (Z1T) Z2 <- matrix (0,nrow = 10, ncol = 6) Z2 [1,1]<-1 Z2 [2,2]<-1 Z2 [3,3]<-1 Z2 [4,1]<-1 Z2 [5,2]<-1 Z2 [6,4]<-1 Z2 [7,1]<-1 Z2 [8,4]<-1 Z2 [9,5]<-1 Z2 [10,6]<-1 print (Z2) Z2T <- t(Z2) print (Z2T) I<-matrix (0, nrow=6, ncol=6) I [1,1]<-1 I [2,2]<-1 I [3,3]<-1 I [4,4]<-1 I [5,5]<-1 I [6,6]<-1 print (I) AI<-matrix(0, nrow=16, ncol=16) AI [1,1]<-2.5 AI [1,2]<-.5 AI [1,3]<-1 AI [1,7:8]<- -1 AI [1,11]<- -1 AI [2,1]<- .5 AI [2,2]<- 2.5 AI [2,4]<- 1 AI [2,7]<- -1 AI [2,10]<- -1 AI [2,13]<- -1 AI [3,1] <- 1 AI [3,3] <- 2 AI [3,8] <- -1 AI [3,11] <- -1 AI [4,2] <- 1 AI [4,4] <- 3 AI [4,5] <- .5 AI [4,8] <- .5 AI [4,9:10] <- -1 AI [4,12:13] <- -1 AI [5,4]<- .5 AI [5,5]<- 1.5 AI [5,9]<- -1 AI [6,6]<- 2.5 AI [6,8]<- .5 AI [6,10]<- .5 AI [6,12]<- .5 AI [6,14:16]<- -1 AI [7,1:2]<- -1 AI [7,7]<- 2 AI [8,1] <- -1 AI [8,3] <- -1 AI [8,4] <- 0.5 AI [8,6] <- .5 AI [8,8] <- 3 AI [8,12] <- -1 AI [8,14] <- -1 AI [9,4:5] <- -1 AI [9,9] <- 2 AI [10,2] <- -1 AI [10,4] <- -1 AI [10,6] <- 0.5 AI [10,10] <- 2.5 AI [10,15] <- -1 AI [11,1] <- -1 AI [11,3] <- -1 AI [11,11] <- 2 AI [12,4] <- -1 AI [12,6] <- .5 AI [12,8] <- -1 AI [12,12] <- 2.5 AI [12,16] <- -1 AI [13,2] <- -1 AI [13,2] <- -1 AI [13,4] <- -1 AI [13,13] <- 2 AI [14:16,6] <- -1 AI [14,8] <- -1 AI [14,14] <- 2 AI [15,10] <- -1 AI [15,15] <- 2 AI [16,12] <- -1 AI [16,16] <- 2 print (AI) S2e <- 6500 P<-matrix(c(2000, 0, 0, 500), nrow=2) print (P) PI<- solve (P) print (PI) k<-PI*S2e print (k) k11<- k [1,1] print (k11) k22<- k [2,2] print (k22) AIk<-AI*k11 print (AIk) Ik<-I*k22 print (Ik) LHS<-rbind(cbind(XT%*%X, XT%*%Z1, XT%*%Z2), cbind(Z1T%*%X, Z1T%*%Z1+AIk, Z1T%*%Z2), cbind(Z2T%*%X, Z2T%*%Z1, Z2T%*%Z2+Ik)) print (LHS) RHS <- rbind(XT%*%Y, Z1T%*%Y, Z2T%*%Y) print (RHS) C<-solve(LHS) print (C) BU <- C%*%RHS print(BU) year<- BU [1:3,1] sex<- BU [4,1] a<- BU [5:20,1] p<- BU [21:26,1] print (year) print (sex) print (a) print (p)