Неточные значения р в Stargazer
Я хочу, чтобы в значении регрессии в Stargazer были те же звезды, что и в «нормальном выходе».
Я производю данные
library("stargazer"); library("lmtest"); library("sandwich")
set.seed(1234)
df <- data.frame(y=1001:1100)
df$x <- c(1:70,-100:-71) + rnorm(100, 0, 74.8)
model <- lm(log(y) ~ x, data=df)
и получить некоторые модельные оценки, где коэффициент на х имеет значение р 0,1023
coeftest(model, vcov = vcovHC(model, type="HC3"))
Я хочу, чтобы эти результаты были в LaTeX. Основываясь на той же функции, я вычисляю непротиворечивые стандартные оценки гетероскедастичности и позволяю stargazer использовать их.
stderr_HC3_model <- sqrt(diag(vcovHC(model, type = "HC3")))
stargazer(model, se=list(stderr_HC3_model))
На выходе звездного наблюдателя есть звезда с коэффициентом, указывающим на значимость, когда альфа = 10%. Я хочу, чтобы Stargazer дал столько же, сколько и самый крутой. (Из-за сравнимости со Stata, где reg L_y x, vce (hc3) дает именно самые лучшие результаты.)
Я поиграл с опциями Stargazer p.auto, t.auto, которые не помогли. Когда я выполняю «stargazer», я не могу просмотреть основной код, как это возможно в других случаях. Что делать?
Ответ Ричардса помог мне. Я указываю шаги, которые я использовал, чтобы выдать более одной регрессии (скажем, ols_a и ols_b).
ses <- list(coeftest(ols_a, vcov = vcovHC(ols_a, type="HC3"))[,2],
coeftest(ols_b, vcov = vcovHC(ols_b, type="HC3"))[,2])
pvals <- list(coeftest(ols_a, vcov = vcovHC(ols_a, type="HC3"))[,4],
coeftest(ols_b, vcov = vcovHC(ols_b, type="HC3"))[,4])
stargazer(ols_a, ols_b, type="text", p=pvals, se=ses)