Please send me a reproducable example, for instance as a variation on: > library(gstat) Loading required package: sp > data(meuse) > coordinates(meuse)=~x+y > v = variogram(zinc~1,meuse,cloud=TRUE) numeric(0) > class(v) [1] "variogramCloud" "data.frame" > as.data.frame(v)$left[1:10] [1] 2 3 3 4 4 4 5 5 5 5
On 03/14/2011 06:43 PM, vetter wrote: > I guess var1 is a variogramCloud. More complete code: > >> var1=variogram(val~1,geoObject,cloud=T) >> var1[1,] > dist gamma dir.hor dir.ver id left right > 1 63834.75 0.08506403 0 0 var1 2 1 >> as.data.frame(var1)$left > NULL > > > str(var1) gives: > > Classes ‘variogramCloud’ and 'data.frame': 9253 obs. of 6 variables: > $ np : num 1 2 65538 3 65539 ... > $ dist : num 63835 97833 158641 251475 312811 ... > $ gamma : num 0.085 2.222 1.438 7.132 5.659 ... > $ dir.hor: num 0 0 0 0 0 0 0 0 0 0 ... > $ dir.ver: num 0 0 0 0 0 0 0 0 0 0 ... > $ id : Factor w/ 1 level "var1": 1 1 1 1 1 1 1 1 1 1 ... > - attr(*, "direct")='data.frame': 1 obs. of 2 variables: > ..$ id : Factor w/ 1 level "var1": 1 > ..$ is.direct: logi TRUE > - attr(*, ".BigInt")= num 65536 > > >> make sure var1 is a variogramCloud, obtained by >> >> variogram(..., cloud=TRUE) >> >> On 03/14/2011 06:27 PM, vetter wrote: >>> It is not working: >>> >>>> as.data.frame(var1)[1,] >>> np dist gamma dir.hor dir.ver id >>> 1 1 63834.75 0.2057185 0 0 var1 >>> >>> and therefore: >>> >>>> as.data.frame(var1)$left >>> NULL >>> >>> >>>> On 03/14/2011 05:47 PM, vetter wrote: >>>>> I computed a variogram cloud using the gstat package. >>>>> >>>>> My results look like this: >>>>> >>>>>> var1[,1] >>>> that was probably a >>>> >>>> var[1,] >>>> >>>>> dist gamma dir.hor dir.ver id left right >>>>> 1 63834.75 0.2057185 0 0 var1 2 1 >>>>> >>>>> My question is: >>>>> How can I extract the values from column "left" and "right" from this >>>>> object? >>>>> >>>> By coercion to a data.frame, and selecting: >>>> >>>> as.data.frame(var1)$left >>>> as.data.frame(var1)$right >>>> >>>>> Thanks Tobias Vetter >>>>> >>>>> _______________________________________________ >>>>> R-sig-Geo mailing list >>>>> R-sig-Geo@r-project.org >>>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >>> _______________________________________________ >>> R-sig-Geo mailing list >>> R-sig-Geo@r-project.org >>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebe...@wwu.de _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo