Dimitris,
Thanks, I shall give this a try as an alternative.
Graham
2009/11/15 Dimitris Rizopoulos :
> 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(
David,
Thanks, its me getting mixed up I actually meant less than or equal to
10. That apart, I guess the code is OK, I just expected, especially
as I increased N that I might have got some means less than 10, but
having gone back to it , I see I need a million iterations before
getting two means
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 cod
On Nov 15, 2009, at 12:12 PM, 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 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 c
5 matches
Mail list logo
