Ошибка: внутренняя ошибка относительно длин векторов коэффициентов в R - PullRequest
Новая
       119

Ошибка: внутренняя ошибка относительно длин векторов коэффициентов в R

0 голосов
/ 07 мая 2020

Я пытаюсь построить график коэффициентов для регрессии 

est_bidders = feols(applied_participants ~ treatment_group::qt('18/02') | fe1 + fe2 + fe3, data)

coefplot(est_bidders)

И получаю странную ошибку: Error in coefplot_prms(object = object, ..., sd = sd, ci_low = ci_low, : Internal error regarding the lengths of vectors of coefficients

I на самом деле не понимаю, что это значит и как я могу это исправить. Когда я запускаю ту же регрессию для других переменных результата, она работает нормально, и в конце я получаю желаемые графики. Вероятно, это как-то связано с самой переменной (в основном она содержит числа от 1 до 86)

UPDATE. Основываясь на комментариях ниже, я использую dput () для воспроизведения: 

> dput(subset[1:15,])

structure(list(applied_participants = c(1L, 2L, 2L, 5L, 3L, 1L, 2L, 14L, 2L, 3L, 1L, 7L, 2L, 1L, 1L), treatment_group = c(1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0), qt = c("16/01", "16/02", "16/02", "16/02", "16/02", "16/02", "16/02", "16/02", "16/02", "16/03", "16/03", "16/03", "16/03", "16/03", "16/03"), good_class = c("27", "17", "22", "20", "43", "10", "45", "27", "43", "26", "22", "26", "26", "26", "26"), customerrequirement_customer_regnum = c("01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073", "01012000073"), year_month = c("2016-02", "2016-04", "2016-04", "2016-04", "2016-04", "2016-04", "2016-05", "2016-06", "2016-06", "2016-07", "2016-07", "2016-08", "2016-08", "2016-09", "2016-09")), row.names = c(NA, -15L), .internal.selfref = <pointer: 0x0000000002551ef0>, class = c("data.table", "data.frame"))

Я использую следующую библиотеку для вычисления коэффициентов регрессии и графика: 

install.packages('fixest')
library(fixest)

И в ситуации с этими точными переменными я использую регрессию: 

est_bidders = feols(applied_participants ~ treatment_group::qt('18/02') | good_class + customerrequirement_customer_regnum + year_month, subset)

Кажется, что регрессия работает нормально, но тогда коэфф-график не работает с указанной ошибкой coefplot(est_bidders)

1 Ответ

1 голос
/ 18 мая 2020

fixest автор здесь. Несколько комментариев:

  • treatment_group::qt('18/02') на самом деле сокращение от interact(treatment_group, qt, ref = "18/02"). Синтаксис variable::interaction очень нестандартен в R, я ввел его, потому что считаю его удобным, но он чрезвычайно сбивает с толку людей, не знакомых с ним (особенно, что в приведенном выше примере переменная name qt также является функцией R ...)
  • Вы обнаружили внутреннюю ошибку. Ошибка, с которой вы столкнулись, не должна никогда возникнуть, поэтому я кое-что не уловил. Не могли бы вы создать проблему на https://github.com/lrberge/fixest/issues с минимальным рабочим примером? Или просто отправьте мне электронное письмо? Поправлю в следующем обновлении пакета.
...