[R-sig-Geo] Problem about ST kriging in gstat

e.s nugget.s range.t nugget.t sill #4999.99935010.7498763 150.00347860.2262177 14.0396906 wireframe(model~spacelag+timelag,variogramSurface(ff,vv),drape =T, ,light.source = c(10,0,10),scales=list(arrows=F), zlab=expression(gamma),main="fitted variogram ") TB.ST = kr

[R-sig-Geo] St kriging

I need matrix form of st kriging. I dont know how creat st matrixes form for ordinary kriging. I need it because I want use some non-separable program function which are don't support in gstat package . I thank if someone introduce me a good reference or R cod. [[alternative HTML ver

[R-sig-Geo] St matrix form

Hi. I am working on a program to fit some non-separable spatio-temporal variogram functions. I have a problem to create variogram matrixes which are need to apply ordinary kriging. But, I dont know how create these matrixes. I need matrixes form of st variogram kriging. I thank if someone intro

[R-sig-Geo] help for optim

Subj: Help for optim Dear all, I want to optimize a function using optim. I have two problems: 1-how can i minimize function which fn value be zero not negative. Because i want optimize parameters which lead to zero fn. 2- how can i change it to maximization? Best, ___

[R-sig-Geo] Help for optim

Dear all, I want to optimize a function using optim. I have two problems: 1-how can i minimize function which fn value be positive because it can be negative? 2- how can i change it to maximization? Best, [[alternative HTML version deleted]] ___

[R-sig-Geo] Help for optim

Dear all, I want to optimize a function using optim. I have two problems: 1-how can i minimize function which fn value be positive because it can be negative? 2- how can i change it to maximization? Best, [[alternative HTML version deleted]] ___

[R-sig-Geo] help for show.vgms

Dear all, How can I add text and change color of variograms in show.vgms()? Best, Saman. -- Saman Monfared Msc, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167 ___ R-sig-Geo

[R-sig-Geo] help for lattice plot

Dear All, How can I plot some different plots in one page in lattice like "par(mfrow=c(2,2))" in graphics package? p1<-bwplot(m~PRED,data=PRED) p2<-bwplot(R~Residual,data=PRED) xyplot(obs~PRED.v|method,data=PRED.v) histogram(~PRED.v|method,data=PRED.v) Thanks, Saman. -- Sa

[R-sig-Geo] help for show.vgms

,as.groups =F) I have some problem: 1-How can I replace my desired name c("OLS","WLS","GLS","ML","REML",,.)instead of model names? 2- How can I create a legend with names: c("OLS","WLS","GLS","ML","REM

[R-sig-Geo] fit.variogram.reml

v.rml model psillrange 1 Nug -1711.462 0.00 2 Sph 3431.108 95241.99 How can I set the best range value? If exist a way to estimate range like( fit.ranges = TRUE) in fit.variogram function? Best, -- Saman Monfared Msc, Department of Statistics, Shiraz University, Shiraz 7145

[R-sig-Geo] Variogram Error

;,header=T) grd = SpatialPixels(SpatialPoints(cbind(grd$x,grd$y))) plot(grd) n =3 library(xts) tgrd = seq(max(index(xx)+31622400*2),max(index(xx)+3*31622400),length=2) pred.grd = STF(grd, tgrd) plot(pred.grd ) cancer.ST = krigeST(Rate~1,xx,pred.grd ,separableModel) stplot(cancer.ST) -- Saman Monfa

[R-sig-Geo] All st-gamma values are 0

;,header=T) grd = SpatialPixels(SpatialPoints(cbind(grd$x,grd$y))) plot(grd) n =3 library(xts) tgrd = seq(max(index(xx)+31622400*2),max(index(xx)+3*31622400),length=2) pred.grd = STF(grd, tgrd) plot(pred.grd ) cancer.ST = krigeST(Rate~1,xx,pred.grd ,separableModel) stplot(cancer.ST) -- Saman

[R-sig-Geo] STSDF construct

Hi, I try to constructa stsdf st class. But I dont know how. I attached my data and cod. my data are cancer rate incidence. Best, Saman. -- Saman Monfared Msc, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167

[R-sig-Geo] Variogram Paradox

