pepsim: Create a set of simulated data

In HegemanLab/ProteinTurnover: Protein Turnover Experiment Analysis

Description

Create a set of simulated data, to check if the change in deviance between two models is really chi-squared.

Usage

pepsim(object, Nobs = 1/pepfittest(object)$sigma2, Nrep = 10)

## S3 method for class 'pepsim'
print(x, ...)

## S3 method for class 'pepsim'
plot(x, n = 30, ...)

Arguments

object	a pepfit object
Nobs	the sample size to be used for the simulation
Nrep	how many simulations to run
x	a pepsim object
...	currently ignored
n	how many simulations to plot

Value

A list of class pepsim containing:

sim	a data frame of simulated relative abundances for each day and count value
Nobs	Nobs, as input
Nrep	Nrep, as input
parent	the original pepfit object

Examples

data(isodata)
a0 <- pepfit(data=isodata, Elements = list(N = 12,
    C = 45, H = 73, O = 15))
a0.sim <- pepsim(a0, Nrep = 50)
plot(plot(a0.sim))

HegemanLab/ProteinTurnover documentation built on May 6, 2019, 11:50 p.m.