हर कोई, गड़बड़ी के लिए खेद है, लेकिन मैं आर के साथ काफी नया हूं, एक महत्वपूर्ण अंतर का सामना कर रहा हूं: मैं रुसिन का एक एनिमेटेड नक्शा बनाना चाहता हूं जिसमें अलग-अलग वर्षों के साथ बेरोजगारी में बदलाव होता है, जैसे। छवि पर आप एक वर्ष के लिए डेटा देख सकते हैं
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")
परिणाम, जिसे मैं प्राप्त करना चाहता हूं, वह यहां एनीमेशन जैसा कुछ है: http://spatial.ly/2011/02/mapping-londons-population-change-2011-2030/ हालांकि, मैंने बहुत कुछ जाना, कई विषयों को पढ़ा में http://stackoverflow.com : निम्नलिखित सहित आर में भूखंड की एक श्रृंखला से एक फिल्म बना रहा है , लेकिन अभी भी सही काम नहीं कर सका।
आपका अग्रिम में ही बहुत धन्यवाद!
मैंने कुछ इस तरह से काम किया है, क्या कोई मुझे बता सकता है कि गलती कहाँ है:
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")
पुन: संपादित करें: क्या कोड के साथ गलत हो रहा है?
—
whuber
चूंकि आपका उदाहरण प्रजनन योग्य नहीं है, इसलिए इसका निवारण करना मुश्किल है। कुछ चीजें 1 से बाहर निकलती हैं) आप एक लूप में स्थानिक परिवर्तन लागू कर रहे हैं, इसलिए आप इसे बार-बार 2 कर रहे हैं) आप "कोशिश" नामक एक ऑब्जेक्ट बना रहे हैं जो कि एक आर फ़ंक्शन 3 भी है) आप वास्तविक कॉलम नामों के माध्यम से पुनरावृति कर सकते हैं। ।, के लिए (i in c ("Var1", "Var2")) जिस तरह से आपके पास वर्तमान में यह कोडित है वह बहुत जटिल है 4) spplot के लिए आपकी कॉल सही नहीं है, आप इसे एक बकवास वेक्टर पास कर रहे हैं।
—
जेफरी इवांस
मुझे इस तरह की समझ न होने के लिए वास्तव में खेद है, लेकिन आर के साथ यह मेरा पहला वास्तविक अनुभव है, मैंने मुख्य प्रश्न में डेटा जोड़ा है, अगर यह आपको परेशान नहीं करेगा तो क्या आप कृपया सुधार के तरीके सुझा सकते हैं जैसे मैं वास्तव में हूं विचारों का हिस्सा है
—
रुविन रफैलोव
