f.quad.plateau: Fitting quadratic plateau model using multiple initial values

In FertBoot: Fertilizer Response Curve Analysis by Bootstrapping Residuals

Description

f.quad.plateau fits quadratic plateau model using multiple initial values. The multiple initial values are randomly sampled in a "cube" of parameter space. More precisely, quadratic plateau model assumes y ~ (a + b * x + c *x^2) * (x <= -0.5*b/c) + (a + -b^2/(4 * c)) * (x > -0.5 * b/c).

Usage

f.quad.plateau(
  d,
  start = list(a = 1, b = 1, c = 1),
  plus_minus = 100,
  n.start = 1000,
  msg = FALSE
)

Arguments

d	data frame with two columns (x and y)
start	initial estimate for non-linear least square (default value: list(a = 1, b = 1, c = 1))
plus_minus	radius of random initial values (default: 100)
n.start	total number of initial points considered (default: 1000)
msg	logical flag whether printing progress

Value

f.quad.plateau returns a list of two components (if converged): nls.summary: summary of the fitted model; nls.model: nls object

Examples

set.seed(3)
x <- rep(1:300, each=4)
a <- 8; b <- 0.05; c <- -1e-4
y <- (a + b * x + c *x^2) * (x <= -0.5*b/c) + (a + -b^2/(4 * c)) * (x > -0.5 * b/c) +
    rnorm(length(x), sd=0.1)
d <- cbind(x,y)

# a converged example:
ans <- f.quad.plateau(d, start=list(a = 7, b = 0.02, c = 1e-5),
    plus_minus=10, n.start=10, msg=FALSE)

summary(ans$nls.model)

FertBoot documentation built on Feb. 20, 2021, 1:06 a.m.