I want to iterate the program until it converge? can anybody help me? I just write it manual, but i need to generalize it so that simulation can stop if the difference between two iteration is less than 0.0001. The following are the program :
rm(list=ls(all=T)) x=rnorm(50,2,0.7) y=rnorm(50,3,0.3) beta=0.80 n=length(x) alpha=mean(y)-beta*mean(x) ###initail value indel=sum((x-mean(x))^2)/(n) inep=sum((y-mean(y))^2)/(n) indel inep inXhat=((x/indel)+(beta/inep)*(y-alpha))/((1/indel)+(beta^2/inep)) del1=sum((x-inXhat)^2)/(n) ep1=sum((y-alpha-beta*inXhat)^2)/(n) del1 ep1 ## 1st iteration Xhat1=((x/del1)+(beta/ep1)*(y-alpha))/((1/del1)+(beta^2/ep1)) del2=sum((x-Xhat1)^2)/(n) ep2=sum((y-alpha-beta*Xhat1)^2)/(n) del2 ep2 abs(del2-del1) abs(ep2-ep1) ## 2nd Iteration Xhat2=((x/del2)+(beta/ep2)*(y-alpha))/((1/del2)+(beta^2/ep2)) del3=sum((x-Xhat2)^2)/(n) ep3=sum((y-alpha-beta*Xhat2)^2)/(n) del3 ep3 abs(del3-del2) abs(ep3-ep2) -- View this message in context: http://r.789695.n4.nabble.com/Printing-the-loop-number-for-each-iteration-tp3683495p4646567.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.