Displaying 1 result from an estimated 1 matches for "prop_".
Did you mean:
prop
2008 Oct 21
1
Simulation of data
...sing value where TRUE/FALSE needed
Thanks in advance!
M?rcio
k<-100
d<-0.00112
z<-rnorm(100, 0, 1)
prop<-rep(0,100)
for (i in 1:100){
if (z[i]>1.75){
prop[i]<-1}
else{prop[i]<-0}}
proportion<-mean(prop)
ecart<-sqrt((proportion*(1-proportion))/k)
while(ecart>=d){
prop_<-0
w<- rnorm(1, 0, 1)
z<-c(z,w)
{if (w>1.75){
prop_<-1}
else{
prop_<-0}
}
prop<-c(prop, prop_)
proportion<-proportion+(prop[k+1]-proportion/(k+1))
ecart<-sqrt((proportion*(1-proportion))/(k+1))
k<-k+1
}
--
View this message in context: http://www.nabble...