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   epistat@gmail.com
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]]


