This seems to explain it.
The definition of
orderis thata[order(a)]is in
increasing order. This works with your example, where the correct
order is the fourth, second, first, then third element.You may have been looking for
rank, which returns the rank of the
elements
R> a <- c(4.1, 3.2, 6.1, 3.1)
R> order(a)
[1] 4 2 1 3
R> rank(a)
[1] 3 2 4 1
soranktells you what order the numbers are in,
ordertells you how to get them in ascending order.
plot(a, rank(a)/length(a))will give a graph of the CDF. To see why
orderis useful, though, tryplot(a, rank(a)/length(a),type="S")
which gives a mess, because the data are not in increasing orderIf you did
oo<-order(a)
plot(a[oo],rank(a[oo])/length(a),type="S")
or simply
oo<-order(a)
plot(a[oo],(1:length(a))/length(a)),type="S")
you get a line graph of the CDF.
I’ll bet you’re thinking of rank.