[R] A simple simulation question
Ivan Calandra
ivan.calandra at u-bourgogne.fr
Thu Jul 12 15:33:11 CEST 2012
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
More information about the R-help
mailing list