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

419595891

7 _) _2 L4 I! y$ V
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？

library(Design)                                                        
reri <- function(datsam)                                                        
, q4 k+ @+ r  \, y+ s7 r! f0 I{                                                        
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
/ v& K+ P  k( C- B1 Hreri <- 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))                                             

1 [$ S) J, R' b1 l# ^7 i————————————我是分割线——————————————————————————————      
8 p3 S4 _5 F  [3 `+ C; r/ n1 r# \
5 P2 ?- R* f( N# g
5 z2 q9 ~% p2 B' J$ M- E/ N
library(Design)
> reri <- function(datsam)
: ]3 ?, M! h8 c) S* G" g; q8 }+ {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}
    Last expression in function is an assignment
2 [# `+ O, u. o; n$ o2 S* }                (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))                                               
/ |4 a) o- M' lWarning messages:
+ P) y' }7 K" U; ]  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
Problem: Object "DM.0813" not found
0 J6 O2 g4 I5 @6 ^, ]Use traceback() to see the call stack

5 P- @( Z0 H6 D4 s8 v( x" [4 O
* \2 Q2 B$ [8 y3 H. {5 T不明白这个错误什么意思，小白求指导……
1 J* g7 w$ P4 h7 G6 F' T
" L/ m5 c9 S$ @1 B
8 A- o; B/ k" h0 w7 z