DINT {RNOmni} | R Documentation |
Applies the rank-based inverse normal transformation (rankNorm
)
to the phenotype y
. Conducts tests of association between the loci in
G
and transformed phenotype, adjusting for the model matrix X
.
DINT(y, G, X = NULL, k = 3/8, test = "Score", simple = FALSE, parallel = FALSE)
y |
Numeric phenotype vector. |
G |
Obs by snp genotype matrix. |
X |
Model matrix of covariates and structure adjustments. Should include an intercept. Omit to perform marginal tests of association. |
k |
Offset applied during rank-normalization. See
|
test |
Either Score or Wald. |
simple |
Return the p-values only? |
parallel |
Logical indicating whether to run in parallel. Must register parallel backend first. |
If simple=T
, returns a vector of p-values, one for each column
of G
. If simple=F
, returns a numeric matrix, including the
Wald or Score statistic, its standard error, the Z-score, and the p-value.
Basic association test BAT
, indirect INT IINT
, omnibus INT OINT
.
## Not run: set.seed(100); # Design matrix X = cbind(1,rnorm(1e3)); # Genotypes G = replicate(1e3,rbinom(n=1e3,size=2,prob=0.25)); storage.mode(G) = "numeric"; # Phenotype y = exp(as.numeric(X%*%c(1,1))+rnorm(1e3)); # Association test p = DINT(y=y,G=G,X=X,simple=T); ## End(Not run)