[R-sig-Geo] variogram using gstat
ONKELINX, Thierry
Thierry.ONKELINX at inbo.be
Wed May 7 13:06:58 CEST 2008
Dear Pujan,
Calculating variograms on large datasets is computationally rather
heavy. So maybe R didn't hang but was still busy calculating the
variogram. Run your script overnight (or over the weekend) and see the
next morning if it finished.
I couldn't find an error in your code at first glance.
HTH,
Thierry
------------------------------------------------------------------------
----
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel biometrie, methodologie en kwaliteitszorg / Section biometrics,
methodology and quality assurance
Gaverstraat 4
9500 Geraardsbergen
Belgium
tel. + 32 54/436 185
Thierry.Onkelinx at inbo.be
www.inbo.be
To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to
say what the experiment died of.
~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data.
~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of
data.
~ John Tukey
-----Oorspronkelijk bericht-----
Van: r-sig-geo-bounces at stat.math.ethz.ch
[mailto:r-sig-geo-bounces at stat.math.ethz.ch] Namens PUJAN RAJ REGMI
Verzonden: woensdag 7 mei 2008 12:22
Aan: r-sig-geo at stat.math.ethz.ch
Onderwerp: [R-sig-Geo] variogram using gstat
Dear list,
I am trying to plot semivariogram suing gstat packages. But as soon as i
run the command variogram the programs hangs!! I don't know what
happened! Actually in my data set i have many rows (96000). My data has
5 columns out of these first 2 are x and y coordinates and the rest 3
are the sample data of different items under study. following is the
command i used to run the program:
A<-read.table(file.choose(new=FALSE), sep = ",",skip=1,header=FALSE)
names(A) <- c("x", "y", "R", "G", "B")
coordinates(A) = ~x+y
variogram(R~1,A,cutoff=800,width=800/15)
My data set has maximum spatial extent of 2000m in x and and around
1600m in y direction. The variation in x and y direction is 12 m ( the
data set is irregular but in grid(12m))
Any help will be greatly appreciated
Thanking you
Pujan
_________________________________________________________________
1.
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
More information about the R-sig-Geo
mailing list