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

419595891

" S( P* E6 ?1 b7 y$ I" i这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
4 O' w* y3 Z  M" z- ~
/ |1 `' M2 f' g7 G0 J: ~library(Design)                                                        
reri <- function(datsam)                                                        
{                                                        
( c6 Q1 M% w. @fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
9 y9 S3 [3 q2 U1 preri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
}                                                        
! }4 x; P: Y/ j$ {; x( xsummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

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

" y4 K* p# q: }; p* \

% o* u$ I, S! T" S: X9 [4 y: s library(Design)
> reri <- function(datsam)
; c! a0 D5 O9 w% R+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
, W; T9 m4 E; }: H4 H+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
8 Y+ p* K( w8 b7 z: p    Last expression in function is an assignment
                (You probably wanted to return the left-hand side)
1 A: g0 ^1 K/ H% H         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:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
) o* G1 j! J: L5 LProblem: Object "DM.0813" not found
Use traceback() to see the call stack


( c8 t+ o9 a* U8 q4 p不明白这个错误什么意思，小白求指导……

- J3 Y: W7 p8 U( ~  g