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

419595891

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

% P5 h9 V* D, a4 [0 W/ `, Glibrary(Design)                                                        
' W; x& z. t/ U) B: Ureri <- function(datsam)                                                        
{                                                        
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
2 r) q6 c$ ]- `5 o! N+ creri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             

————————————我是分割线——————————————————————————————      
7 K' Q2 J  E$ l4 ~
( O- T: c1 E- {+ m
+ u& [- u& [# `0 G
library(Design)
' O; Q; \2 o: W3 `0 e> reri <- function(datsam)
8 K$ Q- w: Y8 E) r+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
& D! d  X- L& x% k8 b+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
0 K7 c5 `) n, g5 H, |$ v    Last expression in function is an assignment
" v% i: S: s" o& {1 q" Q/ D                (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))                                               
' U4 x; W6 Z" _& L% 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
1 P4 a2 f; }" I( r6 E* U% rUse traceback() to see the call stack
$ N4 p+ O( X* q0 u
: n: j+ m' F* W4 _
不明白这个错误什么意思，小白求指导……
) Y- G6 S% C/ Q$ ^
  P+ i4 h6 [5 x