Single_No: Run Single-Trait without Annotation Model

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

Description

Run Single-Trait without Annotation Model

Usage

Single_No(data, dn_col, N_1, pi_init = 0.1, threshold = 1e-06)

Arguments

data	A dataframe with columns: Gene name, mutability, de novo mutation count for traits
dn_col	A character indicating the column name of de novo mutation counts for the trait
N_1	Cohort size for the trait
pi_init	Initial value for risk gene proportion. Default=0.1
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 vector, the second value represents risk gene proportion
beta	Estimated log of gamma for risk genes
Z_mat	Estimated posterior probabilities of genes under different assumptions

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