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