Skip to content

Instantly share code, notes, and snippets.

@tukachev
Created December 7, 2013 01:53
Show Gist options
  • Save tukachev/7836355 to your computer and use it in GitHub Desktop.
Save tukachev/7836355 to your computer and use it in GitHub Desktop.
Младенческая смертность в России (2011-2012) и результаты выборов в госдуму (2011)
# Младенческая смертность в России (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