Dear all,
   In the following code, I was trying to compute each row of the "param"
iteratively based
on the first row.
   This likely is not the best way. Can anyone suggest a simpler way to
improve the code.
   Thanks a lot!
         Hannah


param <- matrix(0, 11, 5)

colnames(param) <- c("p", "q", "r", "2s", "t")

param[1,] <- c(0.5, 0.5, 0.4, 0.5, 0.1)

for (i in 2:11){

param[i,1] <- param[(i-1),3]+param[(i-1),4]/2

param[i,2] <- param[(i-1),4]/2+param[(i-1),5]

param[i,3] <- param[(i-1),1]*(param[(i-1),3]+param[(i-1),4]/2)

param[i,4] <- param[(i-1),1]*(param[(i-1),4]/2+param[(i-1),5])+param[(i-1),2
]*(param[(i-1),3]+param[(i-1),4]/2)

param[i,5] <- param[(i-1),2]*(param[(i-1),4]/2+param[(i-1),5])

}


round(param,3)

	[[alternative HTML version deleted]]

