Description Format Source References See Also Examples
Example of a model structure holding parameter estimates for the model in
mesa.model
using estimate.STmodel
. Estimation
results are also provided for models including spatio-temporal covariates.
A list with elements, see the return description in
estimate.STmodel
.
Contains parametere estimates for the Spatio-Temporal model applied
to monitoring data from the MESA Air project, see
Cohen et.al. (2009) and mesa.data.raw
for details.
M. A. Cohen, S. D. Adar, R. W. Allen, E. Avol, C. L. Curl, T. Gould, D. Hardie, A. Ho, P. Kinney, T. V. Larson, P. D. Sampson, L. Sheppard, K. D. Stukovsky, S. S. Swan, L. S. Liu, J. D. Kaufman. (2009) Approach to Estimating Participant Pollutant Exposures in the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air). Environmental Science & Technology: 43(13), 4687-4693.
estimate.STmodel
for parameter estimation.
createSTmodel
for creation of the originating STmodel
object.
Other example data: MCMC.mesa.model
,
est.cv.mesa
, mesa.data.raw
,
mesa.model
, pred.mesa.model
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 | ##load a model object
data(mesa.model)
##create vector of initial values
dim <- loglikeSTdim(mesa.model)
x.init <- cbind(c( rep(2, dim$nparam.cov-1), 0),
c( rep(c(1,-3), dim$m+1), -3, 0))
rownames(x.init) <- loglikeSTnames(mesa.model, all=FALSE)
## Not run:
##estimate parameters
est.mesa.model <- estimate(mesa.model, x.init, hessian.all=TRUE)
## End(Not run)
##time consuming estimation, load pre-computed results instead
data(est.mesa.model)
#estimation results
print(est.mesa.model)
##compare the estimated parameters for the two starting points
est.mesa.model$summary$par.all
##and values of the likelihood (and convergence info)
est.mesa.model$summary$status
##extract the estimated parameters and approximate uncertainties
x <- coef(est.mesa.model)
##compare estimated parameters
##plot the estimated parameters with uncertainties
par(mfrow=c(1,1),mar=c(13.5,2.5,.5,.5))
with(x, plot(par, ylim=range(c(par-1.96*sd, par+1.96*sd)),
xlab="", xaxt="n"))
with(x, points(par - 1.96*sd, pch=3))
with(x, points(par + 1.96*sd, pch=3))
abline(h=0, col="grey")
##add axis labels
axis(1, 1:length(x$par), rownames(x), las=2)
## Not run:
##example using a few fixed parameters
x.fixed <- coef(est.mesa.model)$par
x.fixed[c(1,2,5:9)] <- NA
est.fix <- estimate(mesa.model, x.init, x.fixed, type="p")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.