The function generates increasing amount of errors to the pedigree structure and calculates inbreeding depression for each simulation. For each level of the error proportion, a number of simulation equal to B are repeated. Errors are generated within the last nbottom generations of pedigree by eliminating fathers.
simul.na.bottom(pedigree, pedigree0, lambda, B, data, model, nbottom)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.
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)
A object of class model used to calculate inbreeding depression.
An integer specifying the number of the last generations where to generate errors.
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