Created
December 7, 2013 01:53
-
-
Save tukachev/7836355 to your computer and use it in GitHub Desktop.
Младенческая смертность в России (2011-2012) и результаты выборов в госдуму (2011)
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
# Младенческая смертность в России (2011-2012) и результаты выборов в госдуму (2011) | |
# по регионам РФ | |
# Исходные данные: | |
# - официальные результаты выборов | |
# http://ru.wikipedia.org/wiki/Выборы_в_Государственную_думу_(2011)#cite_note-87 | |
# - данные Росстата по младенческой смертности в 2012 году | |
# http://www.gks.ru/free_doc/2012/demo/t3-3.xls | |
# Загружаем необходимые пакеты | |
library("psych") | |
# library("Hmisc") | |
# library("survival") | |
# library("splines") | |
# Функция, которая нам понадобится для наведения красоты в корреляционной матрице :) | |
corstars < - function(x){ | |
require(Hmisc) | |
x <- as.matrix(x) | |
R <- rcorr(x)$r | |
p <- rcorr(x)$P | |
mystars <- ifelse(p < .01, "**|", ifelse(p < .05, "* |", " |")) | |
R <- format(round(cbind(rep(-1.11, ncol(x)), R), 3))[,-1] | |
Rnew <- matrix(paste(R, mystars, sep=""), ncol=ncol(x)) | |
diag(Rnew) <- paste(diag(R), " |", sep="") | |
rownames(Rnew) <- colnames(x) | |
colnames(Rnew) <- paste(colnames(x), "|", sep="") | |
Rnew <- as.data.frame(Rnew) | |
return(Rnew) | |
} | |
# Загружаем данные | |
chmortality2 <- read.delim("~/Rprimer/chmortality2.csv") | |
# Исходный файл: https://dl.dropbox.com/u/33111025/r_examples/chmortality2.csv | |
# Если нужны полные данные, в которых присутствует Чукотка, то | |
# качаем файл: https://dl.dropbox.com/u/33111025/r_examples/chmortality.csv | |
# Присоединяем наши данные, чтобы напрямую работать с переменными | |
attach(chmortality2) | |
#Считаем корреляцию | |
corstars(chmortality2[2:11]) | |
# графики 2 по 2 | |
layout(matrix(c(1,2,3,4), 2, 2, byrow=TRUE)) | |
# Диаграмма рассеяния %Голосов за ЕР и младенческая смертность в 2011 году | |
plot(chmort11~er, main='Глоса за ЕР и детская смертность, 2011', ylab='Млад. смертность на 1000 род.', xlab='ЕР, %голосов') | |
abline(lm(chmort11~er), col='blue') | |
# Диаграмма рассеяния %Голосов за ЕР и младенческая смертность в 2012 году | |
plot(chmort12~er, main='Глоса за ЕР и детская смертность, 2012', ylab='Млад. смертность на 1000 род.', xlab='ЕР, %голосов') | |
abline(lm(chmort12~er), col='blue') | |
# Диаграмма рассеяния %Голосов за Яблоко и младенческая смертность в 2011 году | |
plot(chmort11~yabloko, main='Глоса за Яблоко и детская смертность, 2011', ylab='Млад. смертность на 1000 род.', xlab='Яблоко, %голосов') | |
abline(lm(chmort11~yabloko), col='blue') | |
# Диаграмма рассеяния %Голосов за Яблоко и младенческая смертность в 2012 году | |
plot(chmort12~yabloko, main='Глоса за Яблоко и детская смертность, 2012', ylab='Млад. смертность на 1000 род.', xlab='Яблоко, %голосов') | |
abline(lm(chmort12~yabloko), col='blue') |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment