library(ggplot2)
library(ggthemes)
library(babynames)
library(Zelig)
library(ggrepel)
library(HistData)
library(tidyverse)
library(car)
data("Vocab")
options(dplyr.show_progress = FALSE)

Explaining Data

My data is looking into how number of vocabulary words correct in a 10-word test is effected based on number of years of education completed as well as sex. I wanted to use ggplots to see whether or not number of years in education completed had a direct effect on number of vocabulary words correct. The answers come from over 21,000 respondents.

Indepdent Variables

Education and Sex. Education is the number of years in school completed. 12 years completed is equivilent to a high school education. 12 years to 16 years is a college degree. Above that is masters levels education and on. Sex is Male or Female.

Dependent Variables

Vocabulary words correct. Vocabulary words correct means the correct number of vocabulary words in a ten question test.

head(Vocab)

Does Education Effect Number of Words Correct?

g <- ggplot(Vocab, mapping = aes(x = education, y = vocabulary))
g3 <- g + geom_density2d()
g3

g4 <- g + geom_smooth()
g4

g5 <- g1 + geom_smooth()
g5

g6 <- g + geom_point() + stat_smooth(method = "lm")
g6

g8 <- g + geom_point() + stat_quantile()
g8

g9 <- g + geom_point() + stat_summary(fun.data = "mean_cl_boot", colour = "red")
g9

g10 <- g + geom_point() + stat_summary(fun.y = "median", colour = "red", geom = "line")
g10

Do Males or Females get more words correct?

ggplot(Vocab,
  aes(x = factor("Male", "Female"), fill = sex) ) +
  geom_bar() +
  coord_polar(theta = "y") +
  scale_x_discrete("")

ggplot(Vocab, aes(vocabulary, fill = sex ) ) +
  geom_bar()

ggplot(Vocab, aes(vocabulary) ) +
  geom_bar() + facet_grid(sex ~ .)

ggplot(Vocab,
  aes(vocabulary, ..count.. ) ) +
  geom_point(stat = "count", size = 4) + 
  coord_flip()+
  facet_grid( sex ~ . )

Conclusion

From my short analysis and visualizations - we can clearly see that as number of years in education completed increases, so does the number of vocabulary words correct in a 10 word test. We can also see from the visuals above that Females do get more vocabulary words correct then their male counterparts.

