Skip to content

Instantly share code, notes, and snippets.

@chrishanretty

chrishanretty/cses_import.R

Created June 23, 2015 09:01
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save chrishanretty/82bea7e3dd0119cbba5b to your computer and use it in GitHub Desktop.
Save chrishanretty/82bea7e3dd0119cbba5b to your computer and use it in GitHub Desktop.
Download ZIP
Greece turnout analysis
Raw
cses_import.R
dat <- read.fwf("cses4.dat",
widths = c(13, 14, 35, 8, 8, 18, 4, 1, 4, 10, 18, 18, 18, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 4, 2, 2, 4, 1, 1,
6, 1, 2, 2, 4, 3, 3, 2, 4, 1, 2, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 2, 3, 1, 1, 1, 1, 2, 2, 2, 1, 1, 4, 3, 2, 3, 3, 1, 5, 3, 4,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 1, 2, 2, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 1, 1, 2, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1,
1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 4, 3, 16,
16, 16, 17, 17, 17, 18, 17, 18, 16, 16, 16, 16, 16, 16, 16, 17, 17, 17, 16, 16, 17, 17, 16, 16, 17, 17, 17,
16, 16, 16, 16, 17, 16, 16, 17, 16, 15, 16, 16, 16, 16, 17, 17, 3, 3, 16, 16, 17, 17, 18, 16, 16, 16, 3, 16,
16, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 1, 1, 1, 2, 2, 4, 2, 2, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2,
2, 2, 2, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 16, 2, 2, 2, 1, 1, 1, 1, 3, 3,
3, 2, 2, 2, 3, 1, 2, 3, 1, 1, 1, 1, 1, 1, 3, 16, 2, 3, 16, 2, 3, 3, 2, 3, 16, 1, 3, 18, 18, 18, 16, 16, 16, 18, 18,
18, 17, 17, 17, 9, 9, 9, 17, 16, 16, 2, 18, 16, 16, 16, 16, 16, 16, 17, 17, 17, 16, 16, 16, 16, 16, 16, 1, 1,
1, 16, 16, 16, 16),
colClasses = c("character", "character", "character", "numeric", "character", "character", "character",
"numeric", "numeric", "character", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric", "numeric", "numeric", "numeric"),
col.names = c("D1001","D1002_VER","D1002_DOI","D1003","D1004","D1005","D1006","D1007","D1008",
"D1009","D1010_1","D1010_2","D1010_3","D1011_1","D1011_2","D1011_3","D1012_1",
"D1012_2","D1012_3","D1013","D1014_1","D1014_2","D1014_3","D1015","D1016","D1017",
"D1018","D1019","D1020","D1021","D1022","D1023","D1024","D1025","D1026","D1027",
"D1028","D1029","D1030","D2001_M","D2001_Y","D2002","D2003","D2004","D2005","D2006",
"D2007","D2008","D2009","D2010","D2011","D2012","D2013","D2014","D2015","D2016",
"D2017","D2018","D2019","D2020","D2021","D2022","D2023","D2024","D2025","D2026",
"D2027","D2028","D2029","D2030","D2031","D2032","D2033","D2034","D3001_1","D3001_2",
"D3001_3","D3001_4","D3001_5","D3001_6","D3001_7","D3001_8","D3002","D3003_1",
"D3003_2","D3003_3","D3004","D3005_PR_1","D3005_PR_2","D3005_LH","D3005_UH",
"D3006_PR_1","D3006_PR_2","D3006_LH_PL","D3006_LH_DC","D3006_LH_PF",
"D3006_UH_PL","D3006_UH_DC","D3006_UH_PF","D3007_PR_1","D3007_PR_2",
"D3007_LH","D3007_UH","D3008_PR_1","D3008_PR_2","D3008_LH_PL","D3008_LH_DC",
"D3008_UH_PL","D3008_UH_DC","D3009","D3010","D3011_A","D3011_B","D3011_C",
"D3011_D","D3011_E","D3011_F","D3011_G","D3011_H","D3011_I","D3012_A","D3012_B",
"D3012_C","D3012_D","D3012_E","D3012_F","D3012_G","D3012_H","D3012_I","D3013_A",
"D3013_B","D3013_C","D3013_D","D3013_E","D3013_F","D3013_G","D3013_H","D3013_I",
"D3014","D3015_A","D3015_B","D3015_C","D3015_D","D3015_E","D3015_F","D3015_G",
"D3015_H","D3015_I","D3016","D3017","D3018_1","D3018_2","D3018_3","D3018_4",
"D3019","D3020_1","D3020_2","D3020_3","D3020_4","D3020_5","D3020_6","D3021_1",
"D3021_2","D3021_3","D3021_4","D3021_5","D3021_6","D3021_7","D3021_8","D3021_9",
"D3021_10","D3022","D3023_1","D3023_2","D3023_3","D3023_4","D3023_5","D3023_6",
"D3024","D3025_1_A","D3025_2_A","D3025_3_A","D3025_4_A","D3025_1_B","D3025_2_B",
"D3025_3_B","D3025_4_B","D3026","D3027_1","D3027_2","D3027_3","D3027_4","D3028_1",
"D3028_2","D4001","D4002","D4003","D4004_A","D4004_B","D4004_C","D4004_D","D4004_E",
"D4004_F","D4004_G","D4004_H","D4004_I","D4005","D5001_A","D5001_B","D5001_C",
"D5001_D","D5001_E","D5001_F","D5001_G","D5001_H","D5001_I","D5002_A","D5002_B",
"D5002_C","D5002_D","D5002_E","D5002_F","D5002_G","D5002_H","D5002_I","D5003_A",
"D5003_B","D5003_C","D5003_D","D5003_E","D5003_F","D5003_G","D5003_H","D5003_I",
"D5004_A","D5004_B","D5004_C","D5004_D","D5004_E","D5004_F","D5004_G","D5004_H",
"D5004_I","D5005_A","D5005_B","D5005_C","D5005_D","D5005_E","D5005_F","D5005_G",
"D5005_H","D5005_I","D5006_1","D5006_2","D5007","D5008","D5009_A","D5009_B",
"D5009_C","D5009_D","D5009_E","D5009_F","D5009_G","D5009_H","D5009_I","D5010",
"D5011","D5012","D5013_A","D5013_B","D5013_C","D5013_D","D5013_E","D5013_F",
"D5013_G","D5013_H","D5013_I","D5014","D5015","D5016_A","D5016_B","D5016_C",
"D5016_D","D5016_E","D5016_F","D5016_G","D5016_H","D5016_I","D5017_A","D5017_B",
"D5017_C","D5017_D","D5017_E","D5017_F","D5017_G","D5017_H","D5017_I","D5018",
"D5018_A","D5018_B","D5018_C","D5018_D","D5018_E","D5018_F","D5018_G","D5018_H",
"D5018_I","D5019_1","D5019_2","D5019_3","D5019_4","D5019_5","D5020","D5021","D5022",
"D5023_1","D5023_2","D5023_3","D5024_1","D5024_2","D5024_3","D5025","D5026","D5027",
"D5028","D5029","D5030","D5031","D5032","D5033","D5034","D5035","D5036","D5037",
"D5038_1","D5038_2","D5038_3","D5038_4","D5039_1","D5039_2","D5039_3","D5039_4",
"D5040_1","D5040_2","D5040_3","D5040_4","D5041_1","D5041_2","D5041_3","D5041_4",
"D5042_1","D5042_2","D5042_3","D5042_4","D5043_1","D5043_2","D5043_3","D5043_4",
"D5044_1","D5044_2","D5044_3","D5044_4","D5045_1","D5045_2","D5045_3","D5045_4",
"D5046_1","D5046_2","D5046_3","D5046_4","D5047_1","D5047_2","D5047_3","D5047_4",
"D5050_1","D5050_2","D5050_3","D5051_1","D5051_2","D5051_3","D5052","D5054",
"D5055","D5056","D5057","D5058","D5059","D5060","D5061","D5062","D5063","D5064",
"D5065","D5066","D5067","D5068","D5069","D5070","D5071","D5072","D5073","D5074",
"D5075","D5080_1","D5080_2","D5080_3","D5081_1","D5081_2","D5081_3","D5082_1",
"D5082_2","D5082_3","D5083_1","D5083_2","D5083_3","D5084_1","D5084_2","D5084_3",
"D5085_1","D5085_2","D5085_3","D5086","D5088","D5089_1","D5089_2","D5089_3",
"D5091_1","D5091_2","D5091_3","D5092_1","D5092_2","D5092_3","D5093_1","D5093_2",
"D5093_3","D5095","D5096","D5097","D5099","D5100","D5101","D5102","D5103","D5104","D5105"),
n = 33030)
### Create resp age
dat$dob_y <- dat$D2001_Y
dat$dob_y[which(dat$dob_y > 9900)] <- NA
dat$age <- dat$D1008 - dat$dob_y
### Turnout
dat$turned_out <- ((dat$D3005_LH == 1) |
(dat$D3005_UH == 1) |
(dat$D3005_PR_1 == 1) |
(dat$D3005_PR_2 == 1))
dat$turned_out <- as.numeric(dat$turned_out)
dat$turned_out[which(dat$D3005_LH == 9)] <- NA
### country
dat$iso3c <- substr(dat$D1004, 1, 3)
library(ggplot2)
library(plyr)
png(width = 800, height = 550, file = "greece_turnout.png")
ggplot(subset(dat, iso3c == "GRC"), aes(x = age, y = turned_out)) +
geom_point(alpha =0.5, position = position_jitter(width = 0, height = 0.05)) +
stat_smooth(aes(weight = D1010_3)) +
theme_bw()
dev.off()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment