12

TL;DR

R の VGAM パッケージを使用して Tobit 回帰を実行しています。これは、診断できなかったエラーを一貫して表示しているおもちゃのデータセットです。

library(data.table)
library(VGAM)

> sessionInfo()$otherPkgs
$VGAM
Package: VGAM
Version: 0.9-7
Date: 2015-03-06
... <ommitted> ...

reg_data <- structure(list(S = c(1.83271488441825, 0.75411550370994, 0.904938604451928, 
                                 0.75411550370994, 0.75411550370994), H = c(0.6429, 0.7788, 
                                                                            0.6292, 0.8892, 0.2035), W= c(1.52497, 1.1391, 1.59722, 
                                                                                                          1.8406, 1.01865)), .Names = c("S", "H", "W"), class = c("data.table", 
                                                                                                                                                                  "data.frame"), row.names = c(NA, -5L))

minS <- 0.75411550370994
maxS <- 1.83271488441825

m <- vglm(S ~ H, tobit(Upper = maxS, Lower = minS), weights = W, data = reg_data)
Error in lm.wfit(x = cbind(x[!use.i11, ]), y = y[!use.i11, ii], w = w[!use.i11,  : 
  incompatible dimensions

診断の試み

トレースバックあり:

> traceback()
6: stop("incompatible dimensions")
5: lm.wfit(x = cbind(x[!use.i11, ]), y = y[!use.i11, ii], w = w[!use.i11, 
       ii])
4: eval(expr, envir, enclos)
3: eval(slot(family, "initialize"))
2: vglm.fitter(x = x, y = y, w = w, offset = offset, Xm2 = Xm2, 
       Ym2 = Ym2, etastart = etastart, mustart = mustart, coefstart = coefstart, 
       family = family, control = control, constraints = constraints, 
       criterion = control$criterion, extra = extra, qr.arg = qr.arg, 
       Terms = mt, function.name = function.name, ...)
1: vglm(y ~ x, tobit(Upper = maxy, Lower = miny), weights = w, data = X)

のソース コードを表示してlm.wfit、エラーの原因を見つけました。

function (x, y, w, offset = NULL, method = "qr", tol = 1e-07, 
    singular.ok = TRUE, ...) 
{
  <ommitted...>
    if (NROW(y) != n | length(w) != n) 
        stop("incompatible dimensions")
  <ommitted...>
  }

のソース コードで次の情報を見つけましたvglm。

    vglm.fitter <- get(method)
    fit <- vglm.fitter(x = x, y = y, w = w, offset = offset, 
        Xm2 = Xm2, Ym2 = Ym2, etastart = etastart, mustart = mustart, 
        coefstart = coefstart, family = family, control = control, 
        constraints = constraints, criterion = control$criterion, 
        extra = extra, qr.arg = qr.arg, Terms = mt, function.name = function.name, 
        ...)

メソッドのデフォルトはvglm.fitです。

use.i11除外基準が作成された場所、それが何をしているのか、なぜそれが重み、リグレッサー、リグレッサンドの間で競合するディメンションになるのか、私はまだ見つけることができません.

minSandmaxSを 10 桁以下に丸めると正常に実行されることがわかりましたが、これはmaxSが増加したため、最初の観測が右打ち切りではなくなり、minSが増加したため、2 番目、4 番目、5 番目の観測が打ち切られなくなりました。どちらも最尤関数での観測の扱いを変更するため、誤った結果で回帰を汚染していると思われます。

このタイプのエラーが発生する理由を誰かが診断するのを手伝ってくれますか?

4

1 に答える 1

7

パッケージの開発者から、これは実際にはバグであり、プレリリースされたパッケージhereで修正されているという知らせを受けました。これは、次の反復で、または彼の本がリリースされたときに、おそらく CRAN にアップグレードされるでしょう。

于 2015-04-17T12:58:15.117 に答える