pepsimtest: Test a set of simulated data

In HegemanLab/ProteinTurnover: Protein Turnover Experiment Analysis

Description

Fit a separate model to each simulation and find the deviance, to see if the distribution of the deviance is really chi-squared.

Usage

pepsimtest(object, ask = TRUE, n = object$Nrep, model = object$parent)

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

## S3 method for class 'pepsimtest'
plot(x, ylim, ...)

Arguments

object	a pepsim object
ask	Since this may take a while, it will ask for permission before starting. Disable with ask=FALSE.
n	The number of simulations to fit models to.
model	A pepfit object with the model that you want to fit the simulated data to.
x	the pepsimtest object
...	to match the generic print function
ylim	the vertical limits of the plot

Value

A list of class "pepsimtest" containing:

sim.statistic	The deviances from the fitted models for each simulation
obs.test	The test of the original fitted model
sim	The original pepsim object

See Also

pepsim

Examples

## Not run: 
data(isodata)
a0 <- pepfit(data=isodata, Elements = list(N = 12,
    C = 45, H = 73, O = 15)))
a0.sim <- pepsim(a0, Nrep = 50)
a0.test.sim <- pepsimtest(a0.sim, n = 10, ask = FALSE)
a0.test.sim
plot(a0.test.sim)

## End(Not run)

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