[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.



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
tel. + 32 54/436 185
Thierry.Onkelinx at 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
~ 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

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



	[[alternative HTML version deleted]]

R-sig-Geo mailing list
R-sig-Geo at stat.math.ethz.ch

More information about the R-sig-Geo mailing list