[R] R: simulation
Clark Allan
Allan at STATS.uct.ac.za
Fri Mar 4 13:21:27 CET 2005
Hi
thanx for the reply.
i used your code and pasted it into R and ran it a few times. the output
is below. what i want is to get the same output every time the program
is run. is this possible?
another question?
x<-rnorm(100)
y<-rnorm(100)
is x and y independent?
> rnorm(1)
[1] 0.4251004
> rnorm(1)
[1] -0.2386471
> rnorm(1)
[1] 1.058483
> rnorm(1)
[1] 0.8864227
> rnorm(1)
[1] -0.619243
> rnorm(1)
[1] 2.206102
> rnorm(1)
[1] -0.2550270
> rnorm(1)
[1] -1.424495
> rnorm(1)
[1] -0.1443996
> rnorm(1)
[1] 0.2075383
> rnorm(1)
[1] 2.307978
> rnorm(1)
[1] 0.1058024
Dimitris Rizopoulos wrote:
>
> look at ?set.seed, e.g.,
>
> rnorm. <- function(seed, ...){
> set.seed(seed)
> rnorm(...)
> }
> #########
> rnorm.(100, 1, 2)
>
> I hope it helps.
>
> Best,
> Dimitris
>
> ----
> Dimitris Rizopoulos
> Ph.D. Student
> Biostatistical Centre
> School of Public Health
> Catholic University of Leuven
>
> Address: Kapucijnenvoer 35, Leuven, Belgium
> Tel: +32/16/336899
> Fax: +32/16/337015
> Web: http://www.med.kuleuven.ac.be/biostat/
> http://www.student.kuleuven.ac.be/~m0390867/dimitris.htm
>
> ----- Original Message -----
> From: "Clark Allan" <Allan at stats.uct.ac.za>
> To: <r-help at stat.math.ethz.ch>
> Sent: Friday, March 04, 2005 12:07 PM
> Subject: [R] R: simulation
>
> > hi all
> >
> > a simple question
> >
> > i want to run simulations in r. i however want the experiments to be
> > repeated at a later time with exactly the same numbers by other
> > users.
> > can i set the random number seed for rnorm in some way?
> >
> > e.g. is there some arguement that goes with rnorm?
> >
> > please supply an example
> >
> > regards
> > Allan
>
> --------------------------------------------------------------------------------
>
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide!
> > http://www.R-project.org/posting-guide.html
More information about the R-help
mailing list