get_yty: Get transformed statistics: yty

In USCbiostats/hJAM: Hierarchical Joint Analysis of Marginal Summary Statistics

Get transformed statistics: yty

Description

To calculate sufficient statistics based on summary statistics. This yty estimate follows Yang et al. (2012) Nat Gen. Marginal estimates from one SNP will produce one yty estimates. Yang suggests taking the median across all SNPs to obtain a robust estimate. Here we record all yty estimates and output both the median and the entire vector.

To calculate sufficient statistics based on summary statistics. This yty estimate follows Yang et al. (2012) Nat Gen. Marginal estimates from one SNP will produce one yty estimates. Yang suggests taking the median across all SNPs to obtain a robust estimate. Here we record all yty estimates and output both the median and the entire vector.

Usage

get_yty(maf, N_outcome, betas, betas.se)

get_yty(maf, N_outcome, betas, betas.se)

Arguments

maf	A vector of minor allele frequencies
N_outcome	Sample size in the GWAS where we obtained 'betas'
betas	A vector of marginal estimates of effect sizes (betas for continuous outcome; logOR for binary outcome)
betas.se	A vector of the standard errors of marginal effect estimates ('betas').

Value

median of yty estimates across all SNPs; and a vector of all yty estimates

median of yty estimates across all SNPs; and a vector of all yty estimates

USCbiostats/hJAM documentation built on Jan. 26, 2024, 5:27 p.m.