Lab session 0
Introduction to R Practical
Heather Turner
12 September 2016
-
Download
duncan.txt
from the course website. Then, in RStudio, click on the Import Dataset button on the Environment tab in the top-right pane. Select From Local File and select theduncan.txt
file to import and view the data.The data are from the 1971 census of Canada, with the following variables
- occupation
-
named occupation
- education
-
average education in years
- income
-
average income in dollars
- women
-
percentage women in occupation
- prestige
-
prestige score for occupation, from earlier survey
- census
-
Canadian census occupation code
- type
-
type of occupation: bc (blue collar), prof (professional, mangerial and technical) and wc (white collar)
Open a new script file to work in, save it as a
.R
file. (Just give it any file name you want, ending in .R) -
Use
summary
to get a quick summary of the variables. You will see thattype
has missing values. The commandis.na(x)
will return a logical vector indicating whether each element ofx
isNA
(TRUE
) or not (FALSE
). Look up the help for the functionsubset
and find out how to use such a logical vector to obtain the subset of the data for whichtype
isNA
. In summaries involvingtype
these observations will be silently excluded. -
Using
min
, assign the minimum proportion of women to a name. The commandx == a
will return a logical vector indicating whether each element ofx
is equal toa
(TRUE
) or not (FALSE
). Usesubset
to obtain the rows of data where the proportion of women is equal to the minimum value. Repeat the process to obtain rows corresponding to the maximum. -
Create a histogram and then a density plot of the prestige variable and compare the output. Use the code completion tools in RStudio to look at the second argument of
hist
and try modifying. -
Create a boxplot of years of education by occupation type.
-
Create a scatterplot of prestige against income. Look at the help for
log
and then create a plot of prestige against log income with base 10.