grammar {GenABEL} | R Documentation |
Fast approximate score test for association between a trait and genetic polymorphism,
in samples of related individuals. When used with argument "times=1", it is equivalent to
running qtscore
on polygenic residuals from polygenic
.
However, it does not produce correct results with permutations, because the raw
trait values, which are not exchangable, are permuted. Use qtscore
on polygenic residuals when you wnat to have empirical GW significance with GRAMMAR
method.
grammar(h2object,data,snpsubset,idsubset,strata,times=1,quiet=FALSE,bcast=10,clambda=FALSE,propPs=1.0)
h2object |
An object returned by polygenic polygenic mixed model analysis
routine. The sub-objects used are measuredIDs, residualY, h2an$estimates (last element,
total variance, only), and InvSigma.
One can supply grammar with a fake h2object, containing these list elements.
|
data |
An object of gwaa.data-class |
snpsubset |
Index, character or logical vector with subset of SNPs to run analysis on.
If missing, all SNPs from data are used for analysis. |
idsubset |
Index, character or logical vector with subset of IDs to run analysis on.
If missing, all people from data/cc are used for analysis. |
strata |
Stratification variable. If provieded, scores are computed within strata and then added up. |
times |
If more then one, the number of replicas to be used in derivation of empirical genome-wide significance. NOTE: do not use times > 1 unless you are really sure you understand what you are doing! |
quiet |
do not print warning messages |
bcast |
If the argument times > 1, progress is reported once in bcast replicas |
clambda |
If inflation facor Lambda is estimated as lower then one, this parameter controls if the original P1df (clambda=TRUE) to be reported in Pc1df, or the original 1df statistics is to be multiplied onto this "deflation" factor (clambda=FALSE). With GRAMMAR, Lambda is expected ot be less than 1. If a numeric value is provided, it is used as a correction factor. |
propPs |
proportion of non-corrected P-values used to estimate the inflation factor Lambda,
passed directly to the estlambda |
Approximate score test is performed using the formula
$$ σ^4 frac{((G-E[G]) V^{-1} residualY)^2}{(G-E[G]) (G-E[G])} $$
where $sigma^4$ is the square of the residual
variance, $G$ is the vector of genotypes (coded 0, 1, 2) and $E[G]$ is
a vector of (strata-specific) mean genotypic values; $V^{-1}$ is the
InvSigma and $residualY$ are residuals from the trait analysis
with polygenic
procedure.
Compared to score test implemented in mmscore
, grammar
test is faster and computation time grows only linearly with the
number of subjects (with mmscore
this relation is
quadratic). While raw P1df from grammar are not quite correct,
the GC p-values correspond very closely to these from the mmscore
.
Object of class scan.gwaa-class
; only 1 d.f. test is
implemented currently.
Yurii Aulchenko
Chen WM, Abecasis GR (2007) Family-based assoication tests for genome-wide association scans.
grammar
,
qtscore
,
plot.scan.gwaa
,
scan.gwaa-class