Hi > > My algorithm as follows; > y <- c(1,1,1,0,0,1,0,1,0,0) > x <- c(1,0,0,1,1,0,0,1,1,0) > > n <- length(x) > > t <- matrix(cbind(y,x), ncol=2)
Do not use t, it is a function for transposing matrix and after you redefine it you can get nasty surprise in future. tt <- cbind(y,x) is enough > > z = x+y > > for(j in 1:length(x)) { > out <- vector("list", ) > > for(i in 1:10) { > > t.s <- t[sample(n,n,replace=T),] t.s <- tt[sample(n,n,replace=T),] > > y.s <- t.s[,1] > x.s <- t.s[,2] > > z.s <- y.s+x.s > > out[[i]] <- list(ff <- (z.s), finding=any (y.s==y[j])) Here you compare vector y.s with one element of y as y.s is set of (0,1) values y is either 0 or 1, any tests if there is any match so only in rare case where all values in y.s are 0 and y[something] is 1 you get FALSE > kk <- sapply(out, function(x) {x$finding}) finding is (almost) always TRUE therefore kk is TRUE > ff <- out[! kk] > } > > I tried to find the total of the two vectors as statistic by using > bootstrap. Finally, I want to get the values which do not contain the y's > each elemet. In the algorithm ti is referred to "ff". But i get always the > same result ; I do not understand your intention so it is difficult to help. What is total of two vectors? sum? What does it mean "to get values which do not contain y's each element"? Maybe you shall rethink your code and first try to evaluate each line separately to see what it does and if the result is same as you intended. Regards Petr > > ff > list() > > kk > [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE > Because, my "y" vector contains only 2 elements, and probably all of the > bootstrap resamples include "1", or all of resamples include "0". So I can > not find the true matches. Can anyone help me about how to be? > Thanks. > > -- > View this message in context: http://r.789695.n4.nabble.com/indexing- > tp4428210p4428210.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. ______________________________________________ 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.