大师兄的应用回归分析学习笔记(二十二):主成分回归与偏最小二乘(一)

大师兄的应用回归分析学习笔记(二十一):岭回归(三)
大师兄的应用回归分析学习笔记(二十三):主成分回归与偏最小二乘(二)

一、主成分回归

  • W.F.Massy于1965年根据多元统计分析中的主成分分析提出了主成分回归。
  • 主成分回归(principal components regression,PCR)是对普通最小二乘估计的一种改进,他的参数估计是一种有偏估计。
1. 主成分的基本思想
  • 主成分分析(principal components analysis,PCA)首先由Hotelling于1933年提出。
  • PCA是一种降维思想,在随时很少信息的前提下把多个指标利用正交旋转变换转化为几个综合指标的多元统计分析方法。
  • 通常把转化生成的综合指标称为主成分,其中每个主成分都是原始变量的线性组合,且各个主成分之间互不相关。
  • 这样在研究复杂问题时就可以只考虑少数几个主成分且不至于损失太多信息,从而更容易抓住主要矛盾,揭示事物内部变量之间的规律性,同时使问题得到简化,提高分析效率。
2. 主成分的基本性质
  • 性质1:Y的协方差矩阵对角矩阵\Lambda。其中对角线上的值为\lambda_1,\lambda_2,...,\lambda_p
  • 性质2:记\sum=(\delta_{ij})_{p\times p},由\sum^p_{i=1}\lambda_i=\sum^p_{i=1}\delta_{ii}
  • \alpha_k=\frac{\lambda_k}{\lambda_1+\lambda_2+...+\lambda_p}(k=1,2,...,p)为第k个主成分Y_k的方差贡献率
  • \frac{\sum^m_{i=1}\lambda_i}{\sum^p_{i=1}\lambda_i}为主成分Y_1,Y_2,...,Y_m的累计贡献率。
  • 性质3\rho(Y_k,X_i)=\mu_{ki}\sqrt{\lambda_k}/\sqrt{\delta_{ii}}(k,i=1,2,...,p)
  • 式中,第k个主成分Y_k与原始变量X_i的相关系数\rho(Y_k,X_i)称为因子负荷量
  • 因子负荷量是主成分分析解释中非常重要的解释依据,因子负荷量的绝对值大小刻画了该主成分的主要意义及其成因。
  • 性质4\sum^p_{i=1}\rho^2(Y_k,X_i)\delta_{ii}=\lambda_k
  • 性质5\sum^p_{i=1}\rho^2(Y_k,X_i)=\frac{1}{\delta_{ii}}\sum^p_{k=1}\lambda_k\mu^2_{ki}=1
  • X_i与前m个主成分Y_1,Y_2,...,Y_m的全相关系数平方和称为Y_1,Y_2,...,Y_m对原始变量X_i的方差贡献率v_i,即v_i=\frac{1}{\alpha}\sum^m_{k=1}\lambda_k\mu^2_{ki}(i=1,2,...,p)
  • 这一定义说明前m个主成分提取了原始变量x_iv_i的信息,由此可以判断提取的主成分说明原始变量的能力。
3. 案例
  • 以民航客运数据为例:


  • 首先对5个自变量计算主成分:
FACTOR 
  /VARIABLES x1 x2 x3 x4 x5 
  /MISSING LISTWISE 
  /ANALYSIS x1 x2 x3 x4 x5 
  /PRINT INITIAL EXTRACTION 
  /CRITERIA MINEIGEN(1) ITERATE(25) 
  /EXTRACTION PC 
  /ROTATION NOROTATE 
  /SAVE REG(ALL) 
  /METHOD=CORRELATION.
  • 输出中有5个主成分的特征根,最大的是\lambda_1=4.954,最小的是lambda_5=0.001
  • 方差百分比反应主成分所能结束数据变异的比例,也就是包含原始数据的信息比例。
  • 第一个主成分Factor1的方差百分比=99.071%,含有5个原始变量99%以上的信息量,因此本例取一个主成分已经足够了。
  • 可以获得5个主成分得分:


  • 这5个主成分得分每列数据的平均值为0
  • 平方和都是n-1=10
  • 任意两列都是线性无关的
  • 由于第一个主成分Factor1的方差比已经高达99.071%,可以对这一个主成分做主成分分析:


  • 用y对Factor1做普通最小二乘回归,得主成分回归方程:\hat y=20345.15+13639.253Factor1
  • 下面需要把这个表达式还原回用5个原始自变量表达的形式,这需要找出Factor1和5个原始自变量x_1,x_2,x_3,x_4,x_5之间的关系式。
  • 只需要以Factor1为因变量,以x_1,x_2,x_3,x_4,x_5为自变量做回归,得到的回归方程就是所需要的关系式:
  • Factor1 = -1.763 + 8.892 ^{-7} x_1+3.659^{-5}x_2+2.548^{-7}x_3+3.519^{-6}x_4+1.663^{-4}x_5
  • 代会主成分回归方程,得到y对5个原始自变量的主成分回归方程为:\hat y = -3.698.29+0.01212x _1+0.499x_2+0.003475x_3+0.048x_4+2.268x_5
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 230,362评论 6 544
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 99,577评论 3 429
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 178,486评论 0 383
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 63,852评论 1 317
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 72,600评论 6 412
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 55,944评论 1 328
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 43,944评论 3 447
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 43,108评论 0 290
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 49,652评论 1 336
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 41,385评论 3 358
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 43,616评论 1 374
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 39,111评论 5 364
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 44,798评论 3 350
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 35,205评论 0 28
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 36,537评论 1 295
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 52,334评论 3 400
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 48,570评论 2 379

推荐阅读更多精彩内容