我为一个复杂的勘测定义了一个逻辑模型,该模型中的花键变量edadt定义了自由度,但没有定义结点的位置。在运行该模型后,
> log3<-svyglm(compo ~ bs(edadt,degree=1, df = 3), dclus, family = quasibinomial)
我想知道结点的位置在哪里 当我输入
splineKnots(log3)
这是我得到的。
Error in UseMethod("splineKnots") :
no applicable method for 'splineKnots' applied to an object of class "c('svyglm', 'glm', 'lm')"
有什么办法可以知道结点的位置吗?
谢谢你。
该 splineKnots
函数对很多对象不起作用。 我想你得自己写一个方法来解决这个问题。
library(splines)
splineKnots.default <- function(object) attr(object, "knots")
edadt <- rnorm(1000)
splineKnots(bs(edadt,degree=1, df = 3))
注意,这个方法可以在 bs(edadt,degree=1, df = 3)
,不在 log3
.如果你真的想让它工作在 log3
,你必须弄清楚这个对象是什么类,并找出它是否有以下结果 bs(edadt,degree=1, df = 3)
存储在它的任何地方,或公式,所以你可以重建它。