回复 第2楼 的 波波头一头:
看来真是高手在此啊!
还有一问题,在SAS中,我用proc gam来做曲线拟合,在R中也有,SAS中有一个选项是直线加曲线additive 选项for plot,而且直线和曲线的估计都有,R中好象出来只有曲线拟合的图,我把曲线拟合的变量再加到线性变量中,估计值就明显是错的了,
SAS code
ods graphics on;
ods html
gpath="c:\data\image"
style=journal;
proc gam data=slcrp8 plots=components( additive clm unpack);
model crp=spline(age) param(bmxbmi);
title "crp&age";
title "";
run;
ods graphics off;
R code1
CRP.AGE.T<-gam(crp~s(age,bs="cr",fx=TRUE,k=6)+BMXBMI,na.action=na.omit,data=VD)
plot(CRP.AGE.T,,ylim=c(-1.75,1.75), cex.axis=1.35,cex.lab=1.45,mgp=c(2.5,0.75,0),
xlab="Age (yrs)",ylab="Log hS_CRP (mg/dL)",main="")
summary(CRP.AGE.T)
R code2
CRP.AGE.T0<-gam(crp~s(age,bs="cr",fx=TRUE,k=6)+age+BMXBMI,na.action=na.omit,data=VD)
plot(CRP.AGE.T0,,ylim=c(-1.75,1.75), cex.axis=1.35,cex.lab=1.45,mgp=c(2.5,0.75,0),
xlab="Age (yrs)",ylab="Log hS_CRP (mg/dL)",main="")
summary(CRP.AGE.T0)