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

419595891

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

1 O4 n" @" S6 q& j( flibrary(Design)                                                        
. H' ^5 |5 J, Preri <- function(datsam)                                                        
{                                                        
- x" ]4 h1 a. q2 nfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
9 ~# b1 K4 Y/ r* P5 V; {reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
3 F; v/ P7 p6 e3 e2 o0 a}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
: ]& V$ p8 _0 l; K
5 L1 Q% ^& w2 u# Z8 n9 E————————————我是分割线——————————————————————————————      
+ b& s* T  n& {! A$ M  \) n
( z: ?, i& m3 s' n; {

library(Design)
> reri <- function(datsam)
" ~; R6 _' F4 A+ {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}
, M% P: @: l! R* _8 d# D8 ^4 o; E" L    Last expression in function is an assignment
                (You probably wanted to return the left-hand side)
1 v+ B3 `9 B5 ^$ u; [$ s% J         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)
Problem: Object "DM.0813" not found
Use traceback() to see the call stack


不明白这个错误什么意思，小白求指导……
7 H- _9 |+ G  t( P% D* X0 o