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

Description Usage Arguments Value Author(s)

View source: R/flashfm.R

Description

internal function for calcAdjPP for the 2-trait case

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
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 Aug. 9, 2021, 7:36 a.m.