The R Project SVN R

Rev

Rev 84054 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 84054 Rev 84152
Line 1... Line 1...
1
 
1
 
2
R Under development (unstable) (2022-07-07 r82558) -- "Unsuffered Consequences"
2
R Under development (unstable) (2023-03-28 r84086) -- "Unsuffered Consequences"
3
Copyright (C) 2022 The R Foundation for Statistical Computing
3
Copyright (C) 2023 The R Foundation for Statistical Computing
4
Platform: x86_64-pc-linux-gnu (64-bit)
4
Platform: x86_64-pc-linux-gnu (64-bit)
5
 
5
 
6
R is free software and comes with ABSOLUTELY NO WARRANTY.
6
R is free software and comes with ABSOLUTELY NO WARRANTY.
7
You are welcome to redistribute it under certain conditions.
7
You are welcome to redistribute it under certain conditions.
8
Type 'license()' or 'licence()' for distribution details.
8
Type 'license()' or 'licence()' for distribution details.
Line 2174... Line 2174...
2174
Response: cbind(ncases, ncontrols)
2174
Response: cbind(ncases, ncontrols)
2175
 
2175
 
2176
Terms added sequentially (first to last)
2176
Terms added sequentially (first to last)
2177
 
2177
 
2178
 
2178
 
2179
            Df Deviance Resid. Df Resid. Dev
2179
            Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
2180
NULL                           87     367.95
2180
NULL                           87     367.95              
2181
agegp        5  121.045        82     246.91
2181
agegp        5  121.045        82     246.91 < 2.2e-16 ***
2182
tobgp        3   36.639        79     210.27
2182
tobgp        3   36.639        79     210.27 5.485e-08 ***
2183
alcgp        3  127.933        76      82.34
2183
alcgp        3  127.933        76      82.34 < 2.2e-16 ***
2184
tobgp:alcgp  9    5.451        67      76.89
2184
tobgp:alcgp  9    5.451        67      76.89    0.7934    
-
 
2185
---
-
 
2186
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
2185
> ## Try a linear effect of alcohol and tobacco
2187
> ## Try a linear effect of alcohol and tobacco
2186
> model2 <- glm(cbind(ncases, ncontrols) ~ agegp + unclass(tobgp)
2188
> model2 <- glm(cbind(ncases, ncontrols) ~ agegp + unclass(tobgp)
2187
+                                          + unclass(alcgp),
2189
+                                          + unclass(alcgp),
2188
+               data = esoph, family = binomial())
2190
+               data = esoph, family = binomial())
2189
> summary(model2)
2191
> summary(model2)
Line 3588... Line 3590...
3588
> ### * <FOOTER>
3590
> ### * <FOOTER>
3589
> ###
3591
> ###
3590
> cleanEx()
3592
> cleanEx()
3591
> options(digits = 7L)
3593
> options(digits = 7L)
3592
> base::cat("Time elapsed: ", proc.time() - base::get("ptime", pos = 'CheckExEnv'),"\n")
3594
> base::cat("Time elapsed: ", proc.time() - base::get("ptime", pos = 'CheckExEnv'),"\n")
3593
Time elapsed:  5.453 0.122 5.614 0 0 
3595
Time elapsed:  1.134 0.029 1.167 0 0 
3594
> grDevices::dev.off()
3596
> grDevices::dev.off()
3595
null device 
3597
null device 
3596
          1 
3598
          1 
3597
> ###
3599
> ###
3598
> ### Local variables: ***
3600
> ### Local variables: ***