starting values for the optimization of HAPimLDL method

Share:

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

1
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.

Author(s)

S. Dejean, N. Oumouhou, D. Estivals, B. Mangin

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

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.