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

419595891

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

. z' B& r' q4 t3 tlibrary(Design)                                                        
' g7 z; d* e5 r0 H/ f/ R$ Treri <- function(datsam)                                                        
{                                                        
! a% T, R  y' j/ tfitlr <- 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                                                        
" c# a' _' \) P0 x9 A+ j}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

————————————我是分割线——————————————————————————————      
' \+ R! n1 H2 C

0 c. \2 J: ^* k- W  f: p
library(Design)
6 ~+ C2 |" V/ b> reri <- function(datsam)
: ]- Y$ Z" k' w/ ^. P+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
/ j$ i9 o" o2 D2 W$ @  H/ C+ 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
' j) c: n( f% [/ Y. o1 J: 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))                                               
6 X8 g* d8 \3 N6 v+ e( z$ @- SWarning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
Problem: Object "DM.0813" not found
4 u* u$ }: O6 s  _/ m) o8 tUse traceback() to see the call stack


9 z+ D) E. ]1 T- Z" F$ c5 E/ N不明白这个错误什么意思，小白求指导……

# ~/ q9 R& \) r! X, l& }
6 P4 U3 T# W/ T# K3 q. [5 X