twinSIR_profile {surveillance} | R Documentation |
Function to compute estimated and profile likelihood based confidence intervals. Computations might be cumbersome!
## S3 method for class 'twinSIR' profile(fitted, profile, alpha = 0.05, control = list(fnscale = -1, factr = 10, maxit = 100), ...)
fitted |
an object of class |
profile |
a list with elements being numeric vectors of length 4. These vectors must
have the form
|
alpha |
(1-alpha)*100% profile likelihood based confidence
intervals are computed. If |
control |
control object to use in |
... |
unused (argument of the generic). |
list with profile log-likelihood evaluations on the grid and highest likelihood
and Wald confidence intervals. The argument profile
is also returned.
Michael Höhle and Sebastian Meyer
if (surveillance.options("allExamples")) { data("foofit") prof <- profile(foofit, list(c(1,NA,NA,5), c(3,NA,NA,0), c(4, 0.5, 1.1, 10))) prof ## there is also a plot-method for "profile.twinSIR" plot(prof) }