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.

Reply via email to