2017-07-16 3 views
0

lm()の見積もりが同じデータと数式の異なるgeom_smoothと異なるのはなぜですか?具体的には、私のグループ変数 "cat"レベル5の傾きは、lm()の出力で> 0で、geom_smoothでは< 0です(したがって、プロットはサマリーテーブルを反映していません)。Rで異なるgeom_smoothとlm()推定:ベーススロープとインタラクションスロープを一緒に追加するのを忘れました

ここにはthe dataです。 (同様に振る舞うサンプルデータを考え出すよりも簡単。)

モデル:x:cat5用傾斜> 0であるsummary(lm(data=df, y~x*cat))

注意。

プロット:

library(ggplot2) 
plt <- ggplot(df, aes(x=x, y=y, group=cat)) + 
    geom_smooth(method="lm", show.legend=FALSE) + 
    facet_wrap(~cat, nrow=1) + 
    geom_point(aes(color=color) 

(@Pedro Aphaloの答えhere次)geom_smoothの見積もりを取得:

library(ggpmisc)  
my.formula <- y~x 
plt + stat_poly_eq(formula = my.formula, 
      aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
      parse = TRUE) 

お知らせファセット5の傾きが< 0を使用してlm()geom_smooth Areはを正方形や何か別のもの?論文でどのバージョンを報告しますか?可能であれば、私はgeom_smoothのプロットとlm()の要約表を紙で使用することができるように2つを一致させたいと思います。ありがとう!

+0

データ(とモデル)の 'summary'の出力を追加して、データを正しく読み込んだことを確認できますか? – Spacedman

+0

@Spacedman、あなたの出力はあなたが正しくデータを読んだことを確認します。私は以下の質問を解決したので、質問を出力に更新するつもりはありません。 – BonnieM

答えて

3

それはすべて私に見えます。 CAT5の要約行は、次のとおりです。CAT5の相互作用:CAT5のためのスロープを意味

Coefficients: 
      Estimate Std. Error t value Pr(>|t|)  
(Intercept) 1.932248 0.053131 36.368 < 2e-16 *** 
x   -0.006651 0.001962 -3.389 0.000721 *** 
... 
cat5  -1.080554 0.075138 -14.381 < 2e-16 *** 
... 
x:cat5  0.005602 0.002775 2.019 0.043720 * 

は、xの全体的なスローププラスxの傾きである

> -0.006651+0.005602 
[1] -0.001049 

と私は見プロット上

- 0.00105

切片は2、だから、私の知る限り

> 1.932248+(-1.080554) 
[1] 0.851694 

である0.852として示されています物事は同意する。

+0

オハイオ州私のゴシック、@スペースマン。 (自己額を叩く!)私の簡単な間違いをキャッチしてくれてありがとう。私は相互作用勾配にx勾配を加えるのを忘れていました!どうもありがとう!!! – BonnieM

関連する問題