При использовании stargazer для создания таблицы LaTeX в объекте логистической регрессии стандартное поведение должно выводить логические значения для каждой модели. Возможно ли вместо этого получить exp (logit)? То есть, могу ли я получить коэффициенты Odds?
В документации stargazer ниже упоминается "Coef" -аргумент, но я не понимаю, может ли это включить exp (logits).
Coef: список числовых векторов, который заменит значение по умолчанию коэффи- циентов для каждой модели. Имена элементов будут использоваться для соответствия коэффициенты к отдельным ковариатам, и поэтому должны соответствовать ковариационные имена. NULL-вектор указывает, что для данной модели необходимо использовать набор коэффициентов по умолчанию. Напротив, вектор NA означает, что все коэффициенты моделей должны быть пустыми.
Ответ 1
Согласно симбиотическому комментарию в 2014 году, в более поздних версиях '' stargazer '' есть варианты ''. * '' для '' coef '' '' se '' '' t '' '' p '' и '' ci '', позволяющие прямое преобразование этих статистик.
apply.coef a function that will be applied to the coefficients.
apply.se a function that will be applied to the standard errors.
apply.t a function that will be applied to the test statistics.
apply.p a function that will be applied to the p-values.
apply.ci a function that will be applied to the lower and upper bounds of the confidence intervals.
Значение, которое вы можете напрямую использовать...
stargazer(model,
apply.coef = exp,
apply.se = exp)
EDIT: Я заметил, однако, что просто повышение уровня CI не дает того, чего вы ожидаете.
EDIT: вы можете получить правильные CI, используя описанный метод здесь.
Ответ 2
stargazer позволяет вам подменять множество вещей, зависимые метки переменных, ковариационные метки и т.д. Чтобы заменить те, которые вам нужно снабдить векторными метками переменных, это делается для того, чтобы по умолчанию использовать имена строк для публикации, а не имена переменных из R.
Итак, чтобы иметь коэффициенты шансов, вам нужно указать вектор коэффициентов шансов до stargazer. Как вы получаете этот вектор? На самом деле очень легко. Скажем, что ваша модель называется model, тогда ваш код:
coef.vector <- exp(model$coef)
stargazer(model,coef=list(coef.vector))
Если в вашей таблице несколько моделей, список следует развернуть, например. coef=list(coef.vector1,coef.vector2,...), где все векторы в списке будут получены из аналогичного возведения в степень, как указано выше.
Ответ 3
Итак, проблема в том, что вы хотите отобразить отношение (не-лог) шансов, но сохраняйте статистику теста на основе базовой линейной модели. По умолчанию, когда вы используете один из методов "apply", например apply.coef = exp, stargazer будет пересчитывать статистику t и значения p. Мы этого не хотим. Кроме того, стандартные ошибки находятся в базе данных, но мы не можем просто их проинформировать. Мой предпочтительный подход:
- Экспонизировать коэффициенты в stargazer
- отключите авто p и auto t
- отчет (нетрансформированный) t-статистика в таблице вместо стандартных ошибок
В коде это:
stargazer(model, apply.coef=exp, t.auto=F, p.auto=F, report = "vct*")
