
Hi r-users;
 
I have this values:
eign_val <- c(137.810447,3.538721,2.995161,1.685670)
alp    <- 1.6549 ;  lamda <- eign_val
lamda_m <- min(lamda)
 
First I calculated manually:
delta0 <- 1
delta1 <- alp*delta0*(4-lamda_m*(1/lamda[1]+1/lamda[2]+1/lamda[3]+1/lamda[4]))  
delta1
delta2 <- (alp/2)*(delta1*(delta1/alp) + delta0*((1-lamda_m/lamda[1])^2+ (1-lamda_m/lamda[2])^2+(1-lamda_m/lamda[3])^2+(1-lamda_m/lamda[4])^2))
delta2
delta3 <- (alp/3)*(delta2*(delta1/alp) + delta1*((1-lamda_m/lamda[1])^2+ (1-lamda_m/lamda[2])^2+(1-lamda_m/lamda[3])^2+(1-lamda_m/lamda[4])^2)
 + delta0*((1-lamda_m/lamda[1])^3+ (1-lamda_m/lamda[2])^3+(1-lamda_m/lamda[3])^3+(1-lamda_m/lamda[4])^3))
delta3
c(delta1,delta2,delta3)
> c(delta1,delta2,delta3)
[1]  3.224772  6.391966 10.091279
 
Then I wrote this the following code:
 
term <- function(i, lamda, lamda_m)
     {sum(sapply(lamda, function(lamda, lamda_m,i) ((1-lamda_m/lamda)^i), lamda_m, i))}
sm <- sapply(1:N, term,  lamda, lamda_m);sm
#now calculate the deltas
k1 <- 3+1
delta <- rep(1, k1);delta
delta_calc <- function(k1, delta, sm, alp)
{ k <- k1-1
  alp/k*sum(sapply(1:k, function(i, delta, sm, k1) (sm[i]*delta[k1-i]), delta, sm, k1))
}
delta[2:k1] <- sapply(2:k1, delta_calc, delta, sm, alp);delta[2:k1]
 
> delta[2:k1] <- sapply(2:k1, delta_calc, delta, sm, alp);delta[2:k1]
[1] 3.224772 2.804775 2.526796
 
I could not trace why I the answers are different.
I hope somebody can help me with this.
 
Thank you so much for your help.
 


      
	[[alternative HTML version deleted]]

