Default for boot in mantel{ecodist} was Re: Size of subsample in ecodist mantel()

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

Default for boot in mantel{ecodist} was Re: Size of subsample in ecodist mantel()

Nevil Amos
How do I use pboot to set the level of resampling (size of each resample) in mantel()?

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

many thanks

Nevil Amos


        [[alternative HTML version deleted]]

______________________________________________
[hidden email] 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.