ANLY545 HW4

5.2

a

library(vcdExtra)

## Warning: package 'vcdExtra' was built under R version 3.5.1

## Loading required package: vcd

## Warning: package 'vcd' was built under R version 3.5.1

## Loading required package: grid

## Loading required package: gnm

## Warning: package 'gnm' was built under R version 3.5.1

data("AirCrash",package="vcdExtra")
str(AirCrash)

## 'data.frame':    439 obs. of  5 variables:
##  $ Phase     : Factor w/ 5 levels "en route","landing",..: 2 2 2 2 2 2 2 2 2 2 ...
##  $ Cause     : Factor w/ 5 levels "criminal","human error",..: 1 1 1 1 4 4 4 4 4 4 ...
##  $ date      : Date, format: "1993-09-21" "1993-09-22" ...
##  $ Fatalities: int  27 108 125 112 41 19 8 22 14 43 ...
##  $ Year      : int  1993 1993 1996 2002 1993 1993 1994 2000 2002 2004 ...

aircrash.tab <- xtabs(~Phase + Cause, data= AirCrash)
aircrash.tab

##           Cause
## Phase      criminal human error mechanical unknown weather
##   en route       16          63         29      25      24
##   landing         4         114         19      18      55
##   standing        2           0          2       0       0
##   take-off        1          29         24       8       3
##   unknown         0           1          0       1       1

b

mosaic(aircrash.tab, shade=TRUE)

c

mosaic(aircrash.tab, shade = TRUE, labeling_args = list(rot_labels = c (20,20,20)))

5.7

a

data("caith",package="MASS")
caith

##        fair red medium dark black
## blue    326  38    241  110     3
## light   688 116    584  188     4
## medium  343  84    909  412    26
## dark     98  48    403  681    85

library(MASS)
Caith <- as.matrix(caith)
Caith

##        fair red medium dark black
## blue    326  38    241  110     3
## light   688 116    584  188     4
## medium  343  84    909  412    26
## dark     98  48    403  681    85

names(dimnames(Caith))<-c("Eye_Color","Hair_Color")
Caith

##          Hair_Color
## Eye_Color fair red medium dark black
##    blue    326  38    241  110     3
##    light   688 116    584  188     4
##    medium  343  84    909  412    26
##    dark     98  48    403  681    85

b

(Caith.model_of_independence <- loglm(~Eye_Color + Hair_Color, data = Caith, fitted = TRUE))

## Call:
## loglm(formula = ~Eye_Color + Hair_Color, data = Caith, fitted = TRUE)
## 
## Statistics:
##                       X^2 df P(> X^2)
## Likelihood Ratio 1218.314 12        0
## Pearson          1240.039 12        0

c

residuals(Caith.model_of_independence)

##          Hair_Color
## Eye_Color       fair         red    medium       dark     black
##    blue     8.632012 -0.01931269 -2.668179  -5.995870 -3.938806
##    light   11.604525  3.31247206 -1.728753 -12.191903 -6.629923
##    medium  -6.555953 -1.06909312  7.400871  -2.190322 -2.195900
##    dark   -16.185525 -2.76828628 -5.413467  16.277849  8.459859

d

mosaic(Caith.model_of_independence, shade=TRUE)