pepcurve: Fit an exponential curve to fitted pi, r, and alpha values
In HegemanLab/ProteinTurnover: Protein Turnover Experiment Analysis

Description Usage Arguments Value Examples

Fit an exponential curve to fitted pi, r, and alpha values

pepcurve(x, ...)

## S3 method for class 'pepfit'
pepcurve(x, nab = TRUE, fillNA = FALSE, ...)

## Default S3 method:
pepcurve(x, nab = NA, fillNA = FALSE, intercept = c(pi =
  NA, alpha = NA, r = NA), ...)

makeCurve(fit, plim = NA, isotope.method = c("other", "incorporation"), ...)

getCurve(curve)

`x`	A pepfit object, or a data frame with columns TimePoint, and at least one of pi, r, and alpha
`...`	used to pass nab to the appropriate function
`nab`	for pepcurve.default, sets the natural abundance for use in either intercept or the limit of pi; NA finds the best fit instead for pepcurve.pepfit, can be numeric, or TRUE to use the known natural abundance from the pepfit, or FALSE to find the best fit
`fillNA`	should output have NA?
`intercept`	set intercept values for pi, alpha, or r
`fit`	fit from makeFit
`plim`	sent to pepcurve as nab parameter
`isotope.method`	used to set up intercepts correctly for incorporation experiments
`curve`	a fitted curve

A list of class "pepcurve".

The first element has four values. The first three are from the pi fit: the fitted value at day 0, the lower limit, and the halflife. The last is the halflife from the r fit, for which the fitted value at day 0 is set to 1 and the lower limit is set to 0.

Remaining elements contain the pi and r values from the pepfit model and the fitted curve objects for pi and r, I think. This should be checked to make sure it's still true.

data(isodata)
a0 <- pepfit(TimePoint, RelAb, Channel, data=isodata, Elements = list(N = 12,
     C = 45, H = 73, O = 15))
 
c0 <- pepcurve(a0)
c0
plot(c0)