todo mundo, desculpe pela perturbação, mas estou sendo bastante novo com r enfrentou uma dificuldade crucial: quero criar um mapa animado de Russin com mudanças no desemprego com anos diferentes, como. Na imagem você pode ver dados de um ano
require(sp)
require(maptools)
require(RColorBrewer)
require(rgdal)
rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData")
print(load(rus))
unempl <- read.delim2(file="C:\\unempl1.txt", header = TRUE,
sep = ";",quote = "", dec=",", stringsAsFactors=F)
gadm_names <-gadm$NAME_1
total <- length(gadm_names)
pb <- txtProgressBar(min = 0, max = total, style = 3)
order <- vector()
for (i in 1:total){
order[i] <- agrep(gadm_names[i], unempl$region,
max.distance = 0.2)[1]
setTxtProgressBar(pb, i) # update progress bar
}
col_no <- as.factor(as.numeric(cut(unempl$data[order],
c(0,2.5,5,7.5,10,15,100))))
levels(col_no) <- c("<2,5%", "2,5-5%", "5-7,5%",
"7,5-10%", "10-15%", ">15%")
gadm$col_no <- col_no
myPalette<-brewer.pal(6,"Purples")
proj4.str <- CRS("+init=epsg:3413 +lon_0=105")
gadm.prj <- spTransform(gadm, proj4.str)
spplot(gadm.prj, "col_no", col=grey(.9), col.regions=myPalette,
main="Unemployment in Russia by region")
O resultado que desejo obter é algo como animação aqui: http://spatial.ly/2011/02/mapping-londons-population-change-2011-2030/ No entanto, pesquisei bastante no Google, li vários temas em http://stackoverflow.com, incluindo o seguinte: Criando um filme a partir de uma série de gráficos em R , mas ainda não conseguia fazer a coisa certa.
Agradeço antecipadamente!
Eu vim com algo assim, alguém pode me dizer onde está o erro:
require(animation)
require(sp)
require(RColorBrewer)
require(classInt)
require(rgdal)
rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData")
print(load(rus))
unempl1 <- read.delim2(file="C:\\unempl11.txt", header = TRUE,
sep = ";",quote = "", dec=",", stringsAsFactors=F)
unempl2<- read.delim2(file="C:\\unempl12.txt", header = TRUE,
sep = ";",quote = "", dec=",", stringsAsFactors=F)
gadm_names <-gadm$NAME_1
total <- length(gadm_names)
pb <- txtProgressBar(min = 0, max = total, style = 3)
order <- vector()
for (i in 1:total){
order[i] <- agrep(gadm_names[i], unempl1$region,
max.distance = 0.2)[1]
setTxtProgressBar(pb, i) # update progress bar
}
for (l in 1:total){
order[l] <- agrep(gadm_names[l], unempl2$region,
max.distance = 0.2)[1]
setTxtProgressBar(pb, i) # update progress bar
}
col_no_1 <- as.factor(as.numeric(cut(unempl1$data[order],
c(0,2.5,5,7.5,10,15,100))))
col_no_2<- as.factor(as.numeric(cut(unempl2$data[order],
c(0,2.5,5,7.5,10,15,100))))
saveHTML(
for(k in 1:2) {
try<-get(paste("col_no_", k, sep = ""))
levels(try) <- c("<2,5%", "2,5-5%", "5-7,5%",
"7,5-10%", "10-15%", ">15%")
gadm$col_no <- try
myPalette<-brewer.pal(6,"Purples")
proj4.str <- CRS("+init=epsg:3413 +lon_0=105")
gadm.prj <- spTransform(gadm, proj4.str)
spplot(gadm.prj, "col_no", col=grey(.9), col.regions=myPalette,
main="Unemployment in Russia by region")
},img.name = "map", htmlfile = "unrus2.html")
Re A edição: o que está errado com o código?
—
whuber
Como o seu exemplo não é reproduzível, é difícil solucionar problemas. Algumas coisas saltam para fora 1) você está aplicando uma transformação espacial em um loop, então você está fazendo isso repetidamente 2) você está criando um objeto chamado "try", que também é uma função R 3) você pode iterar através dos nomes reais das colunas, ou seja ., para (i em c ("Var1", "Var2")) a maneira como você o codifica atualmente é muito complicada 4) sua chamada ao spplot não está correta, você está passando um vetor sem sentido para ela.
—
Jeffrey Evans
Lamento muito por não ter entendido, mas esta é minha primeira experiência real com R, adicionei os dados na pergunta principal, se não incomodá-lo, por favor, sugere maneiras de melhorar como eu realmente correu ut de idéias
—
Ruvin Rafailov
