Hi Manning, There are two obvious mistakes: - close the brackets for replicate() and - do the things in the correct order
replicate(100, { x=rep(1:10,10) ## first define x y=rnorm(100,x,5) ## and then y because it depends on x plot(y~x) ## then plot, because it depends on x and y abline(lm(y~x)) ## then add the line because it depends on plot() }) ## don't forget the last bracket Now it seems to work! HTH, Ivan -- Ivan CALANDRA Université de Bourgogne UMR CNRS/uB 6282 Biogéosciences 6 Boulevard Gabriel 21000 Dijon, FRANCE +33(0)3.80.39.63.06 ivan.calan...@u-bourgogne.fr http://biogeosciences.u-bourgogne.fr/calandra Le 12/07/12 14:55, Aldous Huxley a écrit : Hi! I would like to post the following question: I was trying to figure out how to do the simulation shown in Fig 10.6 of John Verzani's book 'Using R for Intro Statistics'. It is on page 290, with a description on the previous page. It seems like a simple thing... Just needing to duplicate a procedure. (Perhaps I need to do it with a loop?) This is what I was trying: replicate(100, { abline(lm(y~x)) plot(y~x) y=rnorm(100,x,5) x=rep(1:10,10) } Obviously not correct. Can anyone help me with a simple way to do this? With thanks, Manning ______________________________________________ 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.