2017-04-22 7 views
0

私は2つのロットフィットモデルを1つのプロットにプロットしようとしていますが、信頼区間をプロットするために2番目のロットフィットモデルを得ることはできません。私は、追加locfitモデルをプロットしようとすると、しかしR信頼区間を使って複数のロットフィットモデルをプロットする

plot(1_fit,band="local",type = "l", xlab = "Time", ylab = "Y-Axis",ylim=c(0,22), 
    col = "red",lwd = 5,font=3,main="Local Poly Fit 1",cex.lab=1.5, cex.axis=1.5, 
    cex.main=1.5, cex.sub=1.5) 

:それぞれが、以下のものを使用して信頼区間で自分の罰金にプロットすることができます

1_fit = locfit(Y~Time,data=data_1) 
2_fit = locfit(Y~Time,data=data_2) 

:私は2つのlocfitモデルを作成しました

lines(2_fit,col="blue") 

私はロフフィットラインのみを追加できますが、信頼区間は追加できません。私がやることを試みてきました:

lines(2_fit,band="local",col="blue") 

をしかし、私はこのメッセージなし信頼区間を得る:

Warning message: In plot.xy(xy.coords(x, y), type = type, ...) : "band" is not a graphical parameter

私もlines.locfitを使用してに見えたが、運がなかったしましたRは単に関数lines.locfitを見つけることができないと言っているからです。

私が使用して、同じウィンドウ内で両方のプロットを置くために周りの仕事を持っている:

par(mfrow=c(2,1)) 

しかし、彼らは同じプロット内であった場合、それはプロットがより比較になるだろうと避けたいです。

+0

いつでも 'par(new = TRUE)'を使用して、もう一方の上にプロットすることができます。一致するようにxlim/ylimを設定する必要があります。 '1_fit'はR –

+0

の正式な名前ではありません。代わりに、' predict'を使ってCIを抽出し、直接プロットしてください –

+0

par(new = TRUE)トリックでした!どうもありがとう!ちなみに、1_fitはちょうど一時的な名前であり、実際には自分のコードで使用していませんでしたが、とにかくそれが合法ではないことを知りました。 – user3221962

答えて

0

回答はRichard Telfordによってコメントで提供されました。

plot(1_fit,band="local",type = "l", xlab = "Time", ylab = "Y-Axis",ylim=c(0,22), xlim=c(0,12), col = "red",lwd = 5,font=3,main="Local Poly Fit 1",cex.lab=1.5, cex.axis=1.5, cex.main=1.5, cex.sub=1.5)` 

    par(new = TRUE) 

    plot(2_fit,band="local",type = "l", xlab = "Time", ylab = "Y-Axis",ylim=c(0,22), xlim=c(0,12),col = "blue",lwd = 5,font=3,main="Local Poly Fit 1",cex.lab=1.5, cex.axis=1.5, cex.main=1.5, cex.sub=1.5)` 

私はylimxlimは、ylab、およびxlab等しいだけでなく、メインだったことを確認するために必要な次のコードは、私が必要なものを達成することができました。リチャードのサイドノートも正当な名前ではありませんが、ここではプレースホルダの名前として使用しましたが、渡すための良い知識のようです。

関連する問題