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

419595891

# k) Z- i- |2 r; Z+ J
* k& z( P# z* y2 ^6 m( b这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
  k2 w5 x* x; t4 W
; O, y! e- o/ c' \library(Design)                                                        
reri <- function(datsam)                                                        
' u- Q( V2 D! Q; V! s+ _0 t# G6 @{                                                        
, o, `' n4 N; zfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
% B0 e+ g; t+ w* D- [5 T8 B}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
0 C; c$ @5 p' U9 v3 {
————————————我是分割线——————————————————————————————      
* z6 ~7 u' y4 t0 W


+ V$ U3 c/ N" `4 f library(Design)
> reri <- function(datsam)
+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
8 e! C8 b* k, B6 @0 ]% {8 J+ 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
                (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
; y6 Y+ H3 \2 M8 U  `" h  I> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
Warning messages:
! U6 i  Q+ M. a/ U% R& s  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
Problem: Object "DM.0813" not found
Use traceback() to see the call stack


不明白这个错误什么意思，小白求指导……