diseaseModel | R Documentation |
Create a disease model in form of a disModel
object, for use in e.g.
lod()
.
diseaseModel(model = NULL, chrom = NULL, penetrances = NULL, dfreq = NULL)
model |
An existing
In all of the above, the disease is assumed to be fully penetrant, and with a disease allele frequency of 1e-5. |
chrom |
Either "AUTOSOMAL" or "X". Lower case versions are allowed and will be converted automatically. |
penetrances |
For autosomal models, a numeric of length 3 corresponding
to For X-linked models, a list of two vectors named |
dfreq |
A number in |
An object of class disModel
, which is a list with entries chrom
,
penetrances
and dfreq
.
lod()
# Fully penetrant AD model:
m1 = diseaseModel(model = "AD")
# The above is equivalent to
m2 = diseaseModel(chrom = "Aut", penetrances = c(0,1,1), dfreq = 1e-5)
stopifnot(identical(m1, m2))
# X-linked recessive model:
m3 = diseaseModel(model = "XR", dfreq = 0.01)
# Long version of the above:
m4 = diseaseModel(chrom = "X", penetrances = list(male = c(0,1), female = c(0,0,1)),
dfreq = 0.01)
stopifnot(identical(m3, m4))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.