Displaying 20 results from an estimated 900 matches similar to: "gstat package. "singular" attibute"
2006 Aug 20
3
plot problem
Hello.
I'm pretty much new to R and I'm trying to produce some figures.
It seems to me, that R has some asynchronous way of plotting figures.
When I run this code:
#constructs the semivariogram of SC1929
vgm1 <- variogram(SC1929~1,~U+V,puerto.map$att.data)
# trying to make new plot
dev.set(which=dev.next())
plot(vgm1)
title(main="Semivariogram",font.main=4)
2009 Mar 24
1
Variogram with Gstat
Dear all can you help me? i have this problem: i have a dataset in a text file
in a matrix of 3 columns: x, y, z where x and y are the coordinates and Z are
the mesurements. How can i do a variogram with R?
2009 Sep 08
1
gstat---2 basic plot questions
Hi all--
I'm new to R, statistics and programming, so sorry if this is a really basic
question!
I have plotted a directional variogram, and I want to
a. overlay the omni-directional line over each directional panel
b. display the directional variograms in a single panel with a legend that
associated each line to each degree measurement.
The line I'm using is
2003 Mar 07
1
REML option in gstat
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
2009 Jun 10
2
plot two variograms on a same graph
Hi,
I would know how to plot two variograms on a same graph. I can plot one by one but I would draw both on the same one.
Is it possible? Do i need any special package?
Thanks!
Cordialement
Damien Landais
2005 Nov 08
2
Variogram
Dear All,
Is there 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
2008 Apr 01
1
spatial cross-correlation
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
2008 Aug 05
4
LIDAR Problem in R (THANKS for HELP)
Hi All,
I am a PhD student in forestry science and I am working with LiDAR data set
(huge data set). I am a brand-new in R and geostatistic (SORRY, my
background it?s in forestry) but I wish improve my skill for improve myself.
I wish to develop a methodology to processing a large data-set of points
(typical in LiDAR) but there is a problem with memory. I had created a
subsample data-base but
2006 Apr 25
3
persp plot increasing 'x' and 'y' values expected
hello,
i do the following in order to get an persp-plot
x<-c(2,2,2,2,2,2,3,3,3,3)
y<-c(41,41,83,83,124,166,208,208,208,208)
z<-c(90366,90366,92240,92240,92240,96473,100995,100995,100995,100995)
x<-data$x
y<-data$y
z<-matrix(data$z,length(y),length(x))
persp(x,y,z, col="gray")
but i always get the error message increasing 'x' and 'y' values
expected,
2008 Apr 29
2
Variogram problem
Hello,
I'm french and I have some difficulties in carry out the semiones under R
with an aim to carry out an interpolation by krigeage. My goal is to obtain
a chart of the distribution of precipitations/temperatures in Europe
starting from 73 different stations (and, of course, distributed irregularly
on the chart, where use of the krigeage). Here, I carried out this to test
to obtain
2007 Oct 08
1
semivariogram
I need the most straightforward way to build semivariograms within R
i am currently using s-gems software but i would like to even test R
Thanks
D
[[alternative HTML version deleted]]
2006 Nov 24
2
plot() and Jpeg() increase font size and resolution
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
2010 Nov 17
1
Please, help me with 'mattern' variogram
Hi, R-folks:
I have been tryin many combination of parameter to make Matern variogram to work, but I can't find the available one. I'm near to be crazy.
I tiped:
A?o2003Selg.lf<-likfit(A?o2003Selg,cov.model="matern",ini.cov.pars=c(1.5,14),kappa=2.5,fix.kappa=FALSE,nugget=0.08,lambda=0.008,fix.lambda=FALSE,hessian=TRUE)
the hessian shows:
$hessian
[,1]
2012 Jul 31
2
Spatial regression
Dear list members,
I have a large dataset with point data. Each data point has X, Y, number of
species present and values for one single explanatory variable. Each point
was generated from a polygon grid shapefile and they are at equal
distances.
I conducted GLS analyses introducing variance structures and correlation
structures to improve the model.
I plotted variograms before and after adding
2008 Jul 31
1
anisotropy in vgm model. HELP!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
? 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>
2005 Jun 10
1
In connection with the creation of a grid of point with gstat
Hi ALL,
With the pakage gstat
How can one create a grid of points as in the meuse.grid example, from
measuring sites as in the Meuse example
Thanks in advance!
Alain PONSERO
Conservateur
Réserve Naturelle de la Baie de Saint-Brieuc
site de l'Etoile
22120 HILLION
tel/fax : 02.96.32.31.40
<mailto:aponsero@cabri22.com> aponsero@cabri22.com
2004 Jun 16
0
gstat 0.9-12: cokriging cross validation and class name incompatibilities
I uploaded gstat 0.9-12 to CRAN, which has a few important changes:
1. Cokriging cross validation
Cokriging cross validation is now possible with the function gstat.cv:
you simply pass a multivariable gstat object, and cross validation is done
for the first variable in the object. Optionally, secondary variable
records
at locations coinciding with the validation locations are removed.
2. Class
2004 Jun 16
0
gstat 0.9-12: cokriging cross validation and class name incompatibilities
I uploaded gstat 0.9-12 to CRAN, which has a few important changes:
1. Cokriging cross validation
Cokriging cross validation is now possible with the function gstat.cv:
you simply pass a multivariable gstat object, and cross validation is done
for the first variable in the object. Optionally, secondary variable
records
at locations coinciding with the validation locations are removed.
2. Class
2006 Jun 03
1
default value for cutoff in gstat variogram()
I wonder what is the default value for the argument 'cutoff' when not
specified in the variogram.formula function of gstat. Computing
variogram envelops within gstat, I am comparing the results obtained
with variog in geoR and variogram in gstat, and it took me a while
before understanding that the cutoff default value is not the maximum
distance.
Can Edzer tell us about it?
All the
2017 Sep 24
1
R version 3.3.2, Windows 10: gstat package: Error in fitting a variogram model using 'fit.variogram' function
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