Using summary statistics, calculates joint effect estimates

multibeta(mod, beta1, Gmat, N, ybar, is.snpmat, raf = NULL)

Arguments

mod

joint SNP model with snps separated by "%" e.g. "snp1%snp2"

beta1

named vector single-SNP effect estimates; each effect estimate has name given by SNP id and ids must appear in Gmat columns

Gmat

genotype matrix with SNP columns and individuals rows; could be reference matrix or from sample

N

sample size for trait

ybar

trait mean

is.snpmat

logical taking value TRUE when Gmat is a genotype matrix and FALSE when Gmat is a SNP covariance matrix

raf

named vector of SNP reference allele frequencies (must match SNP coding used for effect estimates); only needed if Gmat is a covariance matrix

Value

joint effect estimates for SNPs in model given by mod