2

重複の可能性：
高レベルのプロット関数、特にplot.rqs（）を使用した複数のプロット

1つの従属変数と1つの独立変数を使用して18分位点回帰を実行しています。$ \ tau = 0.01,0.02、...、0.99$にわたる18の勾配推定値の分布の$6 *3$タイルプロットが必要です。

これまでに私ができたことは次のとおりです。

require(quantreg)

op <- par(xlab="")
y=rnorm(200)
x=rnorm(200)
plot(summary(rq(y~x,tau=1:99/100)),parm=2,main="",mfrow=c(6,3),mar=c(0.5,2,2,2))
par(new=TRUE)
plot(summary(rq(y~x,tau=1:99/100)),parm=2,main="",mfrow=c(6,3),mar=c(0.5,2,2,2))

2番目のプロットは右下まで行きますが、3回目に試すと、グラフウィンドウにも表示されません。

par（mfrow = c（6,3））をplot（）の外に配置しようとしましたが、機能しませんでした。

ここに画像の説明を入力してください

4

2 に答える 2

3

plot.summary.rqsのデフォルトの動作は、パラメーターごとに1つのプロットを使用してタイルプロットを作成することです。これは、見慣れたものかもしれませんが、Rolandが言うように、努力を打ち負かす方法でpar（）を使用しています。 。これはいくつかの方法で克服できます。1つは、関数をやり直すことです。ただし、関数を変更しない場合は、par呼び出し（2つあります）をコメントアウトする補助関数（新しい名前の元の関数のコピー）を作成し、フレームワーク内での後続のプロットに使用できます。

ここに画像の説明を入力してください

于 2012-08-24T16:01:40.047 に答える
1

まず第一に、私はそれらの警告について心配します。

問題はplot.summary.rqs、プロットパラメータを設定およびリセットすることですが、これは良いことではありません。plot.summary.rqs1つの解決策は、この機能なしに基づいて独自のプロット関数を定義することです。

于 2012-08-24T07:46:52.260 に答える