By the way, my "for" function is below, I can't find the mistake
rand.max.t<- function(n){ f<-rep(NA,n) for (i in 1:n) { reassign[i]<-matrix(c(sample(id),data1),203,3) new.data<-reassign[,1] random.cas=reassign[new.data==0,2:3] random.con=reassign[new.data==1,2:3] f<- list(x=random.cas[,1],y=random.cas[,2]) ppregion(xl=min(random.cas[,1])-0.0001,xu=max(random.cas[,1])+0.0001,yl=min(random.cas[,2])-0.0001,yu=max(random.cas[,2])+0.0001) g <- Kfn(f,0.01) random.g <- ((g$y)^2)/pi } random.g } rand.max.t(10) -- View this message in context: http://r.789695.n4.nabble.com/Repeating-tp4614371p4614482.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.