stargazer с выводом таблицы R или для "специальной" сводной статистики

Question

У меня есть категориальная переменная, df$Chain, для уровней 1,2,3 и 4. У меня есть переменная, df$State, для уровней 0, 1. Я хотел бы создать таблицу, показывающую пропорции по штатам. Стандартное отклонение пропорций и t-статистика теста на равенство двух пропорций.

Я могу произвести пропорции через

table.1 <- prop.table(table(dfLong$CHAIN,dfLong$STATE), margin = 2)

который дает

            0         1
  1 0.4430380 0.4108761
  2 0.1518987 0.2054381
  3 0.2151899 0.2477341
  4 0.1898734 0.1359517

и идет в правильном направлении.

Но когда я передаю этот объект Stargazer через

stargazer(table.1, type = 'text')

это складывает пропорции

  0 1  NA  
 -----------
 1 1 0 0.443
 2 2 0 0.152
 3 3 0 0.215
 4 4 0 0.190
 5 1 1 0.411
 6 2 1 0.205
 7 3 1 0.248
 8 4 1 0.136
 -----------

Как я могу получить пропорции в состоянии рядом друг с другом и добавить стандартные ошибки и t-статистику?

Является ли Stargazer правильным инструментом для этого вообще?

Answer 1

Предположим, ваши данные и prop.table является

set.seed(123)
dfLong <- data.frame(CHAIN = sample(1:4, 100, replace = TRUE),
                     STATE = sample(0:1, 100, replace = TRUE))
table.1 <- prop.table(table(dfLong$CHAIN,dfLong$STATE), margin = 2)

Идея состоит в том, чтобы определить class за table.1 как matrix и проинструктировать stargazer распечатать матрицу "как есть":

class(table.1) <- "matrix"
stargazer(table.1, type = "text", summary = FALSE)
# =============
#     0     1  
# -------------
# 1 0.180 0.340
# 2 0.320 0.220
# 3 0.220 0.180
# 4 0.280 0.260
# -------------

Используя этот подход, вы можете cbind другие значения, представляющие интерес (стандартные ошибки и т. д.) для исходной матрицы, и выполните расширенную матрицу через stargazer,