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

419595891

6 L& m: T5 Q( ^
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
0 Z. J  o) K' @, L' ?' c
2 J) V0 t2 f' z1 ]8 q7 ?5 Qlibrary(Design)                                                        
: w  M  G% {, Z" I; c. o0 creri <- function(datsam)                                                        
{                                                        
/ y5 z* f. p4 |2 u( [% yfitlr <- 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                                                        
+ K5 I/ d3 }0 ]6 l* s1 u) Y8 X}                                                        
+ r8 T6 R. G# u7 n% Ysummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
$ v% @7 |8 g6 R& \
————————————我是分割线——————————————————————————————      
8 _; M& q5 r+ V7 B$ v- ]/ u' J/ ~
6 b" r' R* @0 S
+ C7 M8 ?# b$ _7 c' J6 ]* K+ Z* \
library(Design)
5 N9 o9 w: c6 w; M$ P% Z, L> reri <- function(datsam)
+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
3 f2 `4 k0 M* _* s9 k" E: z    Last expression in function is an assignment
7 D9 u+ m: E) X4 }8 E# w                (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
" _- c4 L6 u4 j. l( M  V' a> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
8 \5 a  g+ ]. K" PWarning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
0 U7 _4 }. \% BProblem: Object "DM.0813" not found
3 C6 {3 K* ?* z- V5 a) yUse traceback() to see the call stack


; d/ C; P. V$ o' _, m不明白这个错误什么意思，小白求指导……