[R] Join points with arrows based a TIME variable
Uwe Ligges
ligges at statistik.tu-dortmund.de
Wed Nov 4 18:35:13 CET 2009
zhijie zhang wrote:
> Hi Uwe,
> It works. Why does the following argument generate different results?
> with(data, arrows(x[-length(x)], y[-length(x)], x[-1], y[-1]))
> #correct,multiple joins
In the former you select all elements of the vectors except the first or
last one, respectively.
In the next call you select just one arguemnt of each vector, namely the
first or last one.
Uwe Ligges
> with(data, arrows(x[1], y[1], x[length(x)], y[length(x)])) #wrong,seems to
> be only one join
> From the ?arrows, the second argument should work. What is the problem?
> Thanks a lot.
> #codes:
> x<-c(1:6);y<-c(1:6);time<-c(6:1);data<-cbind(x,y,time);data
> data<-data.frame(data[order(data[,"time"]),])
> with(data, plot(x,y))
> with(data, arrows(x[-length(x)], y[-length(x)], x[-1], y[-1]))
>
> with(data, plot(x,y))
> with(data, arrows(x[1], y[1], x[length(x)], y[length(x)]))
> 2009/11/4 Uwe Ligges <ligges at statistik.tu-dortmund.de>
>
>>
>> zhijie zhang wrote:
>>
>>> Hi,
>>> I have a data set with three variables,X Y and Time. X and Y are the
>>> coordinates of points, i want to join these points according to the Time
>>> sequence using arrows?
>>> Demo Example data:
>>>
>>>> x<-c(1:6)
>>>> y<-c(1:6)
>>>> time<-c(6:1)
>>>> data<-cbind(x,y,time)
>>>> data
>>>>
>>> x y time
>>> [1,] 1 1 6
>>> [2,] 2 2 5
>>> [3,] 3 3 4
>>> [4,] 4 4 3
>>> [5,] 5 5 2
>>> [6,] 6 6 1
>>> I hope to join the six points with points' time=1 as starting point and
>>> points' time=6 as endpoint. So the sequence is time=1,2,3,4,5,6 and join
>>> the
>>> corresponding points with arrows.
>>> Any ideas on it?
>>>
>>
>> data <- data.frame(data[order(data[,"time"]),])
>> with(data, plot(y ~ x))
>> with(data, arrows(x[-length(x)], y[-length(x)], x[-1], y[-1]))
>>
>>
>> Uwe Ligges
>>
>>
>> [[alternative HTML version deleted]]
>>> ______________________________________________
>>>
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>>> http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html>
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>
More information about the R-help
mailing list