sim.fit.mult.add.data.batch: Simulate fitness data under multiplicative and additive...

Description Usage Arguments Details Value Examples

View source: R/simulate_data_fnxs.R

Description

Simulate fitness data under multiplicative and additive models

Usage

1
2
sim.fit.mult.add.data.batch(epi.model, mut.vals, coe.vals, sig.vals, w.wt = 1,
  n.reps.ea = 100, print.status = FALSE, outdir, wts)

Arguments

epi.model

mult/add Epistasis model to simulate under.

mut.vals

Vector of number of mutations to simulate

coe.vals

Vector of stickbreaking coefficients to simulate

sig.vals

Vector of sigma values to simulate

w.wt

Fitness of the wildtype. Default 1.

n.reps.ea

Number of replicates per parametric condition

print.status

TRUE/FALSE. Should loop counters be printed.

outdir

The path to write output files to (see details about file names).

wts

Weight to give mutation on wildtype background vs other backgrounds. Default is c(2,1).

Details

Function contains a loop for combining each value of mut.vals, coe.vals and sig.vals, generating data under the specified model and then fitting it.

Results are written to files; the name of the output files are formed by concatenating the outpath argument to the epi.model argument. Separate files are generated for each number of mutations (because the dimensionality of the output file changes with the number of mutations). The output files are tab-delimited text files with one row per replicate. The first 5 columns provide the parameter values and the rest of the columns give parameter estimates and measures of fit.

Value

Nothing. Instead results are written to output files and deposited in inst/extdata. The files are named by appending the method

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
## Not run: 
sim.fit.mult.add.data.batch("mult",
  c(3,4,5),
  c(0.1, 0.3, 0.5),
  c(0.02, 0.05, 0.08),
  1,
  100,
  print.status=TRUE,
  outdir="~/Desktop",
  c(2,1))
  
## End(Not run)

Stickbreaker documentation built on May 29, 2017, 9:01 a.m.