作者gsuper (统计的巴比伦塔)
看板Statistics
标题[问题] Paired AUC test 的相关性
时间Thu Jun 23 22:00:54 2011
paired AUC test 时
公式为
AUC1 - AUC2
-------------------------------------------------------
( SE(AUC1)^2+SE(AUC2)^2 + 2*r*SE(AUC1)*SE(AUC2) )^0.5
其中的 r 为 r of control 与 r of case 的合成值 (相关系数)
--------------------------------------------------------------
我目前在做的
是组合多重 SNP markers
计算 AUC 和 performance
假设一开始有 4 个 markers
筛选完剩下 2 个
想比较後者是否损失太多 AUC
又因为是对同一组 data 进行筛选
所以选用 paired AUC test
但现在麻烦的是
筛选前後的 markers 是呈 nested relation (大model 完全包覆小model)
所以估计一定会有很强的相关性
(前) Yi = b1*x1 + b2*x2 + b3*x3 + b4*x4 + 残差
(後) Yi = b1*x1 + b2*x2 + 残差
这样是否有办法计算 (b1~b4) versus (b1~b2) 的相关性?
若要用到共变异矩阵
我应该要去查甚麽章节或关键字?
或是有甚麽降维的方法可用?
-------------------------------------------
我是读这两份文件
http://www.graphpad.com/faq/viewfaq.cfm?faq=1191
http://www.anaesthetist.com/mnm/stats/roc/Findex.htm
诚心请教
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.239.247
※ 编辑: gsuper 来自: 140.113.239.247 (06/23 22:01)
※ 编辑: gsuper 来自: 140.113.239.247 (06/23 22:01)
※ 编辑: gsuper 来自: 140.113.239.247 (06/23 22:02)
※ 编辑: gsuper 来自: 140.113.239.247 (06/23 22:07)
※ 编辑: gsuper 来自: 140.113.239.247 (06/23 22:08)
※ 编辑: gsuper 来自: 140.113.239.247 (06/23 22:09)
※ 编辑: gsuper 来自: 140.113.239.247 (06/23 22:10)
1F：→ kingbee:你是做ROC 曲线的比较吗 还是两个回归式的比较? 06/24 16:56
2F：→ kingbee:虽然不太懂你的问题 但想先确认一下你的Yi是什麽型态 06/24 16:56
3F：→ kingbee:连续资料还是二元 06/24 16:57
每个 x1 是三元资料 , 假设有 4组 xi (markers)
Yi(总风险) 就会有 3^4 种可能性 , 也是一种多元分布
e.g.
X1 是 rs671
有 GG GA AA 三种事件
单一风险值分别是 1 , 1 , 2
X2 是 rs799
有 AA AT TT 三种事件
单一风险值分别是 1 , 2 , 4
Yi 就变成 3^2 种组合事件
(事先已确认 X1 与 X2 没有交互作用或共线性)
e.g. log(1) + log(1)
log(1) + log(2)
log(1) + log(4)
log(1) + log(1)
log(1) + log(2)
log(1) + log(4)
log(2) + log(1)
log(2) + log(2)
log(2) + log(4)
※ 编辑: gsuper 来自: 140.113.239.247 (06/24 18:40)
※ 编辑: gsuper 来自: 140.113.239.247 (06/24 18:41)
※ 编辑: gsuper 来自: 218.160.244.59 (06/24 23:22)
※ 编辑: gsuper 来自: 140.113.239.247 (06/27 14:27)