[R] Size of subsample in ecodist mantel()
Sarah Goslee
sarah.goslee at gmail.com
Wed Jul 4 16:51:42 CEST 2012
pboot is the proportion of the sample to select, and so if it's
greater than or equal to 1 you're using the entire sample rather than
a subsample, so of course the limits are equal.
If you look at the usage line in the help for mantel, where defaults
are given, the default for pboot = 0.9.
The details are discussed in the associated JSS paper that is given by
citation("ecodist")
Sarah
On Thu, Jun 28, 2012 at 9:28 PM, nevil amos <nevil.amos at gmail.com> wrote:
> 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.goslee at 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
>>>
--
Sarah Goslee
http://www.functionaldiversity.org
More information about the R-help
mailing list