兩個因素的所有配對的回歸預測表

Question

對於一項作業，我們要求為表中兩個因子變量的所有配對提供預測。 我有兩個因素和一個線性模型。 我想輸出一個表，使得行是一個因子的水平，列是另一因子的水平，表中的值是這些水平的預測值。

我已經很容易地完成了預測，但是顯示該預測會引起問題，因為似乎沒有任何功能可以很好地完成工作。 我以為tapply()可以做到，但我似乎無法將線性模型用作函數？ 有沒有辦法根據我的線性模型的預測結果輸出兩因素表？

newdata = expand.grid(GenericFactor1 = c("1", "2", "3"), GenericFactor2 = c("10","20","30","40"))
pred = predict.lm(linear.model.lm, newdata)
newdata$pred = as.vector(pred)

Answer 1

據我了解，您有兩個因素，並且想要創建一個表，其中行與第一個因素的值相關聯，列與第二個因素的值相關聯，數據是模型的預測用於相關的一對因子值。 outer功能可以很好地處理這種設置：

outer(c(1, 2, 3), c(10, 20, 30, 40), function(x, y) {
  predict.lm(linear.model.lm, newdata=data.frame(GenericFactor1=x, GenericFactor2=y))
})