R中的部分相关性

时间:2013-01-10 04:22:40

标签: r partial correlation

我正在尝试计算R中的部分相关性。我有两个我想要比较的数据集,目前只有一个受控变量。 (这将在未来发生变化)

我已经在网上试图自己解决这个问题,但很难理解我所看过的网站上使用的术语。有人可以解释一下我将如何做这个并提供一个简单的例子吗?

数据采用以下形式:

                Project.Name Bugs.Project Changes.Project Orgs.Project
1     platform_external_svox            4             161            2
3 platform_packages_apps_Nfc           13             223            2
5      platform_system_media           36             307            2
7     platform_external_mtpd            2              30            2
9            platform_bionic           42            1061            4

我希望Bugs.Project和Orgs.Project之间的关联与Changes.Project作为受控变量。我已经下载了ppcor库,因为它看起来具有我需要的功能。但是,我不确定如何使用它。如何将数据添加到矩阵并使用pcor函数?

这就是我一直在尝试的:

y.data <- data.frame(
bpp=c(projRelateBugsOrgs[2]),
opp=c(projRelateBugsOrgs[4]),
cpp=c(projRelateBugsOrgs[3])
)

test <- pcor(y.data)

我刚刚使用了一个我发现的例子,并试图用我的数据代替他们的数据。我不明白我的输出。

看起来像这样:

$estimate
                Bugs.Project Orgs.Project Changes.Project
Bugs.Project       1.0000000    0.3935535       0.9749296
Orgs.Project       0.3935535    1.0000000      -0.1800788
Changes.Project    0.9749296   -0.1800788       1.0000000

$p.value
                Bugs.Project Orgs.Project Changes.Project
Bugs.Project     0.00000e+00  2.09795e-07       0.0000000
Orgs.Project     2.09795e-07  0.00000e+00       0.0264442
Changes.Project  0.00000e+00  2.64442e-02       0.0000000

$statistic
                Bugs.Project Orgs.Project Changes.Project
Bugs.Project        0.000000     5.190442       53.122165
Orgs.Project        5.190442     0.000000       -2.219625
Changes.Project    53.122165    -2.219625        0.000000

$n
[1] 150

$gp
[1] 1

$method
[1] "pearson"

我想我想从$估价表中得到一些东西,但我并不完全确定它给了我什么,

1 个答案:

答案 0 :(得分:10)

在值部分

中读取help('pcor')
  

     

估计两个之间的偏相关系数的矩阵   变量

     

p.value 测试的p值矩阵

     

统计测试统计值

的矩阵      

n 样本数量

     

gn 给定变量的数量

     

方法使用的相关方法

详细信息部分提供

  

详细

     

部分相关是两个变量的相关性   控制第三个或更多其他变量。

为您的结果

$estimate
                Bugs.Project Orgs.Project Changes.Project
Bugs.Project       1.0000000    0.3935535       0.9749296
Orgs.Project       0.3935535    1.0000000      -0.1800788
Changes.Project    0.9749296   -0.1800788       1.0000000

Changes.ProjectOrgs.Project的部分相关性为-0.1800788。这是控制Changes.Project

Orgs.ProjectBugs.Project的相关性

Changes.ProjectBugs.Project的部分相关性为0.9747296。这是控制Changes.Project

Bugs.ProjectOrgs.Project的相关性

Orgs.ProjectBugs.Project的部分相关性为0.3935535。这是控制Orgs.Project

Bugs.ProjectChanges.Project的相关性

您可以从

获得相同的信息(如果您只对第三种情况感兴趣)
pcor.test(y.data$Orgs.Project, y.data$Bugs.Project, y.data$Changes.Project)