s2<-function(datsam,indices){d<-datsam[indices,]}中datsam和 indices是什么意思
附代码:
s1 <- function(datsam,indices){d <- datsam[indices,]}
fit<glm(hypertension~HLS1*lipid_class+GENDER+edu_new+marriage_new+area+index+hhsize+d3kcal1+DM+hhincpc_cpi1+HS_CRP+age_class2,data=chns_final4,family =binomial())
s <- (exp(fit$coef[2]+fit$coef[3]+fit$coef[18])-1)/(exp(fit$coef[2])-1+exp(fit$coef[3])-1)
outs <- boot(data=chns_final4,statistic=s1,R=1000)
print(outs)###输出点估计值
plot(outs)
boot.ci(outs,conf = 0.95,type='perc')


雷达卡


京公网安备 11010802022788号







