Следующий код
PredictNew <- predict (glm.fit, newdata = Predict, X1 =X1, Y1= Y1,
type = "response", se.fit = TRUE)
выдает 3 столбца data.frame
--PredictNew, подогнанные значения, стандартные ошибки и остаточный член масштаба.
Идеально ... Однако, используя модель, оснащенную zeroinfl {pscl}
:
PredictNew <- predict (zeroinfl.fit, newdata = Predict, X1 =X1, Y1= Y1,
type = "response", se.fit = TRUE)
или
PredictNew <- predict (zeroinfl.fit, newdata = Predict, X1 =X1, Y1= Y1,
type = "response", se.fit = TRUE, MC = 2500, conf = .95))
создать только один столбец с вектором значений Я, однако, очень хотел бы иметь стандартные ошибки. Все, что я прочитал, говорит, что они должны быть произведены ..
(Код несколько упрощен, у меня на самом деле четыре переменные и смещение - никаких проб с predict.glm
и se.fit = TRUE
производящих SE).
r
generalized-linear-model
count-data
zero-inflation
KalahariKev
источник
источник
predict()
функцию были встроены стандартные ошибкиzeroinfl()
.Ответы:
Насколько мне известно,
predict
метод для получения результатовzeroinfl
не включает стандартные ошибки. Если ваша цель - построить доверительные интервалы, одной из привлекательных альтернатив является использование начальной загрузки. Я считаю привлекательным, потому что у начальной загрузки есть потенциал, чтобы быть более устойчивым (с потерей эффективности, если все предположения для SE выполнены).Вот примерный код для того, чтобы делать то, что вы хотите. Это не будет работать точно, но, надеюсь, вы сможете внести необходимые исправления.
Я нарисовал этот код на двух написанных мной страницах: одна - параметры начальной загрузки из регрессии Пуассона с
zeroinfl
нулевым раздувом с нулевым раздувом Пуассона, а другая - как получить доверительные интервалы с начальной загрузкой для прогнозируемых значений из отрицательной биномиальной модели с нулевым усечением , Надеемся, что в совокупности это даст вам достаточное количество примеров, чтобы заставить его работать с прогнозируемыми значениями из нулевого раздувания Пуассона. Вы также можете получить некоторые графические идеи :)источник
Error in X.vlm.save %*% coefstart : non-conformable arguments
.