Dear Rusers,
  I'm using geoRglm package to analyze my data, but an error occured when i
do the krige predictions. The information on the packages are:

R version 2.4.1 beta (2006-12-05 r40126)
i386-pc-mingw32

locale:
LC_COLLATE=Chinese_People's Republic of China.936;LC_CTYPE=Chinese_People's
Republic of China.936;LC_MONETARY=Chinese_People's Republic of
China.936;LC_NUMERIC=C;LC_TIME=Chinese_People's
Republic of China.936

attached base packages:
[1] "stats"     "graphics"  "grDevices" "utils"     "datasets"  "methods"
"base"

other attached packages:
     coda   lattice     rgdal   geoRglm      geoR        sp
 "0.10-7" "0.14-16"   "0.5-1"  "0.8-18"  "1.6-11"   "0.9-4"
 My programs are as follows:

*model<-read.csv('d:/deleting/models2.csv',sep=',',header=T)
model[1:3,]
*          lat          long      x
1 30.61688 117.4140    49
2 30.61686 117.4140    79
3 30.61680 117.4139    43
*coordinates(model) <- c("long", "lat")
proj4string(model) <- CRS("+proj=longlat +ellps=WGS84")
model.proj<-spTransform(model, CRS(("+init=epsg:2384")))
model.proj[1:3,]
        *coordinates        x
1 (539697, 3388570)    49
2 (539694, 3388570)    79
3 (539690, 3388570)    43
*model.proj<-as.geodata(model.proj)*

*model<-list(cov.pars=c(1,1),beta=1,family="poisson")
mcmc2.test<-mcmc.control(S.scale=0.2,thin=1)
test2.tune<-glsm.mcmc(model.proj,model=model,mcmc.input=mcmc2.test)  #make
theAcc. rate=0.6*

*mcmc2.test<-mcmc.control(S.scale=0.44,thin=1)
test2.tune<-glsm.mcmc(model.proj,model=model,mcmc.input=mcmc2.test)
#Acc.-rate =  0.605 *

*test2.tune.c<-create.mcmc.coda(test2.tune,mcmc.input=mcmc2.test)
plot(test2.tune.c,density=FALSE,ask=FALSE,auto.layout=FALSE)  #trace plot
autocorr.plot(test2.tune.c,ask=FALSE,auto.layout=FALSE)  #autocorrelation
plot*

*mcmc2<-mcmc.control(S.scale=0.44,thin=100)  #thin=100,in order to reduce
the autocorrelation
test2<-glsm.mcmc(model.proj,model=model,mcmc.input=mcmc2)*

*out2<-output.glm.control(sim.predict=TRUE)
pred.test2<-glsm.krige(test2,locations=model.proj$coordinates,output=out2)*
*#errors
*glsm.krige: Prediction for a generalised linear spatial model
*erro in locations[, 1] : measurement number is wrong*

Where's the problem?
BTW, how should i split the coordinates into two variables,like:
*    *coord1     coord2       x
1  539697    3388570    49
2  539694   3388570    79
3  539690   3388570    43

Any suggestions or help are greatly welcome.


-- 
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]]


