主要内容

partialcorri

调整内部变量的偏相关系数

描述

例子

ρ= partialcorri (yx中变量对之间的样本线性偏相关系数yx,调整中的其余变量x

例子

ρ= partialcorri (yxz中变量对之间的样本线性偏相关系数yx,调整中的其余变量x,先控制两者xy对于z

例子

ρ= partialcorri (___名称,值使用来自前面任何语法的输入参数,返回带有一个或多个名称-值对参数指定的附加选项的样例线性偏相关系数。例如,您可以指定是使用皮尔逊偏相关还是斯皮尔曼偏相关,或者指定如何处理缺失值。

例子

ρpval) = partialcorri (___也返回一个矩阵pvalp-值,用于检验不存在偏相关的假设与存在非零偏相关的单边或双边备选方案。

例子

全部折叠

计算中每对变量的偏相关系数xy输入矩阵,同时控制其余变量的影响x

加载示例数据。

负载carsmall

这些数据包括1970年、1976年和1982年生产的汽车的测量数据。它包括英里/加仑加速度作为绩效衡量标准位移马力,重量作为设计变量。加速度从0加速到60英里每小时需要时间吗加速度对应的是低加速度的车辆。

定义输入矩阵。的y矩阵包括绩效度量,以及x矩阵包括设计变量。

y = (MPG、加速度);x =(位移、马力、重量);

计算相关系数。只包括计算中没有缺失值的行。

ρ= partialcorri (y, x,“行”“完成”
ρ=2×3-0.0537 -0.1520 -0.4856 -0.3994 -0.4008 0.4912

例如,结果表明,在控制了位移和马力的影响后,重量和加速度之间的相关性为0.4912。你可以退回 p -值作为第二个输出,并检查它们,以确认这些相关性是否具有统计学意义。

为了更清晰地显示,请创建一个带有适当变量和行标签的表。

ρ= array2table(ρ,...“VariableNames”, {“位移”“马力”“重量”},...“RowNames”, {“英里”“加速”});disp (的偏相关系数
偏相关系数
disp(ρ)
位移马力重量  ____________ __________ ________ 英里/加仑-0.053684 -0.15199 -0.48563 -0.39941 -0.40075 0.49123的加速度

中对变量之间的偏相关检验xy输入矩阵,同时控制其余变量的影响x加上矩阵中的其他变量z

加载示例数据。

负载carsmall

这些数据包括1970年、1976年和1982年生产的汽车的测量数据。它包括英里/加仑加速度作为绩效衡量标准位移马力,重量作为设计变量。加速度从0加速到60英里每小时需要时间吗加速度对应的是低加速度的车辆。

创建一个新变量逆风,并随机生成数据来表示沿着性能测量路线的平均逆风的概念。

rng (“默认”);%的再现性逆风= (10:-0.2:-9.8)' + 5*randn(100,1);

由于逆风可以影响性能度量,在测试其余变量之间的偏相关性时,控制其影响。

定义输入矩阵。的y矩阵包括绩效度量,以及x矩阵包括设计变量。的z矩阵包含额外的变量,用于计算偏相关时的控制,如逆风。

y = (MPG、加速度);x =(位移、马力、重量);z =阻力;

计算偏相关系数。只包括计算中没有缺失值的行。

(ρ,pval) = partialcorri (y, x, z,“行”“完成”
ρ=2×30.0572 -0.1055 -0.5736 -0.3845 -0.3966 0.4674
pval =2×30.5923 0.3221 0.0000 0.0002 0.0001 0.0000

返回的小 p -value of 0.001 inpval指示,例如,在控制了位移、重量和逆风之后,马力和加速度之间的显著负相关。

为了更清晰地显示,请创建带有适当变量和行标签的表。

ρ= array2table(ρ,...“VariableNames”, {“位移”“马力”“重量”},...“RowNames”, {“英里”“加速”});pval = array2table (pval,...“VariableNames”, {“位移”“马力”“重量”},...“RowNames”, {“英里”“加速”});disp (“考虑逆风的偏相关系数”
偏相关系数,考虑逆风
disp(ρ)
位移马力重量  ____________ __________ ________ 英里/加仑0.057197 -0.10555 -0.57358 -0.38452 -0.39658 0.4674的加速度
disp (“p值,计入逆风”
p值,计入逆风
disp (pval)
排量马力重量____________ __________ __________ MPG 0.59233 0.32212 3.4401e-09加速度0.00018272 0.00010902 3.4091e-06

输入参数

全部折叠

数据矩阵,指定为n——- - - - - -px矩阵。的行x对应观察值,列对应变量。

数据类型:|

数据矩阵,指定为n——- - - - - -py矩阵。的行y对应观察值,列对应变量。

数据类型:|

数据矩阵,指定为n——- - - - - -pz矩阵。的行z对应观察值,列对应变量。

数据类型:|

名称-值参数

指定可选的逗号分隔的对名称,值参数。的名字参数名和价值为对应值。的名字必须出现在引号内。可以以任意顺序指定多个名称和值对参数Name1, Value1,…,的家

例子:“类型”、“枪兵”,“行”,“完成”仅使用不包含缺失值的行中的数据计算斯皮尔曼部分相关性。

要计算的部分关联的类型,指定为由逗号分隔的对组成“类型”,要么皮尔森的“枪兵”皮尔森计算皮尔逊(线性)偏相关。斯皮尔曼计算斯皮尔曼(秩)偏相关。

例子:“类型”,“枪兵”

要在计算中使用的行,指定为组成的逗号分隔对“行”下面是其中之一。

“所有” 使用所有行,不管是否丢失(南)值。
“完成” 只使用没有缺失值的行。
“成对” 使用的每一列中的所有可用值y在计算偏相关系数时p对应于该列的值。的每一列y,将删除相应于中缺失值的行x(和/或z,如果提供)。的列中具有有效值的其余行y,即使在?的其他列中有缺失值y

例子:“行”,“完成”

用于测试的可选假设,指定为逗号分隔的对,由“尾巴”下面是其中之一。

“两个” 检验相关性不为零的替代假设。
“对” 检验相关性大于0的替代假设。
“左” 检验相关性小于0的替代假设。

例子:“尾巴”,“对”

输出参数

全部折叠

样本线性偏相关系数,返回为py——- - - - - -px矩阵。

  • 如果你输入xy矩阵,(j)第th项为样本之间的线性偏相关th列yjth列x,控件的所有列x除了列j

  • 如果你输入xy,z矩阵,(j)第th项为样本之间的线性偏相关th列yjth列x的所有列进行了调整x除了列j,先控制两者xy对于z

p-values,以矩阵形式返回。的每个元素pvalp的对应元素的值ρ.如果pval (i, j)是小的,那么相应的偏相关呢ρ(i, j)在统计上与零有显著差异。

partialcorri计算p-值的线性和排序偏相关使用学生t相关变换的分布。这对于线性偏相关是准确的xz是正常的,但在其他方面是一个大样本近似。

参考文献

斯图尔特,艾伦,K.奥德和S.阿诺德。肯德尔的高级统计理论。第6版,2A卷,第28章,威利,2004年。

罗纳德·A·费雪偏相关系数的分布密特隆3 (1924): 329 - 332

另请参阅

|

介绍了R2013b