ELJAlinear | R Documentation |
A tool for Environment-Wide Association Studies (EnvWAS / EWAS) namely repeated analyses allowing to estimate the relationships between several environmental factors and a health events. This function is especially for linear regressions and allows the addition of adjustment variables.
ELJAlinear(
var,
var_adjust = NULL,
data,
manplot = TRUE,
nbvalmanplot = 100,
Bonferroni = FALSE,
FDR = FALSE,
manplotsign = FALSE
)
var |
A categorical and binary variable. It is generally your outcome. |
var_adjust |
A vector containing the names of the fixed adjustment variables for all the models. |
data |
A dataframe containing all the variables needed for the analysis. |
manplot |
Generate a Manhattan plot of the results of the analysis. |
nbvalmanplot |
The number of variables to include in each Manhattan plot. |
Bonferroni |
Add a dashed bar to the Manhattan plot showing the Bonferroni significance level. |
FDR |
Add a dashed bar to the Manhattan plot showing the False Discovery Rate (Benjamini-Hochberg method) significance threshold. NA if all p-values > FDR corrected p-values. |
manplotsign |
Generates a Manhattan plot with only significant results (p<0.05). |
A Dataframe with results for each variable of the model.
Dunn OJ. Multiple Comparisons Among Means. Journal of the American Statistical Association. 1961;56(293):52โ64. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society: Series B (Methodological). 1995;57(1):289โ300. MLBench ยท Distributed Machine Learning Benchmark. Available from: https://mlbench.github.io/ Smith JW, Everhart JE, Dickson WC, Knowler WC, Johannes RS. Using the ADAP Learning Algorithm to Forecast the Onset of Diabetes Mellitus. Proc Annu Symp Comput Appl Med Care. 1988 Nov 9;261โ5.
### Loading the PIMA dataset contained in the mlbench package
library(mlbench)
data(PimaIndiansDiabetes)
### Using ELJAlinear to perform EWAS analysis
ELJAlinear(var = 'pregnant',data = PimaIndiansDiabetes,manplot = TRUE,
Bonferroni = TRUE,FDR = TRUE, nbvalmanplot = 30, manplotsign = FALSE)
results
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.