如何繪製通過模型估計R

Question

我想用R來做一些建模，我已經開始使用BodyWeight庫，因爲我已經看到一些在線的例子。只是爲了理解和習慣這些命令。

我來我的最終模型，估計和我不知道如何繪製這些估計，但我還沒有在網上見過..

是否有繪製方式用線估計的值，以及每個觀察值的點數？

從哪裏可以找到關於如何做到這一點的信息，我必須自己提取值還是可以說繪製這些模型的估計值？

我只是以R開頭。任何幫助都是值得歡迎的。

謝謝

Answer 1

沒有的功能，只是一個繪製模型的輸出，因爲通常有密謀輸出的許多不同的可能途徑。

看看predict函數適用於您正在使用的任何模型類型（例如，使用lm具有predict.lm函數的線性迴歸）。

然後選擇一個繪圖系統（你可能會想要不同的面板爲不同的飲食水平，所以使用ggplot2或格）。然後看看你能否用文字更清楚地描述你想要情節。然後更新你的問題，如果你卡住了。

---

現在，我們已經確定了哪個數據集使用的是，這裏有一個情節可能：

#Run your model 
model <- lme(weight ~ Time + Diet, BodyWeight, ~ 1 | Rat) 
summary(model) 

#Predict the values 
#predict.lme is a pain because you have to specify which rat 
#you are interested in, but we don't want that 
#manually predicting things instead 
times <- seq.int(0, 65, 0.1) 
mcf <- model$coefficients$fixed 
predicted <- 
    mcf["(Intercept)"] + 
    rep.int(mcf["Time"] * times, nlevels(BodyWeight$Diet)) + 
    rep(c(0, mcf["Diet2"], mcf["Diet3"]), each = length(times)) 
prediction_data <- data.frame(
    weight = predicted, 
    Time = rep.int(times, nlevels(BodyWeight$Diet)), 
    Diet = rep(levels(BodyWeight$Diet), each = length(times)) 
) 

#Draw the plot (using ggplot2) 
(p <- ggplot(BodyWeight, aes(Time, weight, colour = Diet)) + 
    geom_point() + 
    geom_line(data = prediction_data) 
)