mAdd1 {qtlmt} | R Documentation |
Add or drop all possible terms
Description
Add or drop all possible terms to or from a multivariate multiple regression model.
Usage
mAdd1(object, scope, test = c("none","Chisq","F"), k = 0, ...)
mDrop1(object, scope, test = c("none","Chisq","F"), k = 0, ...)
Arguments
object |
initial model. |
scope |
a formula or a vector of variable names, giving lower/upper bound of the model. |
test |
if not "none", the test statistic and p-value will be included in the output. |
k |
penalty on a parameter in AIC. |
... |
additional arguments to |
Value
An object summarizing the differences in fit between the models.
See Also
Examples
data(etrait)
mdf<- data.frame(traits,markers)
## Not run:
mlm<- lm(cbind(T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,T15,T16) ~
m1 + m2 + m3 + m4 + m5, data=mdf)
up<- formula(paste("~", paste("m",1:15,collapse=" + ",sep="")))
oa<- mAdd1(mlm, scope=up, test="F", k=5, data=mdf)
lw<- formula(paste("~ ", paste("m",1:3,collapse=" + ",sep="")))
od<- mDrop1(mlm, scope=lw, test="F", k=5, data=mdf)
## End(Not run)
[Package qtlmt version 0.1-6 Index]