Mọi người, xin lỗi vì đã làm phiền, nhưng tôi còn khá mới với r phải đối mặt với một vấn đề rất quan trọng: Tôi muốn tạo ra một bản đồ hoạt hình của Nga với những thay đổi về tình trạng thất nghiệp với những năm khác nhau. Trên hình ảnh bạn có thể thấy dữ liệu trong một năm
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")Kết quả mà tôi sẵn sàng nhận được là một cái gì đó giống như hoạt hình ở đây: http://spatial.ly/2011/02/mapping-londons-population-change-2011-2030/ Tuy nhiên, tôi đã googled rất nhiều, đọc một số chủ đề trong http://stackoverflow.com bao gồm các mục sau: Tạo phim từ một loạt các cốt truyện trong R , nhưng không thể làm đúng.
cảm ơn bạn trước
Tôi đã đến với một cái gì đó như thế này, bất cứ ai có thể xin vui lòng cho tôi biết đâu là sai lầm:
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 The edit: điều gì đang xảy ra với mã?
—
whuber
Vì ví dụ của bạn không thể tái tạo nên rất khó để khắc phục sự cố. Một vài điều nhảy ra 1) bạn đang áp dụng một biến đổi không gian trong một vòng lặp, vì vậy bạn đang thực hiện nó 2 lần) bạn đang tạo một đối tượng gọi là "thử" cũng là hàm R 3) bạn có thể lặp qua các tên cột thực tế ., for (i in c ("Var1", "Var2")) theo cách mà bạn hiện đang mã hóa rất phức tạp 4) cuộc gọi đến spplot của bạn không chính xác, bạn đang chuyển cho nó một vectơ vô nghĩa.
—
Jeffrey Evans
Tôi thực sự xin lỗi vì không hiểu biết như vậy, nhưng đây là trải nghiệm thực sự đầu tiên của tôi với R, tôi đã thêm dữ liệu vào câu hỏi chính, nếu nó không làm phiền bạn, bạn có thể vui lòng gợi ý các cách để cải thiện vì tôi thực sự chạy ý tưởng
—
Ruvin Rafailov
