回复 第3楼 的 谢益辉:
<br />
#读入数据<br />
testdata <- read.csv("test.csv")<br />
testdata<br />
#ltm包及相关检验<br />
library(ltm)<br />
#临时数据框<br />
testmat <- cbind(height = testdata$height, shoe = testdata$shoe, hand = testdata$hand)<br />
#相关检验<br />
rcor.test(testmat)<br />
以上代码执行后会得到如下结果:
</p>
<p> height shoe hand<br />
height ***** 0.867 0.574<br />
shoe <0.001 ***** 0.516<br />
hand 0.001 0.004 *****</p>
<p>upper diagonal part contains correlation coefficient estimates<br />
lower diagonal part contains corresponding p-values<br />
</p>
我想要的输出形式是:
height和shoe显著相关,相关系数是0.867;
height和hand显著相关,相关系数是0.574;
shoe和hand显著相关,相关系数是0.516。
因为现在变量少还好,如果变量很多,即便是相关分析得到的表格装得下,看起来也不方便,如果可以用文字的形式把显著相关的变量及其系数列出来,似乎更有利于信息的传达,尤其是给非专业人士看的时候。故有此问。希望大家可以指点在下一二。