logitexp: Exponential logit model for two variables
Description
Fit a logit model of v.density on v.fever v.density with
a exponential coefficient for the v.density
Usage
logitexp(v.fever, v.density)
Arguments
v.fever
numeric vector of 0/1 indicating fever or equivalent
v.density
numeric vector of values >= 0 indicating the density
Value
S3 object of class afmodel with 4 components: data, model, coefficients
and the estimated attributable fraction.
Details
logit(v.fever) ~ beta * (v. density ^ tau)
This corresponds to the model 3 describe by
Smith, T., Schellenberg, J.A., Hayes, R., 1994.
Attributable fraction estimates and case definitions for malaria
in endemic areas. Stat Med 13, 2345<U+2013>2358.
# NOT RUN {{
# Get the sample datahead(malaria_df1)
fit <- logitexp(malaria_df1$fever, malaria_df1$density)
fit
senspec(fit, c(1,100,500,1000,2000,4000,8000,16000, 32000,54000,100000))
}
# }