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

419595891

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

/ p- _# g: y" Z3 t3 alibrary(Design)                                                        
3 y* x5 F( k9 T  H( a2 ?0 Q6 oreri <- function(datsam)                                                        
$ M6 L1 _9 Q. ~& ^& t{                                                        
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                                                        
: ~/ ~. s2 G( Y}                                                        
! Z. Z7 `. C) h9 d7 B) bsummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
: K" E$ C: J5 v: p) ^, Z- r' w% R
————————————我是分割线——————————————————————————————      

, L' W8 M: H5 U) F
; H! c0 g7 `4 ]9 v$ b, L. P$ C8 z( N
2 w. g3 @5 N4 l" B' O, p library(Design)
> reri <- function(datsam)
% V! L, y3 l  _+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
& ?3 q, o2 W2 L1 H8 ^) [: _2 r+ 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
/ A+ s/ [: }7 Z8 F6 R' Z                (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))                                               
* `# o/ E# E2 r, O: ]+ Q# Y& }Warning messages:
! R/ j6 Q4 i& R+ G9 F9 T  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
Problem: Object "DM.0813" not found
. H: x( E* g$ qUse traceback() to see the call stack
8 r# ?# G3 z- R. d- Y

# |* S4 @) q4 u: S不明白这个错误什么意思，小白求指导……
: c0 x$ y0 [5 W0 P- Q
5 w0 w$ U+ v" p& m4 B8 n) G
9 u9 n$ l/ m) |( G( M; \