sexbias.test: Test for sex biased dispersal

In hierfstat: Estimation and Tests of Hierarchical F-Statistics

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

Author(s)

Jerome Goudet [e-mail]

References

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

Examples

  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)

hierfstat documentation built on May 6, 2022, 1:05 a.m.