I am performing a MonteCarlo Test for a set of more than 1 million values. Is time consuming to make it run for every single value so I decided to go for a sample size of 10000, 20000, 30000 etc. But, I am having problems re arreging the SaSi.
for(j in 1:SaSi)
run<-100 # run 1000 times (the more often you run, the smoother the output curve becomes, but the longer it takes)
for (i in 1:run) #loop
a<-sample(as.numeric(flux), j, replace=TRUE, prob=NULL) #draw random sample from flux, sample size j
b<-mean(a) #mean value of sample
y[i]<-b # write mean value into empty vector
cv<-sd(y)/mean(y)*100 #calculate Coefficiant of variance for respective sample size
plot(SS, CV, type="l", ylab="CV (%)", xlab= "sample size")
I did it manually like SaSi<-c(100000, 200000, ...), but it does not work. I got this error> Warning message:
In 1:SaSi : numerical expression has 2 elements: only the first used
Any suggestion is very welcome.