[問題] 一個多變量的問題
實驗目的 : 檢定上下兩個矩陣是否顯著差異
____________________
/___/___/___/___/___/|
/___/___/___/___/___/ |
/___/___/___/___/___/ / (右邊這面還有一條線 , 不畫了)
/___/___/___/___/___/ /
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)
推
12/16 01:24, , 1F
12/16 01:24, 1F
→
12/16 03:14, , 2F
12/16 03:14, 2F
→
12/16 09:35, , 3F
12/16 09:35, 3F
推
12/16 10:11, , 4F
12/16 10:11, 4F
→
12/16 10:12, , 5F
12/16 10:12, 5F
→
12/16 10:13, , 6F
12/16 10:13, 6F
→
12/16 10:13, , 7F
12/16 10:13, 7F
問題是
斜率並不相等
斜率大小有 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)
推
12/16 14:01, , 8F
12/16 14:01, 8F
→
12/16 14:01, , 9F
12/16 14:01, 9F
→
12/16 14:01, , 10F
12/16 14:01, 10F
恩恩
我就是想問方法囉
因為我多變量才學個開頭而已
要學到完我看要半年....
可是 Data 過一個月就要教了 (冏)
※ 編輯: gsuper 來自: 140.113.239.247 (12/16 14:48)
→
12/16 15:50, , 11F
12/16 15:50, 11F
→
12/16 15:51, , 12F
12/16 15:51, 12F
→
12/16 16:57, , 13F
12/16 16:57, 13F
→
12/16 16:57, , 14F
12/16 16:57, 14F
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)
→
12/17 09:17, , 15F
12/17 09:17, 15F
→
12/17 09:17, , 16F
12/17 09:17, 16F
→
01/02 15:01,
5年前
, 17F
01/02 15:01, 17F