Description Usage Arguments Details Value Author(s) See Also Examples
The function epitest
takes a simple epidemic "twinstim"
model (one with epidemic = ~1
) and tests if the spatio-temporal
interaction invoked by the epidemic
model component is
statistically significant.
1 2 3 4 5 6 7 8 |
model |
a simple epidemic |
data |
an object of class |
tiles |
(only used by |
method |
one of the following character strings specifying the test method:
|
B |
the number of permutations for the Monte Carlo approach.
The default number is rather low; if computationally feasible,
|
eps.s,eps.t |
arguments for |
fixed |
optional character vector naming parameters to fix at their original
value when re-fitting the |
verbose |
the amount of tracing in the range |
compress |
logical indicating if the |
... |
further arguments for |
object,x |
an object of class |
which |
a character string indicating either the full ( |
teststat |
a character string determining the test statistic to plot, either
|
The test statistic is the reproduction number simpleR0
.
A likelihood ratio test of the supplied epidemic model against
the corresponding endemic-only model is also available.
The null distribution of the test statistic is obtained by a Monte Carlo permutation
approach (via permute.epidataCS
) similar to the
knox
test.
The plot
-method shows a truehist
of
the simulated null distribution together with the observed value.
The coef
-method extracts the parameter estimates from the B
permfits
(by default for the full model which = "m1"
).
a list (inheriting from "htest"
) with the following components:
method |
a character string indicating the type of test performed. |
data.name |
a character string giving the supplied |
statistic |
the observed test statistic. |
parameter |
the (effective) number of permutations used to calculate the p-value (only those with convergent fits are used). |
p.value |
the p-value for the test. For the |
In addition, if method != "LRT"
, the result will have the
following elements:
permfits |
the list of model fits (endemic-only and epidemic)
from the |
permstats |
a data frame with |
The plot
-method invisibly returns NULL
.
The coef
-method returns the B
x length(coef(model))
matrix of parameter estimates.
Sebastian Meyer
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | data("imdepi")
data("imdepifit")
## test for space-time interaction of the B-cases
## assuming spatial interaction to be constant within 50 km
imdepiB50 <- update(subset(imdepi, type == "B"), eps.s = 50)
imdfitB50 <- update(imdepifit, data = imdepiB50,
epidemic = ~1, epilink = "identity",
siaf = NULL, control.siaf = NULL,
start = c("e.(Intercept)" = 1e-6))
## simple likelihood ratio test
epitest(imdfitB50, imdepiB50, method = "LRT")
## permutation test (only few permutations for speed, in parallel)
et <- epitest(imdfitB50, imdepiB50, B = 4 + 25*surveillance.options("allExamples"),
verbose = 2 * (.Platform$OS.type == "unix"),
.seed = 1, .parallel = 2)
et
plot(et)
summary(coef(et))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.