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

419595891

2 m2 c2 ]6 H9 U& b, ~
, t. S- \2 |  k5 z) e这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
5 t+ R' C, g% U3 i. n  [
1 M* L" B4 F( B! @* L) W3 }* wlibrary(Design)                                                        
reri <- function(datsam)                                                        
+ b; p3 E; q" \% s" Y2 e{                                                        
0 h' E" t! N. [* y1 Zfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
}                                                        
5 p2 }/ w1 E* P) d8 Psummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

————————————我是分割线——————————————————————————————      
  a) W+ G2 H, s$ J) M


8 o" l) \  s, e% L' g3 l library(Design)
2 A/ |5 s( F3 m7 E> reri <- function(datsam)
8 _, ~+ T/ X9 G6 ~6 Y' {, j  D* R+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
+ [$ |) B7 L: `4 b  f6 Q+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
; N& {# G5 P% i    Last expression in function is an assignment
3 |. k& G' |' b$ \# S  a                (You probably wanted to return the left-hand side)
1 O; Y% F8 j6 D2 j7 A         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:
2 l# G% k( C4 v  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
" F# S  A8 X' Z  n4 HProblem: Object "DM.0813" not found
Use traceback() to see the call stack

$ W& W$ o- J# {  s3 S
不明白这个错误什么意思，小白求指导……
0 y! e+ D# X/ F; x' d
8 D) ]: s# b7 D- D8 k5 m