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

419595891

  P1 }! n  q$ }+ Z3 r, L
- E. I6 x* E+ K  v这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？

library(Design)                                                        
reri <- function(datsam)                                                        
{                                                        
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
$ h& H+ t. W+ s  q: Z! |% m( A. ]reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
; H! y7 l0 g7 u0 J( K; z9 y}                                                        
" y  H5 m* T- }  R& z) dsummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

————————————我是分割线——————————————————————————————      
. Y* c8 G* E% H0 ^5 z- U. F

4 U7 p7 _0 P; @9 H% K! r+ [
8 E% F/ p$ l, \- e library(Design)
% u( g% u9 d5 A# z1 L- k> reri <- function(datsam)
. |3 u5 H8 C5 ^+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
/ ~5 Q8 g9 u2 C+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
. i* J9 X, ?+ w    Last expression in function is an assignment
* w  L7 e" n$ \& ~6 g; n0 X% c3 K                (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))                                               
Warning messages:
+ B4 Y; d7 s7 j- ?. M  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
Problem: Object "DM.0813" not found
2 g' d; F5 H8 ^9 T+ [$ IUse traceback() to see the call stack


不明白这个错误什么意思，小白求指导……
# u" C8 y8 L6 j0 @5 ^  F, M

; B1 d4 U2 s+ t8 i! W