calcD12: internal function for calcAdjPP for the 2-trait case

In jennasimit/flashfm: flexible and shared information fine-mapping

internal function for calcAdjPP for the 2-trait case

Description

internal function for calcAdjPP for the 2-trait case

Usage

calcD12(
  mod1,
  mod2,
  T1,
  T2,
  beta,
  SSy,
  Sxy,
  xcovo,
  Mx,
  Nqq,
  Vres,
  covY,
  nsnpspermodel
)

Arguments

mod1	vector of model indices for trait 1
mod2	vector of model indices for trait 2
T1	index of trait 1
T2	index of trait 2
beta	list of joint beta estimates for each trait
SSy	matrix of trait cross-products
Sxy	matrix with each column being the cross-product between SNPs and a trait
xcovo	covariance matrix of c("one"=1,X)
Mx	vector of SNP means
Nqq	has N[i,j] = no. with both trait i and trait j measured
Vres	list of variance residuals
covY	covariance matrix of traits
nsnpspermodel	list of number of SNPs per model for each model in mod1, mod2

Value

lbf12-lbf1-lbf2

Author(s)

Jenn Asimit

jennasimit/flashfm documentation built on July 31, 2022, 7:32 p.m.