問題タブ [quantreg]

quantreglapack にリンクしていないため、パッケージのインストールに失敗します:

[をちょきちょきと切る]

[をちょきちょきと切る]

奇妙なことに、新しいユーザーをゼロから作成すると、すべてが完全に機能します (ログに「-llapack」が記録されます)。

いくつかの環境変数がこれに関与していると思われますが、どれを特定できませんでした: LIBRARY_PATH と LD_LIBRARY_PATH を設定解除すると、これはまだ発生します。アイデアをありがとう！

quantregWin7 64bit と Linux Mint Maya 64bit の両方でライブラリを使用していますが、速度に大きな違いがあることに気付きました。

私は Fortran で書かれた Frisch-Newton メソッドを使用しており、.Fortran() 呼び出しの速度のみを測定していました。

長さ1600のベクトルを分析しています。

Windows では約 2.5 秒、Linux では約 22 秒かかります。そして、私は本当に理由を理解していません..

この後、Fortran コードを F2C で C コードに変換し、clapack 関数を使用して動的ライブラリにコンパイルしようとしました。この関数を .C() で呼び出していて、時間効率を再度測定しました。

Windows では 8 秒、Linux では 22.5 秒かかりました。私は本当に混乱しています。

すべてが I7-3770、16GB RAM で行われました。

なぜそれが機能しないのか、誰かが私に説明できますか?

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

Fortran コード:

C コード:

R からの呼び出し:

動的ライブラリのコンパイル:

パッケージを使用してquantreg、R で次の分位点回帰を実行しています。

を介して係数と要約統計量を抽出します。

私が得た結果は、

ただし、R^2/調整済み R^2 を知りたいです。これは、summary()-コマンドが単純な OLS 回帰に対して提供しているように見えますが、分位回帰の場合はそうではありません。

それらを抽出する方法を知っている人はいますか？

パッケージrq()内の関数を使用して分位点回帰を実行しようとしています。quantreg私は自分のデータをセットアップし、モデルなどを作成しました (別のサイトの指示に従います)。下記参照; REとWTPCは連続変数でArmoredあり、以下で定義されます。

summary()ただし、rq()モデルを呼び出すと、これが得られます

これは素晴らしい情報ですが、通常の線形回帰のように標準誤差、p 値などが必要です。ドキュメントと他の例を確認しましたが、これがこれを行う方法のようです。モデル ( QR25) を作成してから呼び出すsummary(QR25)と、エラーと p 値を含む回帰出力が得られます。

これは分位点回帰の最初の試みであるため、以前のステップが欠落しているか、何か他のものを指定する必要がある可能性があります。何かご意見は？これを RStudio バージョン 0.98.501 で実行しています。

quantregパッケージを使用して分位点回帰から係数と上限と下限を抽出したいと思います。ヘルプ ファイルの例を次に示します。

coefficients()係数を取得するために使用できることはわかっています。

しかし、に表示される上限/下限を取得する方法がわかりませんsummary()。を見ましたstr(sf)が、抽出方法がわかりませんでした。

最終的には、タウス、係数、上限/下限をデータフレームに入れてさらに処理したいと考えています。

現在、Roger Koenker のquantregパッケージを使用して、R で分位点回帰モデルを実行しています。5 つのレベルを持つ順序付けされたカテゴリ変数と、それぞれ 7 つ、4 つ、および 2 つのカテゴリを持つ 3 つの順序付けされていないカテゴリ変数 (曜日、季節、夜間かどうか) があります。1 から 5 までの 5 つのレベルである順序付けられたカテゴリ変数を除いて、すべてが最初は文字列としてデータ フレームに格納されます。私の変数が呼び出されるdf_primeことに注意してください。これは、一部の行が削除された元のデータ フレームであるためです。

回帰を実行すると、金曜日がサマリーから完全に削除され、順序付き変数である acuity_id のレベルの 1 つが削除されます。それはおそらくRがダミー変数を作成する方法に関連していますか? 何らかの理由で、seasonWinter の p 値が 1 になっていることに気付きました。

私の出力は次のとおりですが、職場の VM に R があり、コピーして貼り付けることができないため、間隔が異なります。

このquantregパッケージを使用して、R で分位点回帰を計算しています。次のコマンドを使用して QR を呼び出します。

戻り値を見ると、変数 flightHours と flightHoursType の係数はゼロ/NULL です。

追加するse = "nid"と、警告メッセージも表示されます。

警告メッセージ: summary.rq(quantGsReg15, se = "nid") : 7 非正の fis

以前に計算した OLS 回帰は、そうではない可能性が高いことを示しています。さまざまな QR メソッド ( method = "fn") を試しましたが、成功しませんでした。

0.5 分位数に対して同じコマンドを呼び出すと、同様の出力が返されます。ただし、別の警告メッセージが表示されます。

警告メッセージ: rq.fit.br(x, y, tau = tau, ci = TRUE, ...) : 解が一意でない可能性があります

私は何を間違っていますか？私は今30分間グーグルで検索しましたが、成功しませんでした。

通常の OLS 係数とその CI を使用して係数グラフを信頼区間 (CI) でプロットしようとすると、原点を介して回帰を強制するたびにエラーが発生することに気付きました。

したがって、このコードを使用すると (engel は R の分位点回帰の例のデータです):

問題はなく、coefficeint グラフが描画されます。しかし、これを使用すると：

インターセプトが原点を通過するため、問題があります。分位点回帰の最初のコードのように、切片なしでプロットを取得するにはどうすればよいですか。