背景:我有6個人,每個人分配到治療組或對照組(每組3個)。我想比較治療效果,並將“個體”作為混合模型中的隨機效應和“治療”(T或C)作為固定效應。
問題:我發現根據我如何標記我的複製品(個體)得到不同的結果。在實驗中,我將個體組織為1,2,3,在治療組中,4,5,6,在對照組中。這就是我最初在我的混合模型中包含“個體”的方式,但我一直被告知這是不正確的並且使治療效果混淆,因為個體在治療期間沒有重複,並且每個人僅被分配一次治療。相反,我需要在處理的重複結構中具有標稱值(C中的1,2,3和T中的1,2,3)。取決於我如何構建“個體”隨機複製效果:(T-1,2,3和C-1,2,3)與(T-1,2,3和C-4,5,6,)
問題:我不完全理解為什麼治療中的重複結構是正確的,我正在尋求進一步的澄清。
聽起來你有一些不好的建議。 除非您將受試者置於治療範圍內,否則JMP會將數據視為您在對照組和治療組中均有3名受試者進行治療。 實際上,您有6個科目只分配到1組。 這就是為什麼你會遇到不匹配的問題,而且據我所知,你第一次就做對了。
如果使用第二種標記方法,則在處理中嵌入主題ID時應獲得相同的結果。 我附上一個例子來說明我為假研究生成了一些數據,在0和1週時測量了每個治療水平的4名受試者。 我創建了2個主題ID列,以說明我可以這樣做,如果我正確嵌套它可以得到相同的結果。 第一個腳本為每個主題使用唯一ID。 第二個使用為每個治療組重複的名義ID,但我將治療組中的主題嵌套。 我得到完全相同的數字。 這是因為嵌套告訴JMP T組中的Subject 1與C組中的Subject 1不同。
第三個腳本顯示沒有嵌套的錯誤方法。 現在結果是不同的,因為JMP將兩組中的受試者1視為同一個人,並認為我對此人進行了4次測量。 這基本上就是我認為您對數據所做的。
我還想詢問每個人是否有超過1個數據點?例如,在我上傳的數據集中,我在2個時間段測量它們。 如果您只有1個數據點,則主題與殘差混淆(未按照您的建議進行處理)。 如果是這樣的話,那不是一個真正的問題。 JMP只是從REML方差組件表中刪除“Residual”並將其替換為“Subject”。您的治療效果仍然完全可以估計。