proflik.glsm {geoRglm} | R Documentation |
Computes two dimenaional profile likelihood for the
parameters
(phi, nugget.rel) for a model previously
derived using the function likfit.glsm
.
proflik.glsm(mcmc.obj, obj.likfit.glsm, phi.values, nugget.rel.values, messages, ...)
mcmc.obj |
object with the Monte Carlo simulations and corresponding approximating density.
This object should be an output from the function
prepare.likfit.glsm . |
obj.likfit.glsm |
Output file from
likfit.glsm .
|
phi.values |
set of values of the parameter phi for which the
profile likelihood will be computed. |
nugget.rel.values |
set of values of the relative nugget parameter
for which the profile likelihood will be computed.
Only used if obj.likfit.glsm was created with the option fix.nugget = FALSE .
|
messages |
logical. Indicates whether messages should be printed on the screen (or output device)
while the function is running. Note that for this function
additional messages can be obtained by setting the global option
verbose=TRUE |
... |
additional parameters to be passed to the maximization
function. Typically arguments of the type control() which controls the
behavior of the minimization algorithm. For further details, see the documentation
for the minimization function optim . |
An object of the class `"proflik"' which is a list. The element contains values of the pair of parameters and the corresponding value of the profile likelihood.
Ole F. Christensen OleF.Christensen@agrsci.dk,
Paulo J. Ribeiro Jr. Paulo.Ribeiro@est.ufpr.br.
Further information about geoRglm can be found at:
http://www.daimi.au.dk/~olefc/geoRglm.
\lik{likfit.geoRglm}
for the parameter estimation, and
proflik
for the profile likelihood in the Gaussian spatial model.
data(p50) ## Not run: mcmc.5 <- mcmc.control(S.scale = 0.6, thin=20, n.iter=50000, burn.in=1000) model.5 <- list(cov.pars=c(0.6, 0.01), beta=1, family="poisson") outmcmc.5 <- glsm.mcmc(p50, model= model.5, mcmc.input = mcmc.5) mcmcobj.5 <- prepare.likfit.glsm(outmcmc.5) lik.5.sph.nugget <- likfit.glsm(mcmcobj.5, ini.phi = 1, cov.model = "spherical", nugget.rel = 0.385) pr.lik.5.sph.nugget <- proflik.glsm(mcmcobj.5, lik.5.sph.nugget, phi.values = seq(0.5,5,l=10), nugget.rel.values=seq(0.5,5,l=10)) plot(pr.lik.5.sph.nugget) ## End(Not run)