drBootSpline: Function to create a bootstrap sample of splines fitted to...

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

Description

The function passes the numeric vectors conc and test to the function drFitSpline, which fits a smoothed spline to the data and estimates the EC50. Calling drFitSpline several times (specified by the options set in control) creates a bootstrap sample in form of a list of objects of class drFitSpline.

Usage

drBootSpline(conc, test, drID = "undefined", 
             control = grofit.control())

Arguments

conc	Numeric vector, containing concentration information for dose-response. curve fit and EC50 estimation.
test	Numeric vector, containing response values related to concentrations.
drID	Character, identifying the dose-response data.
control	Object of class grofit.control containing a list of options generated by the function grofit.control.

Value

Generates an object of class drBootSpline

raw.conc	Raw data given to the function; equivalent to conc.
raw.test	Raw data given to the function; equivalent to test.
drID	String identifier given to the function; equivalent to drID.
boot.conc	Table of concentration values per column, resulting from each spline fit of the bootstrap.
boot.test	Table of response values per column, resulting from each spline fit of the bootstrap.
boot.drSpline	List containing all drFitSpline objects generated by the call of drFitSpline.
ec50.boot	Vector of estimated EC50 values from each bootstrap entry.
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

drFitSpline, summary.drBootSpline, plot.drBootSpline

Examples

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

Example output

=== Bootstrapping of dose response curve ==========
--- EC 50 -----------------------------------------

Mean  :  24.7465865865866 StDev :  0.113672136805712 
90% CI:  24.5595959215412 90% CI:  24.933577251632
95% CI:  24.5237891984474 95% CI:  24.9693839747258

  meanEC50    sdEC50 ci90EC50.lo ci90EC50.up ci95EC50.lo ci95EC50.up
1 24.74659 0.1136721     24.5596    24.93358    24.52379    24.96938
  meanEC50.orig ci90EC50.orig.lo ci90EC50.orig.up ci95EC50.orig.lo
1      24.74659          24.5596         24.93358         24.52379
  ci95EC50.orig.up
1         24.96938
dev.new(): using pdf(file="Rplots1.pdf")

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