как построить 90% доверительные полосы с помощью locfit в R

Question

Перекрестная публикация от https://stats.stackexchange.com/questions/16811/how-to-plot-90-confidence-bands-with-locfit-in-r.

Я использую пакет locfit для некоторой локальной логистической регрессии в R и его функцию scb для построения одновременных доверительных интервалов.Однако я не вижу способа запросить 90% -ые доверительные интервалы вместо стандартных 95%.Я проверил сами функции, но не вижу, где их можно взломать.

Почему существует это ограничение?Что можно с этим сделать?

Мой код выглядит так:

fit2<-scb(closed_rule ~ bl,deg=1,type=4,xlim=c(0,1),ev=lfgrid(100), family='binomial',data=data,alpha=cbind(0,0.3))

Answer 1

Средство, с помощью которого вы можете указать CI, немного странно, но это можно легко сделать, если вы знаете, как.

См. Функцию crit<-, страница справки которой иллюстрирует ее использование в этом примере.

library(locfit)
data(ethanol)
fit <- locfit(NOx~E,data=ethanol)
crit(fit) <- crit(fit,cov=0.99)
plot(fit,band="local")