6

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

NAS_DELAY=aSEANSON1TO4+bSEANSON2TO4+cSEANSON3TO4+d

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

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

理由がわかりません。

"fit2<-summary(rq(NAS_DELAY ~SEASON1TO4+SEASON2TO4+SEASON3TO4,tau=c(0.1,0.2,0.3,0.4,0.5),data=fddata))
Error in base::backsolve(r, x, k = k, upper.tri = upper.tri, transpose = transpose,  :   singular matrix in 'backsolve'. First zero in diagonal [1]"
In addition: Warning messages:
1: In rq.fit.br(x, y, tau = tau, ...) : Solution may be nonunique
2: In rq.fit.br(x, y, tau = tau, ...) : Solution may be nonunique
3: In rq.fit.br(x, y, tau = tau, ...) : Solution may be nonunique
4: In rq.fit.br(x, y, tau = tau, ...) : Solution may be nonunique
5: In rq.fit.br(x, y, tau = tau, ...) : Solution may be nonunique
6: In summary.rq(xi, ...) : 188771 non-positive fis

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

4

2 に答える 2

0

データセットで構成された因子のため、行列の構成は特異です。

特異点の説明を参照してください:ここ、ここ、およびここ。基本的に、quantreg はデータ行列の反転に依存しており、因子の形式のため、行列は反転できません。

十分なデータがある場合、またはデータにとって意味がある場合、このスレッドは、データに適している場合、いくつかの可能な解決策を示しています。

于 2015-05-07T01:05:46.163 に答える