getP: An internal function for getting empirical p values
In multiDimBio: Multivariate Analysis and Visualization for Biological Data

Description Usage Arguments Value Examples

Simulates p values.

1	getP(ndads, mm, vv, tau2, nperms, nsims, nbins)

`ndads`	a (non-empty) numeric value indicating the number of dads.
`mm`	a (non-empty) numeric value indicating the mean number of offspring per dad per bin (normal dist). mm must be less than vv.
`vv`	a (non-empty) numeric value indicating the variance in offspring per dad per bin (normal dist). vv must be greater than mm.
`tau2`	a (non-empty) numeric value indicating the dad effect (narrow-sense heritability ~ tau2/(tau2+(pi/sqrt(3))^2)).
`nperms`	a (non-empty) numeric value indicating the number of bootstrap permutations to use for caluclating a p value.
`nsims`	a (non-empty) numeric value indicating the number of simulations to run per parameter combination.
`nbins`	a (non-empty) numeric value indicating the number of bins, data are pooled before analysis.

Returns a vector of simulated p values. The list contains:

	ndads <- c(9,18)
	mm <- 4.629634
	vv <- 6.31339
	tau2 <- c(0,0.5)
	nperms <- 2
	nsims <- 2
	nbins <- 3
	getP(ndads, mm, vv, tau2, nperms, nsims, nbins)