vlambda博客
学习文章列表

R语言完成偏相关性分析


今天晚上回家,看了很多偏相关性分析资料。资料里能看懂的就只有这一点儿:“偏相关分析是指当两个变量同时与第三个变量相关时,将第三个变量的影响剔除,只分析另外两个变量之间相关程度的过程,判定指标是相关系数的R值。”其他都看不太懂。 R语言完成偏相关性分析 R语言完成偏相关性分析
另外还有很重要的一点,偏相关系数使用的两个条件:1. 样本间独立性;2.数据服从多元正态分布。
拿大白话来;说“现有变量A,B,C,D,则四个变量同时相关,想分析A,B的相关性,但把C,D影响剔除。”问题的方法。
学习下R语言实现的方法:
#加载ggm包
install.packages("ggm")
library(ggm)
偏相关系数的计算
函数调用格式:pcor(u,s)
其中u是一个数值向量,前两个数值表示要计算相关系数的变量下标,其余数值为条件变量的下标。
偏相关系数显著性检验
函数调用格式:pcor.test(r,q,n)
其中r是由pcor()函数计算得到的偏相关系数,q为要控制的变量数(以数值表示位置),n为  样本大小。
例如:
install.packages("ggm")#安装ggm包;
library(ggm) #读取ggm包;
setwd("") #设置工作目录,双引号内填写工作目录;
x <-read.table("x.txt",sep="\t",header=T) #读取数据
pcor(c(1,4,3,2),cov(x)) #计算偏相关系数
pcor.test(pcor(c(1,4,3,2),cov(x)),2,60)#偏相关系数假设检验

R语言完成偏相关性分析

随便从电脑里找了些数据,试了试没报错,能做出这个数据的偏相关分析结果 。但 其实我自己对此一无所知,将来学了数学也许就懂了。
浑浑噩噩,来日方长,今日晚安。