Re: R and plot.variogram

[willem vervoort <w.vervoort@acss.usyd.edu.au>]

Dear Edzer and Henny
A few more clues, and found the solution. Sorry to have wasted your time 
Edzer. For all you newbies on R and myself.
I did:
> detach("package:geoR")
> plot.variogram(vgm1,model=vgm(1,"Sph",300,1))
Error: couldn't find function "plot.variogram"
# ???????????? I am lost here
> ?plot.variogram
# This brings up the help file for plot.variogram from gstat
> plot(vgm1,model=vgm(1,"Sph",300,1))
# That works (I might almost say of course, I can hear Edzer laugh in 
the background)

Sorry to have wasted your time

Willem

willem vervoort wrote:

> Henny Mehner gave the clue:
> It has to do with the mix-up between geoR and gstat.
>
> > args(plot.variogram)
> > function (x, max.dist, vario.col = "all", scaled = FALSE, var.lines = 
> > FALSE,
> >     envelope.obj = NULL, pts.range.cex, bin.cloud = FALSE, ...)
> > NULL
>
>
> Clearly this is the geoR version of plot.variogram. ?plot.variogram 
> actually tells me that! I am not an expert on R. Is there any way to 
> keep the preferred package in and unload the non-essential one?
>
> Or is there something to do about the confusion in the naming.
>
> Willem
>
> Edzer J. Pebesma wrote:
>
> > Willem, could you retry the commands under the
> > current version of R, 1.8.1?
> > --
> > Edzer
> >
> > Willem Vervoort wrote:
> >
> > > Hi all
> > > I am using R 1.8.0 under Windows XP
> > > I have installed gstat from CRAN
> > > I am having problem ploting the variograms
> > >
> > > > vgm1<-variogram(clay~x+y, loc=~x+y, PC.spdata,witdh=cutoff/50)
> > >
> > >
> > >
> > >
> > >
> > > > tmp<-fit.variogram(vgm1,vgm(1,"Sph",300,1),print.SSE=TRUE)
> > >
> > >
> > >
> > >
> > >
> > > [1] "SSErr:  31366.3780209903"
> > >
> > > All this works
> > >
> > > > plot.variogram(vgm1,model=vgm(1,"Sph",300,1))
> > >
> > >
> > >
> > >
> > >
> > > Error in if (x$output.type == "bin") Ldots$type <- "p" :
> > >       argument is of length zero
> > >
> > > In addition: Warning message:
> > > no finite arguments to max; returning -Inf
> > >
> > > > ?plot.variogram
> > >
> > >
> > >
> > >
> > >
> > > The help text is there!
> > >
> > > > plot.variogram(vgm1)
> > >
> > >
> > >
> > >
> > >
> > > Error in if (x$output.type == "bin") Ldots$type <- "p" :
> > >       argument is of length zero
> > >
> > > In addition: Warning message:
> > > no finite arguments to max; returning -Inf
> > >
> > > Not good, I seem to think it has to do with me. Am I getting the 
> > > class of my data right? The data is pretty simple:
> > >
> > > > PC.spdata[1:10,]
> > >
> > >
> > >
> > >
> > >
> > >   x   y  clay  sand fsand   EC      McN      Tikh
> > >
> > > 1   0 0.5 49.49 34.19 28.27  266 198.3690 126.33059
> > >
> > > 2   0 1.5 47.65 28.75 25.06  911 288.6718 124.15259
> > >
> > > 3   0 2.5 56.88 26.41 18.78 1251 315.1554 121.83433
> > >
> > > 4   0 3.5 61.94 25.00 16.33 1140 302.2132 119.25082
> > >
> > > 5   0 4.5 65.02 21.75 12.62 1188 273.8793 116.30404
> > >
> > > 6   0 5.5 59.79 24.61 13.97 1112 244.8895 112.92241
> > >
> > > 7   0 6.5 61.01 21.72 12.38 1046 220.6816 109.05847
> > >
> > > 8   0 7.5 43.58 52.30  8.54  725 206.2891 104.68629
> > >
> > > 9   0 8.5 51.16 47.36  7.38  713 206.7454  99.79917
> > >
> > > 10 50 0.5 48.10 33.34 27.96  253 157.0218 102.49656
> > >
> > > This is a data.frame. ?spatial.data.frame gives error
> > >
> > > Any clues would be appreciated
> > >
> > > Regards
> > > Willem Vervoort

--
R.W. Vervoort
McCaughey Lecturer Hydrology and Catchment Management
Bldg A03 
The University of Sydney, NSW 2006

ph: 	+61 2 9351 8744 (w)
	+61 2 9713 5625 (h)
fax:	+61 2 9351 5108