VGAM包在R中包含函数AICvlm,这与基本R中的AIC函数不同。下面是一个不同之处的例子:
require(VGAM)
fit = vgam(agaaus ~ s(altitude, df = 6), binomialff, data = hunua)
AICvlm(fit) # I get 386.7385
AIC(fit) # I get 395.2227我对AICvlm上的文档页面感到困惑。具体来说,示例部分只包含AIC函数的示例用法,而没有提到AICvlm。那么,我们这里是否有(a)文档中的错误,或者(b)代码中的错误,或者(c)其他的东西?
发布于 2013-11-01 15:56:30
fit属于"vgam"类,因此AIC()将通过带有签名"vgam"的AIC()的S4方法分派给AICvgam()函数。
R> getMethod(AIC, "vgam")
Method Definition:
function (object, ..., k = 2)
AICvgam(object, ..., k = k)
<environment: namespace:VGAM>
Signatures:
object
target "vgam"
defined "vgam"我怀疑(据称)内部函数(AICvlm())在这种情况下完全有效(即不抛出错误)是由于"vgam"类拥有"vglm"类对象的所有组件(加上一些其他组件)。无论如何,如示例所示,使用AIC()泛型,让VGAM的作者为您整理详细信息。
回答你的具体问题:( c)其他的事情。特别是用户错误。
https://stackoverflow.com/questions/19730003
复制相似问题