如何在R中绘制90％置信度带与locfit

从https://stats.stackexchange.com/questions/16811/how-to-plot-90-confidence-bands-with-locfit-in-r交叉发布。如何在R中绘制90％置信度带与locfit

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

来源

2011-10-10 László

我只瞥了一眼'scb'代码，但似乎95％CIs只是作者写入函数的唯一选项。为什么？可能是因为这对他们来说似乎不是一个重要的功能。你有选择吗？可能你自己写或者使用不同的本地回归工具。 – joran

@joran：我明白。虽然从函数本身的来源来看，我不清楚需要改变的地方到达那里...... –

，通过它可以指定CI的装置是有点古怪，但它可以很容易，一旦你知道该怎么做。

请参阅功能crit<-，其帮助页面说明了此示例中的用法。

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

来源

2011-10-11 03:10:03

谢谢，乔希，这听起来好多了，比我迄今为止（一个天真的应用程序的暴击（）到碰撞的Rb对象）。不过，我不知道这是否适用于scb。这不是一个单一的locfit对象。当我窥视它的来源时，找到locfit电话并不容易，我可能会试图破解它。有什么可做的？ –

嗨拉斯洛。看起来这比我想象的更复杂。使用一个公式调用'locfit（）'''scfit（）''，它又调用'locfit.raw（）'（作为其''lfproc ='参数的默认值。''locfit.raw（）''通过一个编译的C函数（'.C（“slocfit”，...）'执行，下限和上限通过'Z $ L'返回。所以（没有深入挖掘），看起来你可能需要修改C代码以获得不同的配置项（在此包中）。祝您好运！ –

如何在R中绘制90％置信度带与locfit

回答

相关问题