R интервал предсказания для среднего значения нового образца

Question

Учитывая модель регрессии, созданную из одного набора данных, я использовал WinBUGS для построения интервалов прогнозирования (PI) вокруг среднего значения второго набора данных. Я только что открыл функцию «предсказания» в R, но он поставляет PI вокруг каждого прогнозируемого значения во втором наборе данных. Я искал помощь R, здесь и в Сети, и только нашел интервалы для отдельных членов.

Среднее значение этих интервалов явно не совпадает с PI вокруг прогнозируемого среднего значения выборки (и я проверил это по отношению к значению, полученному от WinBUGS).

Как мне получить R, чтобы дать мне PI вокруг среднего?

Answer 1

Раньше была функция R mean.data.frame, но она была устаревшей, а затем удалена. Вы можете получить тот же результат с:

mean.vec <- lapply(na.omit(dfrm), mean) 

Тогда, вероятно, просто:

predict(fit, newdata=data.frame(mean.vec)) 

Я говорю «наверное», потому что вы не представили набор данных, чтобы проверить это с и предоставление таких, на мой взгляд, ваша ответственность , Я не знаю, реплицирует ли этот метод JMP или метод WinBUGS.