Я хочу использовать функцию cluster.bs.glm (т. Е. Пары p-значений начальной загрузки кластера) из пакета clusterSEs в R для нескольких glm объектов.
В частности, я хочу реализовать его с помощью объекта glm.nb, который оценивался с помощью функции glm.nb из пакета MASS, но я столкнулся с ошибкой. Ниже я показываю код, который я использую.
cluster.bs.glm(neg_bin, dat = DATASET, cluster = ~ firm, boot.reps = 1000)
Error in `rownames<-`(`*tmp*`, value = c("(Intercept)", "treat", "post", :
length of 'dimnames' [1] not equal to array extent
Я подозреваю, что оценка тэты отбрасывает это. В конечном счете, я мог бы заставить это работать, используя glm, но сначала мне нужно было бы установить значение тета (не уверен, подходит ли это). Я постараюсь получить рабочий пример, если это будет необходимо. Какие-нибудь мысли? Может быть, удаление тета, чтобы доверительные интервалы могли выводиться на консоль?