Test for sex biased dispersal
Test whether one sex disperses more than the other using the method described in Goudet etal. (2002)
sexbias.test(dat,sex,nperm=NULL,test="mAIc",alternative="two.sided")
dat |
a data frame with n.locs+1 columns and n.inds rows |
sex |
a vector containing the individual's sex |
nperm |
the number of permutation to carry out |
test |
one of "mAIc" (default), "vAIc","FIS" or "FST" |
alternative |
one of "two.sided" (default),"less" or "greater" |
call the function call
res the observation for each sex
statistic the observed statistic for the chosen test
p.value the p-value of the hypothesis
Jerome Goudet jerome.goudet@unil.ch
Goudet J, Perrin N, Waser P (2002) Tests for sex-biased dispersal using bi-parentally inherited genetic markers 11, 1103:1114
data(crocrussula)
sexbias.test(crocrussula$genot,crocrussula$sex)
dat<-qn2.read.fstat(system.file("extdata","qn2_sex.dat",package="hierfstat"))
sexbias.test(dat[[1]],sex=dat[[2]])
## Not run:
sexbias.test(crocrussula$genot,crocrussula$sex,nperm=1000)
sexbias.test(dat[[1]],sex=dat[[2]],nperm=100,test="FST",alternative="greater")
## End(Not run)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.