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

419595891

/ ]! o& ~7 Z" @8 }8 y3 y* E* i6 F; m
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
6 H; o7 K+ X8 r6 k' M
library(Design)                                                        
reri <- function(datsam)                                                        
{                                                        
. F1 c" V/ ], j5 wfitlr <- 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                                                        
4 w% o  t0 X; b/ T2 |9 k}                                                        
) n5 K! `' N5 f2 E! {# Qsummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
: y0 w' x/ K% S1 f
————————————我是分割线——————————————————————————————      
( l% G* V, a1 b) R8 @  z$ z
) p' r% D+ }% ~4 M

library(Design)
> reri <- function(datsam)
+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
2 X+ B8 x1 b3 \  K; N+ 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)
$ f$ i$ k( y8 ?$ d7 T         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))                                               
/ w3 J; N; Q! }, s. S8 l9 E7 S5 kWarning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
/ r1 Q/ w! w- w: X: K7 w0 ?Problem: Object "DM.0813" not found
4 t3 p5 m: f' bUse traceback() to see the call stack
% @8 R, w: u* K% ]8 B

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