Multi_No: Run Multi-Trait without Annotation Model

In JustinaXie/MDATA: Multi-trait for De novo mutation Association Test with Annotations

Description

Run Multi-Trait without Annotation Model

Usage

Multi_No(
  data,
  N_1,
  N_2,
  pi_init = c(0.85, 0.05, 0.05, 0.05),
  threshold = 1e-06
)

Arguments

data	A dataframe with columns: Gene name, mutability, de novo mutation count for trait1, de novo mutation count for trait2
N_1	Cohort size for trait 1
N_2	Cohort size for trait 2
pi_init	Initial value for probabilities of genes under different assumptions. Default=c(0.85,0.05,0.05,0.05)
threshold	Threshold for EM algorithm. Default=1e-6

Value

The estimated model parameters and the posterior probabilities of genes under different assumptions

result	A dataframe that includes estimated posterior probabilities of risk genes for each trait and estimated posterior probability for shared risk gene
pi	Estimated proportion of genes under different assumptions
beta_trait1	Estimated log of gamma for trait 1
beta_trait2	Estimated log of gamma for trait 2
Z_mat	Estimated posterior probabilities of genes under different assumptions

JustinaXie/MDATA documentation built on Dec. 18, 2021, 2:34 a.m.