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

419595891

* m5 A, ?" X6 J! Q1 \; p7 q: a
5 F4 x# [' K* H2 o) M/ K. s2 K这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
2 W* h* D3 ]9 z4 _! l
library(Design)                                                        
reri <- function(datsam)                                                        
: V0 _5 \  s& B- T{                                                        
2 X' c4 c* O- [- J% p+ l9 W9 N0 gfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
, i' Y) W! T$ `7 \& Y2 sreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
6 d1 r8 o8 [5 h4 V}                                                        
+ j. C7 X; b" H, Q, v: esummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

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

# ?3 u. J3 b: F! Z* Y* T, Y
# {" w; q' I0 @. q2 w' D( ^* c
; s0 {- r; f4 w library(Design)
2 u& |" P) H, F5 E  h  w> 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}
7 W3 C; Z7 d$ Q, f& T) T    Last expression in function is an assignment
3 _, W+ {8 m3 O7 z4 |                (You probably wanted to return the left-hand side)
- C( R3 q& O  H         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
% |' I5 e$ i9 x! F3 v  t& }> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
Warning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
) W8 U, T3 @% R; I2 Y6 HProblem: Object "DM.0813" not found
Use traceback() to see the call stack

* F* c. [" o. e* P7 D" P
. l- n9 N( Q5 R( p( T不明白这个错误什么意思，小白求指导……


9 r& S! }  U2 _, z+ M