This function calls the proper simulation function for the parameter of interest with missing paternities.
simul.na(pedigree, pedigree0, lambda, lambda0, B, data,
model, parameter, way, ntop, nbottom, prior, nitt,thin,burnin)A dataset containing the initial pedigree structure. It must have three columns: id, parent1, parent2.
A dataset containing the initial pedigree structure. It must have five columns: id, parent1, parent2, sex, generation.
A vector of real numbers specifying the error proportion to be generated.
A real number specifying the initial error rate.
An integer specifying the number of simulations to be run for each error level.
A dataset containing the phenotypic measurements on the population and the covariates which are included in the model. The trait shoul be named differently than 'trait' (see MCMCglmm)
An object specifying the model to calculate the parameter of interest. It can be a lm or glm for inbreeding depression and a MCMCglmm for heritability.
A string specifying the parameter of interest. It must be 'inbreeding' or 'heritability'.
A string specifying how errors are generated. It must be 'uniform' or 'similar'.
An integer specifying the number of the first generations where to generate errors, when the parameter 'way' is 'top'.
An integer specifying the number of the last generations where to generate errors,when the parameter 'way' is 'bottom'.
Prior distribution for MCMCglmm model
Number of iterations for MCMCglmm model
Thinning interval for MCMCglmm model
Burn in period for MCMCglmm model
For inbreeding a list:
A dataset containing all the simulated values for inbreeding depression
A dataset containing all the simulated standard errors for inbreeding depression
A dataset containing all the simulated p-values for inbreeding depression
A dataset containing all the simulated values for inbreeding coefficient mean
A dataset containing all the simulated values for inbreeding coefficient median
A dataset containing all the simulated values for inbreeding coefficient variance
For heritability a list:
A dataset containing all the simulated values for heritability
A dataset containing all the simulated standard errors for heritability
A dataset containing all the simulated values for additive variance
A dataset containing all the simulated values for environmental variance