我有一个100列的表,我希望运行成对的部分相关,使用pcor.test
包中的ppcor
函数控制第100列的变量。在R中是否存在任何部分相关函数,我可以使用类似rcorr
的返回值,采用整个矩阵的成对相关性,但仅控制一个变量?
答案 0 :(得分:4)
听起来对于n列矩阵,您想要输出前n-1列的成对相关的(n-1)x(n-1)矩阵,控制最后一个(使用{{来自combineLatest:reduce:
包的1}}函数。
您可以使用pcor.test
函数执行此操作,循环遍历每个列并使用ppcor
计算与所有其他列的相关性:
sapply