幫助設計/分析重複測量方案你好,
我正在準備為我的研究方案設計實驗,並試圖確定一個模型來分析我的數據。以下是我的設計設定:
反應變項:心率
隨機因素:主題
固定因素:治療(對照、治療A、治療B)
目標:確定偏離基線心率的統計顯著因素/交互作用
協定:
- 治療前 60 分鐘記錄基線心率
- 在 20 分鐘治療期間記錄 20 分鐘心率
- 記錄治療後 60 分鐘的心率
我很困惑,因為文獻中使用了多種方法,而且網路上關於重複測量資料分析的建議似乎相互矛盾。我已經了解了使用廣義線性混合模型 (GLMM) 相對於重複測量變異數分析的主要好處。我只是對使用我的協議收集的數據的 GLMM 的實際應用有點困惑。
我正在考慮將記錄的心率平均到 5 分鐘的區間中,這將為每個受試者提供 12 個基線、4 個治療期間數據點和 12 個治療後數據點。沒有真正的韻律或理由為什麼我選擇了 5 分鐘。大多數文獻在分析之前都會按前 5 分鐘的區間 (%) 對受試者心率進行標準化,但我厭倦了這樣做,因為它不允許對受試者之間的心率水平進行比較。
現在,我有點陷入困境,不知道如何繼續或理解我擁有什麼樣的模型。
希望熟悉此類協議的人可以提供幫助。
任何幫助深表感謝。
謝謝,
J.P
進階統計建模消費者與市場研究回复:幫助設計/分析重複測量協議謝謝你的幫助!
回复:幫助設計/分析重複測量協議不,不是的。RE和重複結構代表兩種不同的變化來源。但是,對於某些重複結構(如UN),RE不再可估計 - 如果您嘗試這樣做,則會收到錯誤消息。AR(1)具有隨機主題效果。
時間順序對聯合國無關緊要。
回复:幫助設計/分析重複測量協議有兩種治療水平(A =低劑量,B =高劑量)。
基於Cholesterol Stacked.jmp示例,我認為我需要設置數據集和分析:
- 平台 :適合模特
- 個性 :混合模型
- Y,回應 :心率(bpm) - >連續
- 固定效果 :治療。,時區,Reg。研究所。- >所有分類
在重複結構選項卡下,如果我想使用'非結構化':
- 學科 : 學科
- 重複 :Time Reg./Inst。的毗連。- >分類
如果我想使用'AR(1)',請在重複結構選項卡下:
| 子。 | 對待。 | 時區 | 註冊。研究所。 | Time Reg./Inst。的毗連。 | 時間(分鐘) | 心率(bpm) |
| 1 | 一個 | 前處理 | 1 | 預Treatment_1 | 0 | 62 |
| 1 | 一個 | 前處理 | 2 | 預Treatment_2 | 五 | 60 |
| 1 | 一個 | 治療期間 | 1 | 在治療期間_1 | 10 | 64 |
| 1 | 一個 | 治療期間 | 2 | During_Treatment_2 | 15 | 67 |
| 1 | 一個 | 後處理 | 1 | 後Treatment_1 | 20 | 65 |
| 1 | 一個 | 後處理 | 2 | 後Treatment_2 | 25 | 64 |
| 2 | 乙 | 前處理 | 1 | 預Treatment_1 | 0 | 71 |
| 2 | 乙 | 前處理 | 2 | 預Treatment_2 | 五 | 72 |
| 2 | 乙 | 治療期間 | 1 | 在治療期間_1 | 10 | 78 |
| 2 | 乙 | 治療期間 | 2 | 治療期間_2 | 15 | 76 |
| 2 | 乙 | 後處理 | 1 | 後Treatment_1 | 20 | 76 |
| 2 | 乙 | 後處理 | 2 | 後Treatment_2 | 25 | 75 |
| ... | ... | ... | ... | ... | ... | ... |
“重複結構”主題字段是否構成我的“子”。列隨機效果還是需要在“隨機效果”選項卡中指定?此外,使用“非結構化”選項,如果“時間(分鐘)”列未用於該情況,模型如何知道每個時間點的順序?
謝謝!
回复:幫助設計/分析重複測量協議多少治療水平?對於時間組級別,它們應該像pre-1,pre_2,during_1,_2,...等。查看示例數據,膽固醇Stacked.jmp,用於JMP文檔的重複測量示例。該數據的格式和佈局與您的數據類似。
回复:幫助設計/分析重複測量協議謝謝,我一定會看看文字。
總而言之,我的堆疊數據的簡化示例版本如下所示:
| 子。 | 對待。 | 時間組 | 從第一個數據點開始的時間(分鐘) | 心率(bpm) |
| 1 | 一個 | 前處理 | 0 | 62 |
| 1 | 一個 | 前處理 | 五 | 60 |
| 1 | 一個 | 治療期間 | 10 | 64 |
| 1 | 一個 | 治療期間 | 15 | 67 |
| 1 | 一個 | 後處理 | 20 | 65 |
| 1 | 一個 | 後處理 | 25 | 64 |
| 2 | 乙 | 前處理 | 0 | 71 |
| 2 | 乙 | 前處理 | 五 | 72 |
| 2 | 乙 | 治療期間 | 10 | 78 |
| 2 | 乙 | 治療期間 | 15 | 76 |
| 2 | 乙 | 後處理 | 20 | 76 |
| 2 | 乙 | 後處理 | 25 | 75 |
| ... | ... | ... | ... | ... |
回复:幫助設計/分析重複測量協議您需要創建一個時間列(例如,0,10,20,30等)。 與主題ID列一起,兩者唯一地標識每個心率測量值。指定重複結構時,需要這兩個ID列。(注意:對於AR(1),時間列應定義為連續。對於其他人來說,這是絕對的。請參閱“重複協方差結構要求”部分 http://www.jmp.com/support/help/13/Fit_Model_Launch_Window_2.shtml#999895 。
您還需要創建一個單獨的三級治療柱,以用作固定效果。
最後,除了JMP文檔,我強烈推薦這本書 https://www.sas.com/storefront/aux/en/spmixedmodel/59882_toc.pdf 。第5章是關於重複措施的分析。
回复:幫助設計/分析重複測量協議由於30個響應測量中的每一個實際上是每個5分鐘窗口的摘要(平均值),我可以將其濃縮為總結為10分鐘窗口以減少數據點的數量,我想我將研究Exchange和AR1結構。
作為一個更尖銳的問題,我是否需要為“時間組”創建一個列,其中有三個可能的級別:“預處理”,“治療期間”和“治療後”,並將其用作我的混合模特?
我問的原因是因為這些類型的研究通常在施用治療之前具有零時間響應數據點,然後所有數據點被認為是治療後的時間。在我的情況下,我的治療在停止前連續給藥30分鐘,因此還需要定義“治療期間”。
謝謝!
回复:幫助設計/分析重複測量協議1。 您可以將所有測量結果轉換為“時間為零”的差異(或%變化)。是否轉換應部分由您的學習目標指導。
2。重複測量的信息用於估計您將在模型中包含的治療效果以及重複結構。需要注意的是,由於您已經記錄了30個響應點,因此可能會有太多的方差 - 協方差參數來估算UN等結構。 我建議你從Exchangeable或AR(1)開始。
回复:幫助設計/分析重複測量協議你好,健,
謝謝您的回复。在回顧這個討論之前,我剛剛完成了對重複測量教程的JMP Pro混合模型的觀察。他們非常有幫助。
作為後續行動,我有兩個問題:
1)最好將響應變量(心率)保持其原始形式以執行此混合模型分析,或者可以將心率(受試者內部隨時間變化)標準化為從零時間或差異的百分比變化從運行分析之前的零時間開始?我的直覺告訴我,正常化可能不是一個好主意,因為它不允許進行主體間比較,但我很好奇你的意見。
2)在我的情況下,每個受試者在三個有序時間組內及時測量30個響應數據點,即:輸注治療前(12個數據點),輸注治療期間(6個數據點)和輸注後治療(12個數據點)。如何將此信息添加到混合模型平台?
謝謝!
回复:幫助設計/分析重複測量協議鑑於你的反應變量,心率,是連續的,我會用線性混合模型開始分析,這在JMP Pro中可用( 分析=>擬合模型並選擇混合模型如個性 )。數據必須採用高格式,即每個觀察都由“主題”和“時間”列標識。
在JMP混合模型中,可以使用“隨機效果”選項卡指定隨機主題效果,而重複測量之間的相關性由“重複結構”選項卡處理。當然,使用“固定效果”選項卡指定固定效果,例如治療等。
使用JMP Pro 13,有許多不同的方差 - 協方差結構可供選擇。在決定“最佳”結構之前,您可以嘗試使用幾個並評估模型。有關詳細信息,請參閱文檔 http://www.jmp.com/support/help/13/Fit_Model_Launch_Window_2.shtml#999895