如何在 R 中使用 locfit 绘制 90% 置信带答案

【问题标题】：how to plot 90% confidence bands with locfit in R如何在 R 中使用 locfit 绘制 90% 置信带
【发布时间】：2014-10-16 20:57:41
【问题描述】：

我使用 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))

【问题讨论】：

我只看了scb 代码，但似乎 95% CI 只是作者写入函数的唯一选项。为什么？可能是因为这对他们来说似乎不是一个重要的功能。你的选择？可能自己编写或使用不同的本地回归工具。
@joran：我明白了。虽然即使从函数本身的来源来看，我也不清楚需要改变什么才能到达那里......

【解决方案1】：

指定 CI 的方法有点奇怪，但只要知道如何就可以轻松完成。

查看函数crit<-，其帮助页面说明了它在此示例中的用法。

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

【讨论】：

谢谢，乔希，这听起来确实比我目前所听到的要好得多（crit() 对 scb 对象崩溃 R 的天真应用）。不过，我不知道这是否应该适用于 scb。那不是一个单独的 locfit 对象。当我窥视它的源代码时，要找到我可能会尝试破解的 locfit 调用也不是那么容易。有什么办法？
嗨拉斯洛。看起来它比我想象的要复杂。 scb() 使用公式调用 locfit()，而后者又调用 locfit.raw()（作为其 lfproc= 参数的默认值。在 locfit.raw() 中，关键计算由编译后的 C 函数（.C("slocfit", ...)，使用通过Z$L 返回的下限和上限。因此（无需深入挖掘），看起来您可能需要更改 C 代码以获取不同的 CI（在此包中）。祝您好运！