モデルを使って見積もりをする方法R

Question

Rを使ってモデリングをしようとしていますが、私はいくつかの例をオンラインで見たので、BodyWeightライブラリを使い始めました。理解し、コマンドに慣れるだけです。

私は推計で、私の最終的なモデルに来ていると私はこれらの推定値をプロットする方法を思っていたが、私はオンライン何かを見ていない...

は、プロットする方法はあります線での推定値、および各観測値に対する点は？

ここでどのように情報を見つけることができますか、自分で値を抽出する必要がありますか、これらのモデルの推定値をプロットすることは可能ですか？

私はRで始まっています。何か助けを歓迎します。

は出力をプロットする多くの異なる可能な方法は、通常はあるのでただ、モデルの出力をプロット機能がないあなた

Answer 1

ありがとうございます。

使用しているモデルタイプ（たとえば、lmを使用した線形回帰は、predict.lmの関数）を使用して、predict関数を見てください。

次に、プロットシステムを選択します（食事のレベルによってパネルが異なる可能性がありますので、ggplot2またはlatticeを使用してください）。次に、プロットをどのように見たいかを、より明確に記述できるかどうかを確認します。あなたが立ち往生したら、あなたの質問を更新してください。今、私たちは、あなたが使用しているデータセットが同定されてきた

---

、ここでの可能性プロットだ：

#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) 
)