如何从R包mgcv的logLik.gam输出中提取df值?
提取mgcv中logLik.gam返回的自由度(df)值
不用通过截取字符串的方式提取,logLik.gam返回的是一个logLik类对象,自由度值已经作为对象的属性存储了,直接提取属性是更可靠的方法:
先获取对数似然对象:
library(mgcv) # 先拟合一个示例GAM模型 mod <- gam(y ~ s(x), data = your_data) ll_obj <- logLik(mod)提取df值:
直接用attr()函数提取df属性:df_value <- attr(ll_obj, "df")
这样就能得到你需要的有效自由度值,比字符串截取稳定得多——毕竟如果后续logLik的输出格式有变化,字符串截取会直接失效,而属性提取是官方支持的标准方式。
内容的提问来源于stack exchange,提问作者user20536224
相关产品推荐
相关产品推荐

