Description Usage Arguments Value Examples
A multiple linear regression for familial or population data
1 2 3 4 |
formula |
referring to the column(s) in |
phenodata |
a data frame containing columns mentioned in |
genodata |
an object with genotypes to analyze. Several formats are allowed: |
kin |
a square symmetric matrix giving the pairwise kinship coefficients between analyzed
individuals. Under default |
nullmod |
an object containing parameter estimates under the null model. Setting |
regions |
an object assigning regions to be analyzed. This can be: |
sliding.window |
the sliding window size and step. Has no effect if |
mode |
the mode of inheritance: "add", "dom" or "rec" for additive, dominant or recessive mode, respectively. For dominant (recessive) mode genotypes will be recoded as AA = 0, Aa = 1 and aa = 1 (AA = 0, Aa = 0 and aa = 1), where a is a minor allele. Default mode is additive. |
ncores |
number of CPUs for parallel calculations. Default = 1. |
return.time |
a logical value indicating whether the running time should be returned. |
stat |
the statistic to be used to calculate the P values. One of "F" (default), "Chisq", "LRT". |
impute.method |
a method for imputation of missing genotypes. It can be either "mean" (default) or "blue". If "mean" the genotypes will be imputed by the simple mean values. If "blue" the best linear unbiased estimates (BLUEs) of mean genotypes will be calculated taking into account the relationships between individuals [McPeek, et al., 2004, DOI: 10.1111/j.0006-341X.2004.00180.x] and used for imputation. |
write.file |
output file name to write results as they come (sequential mode only). |
... |
other arguments that could be passed to |
A list with values:
results |
a data frame containing P values, numbers of variants and polymorphic variants for each of analyzed regions. |
nullmod |
an object containing the estimates of the null model parameters: heritability (h2), total variance (total.var), estimates of fixed effects of covariates (alpha), the gradient (df), and the total log-likelihood (logLH). |
sample.size |
the sample size after omitting NAs. |
time |
If |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | data(example.data)
## Run MLR with sliding window (default):
out <- MLR(trait ~ age + sex, phenodata, genodata, kin)
## Run MLR with regions defined in snpdata$gene and with
## null model parameters obtained in the first run:
out <- MLR(trait ~ age + sex, phenodata, genodata, kin,
out$nullmod, regions = snpdata$gene)
## Run MLR parallelized on two cores (this will require
## 'foreach' and 'doParallel' R-packages installed and
## cores available):
out <- MLR(trait ~ age + sex, phenodata, genodata, kin,
out$nullmod, ncores = 2)
## Run MLR with genotypes in VCF format:
VCFfileName <- system.file(
"testfiles/1000g.phase1.20110521.CFH.var.anno.vcf.gz",
package = "FREGAT")
geneFile <- system.file("testfiles/refFlat_hg19_6col.txt.gz",
package = "FREGAT")
phe <- data.frame(trait = rnorm(85))
out <- MLR(trait, phe, VCFfileName, geneFile = geneFile,
reg = "CFH", annoType = "Nonsynonymous")
## Run MLR with genotypes in PLINK binary data format:
bedFile <- system.file("testfiles/sample.bed",
package = "FREGAT")
phe <- data.frame(trait = rnorm(120))
out <- MLR(trait, phe, bedFile)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.