如何在R中用locfit绘制90%置信区间_R_Statistics

如何在R中用locfit绘制90%置信区间

r statistics

如何在R中用locfit绘制90%置信区间,r,statistics,R,Statistics,交叉张贴从我使用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=cb

交叉张贴从

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

fit2指定CI的方法有点奇怪，但一旦知道如何指定，就可以很容易地完成
参见函数critI仅浏览了scb
代码，但似乎95%的CI只是作者写入函数的唯一选项。为什么？可能是因为对他们来说，这似乎不是一项重要的功能。你有选择吗？你自己写或者使用不同的本地回归工具，可能。@joran:我理解。尽管即使从函数本身的源代码来看，我也不清楚需要做哪些更改才能达到目的……谢谢，Josh，这听起来确实比我目前所做的（对scb对象crit（）的简单应用）要好得多。不过，我不知道这是否适用于scb。这不是单个locfit对象。当我窥视它的来源时，找到locfit呼叫也不是那么容易，我可能会尝试破解它。怎么办？嗨，拉斯洛。看起来比我想象的要复杂scb（）
带公式调用locfit（）
，后者依次调用locfit.raw（）
（作为其lfproc=
参数的默认值。在locfit.raw（）
中，键计算由编译的C函数（.C（“slocfit”）执行
，上下限通过Z$L返回。因此（无需深入挖掘），您可能需要更改C代码以获得不同的CI（在此包中）。祝您好运！
library(locfit)
data(ethanol)
fit <- locfit(NOx~E,data=ethanol)
crit(fit) <- crit(fit,cov=0.99)
plot(fit,band="local")