我想做一个自变量和一个因变量之间的回归分析（两个变量均为连续型变量，且呈非正态分布），并作限制性立方样条分析，用lm函数后发现没有线性关系，用lrm后发现有非线性相关关系，想问下我这样做对吗？怎样做更合适呢？fit <- lrm(y~rcs(Pb,4),data=data, x=TRUE,y=TRUE)
anova(fit)
p <-round(anova(fit)[,3],3)
ggplot(pred_OR)