ddmodel.Rd
This model is applied on data of a single gene. It will infer effect size for both sample-level variable and positional level functional annotations. We used an EM algorithm to infer parameters.
ddmodel(mut, e, mr, fe, label, ...)
a matrix of mutation status 0 or 1, rows positions, columns are samples.
a vector,phenotype of each sample,
should match the columns of mut
and mr
a matrix, mutation rate of each sample at each mutation (log scale) that is not dependent on sample level factor
a vector, increased mutation rate at each position, depending on e (log scale),
should match the rows of mut
and mr