try this:
tol <- 0.01
mat <- matrix(as.numeric(NA), 1000, 5)
k <- 1
while(any(is.na(mat))){
x <- rnorm(1000, sd = 0.02)
if (abs(mean(x)) < tol) {
mat[, k] <- x
k <- k + 1
}
}
abs(colMeans(mat))
par(mfrow = c(2, 3))
apply(mat, 2, hist)
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/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
http://www.student.kuleuven.be/~m0390867/dimitris.htm
- Original Message -
From: "Ing. Michal Kneifl, Ph.D." <[EMAIL PROTECTED]>
To: "Rhelp"
Sent: Monday, July 23, 2007 4:40 PM
Subject: [R] Help on "looping problem" needed!
>I am wondering if someone could help me out with following problem:
> I have written a for loop which generates a random normal
> distribution
> let us say 1000 times.
> When the restriction is met (mean<0.01), the loop stops, prints
> the mean value and plots a histogram.
>
> for(i in 1:1000) {
> a<-rnorm(1000,0,.2)
> b<-abs(mean(a))
> if(b>.01) next else {print(b);hist(a);break}}
>
> How to reshape the loop when I want to find at least 5 distibutions
> that meet my restriction and save them (assign) under
> names R1R5.
> Could you help me please?
>
> Michael
>
> __
> R-help@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
> and provide commented, minimal, self-contained, reproducible code.
>
Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
__
R-help@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
and provide commented, minimal, self-contained, reproducible code.
