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

419595891

% g1 x! C# Y$ {! D: L7 S& ^8 u+ J
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
0 C1 Z5 M- _" l) E& z' q
library(Design)                                                        
+ S9 h: `) }5 x8 Lreri <- function(datsam)                                                        
& W  s- ?- p: \- c$ _! N- [{                                                        
fitlr <- 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                                                        
}                                                        
3 g& J7 V7 v  u8 R$ O% e2 osummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
1 v3 j0 x( c7 {
! s& Y) c' m6 D5 K————————————我是分割线——————————————————————————————      


0 H& O6 t( V3 c' i* m2 N8 t
  H) G3 i3 H! Q8 ? library(Design)
2 J3 f+ d( e  \: ~> reri <- function(datsam)
) T2 d6 W* y' S. a) @+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
% u' W; Q% i! D2 t2 L+ 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 R: Y$ P1 V$ q                (You probably wanted to return the left-hand side)
  W  T- M) ]# Z+ m+ l& z         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
; l: v& g: y, B( X2 p> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
Warning messages:
+ K$ D' M5 ^9 o8 `& ~  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
9 F% |" ]" M* g7 t' M& R. lProblem: Object "DM.0813" not found
. Q) @4 I+ B  e+ S7 VUse traceback() to see the call stack


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


" X2 z( y7 ?# U, z+ \/ L