当我想用 ns 可视化样条逻辑回归时出现此错误:
Error in quantile.default(x[!outside], knots) : missing values and NaN's not allowed if 'na.rm' is FALSE
data <- read.table("e:/data_logistic.txt", header = TRUE)
dd <- datadist(data)
options(datadist = "dd")
model <- lrm(formula = P_A ~ ns(HEIGHT, 4), data = data, x = TRUE, y = TRUE)
p1 <- ggplot(Predict(model, HEIGHT))
当我使用 rcs 而不是 ns 时,我没有任何问题。我想使用“ns”方法。我该如何解决这个问题?
谢谢。