hierfstat (version 0.5-11)

sexbias.test: Test for sex biased dispersal

Description

Test whether one sex disperses more than the other using the method described in Goudet etal. (2002)

Usage

sexbias.test(dat,sex,nperm=NULL,test="mAIc",alternative="two.sided")

Arguments

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"

Value

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

References

Goudet J, Perrin N, Waser P (2002) Tests for sex-biased dispersal using bi-parentally inherited genetic markers 11, 1103:1114

Examples

# NOT RUN {
  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")
  
# }