Script started on Wed May 14 17:23:25 1997
babcomm% S
S-PLUS : Copyright (c) 1988, 1996 MathSoft, Inc.
S : Copyright AT&T.
Version 3.4 Release 1 for Sun SPARC, SunOS 4.1.3_U1 : 1996 
Working data will be in /biostat/user7/norm/.Data 
> lei3_read.table("lei3.dat",header=T)
attach(lei3,pos=1)
ls()
y_Surv(set,case)
id_rep(1:5,length.out=315)
> >  [1] ".Last.value" "age"         "ageg"        "case"        "est"        
 [6] "gall"        "hyp"         "non"         "ob"          "set"        
> > > fit1_glm(case~est+factor(set),family=binomial,subset=(id<3))
sum1_summary(fit1,correlation=F)
> > 
> sum1$coeffients["est",]
Error in [: No dim attribute for array subset: # NULL
Dumped
> ls()
 [1] ".Last.value" "age"         "ageg"        "case"        "est"        
 [6] "fit1"        "gall"        "hyp"         "id"          "last.dump"  
[11] "non"         "ob"          "set"         "sum1"        "y"          
> attributes(sum1)
$names:
 [1] "call"           "terms"          "coefficients"   "dispersion"    
 [5] "df"             "deviance.resid" "cov.unscaled"   "correlation"   
 [9] "deviance"       "null.deviance"  "iter"           "nas"           

$class:
[1] "summary.glm"

> sum1$coefficients["est",]
    Value Std. Error  t value 
 4.537334  0.8554659 5.303933
> exp(sum1$coefficients["est",1])
[1] 93.44135
> fit2_coxph(y~est+strata(set),subset=(id<3))
sum2_summary(fit2)
> Call:
coxph(formula = y ~ est + strata(set), subset = (id < 3))

  n= 126 

    coef exp(coef) se(coef)    z       p 
est 2.27      9.67    0.606 3.74 0.00018

    exp(coef) exp(-coef) lower .95 upper .95 
est      9.67      0.103      2.94      31.7

Rsquare= 0.176   (max possible= 0.5 )
Likelihood ratio test= 24.4  on 1 df,   p=7.63e-07
Wald test            = 14  on 1 df,   p=0.000183
Efficient score test = 21.1  on 1 df,   p=4.3e-06

> table(est[id<3],case[id<3])
   0  1 
1 33  7
2 30 56
> (56*33)/(30*7)
[1] 8.8
> q()
Warning: Modified object in working data not saved
babcomm% exit
exit

script done on Wed May 14 17:27:13 1997