Младенческая смертность в России (2011-2012) и результаты выборов в госдуму (2011)

childmort.R

# Младенческая смертность в России (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')