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

[vetter]

Problem solved. In the newest R-version it works fine.
Thanks for helping
Tobias Vetter
> strange still NULL:
>
> > data(meuse)
> > coordinates(meuse)=~x+y
> > v = variogram(zinc~1,meuse,cloud=TRUE)
> > as.data.frame(v)$left[1:10]
> NULL
>
> My current Version R version is  2.9.1.
>
>
>> 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
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo