No packages match

iconic - Causal Model Selection with Genetic Instruments and Negative Controls

Provides a model selection workflow for causal inference with genetic instruments and negative controls in observational omics data. The package fits eight estimators of the natural direct and indirect effects (NDE/NIE), diagnoses which are valid for the user's data via instrument-strength checks and negative-control validity screens (A1, A2, A2'), stress-tests them against confounding and pleiotropy violations, and recommends the estimator most likely to be unbiased. A hybrid generative texture model (a torch GAN for sample-level structure and a Gaussian copula for the mediator panel) lets the sensitivity analysis mirror the marginal and joint structure of the user's cohort. The modular API comprises iconic_data() to standardize data, iconic_diagnose() to assess estimator eligibility, iconic_estimate() to fit all eligible estimators, iconic_sensitivity() to sweep instrument exogeneity, iconic_recommend() to rank estimators by identification strength and robustness, and iconic_prospect() to plan future studies. infer_confounding() estimates held-fixed confounding parameters from the data. Helper functions construct exposure instruments (polygenic scores from GWAS summary statistics or PGS Catalog scoring files), mediator instruments (cis-eQTL scans and elastic-net predicted-expression composites), and negative-control panels (principal components of residualized omics matrices), and data can be imported directly from SummarizedExperiment objects. Time-to-event outcomes are supported on the Cox log-hazard-ratio and restricted-mean-survival-time scales. Binary outcomes are supported on the log-odds-ratio (logistic two-stage predictor-substitution outcome stage) and risk-difference (linear probability model) scales. The torch package is required for the generative texture model (sensitivity and prospective analysis); all other functionality works without it.

Last updated

statisticalmethodgeneticsmultiplecomparisonregressiontranscriptomicsrnaseqsurvival

4.18 score