Displaying 1 result from an estimated 1 matches for "sampleelevation".
2008 Aug 07
0
3d kriging et al
...)
options(digits = 10)
channel = odbcDriverConnect("")
sqlTables(channel)
PAHdata <- sqlFetch(channel, "All")
PAHdata
attach(PAHdata)
logTotalPAH = log10(TotalPAH)
PAHdata = cbind(PAHdata, logTotalPAH)
class(PAHdata)
names(PAHdata)
coordinates(PAHdata) = ~Easting + Northing + SampleElevation
class(PAHdata)
summary(PAHdata)
dimensions(PAHdata)
hist(logTotalPAH)
lpah.vgm = variogram(logTotalPAH ~1, PAHdata)
lpah.vgm
plot(lpah.vgm )
lpah.fit = fit.variogram(lpah.vgm, model = vgm( 6, "Sph", 80, 1.5))
lpah.fit
plot(lpah.vgm, lpah.fit)
x = GridTopology(c(534500,3531400, 20), c(2,2,...