Determines a multi-QTL model using a forward regression.
QTL_forward(
mppData = NULL,
trait = 1,
Q.eff,
threshold = 4,
window = 30,
n.cores = 1,
verbose = TRUE
)
Return:
Data.frame
of class QTLlist
with five columns :
1) QTL marker names; 2) chromosomes;
3) interger position indicators on the chromosome;
4) positions in centi-Morgan; and 5) -log10(p-values).
An object of class mppData
Numerical
or character
indicator to specify which
trait of the mppData
object should be used. Default = 1.
Character
vector of possible QTL effects the user want to
test. Elements of Q.eff can be "cr", "par", "anc" or "biall". For details
look at mpp_SIM
.
Numeric
value representing the -log10(p-value) threshold
above which a position can be considered as significant. Default = 4.
Numeric
distance (cM) on the left and the right of a
cofactor position where it is not included in the model. Default = 30.
Numeric
. Specify here the number of cores you like to
use. Default = 1.
Logical
value indicating if the steps of the
forward regression should be printed. Default = TRUE.
Vincent Garin
Forward regression to determine the a multi-QTL model. The function selects successively QTL positions with -log10(pval) above the threshold. Those positions are added as cofactors for following detection run. The procedure stop when no more position has a -log10(pval) above the threshold.
mpp_SIM
,
data(mppData)
QTL <- QTL_forward(mppData = mppData, Q.eff = "par")
Run the code above in your browser using DataLab