LS0tDQp0aXRsZTogIkhvbWV3b3JrIDkgLSBWb2NhYnVsYXJ5Ig0Kb3V0cHV0OiBodG1sX25vdGVib29rDQotLS0NCg0KYGBge3J9DQpsaWJyYXJ5KGdncGxvdDIpDQpsaWJyYXJ5KGdndGhlbWVzKQ0KbGlicmFyeShiYWJ5bmFtZXMpDQpsaWJyYXJ5KFplbGlnKQ0KbGlicmFyeShnZ3JlcGVsKQ0KbGlicmFyeShIaXN0RGF0YSkNCmxpYnJhcnkodGlkeXZlcnNlKQ0KbGlicmFyeShjYXIpDQpkYXRhKCJWb2NhYiIpDQoNCm9wdGlvbnMoZHBseXIuc2hvd19wcm9ncmVzcyA9IEZBTFNFKQ0KYGBgDQoNCiMjIEV4cGxhaW5pbmcgRGF0YQ0KKk15IGRhdGEgaXMgbG9va2luZyBpbnRvIGhvdyBudW1iZXIgb2Ygdm9jYWJ1bGFyeSB3b3JkcyBjb3JyZWN0IGluIGEgMTAtd29yZCB0ZXN0IGlzIGVmZmVjdGVkIGJhc2VkIG9uIG51bWJlciBvZiB5ZWFycyBvZiBlZHVjYXRpb24gY29tcGxldGVkIGFzIHdlbGwgYXMgc2V4LiBJIHdhbnRlZCB0byB1c2UgZ2dwbG90cyB0byBzZWUgd2hldGhlciBvciBub3QgbnVtYmVyIG9mIHllYXJzIGluIGVkdWNhdGlvbiBjb21wbGV0ZWQgaGFkIGEgZGlyZWN0IGVmZmVjdCBvbiBudW1iZXIgb2Ygdm9jYWJ1bGFyeSB3b3JkcyBjb3JyZWN0LiBUaGUgYW5zd2VycyBjb21lIGZyb20gb3ZlciAyMSwwMDAgcmVzcG9uZGVudHMuKg0KDQojIyBJbmRlcGRlbnQgVmFyaWFibGVzDQoqRWR1Y2F0aW9uIGFuZCBTZXguIEVkdWNhdGlvbiBpcyB0aGUgbnVtYmVyIG9mIHllYXJzIGluIHNjaG9vbCBjb21wbGV0ZWQuIDEyIHllYXJzIGNvbXBsZXRlZCBpcyBlcXVpdmlsZW50IHRvIGEgaGlnaCBzY2hvb2wgZWR1Y2F0aW9uLiAxMiB5ZWFycyB0byAxNiB5ZWFycyBpcyBhIGNvbGxlZ2UgZGVncmVlLiBBYm92ZSB0aGF0IGlzIG1hc3RlcnMgbGV2ZWxzIGVkdWNhdGlvbiBhbmQgb24uIFNleCBpcyBNYWxlIG9yIEZlbWFsZS4qDQoNCiMjIERlcGVuZGVudCBWYXJpYWJsZXMNCipWb2NhYnVsYXJ5IHdvcmRzIGNvcnJlY3QuIFZvY2FidWxhcnkgd29yZHMgY29ycmVjdCBtZWFucyB0aGUgY29ycmVjdCBudW1iZXIgb2Ygdm9jYWJ1bGFyeSB3b3JkcyBpbiBhIHRlbiBxdWVzdGlvbiB0ZXN0LioNCg0KYGBge3J9DQpoZWFkKFZvY2FiKQ0KYGBgDQoNCiMjIERvZXMgRWR1Y2F0aW9uIEVmZmVjdCBOdW1iZXIgb2YgV29yZHMgQ29ycmVjdD8NCmBgYHtyfQ0KZyA8LSBnZ3Bsb3QoVm9jYWIsIG1hcHBpbmcgPSBhZXMoeCA9IGVkdWNhdGlvbiwgeSA9IHZvY2FidWxhcnkpKQ0KYGBgDQoNCmBgYHtyfQ0KZzMgPC0gZyArIGdlb21fZGVuc2l0eTJkKCkNCmczDQpgYGANCg0KDQpgYGB7cn0NCmc0IDwtIGcgKyBnZW9tX3Ntb290aCgpDQpnNA0KYGBgDQoNCmBgYHtyfQ0KZzUgPC0gZzEgKyBnZW9tX3Ntb290aCgpDQpnNQ0KYGBgDQoNCmBgYHtyfQ0KZzYgPC0gZyArIGdlb21fcG9pbnQoKSArIHN0YXRfc21vb3RoKG1ldGhvZCA9ICJsbSIpDQpnNg0KYGBgDQoNCmBgYHtyfQ0KZzggPC0gZyArIGdlb21fcG9pbnQoKSArIHN0YXRfcXVhbnRpbGUoKQ0KZzgNCmBgYA0KDQpgYGB7cn0NCmc5IDwtIGcgKyBnZW9tX3BvaW50KCkgKyBzdGF0X3N1bW1hcnkoZnVuLmRhdGEgPSAibWVhbl9jbF9ib290IiwgY29sb3VyID0gInJlZCIpDQpnOQ0KYGBgDQoNCmBgYHtyfQ0KZzEwIDwtIGcgKyBnZW9tX3BvaW50KCkgKyBzdGF0X3N1bW1hcnkoZnVuLnkgPSAibWVkaWFuIiwgY29sb3VyID0gInJlZCIsIGdlb20gPSAibGluZSIpDQpnMTANCmBgYA0KDQojIyBEbyBNYWxlcyBvciBGZW1hbGVzIGdldCBtb3JlIHdvcmRzIGNvcnJlY3Q/DQpgYGB7cn0NCmdncGxvdChWb2NhYiwNCiAgYWVzKHggPSBmYWN0b3IoIk1hbGUiLCAiRmVtYWxlIiksIGZpbGwgPSBzZXgpICkgKw0KICBnZW9tX2JhcigpICsNCiAgY29vcmRfcG9sYXIodGhldGEgPSAieSIpICsNCiAgc2NhbGVfeF9kaXNjcmV0ZSgiIikNCmBgYA0KDQpgYGB7cn0NCmdncGxvdChWb2NhYiwgYWVzKHZvY2FidWxhcnksIGZpbGwgPSBzZXggKSApICsNCiAgZ2VvbV9iYXIoKQ0KYGBgDQoNCmBgYHtyfQ0KZ2dwbG90KFZvY2FiLCBhZXModm9jYWJ1bGFyeSkgKSArDQogIGdlb21fYmFyKCkgKyBmYWNldF9ncmlkKHNleCB+IC4pDQpgYGANCg0KYGBge3J9DQpnZ3Bsb3QoVm9jYWIsDQogIGFlcyh2b2NhYnVsYXJ5LCAuLmNvdW50Li4gKSApICsNCiAgZ2VvbV9wb2ludChzdGF0ID0gImNvdW50Iiwgc2l6ZSA9IDQpICsgDQogIGNvb3JkX2ZsaXAoKSsNCiAgZmFjZXRfZ3JpZCggc2V4IH4gLiApDQpgYGANCg0KIyMgQ29uY2x1c2lvbg0KKkZyb20gbXkgc2hvcnQgYW5hbHlzaXMgYW5kIHZpc3VhbGl6YXRpb25zIC0gd2UgY2FuIGNsZWFybHkgc2VlIHRoYXQgYXMgbnVtYmVyIG9mIHllYXJzIGluIGVkdWNhdGlvbiBjb21wbGV0ZWQgaW5jcmVhc2VzLCBzbyBkb2VzIHRoZSBudW1iZXIgb2Ygdm9jYWJ1bGFyeSB3b3JkcyBjb3JyZWN0IGluIGEgMTAgd29yZCB0ZXN0LiBXZSBjYW4gYWxzbyBzZWUgZnJvbSB0aGUgdmlzdWFscyBhYm92ZSB0aGF0IEZlbWFsZXMgZG8gZ2V0IG1vcmUgdm9jYWJ1bGFyeSB3b3JkcyBjb3JyZWN0IHRoZW4gdGhlaXIgbWFsZSBjb3VudGVycGFydHMuKg0K