2015-04-21 2 views
0

Я пытался извлечь мою сводную статистику из mlogit модели пробегает Zelig с использованием texreg, stargazer и memisc.LATEX Опция Сводной таблицы для Зелиг MLogit Результаты

texreg бросает мне следующую ошибку:

texreg(MLogitRes3) 
Error in (function (classes, fdef, mtable) : unable to find an  
inherited method for function ‘extract’ for signature ‘"vglm"’ 

stargazer кидает мне следующую ошибку:

stargazer(MLogitRes3) 
Error in objects[[i]]$zelig.call : 
$ operator not defined for this S4 class 

memisc (с помощью функции mtable) бросает мне эту последнюю ошибку:

mtable(MLogitRes3) 
Error in UseMethod("getSummary") : 
no applicable method for 'getSummary' applied to an object of class 
"c('vglm', 'vlm', 'vlmsmall')" 

Ни один из эти пакеты поддерживают модели mlogit в пределах zelig? Есть ли у меня варианты экспорта моей сводной статистики в таблицу, используемую в LaTex в другом месте?

ответ

1

Вы пытались воспроизвести это недавно? Я просто проверил код texreg::extract.zelig и это, кажется, есть метод для mlogit:

function (model, include.aic = TRUE, include.bic = TRUE, include.loglik = TRUE, 
    include.deviance = TRUE, include.nobs = TRUE, include.rsquared = TRUE, 
    include.adjrs = TRUE, include.fstatistic = TRUE, ...) 
    ... 
    else if ("mlogit" %in% class(model)) { 
     coefficient.names <- rownames([email protected]) 
     coefficients <- [email protected][, 1] 
     standard.errors <- [email protected][, 2] 
     zval <- [email protected][, 3] 
     significance <- 2 * pnorm(abs(zval), lower.tail = FALSE) 
     gof <- numeric() 
     gof.names <- character() 
     gof.decimal <- logical() 
     if (include.loglik == TRUE) { 
      lik <- logLik(model)[1] 
      gof <- c(gof, lik) 
      gof.names <- c(gof.names, "Log Likelihood") 
      gof.decimal <- c(gof.decimal, TRUE) 
     } 
     if (include.deviance == TRUE) { 
      dev <- deviance(s) 
      if (!is.null(dev)) { 
       gof <- c(gof, dev) 
       gof.names <- c(gof.names, "Deviance") 
       gof.decimal <- c(gof.decimal, TRUE) 
      } 
     } 
     if (include.nobs == TRUE) { 
      n <- nrow(model$data) 
      gof <- c(gof, n) 
      gof.names <- c(gof.names, "Num. obs.") 
      gof.decimal <- c(gof.decimal, FALSE) 
     } 
     tr <- createTexreg(coef.names = coefficient.names, coef = coefficients, 
      se = standard.errors, pvalues = significance, gof.names = gof.names, 
      gof = gof, gof.decimal = gof.decimal) 
     return(tr) 
    } 
    ... 
} 

Во всяком случае, если вы все еще возникают проблемы, вы можете прочитать Section 6 of the texreg article, который дает вам инструкции о том, как определить ваш собственный метод extract для любой модели и/или написать автору пакета (Philip Leifeld), чтобы добавить поддержку.

Смежные вопросы