R语言中svyglm模型样条节点位置查询方法求助

Hey there! I’ve dealt with this exact problem before when working with splines in survey-weighted models—splineKnots() doesn’t play nice with svyglm objects because it doesn’t have a method for that class. But don’t worry, there are two straightforward ways to get those knot positions:

Method 1: Re-generate the Spline Object Directly #

The bs() function calculates knots based on the distribution of your edadt variable, not the survey weights. So you can recreate the spline transformation using the original variable data from your survey object, then extract the knots:

# Extract the edadt variable from your survey design object
edadt_values <- dclus$variables$edadt

# Recreate the bs spline object with the same parameters as your model
spline_obj <- bs(edadt_values, degree = 1, df = 3)

# Pull the knot positions from the spline object's attributes
attr(spline_obj, "knots")

This works because the bs() call in your svyglm uses the exact same logic to set knots when you run the model, so you’ll get identical positions.

Method 2: Extract the Spline Call from the Model Terms #

If you want to avoid manually retyping the bs() parameters, you can pull the exact spline call from your model’s terms and evaluate it:

# Get the terms object from your svyglm model
model_terms <- terms(log3)

# Identify the term that uses bs()
spline_term_index <- which(grepl("bs\\(", labels(model_terms)))

# Extract the bs() call from the terms
spline_call <- model_terms[[spline_term_index + 1]]

# Evaluate the call using the variables from your survey design
spline_obj <- eval(spline_call, envir = dclus$variables)

# Get the knots
attr(spline_obj, "knots")

Why splineKnots() Fails #

The error you’re seeing happens because splineKnots() only has methods for specific object types like lm or spline-specific model objects (like those from the splines package’s standalone functions). Since svyglm objects inherit from glm and lm but don’t have a custom splineKnots method, R doesn’t know how to handle it.

内容的提问来源于stack exchange，提问作者Ángel