请教大神这个S-plus script是否正确？可以改动转换到R运行吗？

419595891

/ v+ j- u5 _% E4 D0 x% [5 i
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？

library(Design)                                                        
reri <- function(datsam)                                                        
/ B. P- `! C" Y0 X+ h+ `& _{                                                        
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
" j: x& V( ^% V4 J4 Yreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

————————————我是分割线——————————————————————————————      



9 `2 x0 A# T' f library(Design)
4 _- ^: [, Z2 Z> reri <- function(datsam)
2 D9 K) ~$ w  T, s1 h( c6 S+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
1 J2 e6 |5 M3 o7 t# S) j% r2 W+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
    Last expression in function is an assignment
                (You probably wanted to return the left-hand side)
         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
Warning messages:
) D1 o; T4 s6 s# z9 i2 G9 V  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
Problem: Object "DM.0813" not found
+ {$ E% f) |) c4 a+ N& G2 [4 fUse traceback() to see the call stack
: F7 O8 ~% h. k  f

不明白这个错误什么意思，小白求指导……

3 ?- J' Q' i0 ]; B. R2 y. t
3 o9 I5 j. c4 Y( f" u; `* j( l9 Z/ q