[英]how to get df from `mgcv::logLik.gam` output
在 R 中,package mgcv
具有logLik.gam
,它給出廣義相加模型的對數似然值以及 df。 但是,我不能只取 df 值。
例如,
dt= gamSim(1,50)
x= dt$x0
y=dt$x1
model= gam(y~s(x))
logLik.gam(model)
[1] 'log Lik.' -1.560092 (df=3.612741)
由此我只能取-1.560092
。 但是要使用 df 的值(即 3.612741),有什么辦法可以從 output 中獲取這個值嗎? 這里的output是一個數字,其實就是一個句子。 那么有沒有什么辦法可以切斷這條 output 線,只取最后一個值呢?
是的。 如果您首先將結果分配給 output 變量名,您就可以使用它。 該 class 的對象有一個打印方法:
res <- logLik.gam(model)
print( res)
'log Lik.' -4.613525 (df=4.322011)
attributes( res)
#-----
$df
[1] 4.322011
$class
[1] "logLik"
attributes( res)$df
#[1] 4.322011
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.