私はいくつかの助けを使うことができます:
abline
を私のプロットに追加する「Control」と呼ばれる2つのレベルの因子変数があり、「consumption_mean」列に対応するデータポイントがあります。
ここで2つをプロットしますが、を追加しようとすると、abline
探している結果がまったく得られません(がcol=
機能していないようです。赤に設定すると、黒になります。
どんな助けでも大歓迎です
赤いデータポイント(Control == 1の場合)を通る赤い回帰直線と、黒いデータポイント(Control == 0の場合)を通る黒い回帰直線が欲しいのですが。
df_mean <- read.table("https://pastebin.com/raw/C2dz6bDT", header = T)
df_mean$Control <- as.factor(df_mean$Control)
plot(df_mean$mean_consumption~df_mean$day, type="n")
points(df_mean$day[df_mean$Control=='1'],df_mean$mean_consumption[df_mean$Control=='1'],col='red',pch=15,cex=.8)
points(df_mean$day[df_mean$Control=='0'],df_mean$mean_consumption[df_mean$Control=='0'],col='black',pch=15,cex=.8)
abline(lm(df_mean$mean_consumption~df_mean$day+(df_mean$Control=='1'),lwd=1,lty=1,col="red"))
abline(lm(df_mean$mean_consumption~df_mean$day+(df_mean$Control=='0'),lwd=1,lty=1,col="black"))
線が赤ではないだけでなく、Control == 1に設定されている最初の回帰は、実際には赤の点にまったく適合しないことに注意してください。
ここにいくつかあります、
1)括弧の1つを移動することで修正できる代わりに、col
およびその他の引数を入力したため、行が赤くなりません。lm
abline
abline(lm(df_mean$mean_consumption~df_mean$day+(df_mean$Control=='1')),lwd=1,lty=1,col="red")
2)異なる「式」を使用してプロットしているため、線とデータポイントが整列していません。コードは次のようになります。
plot(x,y)
abline(lm(y ~ x))
あなたの場合、
x
= df_mean$day[df_mean$Control=='1']
; および
y
=df_mean$mean_consumption[df_mean$Control=='1']
だからあなたabline
は
abline(lm(df_mean$mean_consumption[df_mean$Control=='1']~df_mean$day[df_mean$Control=='1']),lwd=1,lty=1,col="red")
すべての問題を修正すると、目的のプロットが得られます。
この記事はインターネットから収集されたものであり、転載の際にはソースを示してください。
侵害の場合は、連絡してください[email protected]
コメントを追加