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

419595891

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

6 q( `% U9 J# P( i+ I3 }. {2 wlibrary(Design)                                                        
reri <- function(datsam)                                                        
* ?+ S7 w, m1 b* P' y$ K{                                                        
5 B2 u, z" L* t" Y' ^, m! Q0 Nfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
* D- w& f3 O3 d8 mreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
0 {" o% r* m+ w& w. m. r% k}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

1 l+ x" |3 o: u+ j! o$ o" U( g! D————————————我是分割线——————————————————————————————      
  M& R( \, C4 ^8 d4 l, E  U' `' q  G


6 }) D2 R$ W% z! m3 O" Y) N library(Design)
0 l+ R7 a2 h, c4 y! J, ~( c5 Y* r> 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}
    Last expression in function is an assignment
                (You probably wanted to return the left-hand side)
: K! ^. C8 ?( y# i; P         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))                                               
/ y6 b& }- b5 F+ j- l2 rWarning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
Problem: Object "DM.0813" not found
0 I& P# @3 |% |+ ?Use traceback() to see the call stack


不明白这个错误什么意思，小白求指导……
; t' _$ y7 p& d/ w

" L6 q' |, a& `; H