#P37
help(glm) #获得某个函数的帮助信息
?glm #获得某个函数的帮助信息
?'+' #获得某个操作符的帮助信息
example(glm) #查看函数的示例
##
## glm> ## Dobson (1990) Page 93: Randomized Controlled Trial :
## glm> counts <- c(18,17,15,20,10,20,25,13,12)
##
## glm> outcome <- gl(3,1,9)
##
## glm> treatment <- gl(3,3)
##
## glm> print(d.AD <- data.frame(treatment, outcome, counts))
## treatment outcome counts
## 1 1 1 18
## 2 1 2 17
## 3 1 3 15
## 4 2 1 20
## 5 2 2 10
## 6 2 3 20
## 7 3 1 25
## 8 3 2 13
## 9 3 3 12
##
## glm> glm.D93 <- glm(counts ~ outcome + treatment, family = poisson())
##
## glm> anova(glm.D93)
## Analysis of Deviance Table
##
## Model: poisson, link: log
##
## Response: counts
##
## Terms added sequentially (first to last)
##
##
## Df Deviance Resid. Df Resid. Dev
## NULL 8 10.58
## outcome 2 5.45 6 5.13
## treatment 2 0.00 4 5.13
##
## glm> ## No test:
## glm> summary(glm.D93)
##
## Call:
## glm(formula = counts ~ outcome + treatment, family = poisson())
##
## Deviance Residuals:
## 1 2 3 4 5 6 7 8
## -0.6712 0.9627 -0.1696 -0.2200 -0.9555 1.0494 0.8472 -0.0917
## 9
## -0.9666
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 3.04e+00 1.71e-01 17.81 <2e-16 ***
## outcome2 -4.54e-01 2.02e-01 -2.25 0.025 *
## outcome3 -2.93e-01 1.93e-01 -1.52 0.128
## treatment2 1.34e-15 2.00e-01 0.00 1.000
## treatment3 1.42e-15 2.00e-01 0.00 1.000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for poisson family taken to be 1)
##
## Null deviance: 10.5814 on 8 degrees of freedom
## Residual deviance: 5.1291 on 4 degrees of freedom
## AIC: 56.76
##
## Number of Fisher Scoring iterations: 4
##
##
## glm> ## End(No test)
## glm>
## glm> ## No test:
## glm> ## an example with offsets from Venables & Ripley (2002, p.189)
## glm> utils::data(anorexia, package = "MASS")
##
## glm> anorex.1 <- glm(Postwt ~ Prewt + Treat + offset(Prewt),
## glm+ family = gaussian, data = anorexia)
##
## glm> summary(anorex.1)
##
## Call:
## glm(formula = Postwt ~ Prewt + Treat + offset(Prewt), family = gaussian,
## data = anorexia)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -14.108 -4.277 -0.548 5.484 15.292
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 49.771 13.391 3.72 0.00041 ***
## Prewt -0.566 0.161 -3.51 0.00080 ***
## TreatCont -4.097 1.893 -2.16 0.03400 *
## TreatFT 4.563 2.133 2.14 0.03604 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 48.7)
##
## Null deviance: 4525.4 on 71 degrees of freedom
## Residual deviance: 3311.3 on 68 degrees of freedom
## AIC: 490
##
## Number of Fisher Scoring iterations: 2
##
##
## glm> ## End(No test)
## glm>
## glm> # A Gamma example, from McCullagh & Nelder (1989, pp. 300-2)
## glm> clotting <- data.frame(
## glm+ u = c(5,10,15,20,30,40,60,80,100),
## glm+ lot1 = c(118,58,42,35,27,25,21,19,18),
## glm+ lot2 = c(69,35,26,21,18,16,13,12,12))
##
## glm> summary(glm(lot1 ~ log(u), data = clotting, family = Gamma))
##
## Call:
## glm(formula = lot1 ~ log(u), family = Gamma, data = clotting)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.0401 -0.0376 -0.0264 0.0290 0.0864
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.016554 0.000928 -17.9 4.3e-07 ***
## log(u) 0.015343 0.000415 37.0 2.8e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for Gamma family taken to be 0.002446)
##
## Null deviance: 3.51283 on 8 degrees of freedom
## Residual deviance: 0.01673 on 7 degrees of freedom
## AIC: 37.99
##
## Number of Fisher Scoring iterations: 3
##
##
## glm> summary(glm(lot2 ~ log(u), data = clotting, family = Gamma))
##
## Call:
## glm(formula = lot2 ~ log(u), family = Gamma, data = clotting)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.0557 -0.0293 0.0103 0.0171 0.0637
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.023908 0.001326 -18.0 4.0e-07 ***
## log(u) 0.023599 0.000577 40.9 1.4e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for Gamma family taken to be 0.001813)
##
## Null deviance: 3.118557 on 8 degrees of freedom
## Residual deviance: 0.012672 on 7 degrees of freedom
## AIC: 27.03
##
## Number of Fisher Scoring iterations: 3
##
##
## glm> ## Not run:
## glm> ##D ## for an example of the use of a terms object as a formula
## glm> ##D demo(glm.vr)
## glm> ## End(Not run)
## glm>
## glm>
help.search("regression") #获得某个主题的帮助信息
??regression #获得某个主题的帮助信息
library(help="grDevices") #获得某个包的帮助信息