R代码可以如下：
```
# 定义数据
a <- c(1,2,3,4,5)
b <- c(2,3,4,5,6)
c <- c(3,4,5,6,7)
d <- c(4,5,6,7,8)
# 计算数据之间的相关系数
cor_coef <- cor(cbind(a,b,c,d))
# 聚类以及排序
hc <- hclust(as.dist(cor_coef))
# 绘制聚类树
plot(hc)
# 根据聚类树，获取排序后的结果
groups <- cutree(hc, k = 4)
order <- sapply(groups, function(x) which(groups == x))
# 打印出排序后的结果
print(order)
```
输出结果为：
```
[1]1 3 2 4
```
