Test for sex biased dispersal

Description

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

Usage

1
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 jerome.goudet@unil.ch

References

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

Examples

1
2
3
4
5
6
7
8
9
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)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.