Displaying 1 result from an estimated 1 matches for "outmcmc".
Did you mean:
out_cmd
2012 Oct 04
1
geoRglm with factor variable as covariable
...re covar.col is a factor variable (years in this case 90-91-92-93)).
Then I run the model as follow:
/
model.5 = list(cov.pars=c(1,1), cov.model='exponential', beta=1,
family="poisson")
mcmc.5 = mcmc.control(S.scale = 0.25, n.iter = 30000, burn.in=50000, thin =
100) #trial error
outmcmc.5 = glsm.mcmc(geoData9093, model= model.5, mcmc.input = mcmc.5)
mcmcobj.5 = prepare.likfit.glsm(outmcmc.5)
lik.5 = likfit.glsm(mcmcobj.5, ini.phi = 0.3, fix.nugget.rel = F)/
And the summary of lik.5 is:
likfit.glsm: estimated model parameters:
beta sigmasq phi tausq.rel
&qu...