Thanks Sarah, It is not clear to me exactly how I set this value. if I enter a value for pboot then the ulim==llim
> X<-dist(1:100)> Y<-dist(1:100+50*rnorm(100))> length(X)[1] 4950> > print(mantel(X~Y,nperm=1000,nboot=1000,pboot=10)) mantelr pval1 > pval2 pval3 llim.2.5% ulim.97.5% 0.1396906 0.0010000 1.0000000 0.0010000 0.1396906 0.1396906 if I do not set a value for pboot then as expected ulim>llim > print(mantel(X~Y,nperm=1000,nboot=1000)) mantelr pval1 pval2 > pval3 llim.2.5% ulim.97.5% 0.1396906 0.0010000 1.0000000 0.0010000 0.1005011 0.1805416 What is the default for pboot? this does not appear to be dealt with in help for mantel On Thu, Jun 28, 2012 at 8:35 PM, Sarah Goslee <sarah.gos...@gmail.com>wrote: > You can set it using the pboot argument. > > Sarah > > > On Thursday, June 28, 2012, nevil amos wrote: > >> What is the size of the boostrapped subsample in ecodist mantel() >> >> thanks >> >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> 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. >> > > > -- > Sarah Goslee > http://www.stringpage.com > http://www.sarahgoslee.com > http://www.functionaldiversity.org > [[alternative HTML version deleted]] ______________________________________________ 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.
