Dear friends,
I have learned to use geoRglm package for some time, but met an urgent
problem which can't be solved by me. I think it should be some wrong
understanding on
the arguments in geoRglm package. Could anybody help me to correct the
following geoRglm programs?
My dataset includes 1 responsible and 6 explanatory varaibles, the
following are my programs:
a<-read.csv('d:/deleting/models.csv',header=T,sep=',')
geodata<-as.geodata(a,coords.col=c(9,8),data.col=11,covar.col=c(2,4:7,10),
na.action="none",borders=T) #data.col define the responsbile
var,covar.colspecify 6 explanatory varaibles
model2<-list(trend=trend.spatial
(~covariate[,1]+covariate[,2]+covariate[,3]+covariate[,4]+covariate[,5]+covariate[,6]),
cov.pars=c(1,1,1,1,1,1),beta=58,family="poisson") #covariate[,*]specify 6
explanatory varaibles
mcmc2.test2<-mcmc.control(S.scale=0.5,thin=1)
test2.tune<-glsm.mcmc(geodata=geodata,model=model2,mcmc.input=mcmc2.test2)
Error is: glsm.mcmc(geodata = geodata, model = model2, mcmc.input =
mcmc2.test2) :
* size of beta is incompatible with trend specified
*Maybe there are wrong specify both in cov.pars=c(1,1,1,1,1,1) and beta=58
(model2).
I have try to use it to analyze my data for a long time, but failed.
Any comments or help are greatly appreciate.
--
With Kind Regards,
oooO:::::::::
(..):::::::::
:\.(:::Oooo::
::\_)::(..)::
:::::::)./:::
::::::(_/::::
:::::::::::::
[***********************************************************************]
Zhi Jie,Zhang ,PHD
Tel:86-21-54237149 [EMAIL PROTECTED]
Dept. of Epidemiology,school of public health,Fudan University
Address:No. 138 Yi Xue Yuan Road,Shanghai,China
Postcode:200032
[***********************************************************************]
oooO:::::::::
(..):::::::::
:\.(:::Oooo::
::\_)::(..)::
:::::::)./:::
::::::(_/::::
:::::::::::::
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
