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

419595891

" e  e7 L% Q6 ]9 {+ d3 e# G5 y& X
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
. J* X4 A6 H& v; B& X9 }: g' y
3 i! |: T, N8 ^: {& `% Alibrary(Design)                                                        
reri <- function(datsam)                                                        
{                                                        
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                                                        
}                                                        
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
( V6 @) E$ l6 U
* f: v( h/ b0 l2 M8 _7 V9 w————————————我是分割线——————————————————————————————      



* l+ X) o' M. n( o$ C+ V8 ^ library(Design)
* Y- |8 V, p# j$ o- S: E> reri <- function(datsam)
3 g" ~4 B/ z5 X+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
) ~( b9 M. F8 B9 {; u: X+ 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
                (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
# Q6 _; T- c/ ]' s> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
% t' k/ c% B( @" W' eWarning messages:
- }8 n% g% n" g  o; ?; C; x/ h  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
! Q, d' n, F4 V- t- bProblem: Object "DM.0813" not found
. a/ U' {+ S. L* h( s7 z, ?0 eUse traceback() to see the call stack
" w( v3 b, t5 a# ?; v+ I: S- V
) `5 k: l2 C+ o
. X, N" S* T8 R# u. g7 A! U* J不明白这个错误什么意思，小白求指导……

6 ~- R% n! q. \8 b* B& S