Learn R Programming

multic (version 0.4.3.1)

addGE: Assess combinations of univariate multic objects

Description

Determine whether there is any evidence that running a multivariate multic model will significantly increase the evidence of a genetic effect.

Usage

addGE(multic.objs, combine=2, plotit=FALSE, ibd.dist, statistic=c("lrt", "wald"), legend=TRUE, ylim=NULL, ...)

Arguments

multic.objs
A list of 1-trait multic objects.
combine
Indicate how many traits should be examined together. The program will then look at all N traits choose 'combine'.
plotit
Logical, default=FALSE. If TRUE, a LOD plot is generated with a separate line for each combination of traits.
ibd.dist
The default is to use the distances from the first multic object. This options allows the user to provide a different set of distances.
statistic
Character, default="lrt". This determines whether the Wald statistic (MG/SE) or the LRT is used when combining the traits.
legend
Logical, default=TRUE. If TRUE and if plotit=TRUE then a legend is automatically provided.
ylim
The extremes of the range of the y-axis to pass to the plot function.
...
Allows for graphical parameters to be passed to the plot function (only applicable when plotit=TRUE).

Value

A data frame is returned if the function is assigned to new object. Included are the various combinations (listed in order 1-N), the Chi-square statistic, the p-value, the distance, and the LOD score.

Side Effects

If plotit=T, a plot is generated on the current graphical device.

References

M. de Andrade, C. Olswold, J.P. Slusser, L.A.Tordsen, E.J. Atkinson, K.G. Rabe, and S.L.Slager. Identification of genes involved in alcohol consumption and cigarette smoking. BMC Genetics, 6:S112, 2005.

See Also

multic, gene.eff

Examples

Run this code
## Not run: 
# add2 <- addGE(list(bmi10, dia10, sys10), combine = 2, plotit = T, ylim=c(0,8), legend=F)
# add3 <- addGE(list(bmi10, dia10, sys10), combine = 3, plotit = F)
# lines(add3$cM, add3$lod, col=4, lwd=2, lty=4)
# key(corner=c(0,1), lines=list(lwd=2, col=1:4, lty=1:4),
#     text=list(c('BMI-Dia','BMI-Sys','Dia-Sys','BMI-Dia-Sys'), col=1:4))
# ## End(Not run)

Run the code above in your browser using DataLab