Load in Some Useful R Packages

library(tidyverse)
library(psych)
library(haven)

Load in the Data

gss <- read_dta("descriptive_gss.dta")

glimpse(gss)
Rows: 2,765
Columns: 16
$ id       <dbl> 2331, 2003, 1221, 2051, 2465, 546, 1291, 732, 303, 2700, 855, 623, 159, 886, 2036, 27…
$ hrs1     <dbl+lbl> NA, NA, NA, NA, 50, 60, 40, 25, NA, 40, 64, 45, 60, 85, NA,  1, NA, 45, NA, 50, 4…
$ marital  <dbl+lbl> 1, 3, 1, 1, 1, 1, 5, 2, 1, 1, 1, 1, 1, 1, 2, 5, 1, 1, 5, 5, 3, 3, 5, 1, 5, 5, 1, …
$ childs   <dbl+lbl> 3, 8, 3, 2, 0, 0, 0, 3, 3, 2, 3, 3, 2, 3, 6, 0, 4, 2, 0, 0, 3, 3, 0, 2, 0, 0, 3, …
$ age      <dbl+lbl> 71, 69, 40, 60, 31, 37, 23, 86, 70, 42, 41, 30, 43, 48, 70, 24, 60, 36, 25, 44, 6…
$ educ     <dbl+lbl> 18, 11, 19, 13, 11, 19, 11, 11, 13, 12, 12, 12, 13, 20, 12, 16, 16, 16, 19, 12, 1…
$ sex      <dbl+lbl> 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, …
$ polviews <dbl+lbl>  4,  4, NA,  5,  4,  2,  1,  4,  4,  3,  6, NA,  5,  5, NA, NA, NA, NA,  4,  4,  …
$ wwwhr    <dbl+lbl> NA, NA,  7,  1,  0,  3, NA, NA, NA,  3,  5, 20,  1,  2, NA,  2,  2,  5,  5, NA, N…
$ trustpeo <dbl+lbl>  4,  1, NA,  1,  2,  2, NA,  1,  4,  1,  2, NA,  5,  2, NA, NA, NA, NA,  2,  2, N…
$ wantbest <dbl+lbl>  2,  4, NA,  2,  2,  4, NA,  4,  1,  2,  2, NA,  1,  3, NA, NA, NA, NA,  2,  2, N…
$ advantge <dbl+lbl>  3,  2, NA,  1,  2,  2, NA,  2,  2,  2,  1, NA,  4,  4, NA, NA, NA, NA,  4,  2, N…
$ goodlife <dbl+lbl>  4, NA, NA, NA, NA,  1,  2,  2, NA, NA, NA, NA,  3, NA, NA, NA, NA, NA,  3, NA,  …
$ deckids  <dbl+lbl> NA, NA,  4, NA, NA, NA, NA, NA, NA, NA, NA,  3, NA, NA, NA,  3,  4,  4, NA, NA, N…
$ strsswrk <dbl+lbl> NA, NA,  5, NA, NA, NA, NA, NA, NA, NA, NA,  2, NA, NA,  3,  2,  3,  2, NA, NA, N…
$ satjob7  <dbl+lbl> NA, NA,  3, NA, NA, NA, NA, NA, NA, NA, NA,  3, NA, NA, NA,  2, NA,  2, NA, NA, N…
str(gss$satjob7)
 dbl+lbl [1:2765] NA, NA,  3, NA, NA, NA, NA, NA, NA, NA, NA,  3, NA, NA, NA,  2, NA,  2, NA, NA, NA...
 @ label       : chr "job satisfaction in general"
 @ format.stata: chr "%37.0g"
 @ labels      : Named num [1:10] 0 1 2 3 4 5 6 7 8 9
  ..- attr(*, "names")= chr [1:10] "nap" "completely satisfied" "very satisfied" "fairly satisfied" ...

Describe the Dataset

describe(gss, fast = TRUE)

Create a Quick Histogram

hist(gss$educ)

Another Histogram

hist(gss$sex, freq = FALSE)

LS0tCnRpdGxlOiAiV2VlayAxIEhhbmRvdXQgLSBJbnRybyB0byBSIGFuZCBSU3R1ZGlvIgphdXRob3I6ICJEci4gQiIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKIyBMb2FkIGluIFNvbWUgVXNlZnVsIFIgUGFja2FnZXMKCmBgYHtyfQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShwc3ljaCkKbGlicmFyeShoYXZlbikKCmBgYAoKIyBMb2FkIGluIHRoZSBEYXRhCmBgYHtyfQpnc3MgPC0gcmVhZF9kdGEoImRlc2NyaXB0aXZlX2dzcy5kdGEiKQoKZ2xpbXBzZShnc3MpCmBgYAoKYGBge3J9CnN0cihnc3Mkc2F0am9iNykKYGBgCgojIERlc2NyaWJlIHRoZSBEYXRhc2V0CmBgYHtyfQpkZXNjcmliZShnc3MsIGZhc3QgPSBUUlVFKQpgYGAKCiMgQ3JlYXRlIGEgUXVpY2sgSGlzdG9ncmFtCmBgYHtyfQpoaXN0KGdzcyRlZHVjKQpgYGAKCiMgQW5vdGhlciBIaXN0b2dyYW0KYGBge3J9Cmhpc3QoZ3NzJHNleCwgZnJlcSA9IEZBTFNFKQoKYGBgCgoKCgoKCg==