Pertama

ftn <- function(x,yij,i,j){
  x <- matrix(c(10,15,17,22,23,yij),ncol=3,byrow=TRUE)
  mu <- mean(x)
  alfai <- mean(x[i,])- mu
  betaj <- mean(x[,j])- mu
  yij <- mu + alfai + betaj
  return(yij)
}
emrancob <- function(ftn, x,yij,i, j, tol = 1e-9, max.iter = 100) {
  yold <- yij
  ynew <- ftn(x,yold,i,j)
  iter <- 1
  cat("At iteration 1 value of x is:", ynew, "\n")
  
  while ((abs(ynew-yold) > tol) && (iter < max.iter)) {
    yold <- ynew;
    ynew <- ftn(x,yold,i,j);
    iter <- iter + 1
    cat("At iteration", iter, "value of x is:", ynew, "\n")
  }
  
  if (abs(ynew-yold) > tol) {
    cat("Algorithm failed to converge\n")
    return(NULL)
  } else {
    cat("Algorithm converged\n")
    return(ynew)
  }
}
emrancob(ftn, x,0,2,3, tol = 1e-4, max.iter = 100)

## At iteration 1 value of x is: 9 
## At iteration 2 value of x is: 15 
## At iteration 3 value of x is: 19 
## At iteration 4 value of x is: 21.66667 
## At iteration 5 value of x is: 23.44444 
## At iteration 6 value of x is: 24.62963 
## At iteration 7 value of x is: 25.41975 
## At iteration 8 value of x is: 25.9465 
## At iteration 9 value of x is: 26.29767 
## At iteration 10 value of x is: 26.53178 
## At iteration 11 value of x is: 26.68785 
## At iteration 12 value of x is: 26.7919 
## At iteration 13 value of x is: 26.86127 
## At iteration 14 value of x is: 26.90751 
## At iteration 15 value of x is: 26.93834 
## At iteration 16 value of x is: 26.95889 
## At iteration 17 value of x is: 26.9726 
## At iteration 18 value of x is: 26.98173 
## At iteration 19 value of x is: 26.98782 
## At iteration 20 value of x is: 26.99188 
## At iteration 21 value of x is: 26.99459 
## At iteration 22 value of x is: 26.99639 
## At iteration 23 value of x is: 26.99759 
## At iteration 24 value of x is: 26.9984 
## At iteration 25 value of x is: 26.99893 
## At iteration 26 value of x is: 26.99929 
## At iteration 27 value of x is: 26.99952 
## At iteration 28 value of x is: 26.99968 
## At iteration 29 value of x is: 26.99979 
## At iteration 30 value of x is: 26.99986 
## Algorithm converged

## [1] 26.99986

Pertama Rumus Di Anstat

\[Y_{ij} = \frac {(rB_i+tT_j-G)}{(r-1)(t-1)}\]

t = 2
r = 3
T2 = 22 + 23
B3 = 17
G = 10 + 15 + 17 + 22 + 23
Y23 = (r*B3 + t*T2 - G)/((r-1)*(t-1))
Y23

## [1] 27

KEDUA

ftn <- function(x,yij,i,j){
  x <- matrix(c(12,yij,18,22,23,24),ncol=3,byrow=TRUE)
  mu <- mean(x)
  alfai <- mean(x[i,])- mu
  betaj <- mean(x[,j])- mu
  yij <- mu + alfai + betaj
  return(yij)
}
emrancob(ftn, x,14,1,2, tol = 1e-4, max.iter = 100)

## At iteration 1 value of x is: 14.33333 
## At iteration 2 value of x is: 14.55556 
## At iteration 3 value of x is: 14.7037 
## At iteration 4 value of x is: 14.80247 
## At iteration 5 value of x is: 14.86831 
## At iteration 6 value of x is: 14.91221 
## At iteration 7 value of x is: 14.94147 
## At iteration 8 value of x is: 14.96098 
## At iteration 9 value of x is: 14.97399 
## At iteration 10 value of x is: 14.98266 
## At iteration 11 value of x is: 14.98844 
## At iteration 12 value of x is: 14.99229 
## At iteration 13 value of x is: 14.99486 
## At iteration 14 value of x is: 14.99657 
## At iteration 15 value of x is: 14.99772 
## At iteration 16 value of x is: 14.99848 
## At iteration 17 value of x is: 14.99899 
## At iteration 18 value of x is: 14.99932 
## At iteration 19 value of x is: 14.99955 
## At iteration 20 value of x is: 14.9997 
## At iteration 21 value of x is: 14.9998 
## At iteration 22 value of x is: 14.99987 
## Algorithm converged

## [1] 14.99987

KETIGA

# pendugaan data tersembunyi (sebaran multinomial)
x3<- 18
x4<- 20
x5<- 34
x2d1<- 0
x2d2<- 0
theta0<- 2
err<- 10
iterasi<- 1
while(err>10^-4){
  theta1<- (x2d1+x5)/(x2d1+x3+x4+x5)
  x1d2<- 125*(0.5/(0.5+0.25*theta1))
  x2d2<- 125*((0.25*theta1)/(0.5+0.25*theta1))
  err<- abs(theta1-theta0)
  theta0<- theta1
  x2d1<-x2d2
  iterasi<- iterasi + 1
}
x1d2   # x1

## [1] 95.17232

x2d2   # x2

## [1] 29.82768

theta1

## [1] 0.6268142

iterasi

## [1] 7