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

419595891

$ d, k, Z/ s  M+ e& t+ N9 R
+ u$ Q( n2 w# V: [# M' w7 e这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
  i6 ^$ Y4 e5 e7 B
library(Design)                                                        
, o9 Q' C+ V; \reri <- function(datsam)                                                        
{                                                        
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
+ {4 c: E/ W# l0 G' o' K0 Oreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
& A3 e% Y9 q6 B6 l4 {}                                                        
9 z6 c+ O/ z1 Z: ]# Vsummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
0 L1 ^6 e. p0 c2 L6 g* _7 `" Z0 }
————————————我是分割线——————————————————————————————      
  F  S3 [6 \" z# q/ |. m+ j+ X

* K+ U7 g# Y# x: Y6 b1 l( e! i0 Y
+ ~( }; K$ _( V9 A& i; z9 z library(Design)
> reri <- function(datsam)
+ {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}
# S( |% q7 g$ s) I! b  M  y  g8 D/ N2 v    Last expression in function is an assignment
                (You probably wanted to return the left-hand side)
         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:
* y; l" O% n* B$ U  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
4 a0 m4 ?4 L  z3 L! R' {! ?
/ g) U6 f7 I0 i
不明白这个错误什么意思，小白求指导……

# Y# e$ X6 _* N6 Q4 L