Dear Paulo Ribeiro, Sorry to send you again my previous e-mail. I am keen to understand why I am not able to use krige.glsm() with my rectangular grid.
I hope that I am not disturbing you Best regards, Hocine PhD Student UCC Cork, Ireland --- En date de : Mer 16.6.10, Lho Lho <manguefl...@yahoo.fr> a écrit : > De: Lho Lho <manguefl...@yahoo.fr> > Objet: glsm.krige in geoRglm > À: "Paulo J Ribeiro Jr" <paulo...@c3sl.ufpr.br> > Date: Mercredi 16 juin 2010, 11h22 > Dear Paulo Ribeiro > > Concerning the function krige.glsm(), I copied and past the > example code you have in geoRglm package for binomial data > (b50). > > The following code (from the package) is: > > test2 <- glsm.krige(outmcmc.5, > locations=matrix(c(0.15,0.15,0.005,0.05),2,2)) > > My question is: > > Why when I try to put a rectangular grid in the locations > argument above (instead of your matrix) > > rectangleb50.grid<-pred_grid(b50$coords[,1:2], by=0.3) > > test2 <- glsm.krige(outmcmc.5, > locations=rectangleb50.grid) > > I have these 50 warnings() from R: > > "In asympvar(kpl.result$predict, messages = FALSE) : > value of argument lag.max is not suffiently long" > > Actually, I am working with counts data and I can perform > glsm.krige only with locations=geodata$coords[,1:2]. > > If I try to put my rectangular.grid (which I built > according to my data) in the locations argument I have the > same warning messages. > > My purpose is to display a image() using the glsm.krige > output > > Please, would you like to provide me an explanation? > > Best regards > > lho > > > > > > > > > --- En date de : Lun 24.5.10, Lho Lho <manguefl...@yahoo.fr> > a écrit : > > > De: Lho Lho <manguefl...@yahoo.fr> > > Objet: Re: Help with proflik() 2 > > À: "Paulo J Ribeiro Jr" <paulo...@c3sl.ufpr.br> > > Date: Lundi 24 mai 2010, 17h18 > > Hi Paulo! > > > > Thank you for your help and sorry for the delay of > my > > answer I was away from my office. > > > > It is still not working neither (*100) nor (*1000). > > > > However, I tried in the likfit() function to choose > > trend="1st" instead of trend= "2nd" and then, the > proflik() > > function worked. > > > > I hope that the problem raised here is only due to > me. > > > > Actually, the quadratic trend surface (trend="2nd") > makes > > the data "bell" shaped according to the > plot.geodata() > > display. > > Then, the geostatistical results I could draw from it > could > > be of interest if compared with geoRglm results > > (provided I assume that my burrows counts is poisson > > distributed with units.m="m2"). > > > > I am going to try again with an other set of data (and > I > > cross my finger!!) > > > > Please, let me know if you find out a solution!!! > > > > Hocine > > > > > > > > > > > > > > > > --- En date de : Dim 23.5.10, Paulo J Ribeiro Jr > <paulo...@c3sl.ufpr.br> > > a écrit : > > > > > De: Paulo J Ribeiro Jr <paulo...@c3sl.ufpr.br> > > > Objet: Re: Help with proflik() > > > À: "Lho Lho" <manguefl...@yahoo.fr> > > > Date: Dimanche 23 mai 2010, 21h04 > > > Hi > > > > > > Yhanks for reporting and sending detailed > > information. > > > I have not yet inspected/run the files and > commands > > > but I suspect that the probelem may be due to > > numerical > > > instability > > > due to the different orders of magnitude between > your > > phi > > > (range) > > > parameters and the others. > > > > > > As a quick try you could multiply all the > coordinates > > by a > > > constant (100 > > > or 1000), changind range parameters (phi) values > > > accordingly and see if > > > the problem persists > > > > > > Please let me know > > > > > > best > > > P.J. > > > > > > > > > Em Dom, 2010-05-23 às 17:37 +0000, Lho Lho > escreveu: > > > > Dear Paulo Justiniano Ribeiro > > > > > > > > After performing the proflik() function I > had > > this > > > message from R: > > > > > > > > ##proflik:computing profile likelihood for > the > > sill > > > > Error in solve.default(xix, xiy) : > > > > system is computationally singular: > > > reciprocal condition number=5.05981e-19## > > > > > > > > I guess it is related to "matrice" with the > term > > > solve.default(??) > > > > > > > > Unfortunately, I was not able to overcome > this > > > problem. > > > > > > > > Would you like to help me? > > > > > > > > > > > > Please find attached my data in text format > > (Apple). > > > > > > > > My R code was : > > > > > > > > > > > > > > ini.values=expand.grid(x=seq(500,1300,by=25),y=seq(0,0.40,by=0.01)) > > > > inivalues<-as.matrix(ini.values) > > > > > > > > > > > > > > Burrows.likfit<-likfit(Burrows.geo,trend="2nd",fix.lambda=TRUE,lambda=1,ini.cov.pars=ini.values,fix.nugget=FALSE,nugget=0,lik.method="ML",cov.model="mat") > > > > > > > > > > > > > > profile<-proflik(Burrows.likfit,geodata=Burrows.geo,sill.values=seq(1210,1213,l=4),range.values=seq(0.026,0.029,l=4),nugget.values=seq(31,34,l=4), > > > uni.only=T)) > > > > > > > > > > > > Hocine BENCHIKH > > > > PhD student UCC > > > > Ireland > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
