[R] A simple simulation question

[Ivan Calandra]

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.calandra at 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