嵌套和交叉R的anova模型的代碼

我要構造/測試的模型是：因變量=因子A +因子B +因子C因子A和C之間的相互作用因子B和C +因子之間的相互作用B嵌套在A因子中嵌套和交叉R的anova模型的代碼

我在網上遇到的一個示例在文件「ANOVA：高級設計」（http://web.grinnell.edu/individuals/kuipers/stat2labs/Handouts/DOE%20Advancede.pdf）（感謝作者在線共享此文件）中進行了描述。在這個文件中，分割繪圖/重複測量設計（幻燈片9-10）中描述的例子與我的情況類似。這裏，因子A是品牌，因子B是盒子，因子C是溫度。如果我們假設（1）盒子是一個固定的效果（即這3個盒子代表所有可能的盒子等級），（2）每個盒子中的所有袋子都被分配到一個溫度，以及（3）兩個溫度級別（例如，有四個溫度級別，10,20,30,40），並且每個箱子內分配給特定溫度的袋子數量是隨機確定的（即，分配給不同溫度的袋子數量不相等這可能是因爲某些箱子裏沒有袋子被分配到一定的溫度），那麼這個例子幾乎和我想描述的一樣。另外，我的設計不平衡。

我想測試哪些因素以及這些因素如何影響因變量。假設是3路（以爆米花，品牌，溫度，盒子爲例）anova的假設。在爆米花的例子中，零假設可能是：品牌，溫度和/或盒子不會影響彈出的內核％。備選假設與null相反。另外，在我的情況下，可能盒子也可能是一個隨機效果，就像盒子一樣，但我想將這兩種情況考慮在內（盒子是固定的和隨機的效果）。

解決此問題的適當方法是什麼？

謝謝。

來源

2016-01-04 D. Dai

這是一個非常普遍的問題，並丟失了一些重要信息。我會建議首先確定統計數據，然後代碼將很容易搜索。雖然您發佈了有關線性預測變量的一些信息，但要指定模型，您需要考慮需要實現哪個鏈接函數，這取決於依賴數據的結構。當提問時，它也有助於發佈示例代碼和數據。 –

感謝您的建議。我會嘗試編輯我的問題。 –

我編輯了這個問題，並提供了一個我在網上找到的例子。 –

我不是100％肯定我們同意對術語，但我會採取了一槍......

你說要

因子A + B因子+因子C +互動嵌套倍之內的因素B和C + B因子之間因素A和C +互動之間的

要注意的主要事情是，「B嵌套內的」等同，因爲我是這個世界至少熟悉，以「包括A的主要作用和A之間的相互作用nd B，但不是B的主要作用「（即， ~A/B == ~A+A:B。但那你說你做想要B因子的主效應，所以這個看起來有點奇怪。下面您的規格完全相同會給

~ A + B + C + A:C + B:C + A/B

但這相當於

~ A + B + C + A:C + B:C + A + A:B

[R自動丟棄多餘A項，所以這也相當於

~ A + B + C + A:C + B:C + A:B

但由於這是基本上是主效應加上所有雙向交互，你也可以把它寫成

~(A+B+C)^2

因爲多餘的條款將被丟棄，你能等效於許多不同的方式寫：~A*B+A*C+B*C（A*B相當於A+B+A:B）或~A*C+B*C+A/B ...如果你想檢查什麼R在實際生產中，可以使用colnames(model.matrix(my_formula,my_data))。

這一切都是假設我們在lm() /固定效果方面的工作...

來源

2016-01-04 19:18:27

非常感謝您的回答。這對我來說是非常有幫助的！關於因素B，它不一定包含在模型中，因爲它的影響可能實際上是由因子C引起的。同樣，因素B和C似乎是相互影響的。 –

@Ben，你這篇文章對我來說聽起來是錯誤的：「包括A的主要作用和A與B之間的相互作用，但不是A的主要作用」，你的意思是「包括A的主要作用和A和B之間的相互作用，但不是主要影響** B **「 – Mattma

謝謝，修正。 –

嵌套和交叉R的anova模型的代碼

回答

相關問題