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

419595891

! e; X! B/ p& D9 v% G% G
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
+ p) D; W/ ?7 S; k% \$ k5 g# [
library(Design)                                                        
7 O' g& g# l5 j& T& \reri <- function(datsam)                                                        
{                                                        
5 }% g, @, s7 W+ _( T4 dfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
1 ^7 @" m+ j) \# k. ireri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
2 k4 h: ?# b& L* r2 N8 g}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
' M3 k0 e2 ^5 K1 P1 {  v' O
————————————我是分割线——————————————————————————————      

/ W# D; F7 o  h" x
: o- V0 G" b6 X/ ~  n0 ~
8 s) D( ]0 @6 M2 o& t# h library(Design)
' q' S% K1 c! [) S  _  p> reri <- function(datsam)
! l2 H' f: U. x% O/ k+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
3 X$ J7 `& I- j0 w& \+ 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
8 M. I$ D. e8 e; x- F: O7 q' ~, j- |5 Z                (You probably wanted to return the left-hand side)
  I( ^$ ?0 M% A1 j9 l$ l         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
/ c1 r* o- y* i, i1 {9 F& d7 E. Y> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
" `! e: k. v/ `# rWarning 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
6 q1 ?4 a- x  F' w3 Z3 W4 m

不明白这个错误什么意思，小白求指导……
( ?8 \& Y& d- t' r  m3 N