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
[[alternative HTML version deleted]]
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
I think that you have the statements out of order and I know that you are
lacking the last ) to match the replicate(
Try this"
replicate(100, {
x=rep(1:10,10)
y=rnorm(100,x,5)
plot(y~x)
abline(lm(y~x))
})
John Kane
Kingston ON Canada
> -----Original Message-----
> From: aldoushuxley11 at gmail.com
> Sent: Thu, 12 Jul 2012 13:55:04 +0100
> To: r-help at r-project.org
> Subject: [R] A simple simulation question
>
> 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
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at 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.
____________________________________________________________
FREE 3D MARINE AQUARIUM SCREENSAVER - Watch dolphins, sharks & orcas on your
desktop!