Statistics 板


LINE

※ 引述《iphone2003 (307)》之铭言: : 大家好 : 我最近在练习把各种回归提到的东西实作出来 : 使用的是python : 然後使用statsmodels来验证我的结果 : 那我正在写WLS : 作法是令一个矩阵W,其中它的各个对角线元素为weighted值 : 然後把原本的矩阵X,转变为Xw= W^0.5 * X : 向量y也做一样的转换:yw = W^0.5 *y : 然後利用Xw和yw去做一般的OLS计算 : 所以这时候算的系数向量会是 : b = (Xw' * Xw)^(-1) * Xw' * yw : = (X' * W * X)^(-1) * X' * W * y : 那我现在遇到的问题是 : 我在计算SST SSR SSE时 : 一开始想说应该和一般的计算方法一样 : 就是像是SSE=(y-yhat)' * (y-yhat) : 其中yhat是X * b : 结果算出来好像怪怪的 : 後来改成这样算 : SSE = (yw-ywhat)' * (yw-ywhat) : 其中ywhat是Xw * b : 这样结果就对了 : 但是我不懂为什麽这里要改用yw去算? : 另外同样的改法,用在SST和SSR却失败 : 算出来的答案不太对 : 因此想问,是否WLS中,计算SSR、SSE和SST的方法和一般OLS的不太一样呢? : 是的话,又是为什麽要这样计算呢? : 谢谢! : (这里我的SSR SSE SST分别代表回归变异量、误差变异量和总变异量) WLS的optimization problem是 min (y-XB)^T W (y-XB) 这个量是什麽意思? 你一般在做optimization就是在最小化SSE 所以SSE = (y-XB)^T W (y-XB) = sum_i w_i (y_i-X_i B)^2 自己拆一下就跟你那个式子一样了 那SST是什麽? SST是总变异量,在OLS是 sum_i (y_i-y bar)^2 WLS因为多了weight 所以变成 sum_i w_i (y_i - y bar)^2 写成矩阵就是 (y-y bar)^T W (y-y bar) SST = SSR + SSE,那SSR就自己推导吧 我看你现在只是在练习把演算法做出来而已 不是真的在理解它到底为何这样做,跟背後的概念 这样是无助於你学习的 --



※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 118.163.170.73
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1552271897.A.624.html
1F:推 iphone2003: 感谢你!这样看就懂了,本来把W并进Xy所以搞不太清楚 03/11 16:30
2F:→ iphone2003: ,不过分开来看就能理解了 03/11 16:30
※ 编辑: celestialgod (119.14.59.166), 03/11/2019 21:06:26







like.gif 您可能会有兴趣的文章
icon.png[问题/行为] 猫晚上进房间会不会有憋尿问题
icon.pngRe: [闲聊] 选了错误的女孩成为魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一张
icon.png[心得] EMS高领长版毛衣.墨小楼MC1002
icon.png[分享] 丹龙隔热纸GE55+33+22
icon.png[问题] 清洗洗衣机
icon.png[寻物] 窗台下的空间
icon.png[闲聊] 双极の女神1 木魔爵
icon.png[售车] 新竹 1997 march 1297cc 白色 四门
icon.png[讨论] 能从照片感受到摄影者心情吗
icon.png[狂贺] 贺贺贺贺 贺!岛村卯月!总选举NO.1
icon.png[难过] 羡慕白皮肤的女生
icon.png阅读文章
icon.png[黑特]
icon.png[问题] SBK S1安装於安全帽位置
icon.png[分享] 旧woo100绝版开箱!!
icon.pngRe: [无言] 关於小包卫生纸
icon.png[开箱] E5-2683V3 RX480Strix 快睿C1 简单测试
icon.png[心得] 苍の海贼龙 地狱 执行者16PT
icon.png[售车] 1999年Virage iO 1.8EXi
icon.png[心得] 挑战33 LV10 狮子座pt solo
icon.png[闲聊] 手把手教你不被桶之新手主购教学
icon.png[分享] Civic Type R 量产版官方照无预警流出
icon.png[售车] Golf 4 2.0 银色 自排
icon.png[出售] Graco提篮汽座(有底座)2000元诚可议
icon.png[问题] 请问补牙材质掉了还能再补吗?(台中半年内
icon.png[问题] 44th 单曲 生写竟然都给重复的啊啊!
icon.png[心得] 华南红卡/icash 核卡
icon.png[问题] 拔牙矫正这样正常吗
icon.png[赠送] 老莫高业 初业 102年版
icon.png[情报] 三大行动支付 本季掀战火
icon.png[宝宝] 博客来Amos水蜡笔5/1特价五折
icon.pngRe: [心得] 新鲜人一些面试分享
icon.png[心得] 苍の海贼龙 地狱 麒麟25PT
icon.pngRe: [闲聊] (君の名は。雷慎入) 君名二创漫画翻译
icon.pngRe: [闲聊] OGN中场影片:失踪人口局 (英文字幕)
icon.png[问题] 台湾大哥大4G讯号差
icon.png[出售] [全国]全新千寻侘草LED灯, 水草

请输入看板名称,例如:Tech_Job站内搜寻

TOP