Learn R Programming

multic (version 0.2.1)

multic: Create a multic object

Description

Calculate the polygenic and major gene models for quantitative trait linkage analysis using variance components approach.

Usage

multic(formula, data = sys.parent(), famid, id, dadid, momid, sex,
       mloci.out = NULL, share.out = "kinship", longitudinal = FALSE,
       subset = NULL, ascertainment = NULL,
       control = multic.control(...), ...)

Arguments

formula
a formula object, with the traits on the left of a ~ (tilde) operator and the covariates, separated by + operators, on the right. The traits may be a single numeric vector or a matrix. Commonly, traits are aggragated together
famid
integer, numeric, or character vector specifying each individual's family identifier. Members of the same family must have the same famid and each family must have a unique famid. Any missing data will result in an error messa
id
integer, numeric, or character vector specifying each individual's identifier. Members of the same family must have a unique id within the family. id does not have to be universally unique among all individuals. Any missing
dadid
integer, numeric, or character vector specifying each individual's father identifier. This father identifier must have the same famid as the individual. Any missing data will result in an error message and the termination of multic
momid
integer, numeric, or character vector specifying each individual's mother identifier. This mother identifier must have the same famid as the individual. Any missing data will result in an error message and the termination of multic
sex
integer, numeric, or character vector specifying each individual's sex. Acceptable forms of sex-coding are "M", "m", or 1 for male and "F", "f", or 2 for female. Any missing da
data
data.frame in which to interpret the variables named in formula, famid, id, dadid, momid, sex, subset, and ascertainment arguments. If
mloci.out
a character value specifying a path to an mloci.out (or similarly formatted) file. Specifying a non-empty mloci.out file will allow multic to calculate sporadic models using the ibd (identity by decent) information in the mloci.out file. Due to the gene
share.out
a character value specifying a path to a share.out (or similarly formatted) file. This file contains the amount of genetic material shared between each family member pair based on family structure only. It also contains boolean values to indicate whethe
longitudinal
logical flag: if TRUE, then fomula will be interpreted as a longitudinal model. In this case, the fomula argument requires special formatting as described. The number of traits on the left side of the ~ (tilde) is the number o
subset
a logical vector specifying which subset of the rows in data to use in the fit.
ascertainment
vector specifying each individual's ascertainment (effected) status. Acceptable forms of ascertainment are T, TRUE, or 1 for a proband (effected) and F, FALSE, or 0 for a non-p
control
list of iteration and algorithmic constants. See multic.control for their names and default values. These can also be given directly as arguments to multic itself, instead of through a multic.control ob
...
further arguments passed to multic.control to alter multic's default behavior.

Value

  • an object of class "multic". See multic.object for more details.

Side Effects

Many temporary files are created during multic's execution. These files are deleted afterwards (by default). If they are not deleted (due to a crash or some other unexpected action), use the included function clean() to delete them. Also, multic copies, gunzip's, and removes the copies of share.out and mloci.out (if specified).

References

Amos, C. I. (1994). "Robust variance-components approach for assessing genetic linkage in pedigrees." American Journal of Human Genetics 54(3): 535-543.

Almasy, L. and J. Blangero (1998). "Multipoint quantitative-trait linkage analysis in general pedigrees." American Journal of Human Genetics 62(5): 1198-1211.

Details

See the technical report.

See Also

multic.object, multic.control, phi2share, solar2mloci, solar2multic, sw2mloci

Examples

Run this code
# Call multic with a univariate formula with two covariates and no
# markers (no mloci.out argument).
fit.ibd.uni <- multic(k.trig ~ sex.x + agexam,
                      data = ped.phen.data,
                      famid, id, fa, mo, sex.x,
                      share.out = 'multicInput/share.out')

# Call multic with a bivariate formula with three covariates, no
# markers (no mloci.out argument), and calculate the family log
# likelihoods.
fit.ibd.bi <- multic(cbind(k.trig, k.chol) ~ sex.x + agexam + agexam2,
                     data = ped.phen.data,
                     famid, id, fa, mo, sex.x,
                     share.out = 'multicInput/share.out',
                     calc.fam.log.liks = TRUE)

# Call multic with a longitudinal formula with six covariates letting
# the kinship library calculate the share.out argument.
long.fit <- multic(cbind(sbpA, sbpB, sbpC) ~
                     sexA + ageA + bmiA + generA + ageAg + smkA +
                     sexB + ageB + bmiB + generB + ageBg + smkB +
                     sexC + ageC + bmiC + generC + ageCg + smkC,
                   data = long.data,
                   famid, id, dadid, momid, sex,
                   longitudinal = TRUE)

Run the code above in your browser using DataLab