从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))
答案 0 :(得分:2)
您可以指定CI的方法有点奇怪,但一旦您知道如何就可以轻松完成。
请参阅函数crit<-
,其帮助页面说明了此示例的用法。
library(locfit)
data(ethanol)
fit <- locfit(NOx~E,data=ethanol)
crit(fit) <- crit(fit,cov=0.99)
plot(fit,band="local")