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

419595891

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

library(Design)                                                        
$ Z% v0 u( r. u  z' h- yreri <- function(datsam)                                                        
1 {3 ?) `& v# W( k  T- }) ]{                                                        
: D  d% p# Z/ g1 E9 Tfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
# A/ e: Z9 e# Q7 y& u! ^. Breri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

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

' S  @0 I# U" g; h$ u
& u7 U! }7 r# k3 ~" [$ J/ t
library(Design)
( P7 ?+ B" ?2 ~& p. |> 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)
         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
4 \  P3 ?! d9 X% U> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
0 {- ~" B6 I3 p% J# ]+ \, L' m6 pWarning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
, ~6 r6 w  c6 a) t. _3 k5 a; q$ PProblem: Object "DM.0813" not found
Use traceback() to see the call stack

$ t. W. U1 M2 ]' m
+ ~" g! w+ s- m* a* |3 U6 Z' Q& z不明白这个错误什么意思，小白求指导……
* P; P% Z* R  e  u& B

, i3 X0 [6 e! l8 l2 Z" A