問題タブ [quantreg]

Rのパッケージを使用しquantregて分位点回帰モデルを計算しました。このモデルでは、従属変数 (Y) はNAS_DELAYであり、独立変数 (X) はSEANSON1TO4、SEANSON2TO4、SEANSON3TO4です。
モデルは次のとおりです。

, SEANSON1TO4,はダミー変数でSEANSON2TO4、SEANSON3TO40 または 1 です。R を使用して、切片やその他の回帰係数を計算しましたが、結果は

"rq.fit.br(x,y,tau=tau,....)singular design matrix のエラー ;さらに: 警告メッセージ 1:in summary.rq(xi,....):278951 non-positivrフィス」。

理由がわかりません。

私は何を間違っていますか？

nlrq()私は（quantregパッケージから）非線形分位点回帰モデルを当てはめようとしています。しかし、関数が開始値に非常に敏感であることがわかりました。そこでセルフスターターSSexp()を使ってみました。セルフスターター関数が作成したことを読んだことnls()も、それで機能します。しかし、常に次のエラーメッセージが表示されます。

私の入力コードは次のとおりです。

getInitialセルフスターター機能を作成するために何らかの形で接続されている漠然とした考えがあります。しかし、率直に言って、どこに問題があるのか​​ わかりません。誰でもこれに関して私を助けることができますか?

追加情報

R サマリー オブジェクトから複数のプロットの各プロットのタイトルを変更しようとしています。私は五分位回帰を行い、係数をプロットしたいと思います。私は次のことをします

plot(qr.1) は、x1 と x2 の係数のプロットを提供します。ただし、各プロットのタイトルは x1 と x2 です。タイトルを x1 から変数名 (例: 一人当たりの収入) に変更したいと思います。いくつかの標準的なアプローチを試みましたが、失敗しました。これは、集計オブジェクトをプロットするためだと思います。要約オブジェクトから複数のプロットにタイトルを追加する方法を知っていますか?

ありがとう

勇道

そのため、Rのパッケージを使用してquantreg分位点回帰分析を実行し、予測子の効果が結果の分布全体でどのように変化するかをテストしています。

次にanova.rq、すべてのモデルで Wald 検定を実行し、共変量ごとに pvalue を出力して、結果の分布全体で各共変量の効果が大きく異なるかどうかを示します。

それはうまくいきます。ただし、私の特定のデータ (および一般的には?) については、推定値とその誤差をブートストラップすることが望ましいです。上記のコードを少し変更して実行します。

ここで行き詰まります。はquantreg現在、ブーストラップ推定値に対して anova (Wald) 検定を実行できません。パッケージの情報ファイルにはquantreg、ブーストストラップの方法について「anova.rq で使用するメソッドの拡張を行う必要がある」と具体的に記載されています。

anova.rq メソッドの詳細を見てみます。ブートストラップ時に分位モデルに存在しない 2 つのコンポーネントが必要であることがわかります。

1) Hinv(逆ヘッセ行列)。パッケージ情報ファイルには、「推定された共分散行列をそのサンドイッチ構成部分に分割する方法がないことに注意してくださいse = "boot"」と具体的に記載されています。

2)J情報ファイルによると、これは "と の場合に返される勾配行列のスケーリングされていない外積です。 フーバー サンドイッチはです。コンポーネントに関しては、 の場合はコンポーネントがありません。 (フーバー サンドイッチを作成するには、追加する必要があることに注意してください。マヨネーズはcov=TRUEse != "iid"cov = tau (1-tau) Hinv %*% J %*% HinvHinvJse == "boot"tau (1-tau)自分で。） 」

Hinvブートストラップされた推定値から計算または推定できますJか? そうでない場合、続行する最善の方法は何ですか？

これに関するヘルプは大歓迎です。ここに質問を投稿するのは初めてのタイミングですが、過去に他の人の質問への回答から大きな恩恵を受けました.