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

419595891

4 s9 c; w' K4 E1 g6 N" e$ ~5 J- D
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？

library(Design)                                                        
reri <- function(datsam)                                                        
{                                                        
5 |" S+ ]$ \% i7 F6 Xfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
9 B: W  _" J% Q$ T+ d' z3 o: Rreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
! W& j4 J( O# e- K}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
6 l) n9 A& _, Y7 l- k" b
————————————我是分割线——————————————————————————————      
  ^: n. q: d! a) i* o! z* f* S


library(Design)
> reri <- function(datsam)
8 O& W  H+ ^+ b0 g: I5 q+ {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}
4 ]7 z* x0 H* Z) s. l% |5 z5 D+ \    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
1 Z# N9 f" o6 M/ g+ f> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
1 h6 n" g# a8 S' l0 ~: xWarning messages:
' S! ^, Y5 W4 W, l% U$ j, l  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
4 b. p  T" D) Q3 \Problem: Object "DM.0813" not found
% k. |- @3 _8 t) |Use traceback() to see the call stack
5 Q1 u4 x8 ?, O2 G5 d2 N

) e. G% h+ h6 {6 N1 x# A不明白这个错误什么意思，小白求指导……
& b1 d9 b8 L5 X) K

6 I( p0 O% }$ Q9 _, O; C