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

419595891

1 A$ i% V6 N) F0 @: q1 m. y
3 o1 o+ B2 @. n: }' s& F这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
/ k* C9 w6 u+ R& Y7 }
library(Design)                                                        
reri <- function(datsam)                                                        
{                                                        
fitlr <- 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                                                        
: q. C0 z" I( P4 O" v  l}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
1 T# D2 Z6 y4 N. x- W2 H1 R; t
! R$ Y( \) y; B) N( r( `9 }9 L————————————我是分割线——————————————————————————————      
4 ?% r' G9 t" \1 D
$ R* y1 \0 `* P. q$ ]' o6 j
7 Y/ V% Q4 {! d0 l! O
# |  J) h, @, R  I3 X! w library(Design)
> reri <- function(datsam)
* G6 p, _, G5 T+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
! I$ v. j  P0 O) O, m0 v+ 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
1 f. V$ O  b8 i9 t: s                (You probably wanted to return the left-hand side)
+ F" \/ d5 k9 @9 S5 Z         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))                                               
, b! j: _1 H# z3 U7 LWarning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
! r* _8 T; `! W3 GProblem: Object "DM.0813" not found
! j  v( L/ S8 T  l. s& oUse traceback() to see the call stack


不明白这个错误什么意思，小白求指导……

- m% i' M" i) q! _) O6 u
/ z7 o. j3 J; z3 R