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

419595891

" F3 f2 t% ^5 m1 }  v9 l3 J# e6 ?
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？

4 v  n( ?, Q' M' qlibrary(Design)                                                        
reri <- function(datsam)                                                        
{                                                        
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
& f' N; [8 Z4 H5 |* u) p. Jreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
/ W( G  j# c/ \7 H; H' }6 z}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
" c; c) k3 G& o: i6 s* V
————————————我是分割线——————————————————————————————      



library(Design)
# M( D- |! ]$ T: a' a> 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
0 W) m4 M1 Q% j                (You probably wanted to return the left-hand side)
+ _- z6 ^* J. W         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
5 J  M6 j+ Q( T7 b! @* W0 g> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
( ~4 X* n  c& H+ ]4 H1 kWarning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
( `: K# Z) M9 {  TProblem: Object "DM.0813" not found
Use traceback() to see the call stack

0 I% O- v; W3 w% ?6 a
不明白这个错误什么意思，小白求指导……


  O. j  a8 C: q$ M! H# n5 s