我试图计算杠杆和厨师的距离为我开发的模型使用nlme，但我得到的错误说
cooks.distance.lme（模型）中的错误：未为“nlme”对象实现
如何计算nlme中的Cook距离？

formula = log_Lobs ~ log(150*((1 + ((150/Lt_1)^(1/exp(p))-1)*exp(-exp(k)*td/365))^(-exp(p))))
model <- do.call(nlme,
                    list(formula,
                         fixed = c(p ~ 1, k ~ 1 + season2),
                         random = k ~ 1 | id,
                         data = data_select,
                         start = list(fixed = c(p, k)),
                         na.action = na.exclude,
                         control=list(maxIter=1e6, msMaxIter = 1e6, msVerbose = TRUE)
                    ))

#Calculate leverage
lev<-hat(model.matrix(model))

字符集
错误（x$项%）||% attr（x，“terms”）%||% stop（“no terms component nor attribute”）：无术语组件或属性

#Calculate Cook's Distance
cd<-cooks.distance(model)

型
cooks.distance.lme（模型）中的错误：未为“nlme”对象实现