今天的内容应该算是比较简单的，如果看着看着觉得复杂了，那绝对是被我乱整给整复杂的

咱提个简单明了的问题先，知道协方差矩阵不？或者相关系数矩阵也行，反正这俩长得差不多

R里头提供的cov和cor可以直接产生协方差或相关系数的矩阵，但有点了解的亲们应该知道，这种显示方式从表格形式来说是一种宽格式，而且对角线两边的数是重复的，也就是说它虽然显示得好看，但要后续做提取啥的还是有点小繁琐的，尤其是当我们希望提取出数值的同时还要带上相关的标题，那更是格式转着转着，一不小心就找不着了

上例子：

从cov和cor的帮助文件里，咱可以看到一点端倪，那就是它其实支持两种输入方式，一种是像上面这样把整个二维表放进去作为x，这时计算出来的是协方差矩阵，还有一种是同时输入x和y，但两者都是向量，这样计算结果就只有一个数值

如果像上边那样需要3列数据间的协方差，但又不想用宽格式来显示的话，那咱需要的表至少得包含两部分：

1） 由三列数据的列标题两两组合得到的信息列

2） 由标题组合计算出来的协方差

先上手做第一部分，3列标题两两组合，意味着是3个中抽取两个，用函数choose可以先预判下这个组合的长度（这个不是必要步骤，可以省略），而用函数combn可以直接生成组合后的矩阵

Index <- names(Table1[,c(-1,-3)]) #抽取出3列的标题

COVlist <- combn(x = Index,m = 2) #两两组合，抽取数为2

输出的形式行列方向是反的，所以后边做的时候都加了转置函数t()

接下来写得有点长，做完才发觉自己不知不觉的就挖了个坑把自己埋了

#按列名称从Table1中提取数据的function

Pickout <- function(x) {Table1[,colnames(Table1) == x]}

Inputlist <- apply(t(COVlist),c(1,2),Pickout) #把Table1提取出来的6列组成3+3的一个list

COVwithPick <- function (x) {cov(Inputlist[,,1][,x],Inputlist[,,2][,x])}

#做个提取指定列的function计算协方差cov

Value <- apply(data.frame(c(1:3)),1,COVwithPick)

#把标题和协方差计算结果合并到一起

COVresult <- data.frame(t(COVlist),Value)

做倒是做出结果来了，3列两两组合出三个协方差，就是步骤嘛，看着有点费事

啊咧~~~早知道费这劲我还用啥组合数列呀，唉~~~~