dd <- datadist(data) options(datadist='dd')
fit1 <- ols(salary ~ rcs(salbegin,4)+eduy,data=data)
anova <- anova(fit1)anova
beta<-Predict(fit1,salbegin)ggplot(beta,anova=anova,pval=T)
ggplot(anova=anova,pval=T)+ #画曲线 geom_line(data=beta,aes(salbegin,yhat),linetype=1,size=1,alpha = 0.9,colour="red")+ #画置信区间 geom_ribbon(data=beta,aes(salbegin,ymin = lower, ymax = upper),alpha = 0.3,fill="red")+ #x轴任意刻度:增加一条竖线 geom_vline(aes(xintercept=13000), colour="#BB0000", linetype="dashed")+ theme_classic()
联系客服