Q_unibichrom: Calculates Q-matrix for univariate profile likelihoods from...

Description Usage Arguments Value

Description

Calculates Q-matrix for univariate profile likelihoods from BiChroM model

Usage

1
Q_unibichrom(log.pars, log.theta0, size, param.name)

Arguments

log.pars

vector of size 9 indicating values in ln of nuisance parameters (excluding one of interest for the profile likelihood)

log.theta0

vector of size 1 indicating values parameters of interest in BiChroM. Options are (lambda0, lambda1, mu0, mu1, rho0, rho1, q01,q10, e01, e10)

size

Maximum number of chromosomes in the sample (recommended no more than 50, states larger than that should be coded as 51)

param.name

vector of size 1 indicating the names of parameters of interest. E.g. "rho0"

Value

Q a sparse matrix of size 2*(size+1)


fmichonneau/chromploid documentation built on May 16, 2019, 1:43 p.m.