initialize_chains: A Gibbs Sampler for genotyping copy number in trios

In githubmpc/marimba2: Gibbs sampler for genotyping CNPs

Description

This function allows you to obtain a posterior probability distribution for copy number genotyping on Illumina arrays for an fixed genetic model where tau1=1, tau2=0.5 and tau3=0

Usage

initialize_chains(states, N, K, S)

Arguments

S	The number of iterations for the Gibbs Sampler (suggest default 10000).
y	Data$response where data is as generated by the simulate_data function.
mu	Defaults to -.02.
xi	Defaults to sqrt(200).
nu	Defaults to 3.
lambda	Defaults to 0.05.
alpha	Alpha1 and alpha2 and psi refer to HWE parameters as specified in Cardin et al (alpha=w, psi=lambda).

githubmpc/marimba2 documentation built on May 17, 2019, 9:11 a.m.