search for: vgm

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080731/d5bf0e31/attachment.pl>

...<- sqlFetch(channel, "All") PAHdata attach(PAHdata) logTotalPAH = log10(TotalPAH) PAHdata = cbind(PAHdata, logTotalPAH) class(PAHdata) names(PAHdata) coordinates(PAHdata) = ~Easting + Northing + SampleElevation class(PAHdata) summary(PAHdata) dimensions(PAHdata) hist(logTotalPAH) lpah.vgm = variogram(logTotalPAH ~1, PAHdata) lpah.vgm plot(lpah.vgm ) lpah.fit = fit.variogram(lpah.vgm, model = vgm( 6, "Sph", 80, 1.5)) lpah.fit plot(lpah.vgm, lpah.fit) x = GridTopology(c(534500,3531400, 20), c(2,2, 0.5), c(176,126, 31)) class(x) coordinatevalues(x) x = SpatialGrid(x) class(x)...

...y.dat)=~X+Y pred.dat<-read.table("C:/temp/Predict_Location_XY.txt", header = FALSE) names(pred.dat)<-c("x","y","p") my.pred.dat<-data.frame(pred.dat) coordinates(my.pred.dat)<-~x+y gridded(my.pred.dat) = TRUE #Set up 2 exponential variograms lzm.vgm.exp.noNug<-vgm(1.1,"Exp",2000,0) lzm.vgm.exp.Nug<-vgm(0.7,"Exp",3000,0.4) #Krige the exponential variograms lzm.krig.exp.noNug<-krige(Kh_Log10~1,my.dat,my.pred.dat,model=lzm.vgm.exp.noNug) lzm.krig.exp.Nug<-krige(Kh_Log10~1,my.dat,my.pred.dat,model=lzm.vgm.exp.Nug)...

Dear Members, I am trying to fit a variogram model using fit.variogram function from the gstat package. The figure showing my experimental variogram can be seen here: https://i.stack.imgur.com/UZXw4.png My code line for this operation is: > c2.vgm.fit<-fit.variogram(c2.vgm.exp,vgm(nugget=0, psill=400,model="Exp",range =40000),fit.method = 7) The system throws following error message: Error in switch(model, exponential = fit.exponential(v.object, c0 = nugget, : EXPR must be a length 1 vector Any help to know the cause of t...