k2<-krige.cv(log(zinc)~1,meuse,v.ff) mean(k1$residual) mean(k2$residual) mean(k1$residual^2) mean(k2$residual^2) cor(k1$var1.pred,k1$observed) cor(k2$var1.pred,k2$observed) cor(k1$var1.pred,k1$residual) cor(k2$var1.pred,k2$residual) Best, Saman. -- Saman Monfared Msc, Department of Statistics,

[R-sig-Geo] error in STSDF class

",1:length(values), sep = "_") mydata = data.frame(values = signif(values, 3), ID=IDs) library(spacetime) stfdf = STSDF(sp,time, data = mydata) order(time, t) == t is not all TRUE Thanks.-- Saman Monfared Msc, Department of Statistics, Shiraz University,

[R-sig-Geo] NA in ST-variogram

pars.M$par vv$model = ExpVgmMetric(pars.M$par, vv$spacelag, vv$timelag, nugget = 0) plot(vv) but it's sill is negative and in plot part model only has a single value and dose not change by changing spatial and temporal lags. Best, Saman. -- Saman Monfared Msc, Department of Stat

[R-sig-Geo] spatio-temporal data

rstand how they constructed. for example can I input my data as a below then construct a st class?how? locations x ytime point1cord1 cord2 time1 point1cord1 cord2 time2 .. . . .. . . Best. -- Saman Monfared Msc, Dep

[R-sig-Geo] help for spdep

? I couldn't find a good refrences to understand tese model!!! How can I predict this model for new data?? How can I cross validate results of CAR and SAR?? -- Saman Monfared Msc, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 53051

[R-sig-Geo] Help for spdep package

patulm for new data?? Thanks. -- Saman Monfared Msc, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167 ___ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/list

[R-sig-Geo] Help for spdep package

? I couldn't find a good refrences to understand tese model!!! How can I predict this model for new data?? How can I cross validate results of CAR and SAR?? -- Saman Monfared Msc, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 53051

[R-sig-Geo] variogram models in gstat

Dear all, How can I see mathematical formulation of variogram models in gstat? Some of them are in gstat manual but there aren't others like ("Exc,"Ste","Mat","Hol"). Thanks. -- Saman Monfared Msc, Department of Statistics, Shiraz University, Shiraz 714

[R-sig-Geo] formulas of variograms in gstat

Dear all, How can I see mathematical formulation of variogram models in gstat? Some of them are in gstat manual but there aren't others like ("Exc,"Ste","Mat","Hol"). Thanks. -- Saman Monfared Msc, Department of Statistics, Shiraz University, Shiraz 714

[R-sig-Geo] spatio-temporal kriging

atio-temporal grid for prediction in new space time locations? 2.How can I perform krigimg methods? indeed I don't know how calculate weight and covariance matrix for prediction. thanks for your help. -- Saman Monfared Msc, Department of Statistics, Shiraz University,

[R-sig-Geo] help for fit.vaiogram.gls

del (**) is my answer?? I have changed the start values several times but The warning remain. Thanks. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167

[R-sig-Geo] calculate SSerror in fit.variogram

Hi. I have used fit.variogram in gstat. How can I calculate SSerror between the sample variogram and the fitted variogram model in different estimation method of fit.method=2-7. Thanks. Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar

[R-sig-Geo] transform coordinates

Hi. How can I transform a coordinate system like (N26 5 21, E39 4 76) to utm system?? THanks. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167 ___ R

[R-sig-Geo] nlm with condition

Dear all. I try to minimize a function using nlm. my function has three parameter that all of them are positive. How can I do this condition in parameter estimation with nlm. Thanks. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email

[R-sig-Geo] optimization with nlm

n I optimize parameters of a model? How can I use nlm? I read an example in help of nlm but I couldnt understand it. Help me. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305

[R-sig-Geo] Fwd: help in CompRandFld package

-- Forwarded message -- From: Saman Monfared Date: Tue, 1 May 2012 12:31:30 +0430 Subject: help in CompRandFld package To: r-sig-geo Dear All. I Try to use CompRandFld package for Spatio-temporal modeling of Irish wind data. I write Irish wind data in txt format. I calculate

[R-sig-Geo] help in CompRandFld package

the leading minor of order 11 is not positive definite. Please help me. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167 ___ R-sig-Geo mailing list R-sig

[R-sig-Geo] spatio-temporal

e_t=1,sill=sill))$data the error is : Error: could not find function "RFsim" please guide me about this. Thanks. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167 _

[R-sig-Geo] help for emperical spatio-temporal

Dear All. I need help for calculate empirical spatio-temporal variogram or covariance. for example Irish wind data. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167

[R-sig-Geo] help to analisys spatio-temporal

Hi. Are there packages in R that can perform spatio-temporal analysis such as Testing stationary and sepreability, fitting covariance function(Sepreable and non-sepreable) and kriging?? Please guid me. Thanks -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz

[R-sig-Geo] spatio-temporal

Hi. Is there packages in R that can perform spatio-temporal analysis such as Testing stationary and sepreability, fitting covariance function(Sepreable and non-sepreable) and kriging?? Please guid me. Thanks -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454

[R-sig-Geo] extract value of polygone

Hi. I try to extract values of temperature in a topography map. It is like a contour plot. I need the coordinates and temperature of some points. how can I perform it?? thanks. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar

[R-sig-Geo] extract value of polygone

Hi. I try to extract values of temperature in a topography map. It is like a contour plot. I need the coordinates and temperature of some points. how can I perform it?? thanks. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar

[R-sig-Geo] difrrence between regression kriging&residual kriging

-- Forwarded message -- From: Saman Monfared Date: Sat, 14 Jan 2012 18:25:51 +0330 Subject: difrrence between regression kriging&residual kriging To: r-sig-geo Hi I'm very confused about the concepts of regression kriging ,kriging with external drift (KED) and

[R-sig-Geo] difrrence between regression kriging&residual kriging

iging??for example(residual cokriging-regression cokriging,...) Which packages in R can perform them?? -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167 ___

[R-sig-Geo] how can I add back transformation for ploting

or transformed data how can I add back the transformation for ploting by image(pp,1) and countour(pp,1,add=T)?? Thanks. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167 _

[R-sig-Geo] problem in cokriging plot

Hi. I have some problems in prediction of cokriging 1)- I don't know how make a grid for prediction, I have a polygone and I select 3000 points into this polygone. I want make grid like "meuse.grid" nc.sids <- readShapePoly("City Boundary.shp") pts <- spsample(nc.sids, n=3000, type="regular") grd

[R-sig-Geo] fit variogram base on cross validation

0 E[2] Lin 123.96281576 75017 T[1] Nug 4.41847897 0 T[2] Lin 12.22390243 75017 B.E[1] Nug -1.26696750 0 B.E[2] Lin 1.26071039 75017 B.T[1] Nug -1.25868233 0 B.T[2] Lin -2.4909 75017 E.T[1]Nug 2.21925267 0 E.T[2]Lin -23.10121223 75017

[R-sig-Geo] Arbitrary cross variogram

an R but I don't know how fit this cross variogram in R. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email: samanmonfar...@gmail.com Tel: +98 917 5305167 [[alternative HTML version deleted]] __

[R-sig-Geo] cross variogram

t;T", form = T~ 1,data=dd) vm <- variogram(g) plot(vm) vm.fit <- fit.lmc(vm, g, model=m,fit.ranges =T) plot(vm, vm.fit,col="black",type="b",lwd=3,main="Fig1") cv.c2 <- gstat.cv(g) summary(cv.c2$residual) sqrt(mean(cv.c2$residual^2)) mean(cv.c2$residual

[R-sig-Geo] How can i chose number of lags

Dear All. How can I chose number of lags in variogram function (gstat) for example I want calculate empirical variogram in 12 lags.it is 15 by default in gstat package. ٌWhat should I do?? -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran Email

[R-sig-Geo] Cross variogram fitting

ot;black") but when I use predict.gstat(g.fit, dd) I faced with error ""non-positive definite coefficient matrix in structure 1non-positive definite coefficient matrix in structure 2Warning: No Intrinsic Correlation or Linear Model of Coregionalization found"" -- Sa

[R-sig-Geo] Please help me

Dear All. Hi. I want perform a cokriging for three co-variable with gstat package but I don't know how fit the best cross-variogram model that have a minimum error. Please guide me. Thanks. -- Saman Monfared Msc Student, Department of Statistics, Shiraz University, Shiraz 71454, Iran

Re: [R-sig-Geo] Getting vector index (i.e. layer number) for each cell in raster brick

;> >> Although this *appears* to work, it seems to be unnecessarily >> complicated for something that I am sure is done in a only a few >> lines. >> >> I will greatly appreciate any pointers for how to go about solving >> this problem. >> >> Many th