[R-sig-Geo] Kriging

[Thomas.Adams at noaa.gov]

Maybe something like this:

http://www.ars.usda.gov/sp2UserFiles/ad_hoc/12000000SpatialWorkshop/01VinyardOverview.pdf



----- Original Message -----
From: milton ruser <milton.ruser at gmail.com>
Date: Tuesday, August 25, 2009 7:50 pm
Subject: Re: [R-sig-Geo] Kriging

> Hi Kabeli,
> 
> I never saw Brian Vinyard slides 40! :-)
> It is accessible on a internet site?
> 
> bests
> milton
> 
> On Tue, Aug 25, 2009 at 6:56 PM, KABELI MEFANE 
> <kabelimefane at yahoo.co.uk>wrote:
> >   Hi all
> >
> > Please help and correct me, to predict Z(S0) at (0.5,0.5) given 
> Z(S1) = 3
> > at (0,0), Z(S2) = 5 at(0,1), Z(S3) = 6(1,0) and Z(S4) = 4(1,1). Let
> >
> > ã (h)=h^2  h<1
> >
> >        = 1     h>=1 using geo.
> >
> >
> >
> > I did this: coords<-matrix(c(0,0,1,1,0,1,0,1), nrow=4, ncol=2)
> > > data<-c(3,5,6,4)
> > > coordata<-data.frame(coords=coords,data=data)
> > > geo<-as.geodata(coordata)
> > > plot(geo$coords, xlim = c(0, 1.2), ylim = c(0,1.2),
> > +     xlab = "X-axis", ylab = "Y-axis")
> > > locations<- matrix(c(0.5, 0.5), ncol = 2)
> > > text(locations, as.character(0), col = "red")
> > > polygon(x = c(0, 1, 1, 0), y = c(0, 0, 1, 1),lty = 2)
> > > szero<- krige.conv(geo, locations=locations,
> > +         krige=krige.control(cov.pars=c(1,1)))
> >
> >
> >
> > how do i show this points, lines and values on the gragh like it 
> was done
> > an erview of spatial statistics by Brian Vinyard slides 40 or 
> something> similar?.
> >
> >
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > R-sig-Geo at stat.math.ethz.ch
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
> >
> 
> 	[[alternative HTML version deleted]]
> 
> 
-------------- next part --------------
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo