[R-sig-Geo] How to get point pairs form variogram cloud in gstat package

[Edzer Pebesma]

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 at r-project.org
>>>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at 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.pebesma at wwu.de