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

419595891

! I% A! Y6 M( a- h0 V$ S# d
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？

library(Design)                                                        
reri <- function(datsam)                                                        
) ^, Y+ v9 j0 p4 _  C0 D{                                                        
& l1 A3 F  Q$ U# Afitlr <- 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                                                        
$ G6 ^" b# A4 e9 `9 t8 N}                                                        
: l* a& H' W- s( A* I* g1 J# psummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
/ Q+ ]! {# |) p" Z- E
* R0 C3 Z  V: ~————————————我是分割线——————————————————————————————      

6 N. d( f, r8 U! J- J

library(Design)
4 l+ O  Y+ D+ d> reri <- function(datsam)
+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
/ I+ U7 q8 J* |4 z% g; f+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
( U1 t' {+ x& k0 O    Last expression in function is an assignment
                (You probably wanted to return the left-hand side)
6 b; O) Z* `  k! M# d/ t         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))                                               
+ j7 f( w; f: O2 ^  uWarning messages:
+ B# X& H  F* O& s3 L  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
1 _: [1 Q" m  N, BProblem: Object "DM.0813" not found
Use traceback() to see the call stack
" ^' Z0 ?3 p- g9 ~' j8 G% w: ]

3 N* [' F+ g! X不明白这个错误什么意思，小白求指导……
. p( l/ t3 }) c- W+ q
! l/ ~. A" E2 w* Y8 s4 C2 X# t