Example twinstim Fit for the imdepi Data
data("imdepifit")an object of class "twinstim"
common methods for "twinstim" fits,
exemplified using imdepifit, e.g.,
summary.twinstim, plot.twinstim,
and simulate.twinstim
data("imdepi", "imdepifit")
## Not run: 
## reproduce "imdepifit"
myimdepifit <- twinstim(
    endemic = addSeason2formula(~ offset(log(popdensity)) + I(start/365-3.5),
                                S = 1, period = 365, timevar = "start"),
    epidemic = ~ type + agegrp, siaf = siaf.gaussian(),
    data = imdepi, subset = !is.na(agegrp),
    optim.args = list(control = list(reltol = sqrt(.Machine$double.eps))),
    ## the historical default for reltol is 1e-6, which is rather large
    model = FALSE, cumCIF = FALSE
)
stopifnot(all.equal(imdepifit, myimdepifit))
## End(Not run)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.