Dimitris, Thanks, I shall give this a try as an alternative.
Graham 2009/11/15 Dimitris Rizopoulos <d.rizopou...@erasmusmc.nl>: > try the following: > > LL <- c(12.5,17,12,11.5,9.5,15.5,16,14) > n <- length(LL) > N <- 1000 > threshold <- 10 > > smpls <- sample(LL, N*n, replace = TRUE) > dim(smpls) <- c(n, N) > cnt <- sum(colMeans(smpls) > threshold) > cnt > > > I hope it helps. > > Best, > Dimitris > > > Graham Smith wrote: >> >> I am trying to modify some code from Good 2005. >> >> I am trying to resample the mean of 8 values and then count how many >> times the resampled mean is greater than 10. But my count of means >> above 10 is coming out as zero, which I know isn't correct. >> >> I would appreciate it if someone could look at the code below and tell >> me what I am doing wrong. >> >> Many thanks, >> >> Graham >> >>> LL<- c(12.5,17,12,11.5,9.5,15.5,16,14) >>> N<-1000 >>> n<-length(LL) >>> threshold<-10 >>> cnt<-0 >>> for(i in 1:N){ >> >> + LLb <- sample (LL, n, replace=TRUE) >> + if (mean(LLb)<=threshold) cnt<-cnt+1 >> + } >>> >>> cnt >> >> [1] 0 >> >> ______________________________________________ >> 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. >> > > -- > Dimitris Rizopoulos > Assistant Professor > Department of Biostatistics > Erasmus University Medical Center > > Address: PO Box 2040, 3000 CA Rotterdam, the Netherlands > Tel: +31/(0)10/7043478 > Fax: +31/(0)10/7043014 > ______________________________________________ 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.