重複分析<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />大家好,
我想根據附帶的示例處理表中包含的數據。該表包含100列+一列作為控制措施,並包含30到80行。
目的是以迭代的方式應用統計測試,從對照與列51,對照與列52,對照與列53等開始,以便檢測何時出現顯著差異,在連續列上存在多長時間以及當意義喪失。
例如,可以從列63到列78檢測差異,然後從c79到c86丟失,然後再次出現在c87到c98上。
我不知道如何編寫一個腳本,允許自動分析這種類型的表,給出一個報告差異發生的時刻。
有什麼建議嗎?
Re:重複分析這是交換建議的建設性方式!謝謝你,凱文。有趣的是,我們可以與R和JMP來回走動。我一定會仔細研究一下變更點程序。
Re:重複分析嗨,保羅!
沒有更多的信息,我所做的不僅僅是在黑暗中拍攝。
但是,如果我懷疑您正在搜索順序收集的數據以試圖辨別數據生成過程的變化,那麼變換點方法可能更合理。
有許多方法可以檢測變化點,並且有許多參考資料可以追溯到很多年。
CRAN有一個變點在R中實現幾個最近研究過的方法的包.JMP有一些簡潔的R集成,你可以在JMP數據集上執行R代碼並將結果返回到JMP中。 嘗試R.Killick,P.Fearnhead&IAEckley(2012)Optimal Detection of Changepoints with Linear Computational Cost,Journal of the American Statistical Association,107:500,1590-1598中引用的Pruned Exact Linear Time(PELT)方法。
Re:重複分析謝謝凱文,任何建議表示讚賞。
根據我的理解,保羅沒有選擇任何測試,但更多的是要求如何在JMP中運行一系列雙樣本測試(一個樣本總是相同),給出他的原始數據集。我的觀點是在這種情況下不需要編寫腳本。
再次感謝您的貢獻。
Re:重複分析謝謝你的留言Kevin,
我不想打擾傑出的統計學家,無論他們是活著還是死了。
一世將會滿意至知道你的建議是什麼。
Re:重複分析有人應該提一下,你們所採取的方法可能會導致相當數量的死去的傑出統治者像他們各自的墳墓中的車床一樣旋轉。
我無法準確地看出你想從你的描述中做出什麼決定,但這聽起來更像是對我的變換點檢測。 如果是這樣,從任何雙樣本測試中分析許多順序聚合的p值可能被認為是最好的形式。
Re:重複分析嗨保羅,
謝謝,這幫了很多忙。
似乎腳本在這裡不是強制性的。您可以執行以下操作:
1。表>堆棧:獲取一個新數據表,其中包含來自控制列的數據和堆疊的所有測量列。這將創建一個包含列名稱的列標籤,以及包含所有數據的第二列。
附件顯示我的結果。
2。分析>適合Y乘X,X =標籤,Y =數據>熱點>比較均值>使用控制,Dunnett:選擇您的控制列,您將獲得所有成對測試。您可以右鍵單擊LSD閾值矩陣下的表,例如對其進行排序或將其轉換為數據表(可能用於在Graph Builder中繪製p值圖形)。
對於Dunnett的測試也見比較手段 。
希望有所幫助,
沃爾克
Re:重複分析嗨Volker,謝謝你照顧我的問題,
這是一個已經分析過的數據表,其中Ctrl / test比較顯示出顯著差異(在測試中評估,但Wilcoxon更合適)。
下面的每一行給出了Ctrl數據(控制列)和給定測試列(編號從100到200)之間的比較
Col nbr P值
...
160 0.01
161 0,0022
162 ns
163 ns
164 0,03
165 ns
166 0,008
167 0,04
168 0,04
169 ns
170 0,004
171 ns
172 ns
173 ns
174 0,04
175 ns
176 <0,0001
177 0,008
178 ns
179 0,004
180 0,02
181 0,003
182 0,004
183 0,009
184 0,04
185 0,005
186 ns
187 0,006
188 0,03
189 0,02
190 0,0002
191 ns
192 0,001
等等...
Re:重複分析嗨保羅,
這是我從聊天中學到的東西:
差異的定義:Wilcoxon測試將Ctrl列中的值與其他100列中的值進行比較(重複比較)。
示例表僅用於說明數據結構。上述行為不適用於此表。
目的是對每種分析表的差異的起始和偏移時間進行排序。
如果您可以提供另一個表格以及一些分析結果,您希望看到它們會有很大幫助。
BR,Volker