matthewfindleych2mcom r users working sp gstat packages past couple days effort analyze set 200 soil samples collected various eastings northings depths containing wide range measured hydrocarbon concentrations thus far ive managed import data logtransform concentrations assign coordinates generate fit variogram model krige data 3d grid ive come hard limits r proficiency im looking helpful suggestions visualize 3d grid expected values subset data 2d slices grid xy xz planes relatively simple grid math operations removing coordinates higher elevations ground surface ive made 2d grid ground surface values reverse logtransformation ultimately estimate total mass hydrocarbons soil heres ive gathered supporting documentation need cast 3d expected concentration values sps spatialpixelsdataframe class order make irregular subsets rectangular grid subsetting grid probably involves use heres code far libraryrodbc librarygstat optionsdigits 10 channel odbcdriverconnect sqltableschannel pahdata sqlfetchchannel pahdata attachpahdata logtotalpah log10totalpah pahdata cbindpahdata logtotalpah classpahdata namespahdata coordinatespahdata easting northing sampleelevation classpahdata summarypahdata dimensionspahdata histlogtotalpah lpahvgm variogramlogtotalpah 1 pahdata lpahvgm plotlpahvgm lpahfit fitvariogramlpahvgm model vgm 6 sph 80 15 lpahfit plotlpahvgm lpahfit x gridtopologyc5345003531400 20 c22 05 c176126 31 classx coordinatevaluesx x spatialgridx classx grd asx spatialpixels lpahkriged krigelogtotalpah1 pahdata x model lpahfit classlpahkriged lpahkriged aslpahkriged spatialpixelsdataframe help would appreciated using r 270 windows xp desktop computer happily share scrubbed anonomized data requested matt findley ch2m hill corvallis oregon alternative html version deleted ______________________________________________ rhelprprojectorg mailing list httpsstatethzchmailmanlistinforhelp please read posting guide httpwwwrprojectorgpostingguidehtml provide commented minimal selfcontained reproducible code thu 07 aug 2008 140135 0600 r 3d kriging et al