simul.na.bottom.herit: Simulation phase for P-SIMEX on heritability with missing...

Description Usage Arguments Value Author(s)

Description

The function generates increasing amount of errors to the pedigree structure and calculates heritability 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.

Usage

1
2
simul.na.bottom.herit(pedigree, pedigree0, lambda, B, 
data, model, nbottom, prior, nitt,thin,burnin)

Arguments

pedigree

A dataset containing the initial pedigree structure. It must have three columns: id, parent1, parent2.

pedigree0

A dataset containing the initial pedigree structure.

lambda

A vector of real numbers specifying the error proportion to be generated.

B

An integer specifying the number of simulations to be run for each error level.

data

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)

model

A object of class MCMCglmm used to calculate heritability.

nbottom

An integer specifying the number of the last generations where to generate errors.

prior

Prior distribution for MCMCglmm model

nitt

Number of iterations for MCMCglmm model

thin

Thinning interval for MCMCglmm model

burnin

Burn in period for MCMCglmm model

Value

A list:

h

A dataset containing all the simulated values for heritability

se_h

A dataset containing all the simulated standard errors for heritability

VA

A dataset containing all the simulated values for additive variance

VE

A dataset containing all the simulated values for environmental variance

Author(s)

Erica Ponzi


PSIMEX documentation built on May 2, 2019, 3:34 p.m.