当使用R中的b作为变量时,从glm中提取结点

Question

我使用以下内容拟合模型:

mymodel <- glm(LS ~ bs(LA, df = 8) + bs(IN, df = 7),
               family = binomial, data = mydata, na.action = na.omit)

没问题,我现在有模型适合我试图提取使用的结点.我跟着一篇关于使用attr和提取结点的帖子str.那是一个只是样条的模型.我认为这个结在某种程度上就在结构中

我打电话str(mymodel$terms)有..-attr(*, "variables").我无法继续前进,attr但我相对肯定这基本上是我需要做的.任何获得结的指导都表示赞赏.

Answer 1

你走在正确的轨道上。mymodel$terms包含有关模型中的项的信息，并且attr(mymodel$terms, "predvars")是一个语言对象，它是具有节点计算值的预测变量列表。

为了让他们出去：

x <- attr(mymodel$terms, "predvars")
x[[2]] # bs(LA, degree=3, knots=<vector>, Boundary.knots=<vector>, intercept=FALSE)
x[[2]]$knots
x[[2]]$Boundary.knots