[R-sig-Geo] variogram using gstat

[ONKELINX, Thierry]

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