glmnet_pca
, scores_glmnet
, anova2
, and anova_glmnet
.
Functions are well documented.anova2
and anova_glmnet
now work when genotype matrix X
has rownames with special characters.
Before, for example if a locus name had the common format "chr:pos" where chr and pos are both integers, these functions would fail with this error: "Error in terms.formula(formula, data = data) : invalid model formula in ExtractVars".
Now all variables are quoted with backticks, which avoids this problem.
Tested to work when names have math operations such as "a+b" (treated as a single variable name now).
Additionally, an informative error is thrown if a locus name actually contains a backtick (only special character we can't handle).anova_glmnet
now returns -log10 p-values as a sparse matrix.
Pros: saves memory and runtime.
Cons: had to change to return -log10 p-values instead of p-values (to retain sparsity, so that untested variants are all zero)glmnet_pca
how supports cross validation (a wrapper around cv.glmnet
if cv = TRUE
).anova_glmnet_single
, for calculating scores on a single model/lambda (by default, approximately the best in terms of a precalculated cross-validation).
This is in contrast to the previous function anova_glmnet
, which calculates scores for all models/lambdas, which is much slower and not generally recommended anymore.anova_single
for scoring a model specified by locus indexes only.NEWS.md
slightly to improve its automatic parsing.LazyData: true
(to avoid a new "NOTE" on CRAN).URL
and BugReports
)Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.