gcBootSpline: Function to generate a bootstrap sample from a dataset of...

In grofit: The package was developed to fit fit many growth curves obtained under different conditions.

Description

Bootstrapping of a growth curve.

Usage

gcBootSpline(time, data, gcID = "undefined",
             control = grofit.control())

Arguments

time	Numeric vector containing the data for x-axes.
data	Numeric vector giving the growth values belonging to each element of time.
gcID	Vector (of any length) identifying the growth curve data.
control	Object of class grofit.control containing a list of options generated by the function grofit.control.

Details

A bootstrap sample (with replacement) is generated from the original data. For each of the bootstrap samples, the characteristic parameters λ, μ, A and integral are estimated.

Value

Generates an object of class gcBootSpline

raw.time	Raw data given to the function; equivalent to time.
raw.data	Raw data given to the function; equivalent to data.
gcID	Identifier, given to the function as gcID.
boot.time	Table of concentration values per column, resulting from each spline fit of the bootstrap.
boot.data	Table of response values per column, resulting from each spline fit of the bootstrap.
boot.gcSpline	List containing all gcFitSpline objects generated by the call of gcFitSpline.
lambda	Vector containing all the estimates from the bootstrap sample for the lag-phase.
mu	Vector containing all the estimates from the bootstrap sample for the maximum slope.
A	Vector containing all the estimates from the bootstrap sample for the maximum growth value.
integral	Vector containing all the estimates from the bootstrap sample for the integral.
bootFlag	Logical, indicating an empty bootstrap sample.
control	Object of class grofit.control containing list of options passed to the function as control.

See Also

gcFitModel, gcFitSpline, gcFit, summary.gcBootSpline, plot.gcBootSpline

Examples

x <- 1:30
y <- 1/(1+exp(0.5*(15-x)))+rnorm(30)/20
TestRun <- gcBootSpline(x,y,"ID",grofit.control(nboot.gc=50))
print(summary(TestRun))
plot(TestRun)

Example output

      mu.bt lambda.bt     A.bt integral.bt   stdmu.bt stdlambda.bt    stdA.bt
1 0.1707368  11.25355 1.037294    13.93107 0.02660015     2.699683 0.02421786
  stdintegral.bt ci90.mu.bt.lo ci90.mu.bt.up ci90.lambda.bt.lo
1      0.8351957     0.1269795      0.214494          6.812566
  ci90.lambda.bt.up ci90.A.bt.lo ci90.A.bt.up ci90.integral.bt.lo
1          15.69452    0.9974556     1.077132            12.55717
  ci90.integral.bt.up ci95.mu.bt.lo ci95.mu.bt.up ci95.lambda.bt.lo
1            15.30497     0.1186005      0.222873          5.962166
  ci95.lambda.bt.up ci95.A.bt.lo ci95.A.bt.up ci95.integral.bt.lo
1          16.54492    0.9898269     1.084761            12.29409
  ci95.integral.bt.up
1            15.56805

grofit documentation built on May 30, 2017, 4:08 a.m.