ABneutralSOMA: Model with no selection (outneutral)

Description Usage Arguments Value Examples

View source: R/ABneutralSOMA.R

Description

This model assumes that somatically heritable gains and losses in cytosine methylation are selectively neutral.

Usage

1
ABneutralSOMA(pedigree.data, p0uu, eqp, eqp.weight, Nstarts, out.dir, out.name)

Arguments

pedigree.data

pedigree data.

p0uu

initial proportion of unmethylated cytosines.

eqp

equilibrium proportion of unmethylated cytosines.

eqp.weight

weight assigned to equilibrium function.

Nstarts

iterations for non linear LSQ optimization.

out.dir

output directory.

out.name

output file name.

Value

ABneutralSoma RData file.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
#Get some toy data
inFile <- readRDS(system.file("extdata/soma/","outputSoma.rds", package="AlphaBeta"))
pedigree <- inFile$Pdata
p0uu_in <- inFile$tmpp0
eqp.weight <- 0.001
Nstarts <- 2
out.name <- "ABneutralSOMA_CG_estimates"
out <- ABneutralSOMA(pedigree.data = pedigree,
                  p0uu=p0uu_in,
                  eqp=p0uu_in,
                  eqp.weight=eqp.weight,
                  Nstarts=Nstarts,
                  out.dir=getwd(),
                  out.name=out.name)

summary(out)

AlphaBeta documentation built on Nov. 8, 2020, 6:30 p.m.