lognormal {drcte}

Log-normal distribution of times to an event

Description

This function provides the truncated log-normal cumulative distribution of event times, to be used to fit time-to-event methods to the time-course of seed germination for a seed lot.

Usage

lognormal(fixed = c(NA, NA, NA), names = c("b", "d", "e"))
lognormalSurv(fixed = c(NA, NA, NA), names = c("b", "d", "e"))

Arguments

fixed

numeric vector. Specifies which parameters are fixed and at what value they are fixed. NAs for parameter that are not fixed. It is often used to fix the 'd' parameter to 1, which provides the usual log-logistic distribution, with no ungerminated fraction.

names

a vector of character strings giving the names of the parameters. The default is reasonable.

Details

The log-normal distribution of germination times is parameterised as:

P(t) = d * pnorm(b*(log(t + 0.000001) - log(d)))

where 't' is the time and 'P(t)' is the proportion of germinated seeds at time = t. LogNormalSurv is parameterised as usual in survival analysis:

P(t) = d * pnorm(b*(log(t + 0.000001) - d))

i.e. with 'e' on a log-scale and a log-link on 'b'. With no cured fraction (d = 1), logNormalSurv gives same results as 'survereg' in the 'survival package', with 'dist = "lognormal"'.

Value

The value returned is a list containing the nonlinear function, the self starter function and the parameter names.

Author(s)

Andrea Onofri

References

Kleinbaum, D.G., Klein, M., 2005. Survival analysis. Springer Science, USA.

Examples

library(drcte)
data(chickweed)
modb <- drmte(count ~ start + end,
               data=chickweed,
               fct = lognormal())