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

419595891

这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
4 G& ?0 s/ N( Z4 G: Q5 j
# V" `6 w, i, U7 t! c& Vlibrary(Design)                                                        
reri <- function(datsam)                                                        
{                                                        
7 e9 G2 z2 R. m6 efitlr <- 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                                                        
- k! R( P% G0 y' d& w}                                                        
' E+ ^: i4 Z' r; esummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

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


0 f. l$ b6 k. e) i0 d1 I  D
library(Design)
5 z" C( P2 W4 P# }1 D> 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}
9 U2 |. O+ V# K2 x    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
> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
Warning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
2 f, u7 B' K. sProblem: Object "DM.0813" not found
8 h7 b) [; U5 q# O/ ^) F2 r% `Use traceback() to see the call stack
# G4 P$ v2 A  X4 P# j
( r4 _. H9 E! l, Q7 e
& S& G3 q3 J% i2 b. o( p, v  u不明白这个错误什么意思，小白求指导……

6 z& }) o5 ?' w# X/ k% @
: g4 c' H9 R; v8 l3 W* F- W