作者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