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

419595891

/ k/ }6 B* k2 M( k' Z5 [这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
7 u" u) `7 G! r9 m; o
library(Design)                                                        
reri <- function(datsam)                                                        
% H" |) G3 ^5 ~1 m0 R! j+ N) F{                                                        
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
% b8 A; T* `* I; n( L# Xreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
}                                                        
% {9 C) P+ r. [/ i( X. n+ osummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

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


3 ^; a" S  @3 m8 T! ~3 W
library(Design)
, S6 M. p" K: x9 l/ {> reri <- function(datsam)
: a, z  P* r9 _/ v' e  m+ {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
- l  J8 n3 W+ Y) K" z5 W' Q0 X                (You probably wanted to return the left-hand side)
% f3 Z# X( i' l) n' p% D6 e         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))                                               
) @' _, ?* Y; dWarning messages:
% m  j9 G* S6 V. U7 r  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
3 b, U8 ?# v' W- Q7 Z4 z2 U

+ s) w  D/ T: E* @. Y不明白这个错误什么意思，小白求指导……

5 |+ ]1 ?9 l& W9 V6 U