Learn R Programming

HAPim (version 1.3)

depart.LDL: starting values for the optimization of HAPimLDL method

Description

The function calculates the starting value of the error variance and the starting value of the QTL effect for the optimization of HAPimLDL method. It can be viewed as an internal function. The user does not have to call it by himself.

Usage

depart.LDL(moyenne.pere, perf, CD, PLA, desc.pere)

Arguments

moyenne.pere
results provided by moyenne.pere() function, mean of half-sib family performances.
perf
numeric vector of length=number of individuals which contains the performances of individuals.
CD
numeric vector of length=number of individuals which contains the CD of individuals. var(perf$_i$)=s/CD$^2_i$
PLA
numeric vector (number of individuals) which contains transmission probabilities at a single test position.
desc.pere
results provided by descendant.pere() function, numeric matrix (number of sires x 2) which gives for each sire, the first and last indexes of its progeny.

Value

The returned value is a numeric vector of length=2 which contains estimates of the error variance and the Q allele effect.

References

publication to be submitted: C. Cierco-Ayrolles, S. Dejean, A. Legarra, H. Gilbert, T. Druet, F. Ytournel, D. Estivals, N. Oumouhou and B. Mangin. Combining linkage analysis and linkage disequilibrium for QTL fine mapping in animal pedigrees.

See Also

moyenne.pere, descendant.pere