Assignment 1

This is a sample R script, formatted so that you can use File -> Compile Report or Ctrl+Shift+K to run analyses and turn it into an HTML document, or a Word document. The general rule is that everything is just R code, but comments that begin with #' are treated as text in the output.

I use section headers, #' ## Section to break the script into sections and provide other explanation.

Load packages

Generally, I find it useful to load all packages I’m using at the beginning of a script.

library(vcdExtra)   # for data and graphics

Load data

Illustrate the Abortion data. For details: help(Abortion, package = "vcdExtra").

data(Abortion, package = "vcdExtra") 

Look at what’s there. The Abortion table has 1100 observations.

str(Abortion)

##  'table' num [1:2, 1:2, 1:2] 171 152 138 167 79 148 112 133
##  - attr(*, "dimnames")=List of 3
##   ..$ Sex             : chr [1:2] "Female" "Male"
##   ..$ Status          : chr [1:2] "Lo" "Hi"
##   ..$ Support_Abortion: chr [1:2] "Yes" "No"

Print the table nicely

ftable() flattens the table to a two-way display, with some variables assigned to the rows and others assigned to the columns.

ftable(Abortion)

##               Support_Abortion Yes  No
## Sex    Status                         
## Female Lo                      171  79
##        Hi                      138 112
## Male   Lo                      152 148
##        Hi                      167 133

Make a plot

I want to make a fourfold plot of Status by Support_Abortion each level of Sex

fourfold(aperm(Abortion, 3:1))