作者gsuper (绿色苏打心)
看板Statistics
标题[问题] 一个多变量的问题
时间Tue Dec 15 21:20:05 2009
实验目的 : 检定上下两个矩阵是否显着差异
____________________
/___/___/___/___/___/|
/___/___/___/___/___/ |
/___/___/___/___/___/ / (右边这面还有一条线 , 不画了)
/___/___/___/___/___/ /
control |___|___|___|___|___| /
experiment |___|___|___|___|___|/
Y轴 : control组和实验组 , 上层与下层为 paired data
Z轴 : 每次实验做一个向量 (左到右) , i.e. , 2*4 = 8次实验
X轴 : 每个格子代表浓度
真实 data 是 16*2*22 的 array space
但有一个现象
就是与X轴平行的向量
会有一个斜率
但所有向量的斜率又并非相等
所以图型变成以下
____
___/___/| 总之是像阶梯状
___/__/___/ |
/__/__/___/ /
/__/__/___/ /
/__/___|___| /
/___|___|___|/
|___|___|
|___| y轴和z轴都是常态分布
最主要的问题在於 x 轴斜率
请问是否有对应的统计方法
或是有什麽非简单平均数可以用来降维?
--
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.239.247
※ 编辑: gsuper 来自: 140.113.239.247 (12/15 21:20)
※ 编辑: gsuper 来自: 140.113.239.247 (12/15 21:22)
※ 编辑: gsuper 来自: 140.113.239.247 (12/15 21:26)
※ 编辑: gsuper 来自: 140.113.239.247 (12/15 21:26)
※ 编辑: gsuper 来自: 140.113.239.247 (12/15 21:27)
※ 编辑: gsuper 来自: 140.113.239.247 (12/15 21:36)
※ 编辑: gsuper 来自: 140.113.239.247 (12/15 21:36)
1F：推 west1996:图画的真好XD 12/16 01:24
2F：→ gsuper:画很久....不过远远比不上那个画30面体的神人 12/16 03:14
3F：→ kingbee:看不懂得你的问题 感觉是变异数分析须要先跑主成分 12/16 09:35
4F：推 west1996:看得不是很懂+1 不过既然是要比较control和experiment 12/16 10:11
5F：→ west1996:有没有显着差异，又已经是paired data，那x轴方向有一个 12/16 10:12
6F：→ west1996:斜率在有差吗？paired data应该会自动吸收掉那个影响不是 12/16 10:13
7F：→ west1996:吗？还是我误会了？ 12/16 10:13
问题是
斜率并不相等
斜率大小有 signal dependent 的性质
比方说 :
实验组侦测到有高表现
可能为 c( 1000,2000,3000,4000,5000,6000 ) 对index作图 , 呈现还算不错的线性
而控制组本身是无表现
可能为 c( 2 , 5 , 10 , 20 , 50 , 10 ) 100 以下视为 0 , 斜率为零
又或是控制组为低表现
可能为 c( 100 , 150 , 200 , 250 , 300 , 350 ) 与实验组的斜率不同
又或是控制组为中表现
可能为 c( 500,600,700,800,900,1000 ) 与实验组的斜率不同
---------------------------------------------------------------------
在上面几种状况
当数值都是几千在跑的时候
斜率较大
而在中表现和低表现时
斜率渐趋平缓 , 且可信度降低
-------------------------------------------------------------------
假设控制组是 C 矩阵
实验组是 E 矩阵
则相减以後是 D 矩阵 (differencial matrix)
D 矩阵的左到右向量 , 会有以下几种状况
1. C >>>>> E , i.e. , E = 0向量 , 则 D 的斜率为 C
2. C >> E
3. C > E D的斜率渐渐趋向平缓
4. C = E D的斜率完全被抵销 , 呈现常态分布
-------------------------------------------------------------------
※ 编辑: gsuper 来自: 140.113.177.165 (12/16 13:57)
8F：推 west1996:所以现在的问题在於如果不是常态就不能用原本想用的工具 12/16 14:01
9F：→ west1996:来分析，是这样吗？或许可以考虑无母数的方法避开母体假 12/16 14:01
10F：→ west1996:设 12/16 14:01
恩恩
我就是想问方法罗
因为我多变量才学个开头而已
要学到完我看要半年....
可是 Data 过一个月就要教了 (冏)
※ 编辑: gsuper 来自: 140.113.239.247 (12/16 14:48)
11F：→ andrew43:这和多变量有关系吗? 12/16 15:50
12F：→ andrew43:我看不出来哪里有多变量囧 12/16 15:51
13F：→ kingbee:向量 斜率 指的是什麽? 可以举的数据出来吗 12/16 16:57
14F：→ kingbee:实验数据列出来 12/16 16:57
data 转过 log2
因为其为 log normal distribution
C矩阵
N代表 Normal cell (
Z轴 )
上方 81~95 共8个 probes (
X轴 ) , 16个太长就不贴了
81 83 85 87 89 91 93 95
056N.CEL 13.52 13.56 13.53 13.51 13.54 13.54 13.5 13.5
064N.CEL 13.13 13.21 13.16 13.13 13.13 13.19 13.17 13.2
070N.CEL 13.01 13.01 13 12.97 12.97 13.08 13.02 13.01
095N.CEL 13.39 13.4 13.4 13.42 13.42 13.38 13.45 13.47
097N.CEL 13.02 13.04 13.06 13.04 13.04 13.04 13.05 13.07
104N.CEL 13.68 13.68 13.7 13.66 13.68 13.69 13.68 13.66
107N.CEL 12.94 12.95 12.96 12.94 12.95 12.94 12.93 12.94
122N.CEL 12.89 12.84 12.91 12.86 12.89 12.89 12.88 12.87
127N.CEL 12.8 12.78 12.78 12.81 12.83 12.82 12.8 12.85
134N2.CEL 13.07 13.07 13.06 13.07 13.06 13.12 13.11 13.14
156N.CEL 13.04 13.07 13.08 13.07 13.08 13.07 13.05 13.07
159N.CEL 13.18 13.2 13.22 13.16 13.15 13.22 13.16 13.21
164N.CEL 12.83 12.87 12.89 12.87 12.87 12.89 12.89 12.86
181N.CEL 12.68 12.69 12.66 12.69 12.68 12.75 12.66 12.67
184N.CEL 12.79 12.75 12.72 12.81 12.78 12.79 12.76 12.75
196N.CEL 13.1 13.04 13.02 13.02 13.08 13.02 13.03 12.99
208N.CEL 13.3 13.31 13.34 13.31 13.32 13.33 13.32 13.32
209N.CEL 13.72 13.7 13.77 13.73 13.75 13.82 13.79 13.78
210N.CEL 13.44 13.43 13.4 13.4 13.43 13.45 13.46 13.42
212N.CEL 13.32 13.3 13.32 13.3 13.29 13.34 13.32 13.31
226N.CEL 13.72 13.71 13.75 13.7 13.74 13.71 13.72 13.69
234N.CEL 13.46 13.41 13.43 13.42 13.48 13.47 13.47 13.45
E矩阵
T 代表 tumor
81 83 85 87 89 91 93 95
056T.CEL 13.17 13.2 13.26 13.22 13.24 13.25 13.23 13.23
064T.CEL 13.29 13.27 13.31 13.25 13.25 13.32 13.31 13.26
070T.CEL 13.37 13.3 13.33 13.33 13.31 13.29 13.4 13.34
095T.CEL 13.08 13.08 13.04 13.03 13.03 13.06 13.02 13.09
097T.CEL 12.76 12.79 12.76 12.8 12.76 12.8 12.83 12.8
104T.CEL 13.2 13.16 13.18 13.12 13.18 13.19 13.21 13.2
107T.CEL 12.86 12.81 12.82 12.83 12.84 12.87 12.83 12.8
122T.CEL 12.83 12.81 12.81 12.81 12.85 12.82 12.82 12.84
127T.CEL 12.91 12.93 12.91 12.88 12.87 12.9 12.87 12.87
134T2.CEL 13.1 13.12 13.11 13.09 13.09 13.12 13.15 13.14
156T.CEL 13.23 13.25 13.23 13.25 13.26 13.26 13.26 13.22
159T.CEL 13.4 13.39 13.36 13.39 13.38 13.39 13.41 13.4
164T.CEL 12.71 12.71 12.7 12.7 12.71 12.8 12.72 12.71
181T.CEL 12.84 12.88 12.85 12.84 12.83 12.87 12.86 12.86
184T.CEL 12.83 12.8 12.77 12.77 12.9 12.81 12.77 12.73
196T.CEL 13.04 13.05 13.06 13.05 13.01 13.05 12.99 13.04
208T.CEL 13.6 13.63 13.65 13.62 13.64 13.66 13.65 13.61
209T.CEL 13.79 13.79 13.79 13.79 13.81 13.79 13.81 13.77
210T.CEL 13.22 13.23 13.22 13.22 13.23 13.28 13.26 13.21
212T.CEL 13.19 13.2 13.18 13.07 13.08 13.14 13.1 13.11
226T.CEL 13.61 13.61 13.68 13.65 13.68 13.71 13.7 13.68
234T.CEL 13.8 13.77 13.8 13.77 13.77 13.82 13.78 13.75
C 和 D矩阵黏起来 (
Y轴 )
※ 编辑: gsuper 来自: 140.113.239.247 (12/16 17:57)
※ 编辑: gsuper 来自: 140.113.239.247 (12/16 17:59)
15F：→ kingbee:跟斜率、向量有什麽关系? 建议去看本实验设计 (沈明来) 12/17 09:17
16F：→ kingbee:你的问题太深奥了....我有看没有懂 orry 12/17 09:17