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

419595891

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

) |  c( W1 w' i+ J. Dlibrary(Design)                                                        
reri <- function(datsam)                                                        
$ ]( }: F* h& _) d6 {7 N1 k- z{                                                        
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
9 f6 C/ v/ l* Y4 b5 F9 e1 _. J6 h8 G% t: Nreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
0 w; F- A) Z8 E! k" }5 q. w7 z- \' U; b}                                                        
  v/ J# O6 A. D! X- osummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
3 O: p$ b$ J! C3 Z7 _0 p4 J
2 B/ m6 l/ g" O( }) X  o. I- I  N6 U————————————我是分割线——————————————————————————————      

3 c  q) d5 x. I1 K2 S8 G+ J! ~

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}
5 {! `8 d3 b: N5 J( j/ M* @) H    Last expression in function is an assignment
& K, v/ J$ ^/ f7 z( l2 }6 s' K' B# Y                (You probably wanted to return the left-hand side)
% j. p3 `9 w: o9 {1 O6 ]         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
' N% @/ Q2 q: x3 p# q> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
* g  S! j$ y+ J6 a$ q: WWarning messages:
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
7 T" B' [* |" v8 e0 XProblem: Object "DM.0813" not found
Use traceback() to see the call stack
! B1 e1 _3 S0 m) g" e2 {: `

% r4 d7 b' ?8 l不明白这个错误什么意思，小白求指导……