...es: "solve.c", line 88: singular matrix in function Usolve() "lufactor.c", line 208: singular matrix in function m_inverse() Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, : m_inverse In addition: Warning message: In fit.variogram(vg.aux, model = vgm(psill = 0.1, model = "Gau", : Warning: singular model in variogram fit Mi función del variograma es : v.fit <- fit.variogram(vg.aux, model=vgm(psill=0.15, model='Gau', range=5000, nugget=0.05)) ¿Alguien me podría ayudar? Gra...

Dear list, I am having troubles increasing the fontize when plotting a variogram{gstat} and its model (vgm) with plot and using jpeg(). Also the resolution in the jpeg call does not work. I am using R2.4 on Mandriva 10.2 linux. I can change fontsize with cex.axis in a normal plot, so I presume it has to do with plotting the variogram model. Any help on how to increase the font size and resolution would...

Hi, Is there any way that I can supress error messages so that they don't stop for loops running? I'm using the gstat package and have created a variogram model for the double exponential model: "dexp <- vgm( .5, "Exp", 10,.1, add.to = vgm(.5, "Exp", 100))" and fitted it using "dexp.fit <- fit.variogram(temp2, dexp, fit.method=6, debug.level=0)" Occasionally I get the following error message and the model doesn't fit "Error in if (direct[dire...

...t I get incomplete line plots when I plot the output of variog4, compared to the plot for the output of variog. I'm using the same geodata for each plot. I've attached 2 graphs to illustrate. The first is a plot of a 90- degree directional variogram created using variog: > trees90.vgm<-variog(trees.geo,max.dist=100,direction=pi/2) > plot(trees90.vgm) The second is a plot of a multidirectional variogram created using variog4: > trees.vgm4<-variog4(trees.geo,max.dist=100) > plot(trees.vgm4) The graphs are also available at: www.tiem.utk.edu/~mmfuller/R_plots....

...tion Usolve() > > > > "lufactor.c", line 208: singular matrix in function m_inverse() > > Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, : > > m_inverse > > In addition: Warning message: > > In fit.variogram(vg.aux, model = vgm(psill = 0.1, model = "Gau", : > > Warning: singular model in variogram fit > > > > > > Mi funci?n del variograma es : > > v.fit <- fit.variogram(vg.aux, model=vgm(psill=0.15, model='Gau', range=5000, > >...

...proj4string(grid)=CRS("+init=epsg:28992") proj4string(data_c)=CRS("+init=epsg:28992") g <- gstat(id = "Cu", formula =Cu~1,data = data_c) g <- gstat(g, id = "S", formula = S~1,data =data_c) g <- gstat(g, id = c("Cu", "S"), model = vgm(cov(data_c$Cu,data_c$S) , "Sph", 200, 0)) g<fit.lmc(variogram(g),g,model=vgm(psill=cov(data_c$Cu,data_c$S),model="Sph",range=200,nugget=0)) k.c <- predict.gstat(g, grid,block = c(30,30,15)) Thanks and Regards, Nishani Musafer

...h xyplot in lattice; you'll find where it gets it's defaults by library(lattice) trellis.par.get() ?trellis.par.set -- Edzer Alexander.Herr at csiro.au wrote: > > Dear list, > > I am having troubles increasing the fontize when plotting a > variogram{gstat} and its model (vgm) with plot and using jpeg(). Also > the resolution in the jpeg call does not work. I am using R2.4 on > Mandriva 10.2 linux. > > I can change fontsize with cex.axis in a normal plot, so I presume it > has to do with plotting the variogram model. Any help on how to > increase t...

...times that the passwords in the Icecast config file & oddcast config box are both the same! [2002-10-20 17:21:30] INFO main/D:\Ed\xiph\icecast\src\main.c icecast server started [2002-10-20 17:23:26] INFO connection/D:\Ed\xiph\icecast\src\connection.c Source logging in at mountpoint "/vgm.ogg" **** [2002-10-20 17:23:26] INFO connection/D:\Ed\xiph\icecast\src\connection.c Source (/vgm.ogg) attempted to login with invalid or missing password **** Any help? - Moguta www.ffradio.net >===== Original Message From "Dave St John" <dstjohn@mediacast1.com> =====...

...e anybody has the experience in using variogram(gstat) ? Please kindly give me some hints about the results. I used variogram() to build a semivariogram plot as: tr.var=variogram(Incr~1,loc=~X+Y,data=TRI2TU,width=5) then fir the variogram to get the parameters as: v.fit = fit.variogram(tr.var,vgm(0.5,"Exp",300,1)) v.fit model psill range 1 Nug 1.484879 0.00000 2 Exp 3.476700 29.70914 This is the output of v.fit. Can anybody help me write the exponential formula for this variogram? I have the problem in understanding the result. Thanks! Leaf

Hi; I cannot find in the R html documentation a way to evaluate cross-correlation in 2D data sets. I would like to evaluate cross-correlation in a series of moving windows between two maps. i,e, specify several windows inside the complete 2D spatial matrixes and for each one ofthese windows evaluate the 2D cross-correlation (commonly conducted in the spectral domain). Thanks in advance and best

...or radius 45? minus radius) 8. Fig 8- semi-variogram with anisotropy (0?, 45?, 90?, 135?). Best shape is from 135? 9. Fig 9- semi-variogram fit with Gaussian Model. R code is: > v = variogram(Z~X+Y, subground, cutoff=1800, width=200, alpha=c(135)) > v.fit = fit.variogram(v, vgm(psill = 1, model="Gau", range=1800, nugget= 0, anis=c(135, 0.5))) R code: testground2 <- read.table(file="c:/work_LIDAR_USA/R_kriging_new_set/ground_26841492694149_x yz.txt", header=T) class (testground2) coordinates (testground2)=~X+Y # this makes testground a Spat...

Hi, please help!! I've been trying to fit variogram models using the REML method in the gstat package. Every time the Windows GUI crashes. For example library(gstat) data(meuse) x <- variogram(zinc ~ 1, ~x + y, meuse) v <- vgm(140000, "Sph", 800, nug = 10000) plot(x, model = fit.variogram(x, model = v, fit.method=5)) Other fit methods are non problematic (eg. fit.method=7 for WLS or fit.method=1 for OLS) I've tried the code on Windows XP, 98 and NT and all fail, how about non-Windows platforms? Any sugg...

Hello everybody, I want to run the following R-code, but it does'nt work. > res.vgm = variogram(ertrag ~ rep, # !!! + loc = ~ xpos + ypos, # !!! + width = start.range / 10, + data = d) Error in vector("double", length) : element 2 is empty; the part of the args list of '.Internal' being evaluated was:...

...e.offset = cc,cellsize = cs, cells.dim = cd) getClass("SpatialGrid") p4s <- CRS(proj4string(rsa2)) x2_SG <- SpatialGrid(rsa2_grd, proj4string = p4s) x2_SP <- SpatialPoints(cbind(x2$X, x2$Y)) v <- variogram(log1p(tsport_ace) ~ 1, x2, cutoff=100, width=9) te<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0)) y <- krige(tsport_ace~1, x2, x2_SG, model = ve.fit) spplot(y, 1, col.regions = bpy.colors(100), sp.layout = list("sp.lines",as(rsa2, "SpatialLines"),no.clip = TRUE)) I'm looking forward to your response. Thanks. Best regards, Pearl del...

...validation. Right now I am having difficulty using the cross validation of gstat. Below are my commands with the tsport_ace as the variable: nfold <- 3 part <- sample(1:nfold, 69, replace = TRUE) sel <- (part != 1) m.model <- x2[sel, ] m.valid <- x2[-sel, ] t<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0)) cv69 <- krige.cv(tsport_ace ~ 1, x2, t, nfold = nrow(x2)) The last line gives an error saying: Error in SpatialPointsDataFrame(coordinates(data), data.frame(matrix(as.numeric(NA), : row.names of data and coords do not match I don't know what is wrong....

Hello, I'm using the gstat package within R for an automated procedure that uses ordinary kriging. I can see that there is a logical ("singular") atrtibute of some adjusted model semivariograms: .- attr(*, "singular")= logi TRUE I cannot find documentation about the exact meaning and the implications of this attribute, and I dont know anything about the inner calculations