Description Usage Arguments Value
Performs a genome scan on heteroscedastic data.
1 |
G |
A s by p matrix of genotypes, where s is the number of strains and p is the number of snps to be tested. This matrix can have missing values, and each SNP should be coded as 0, 1, 0.5 or NA |
y |
A N length vector of phenotypes for each individual organism, where N is the total number of individuals |
strains |
A s by N incidence matrix that maps every individual to a strain |
X |
A s by q matrix of covariates (optional) |
K |
A s by s genomic relationship matrix. Will be calculated if unspecified. |
weights |
A string specifying the weights to be used. The following are permitted: "none", "samplevars", "limma", "counts", and "user" |
user_weights |
A s length vector of weights for each strain, used if weights = "user" |
A list containing:
pvalue: A p length vector of p-values for every snp
ML0: A p length vector of log maximum likelihood under the null hypothesis
ML1: A p length vector of log maximum likelihood under the alternative hypothesis
beta: A p length vector of regression parameter estimate for slope under the alternative hypothesis
s2: A p length vector of the sum of variance components under the null hypothesis
h2: A p length vector of the heritability estimate under the null hypothesis
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.