这个图怎么做的 求指点

这个图怎么做的 求指点

A piece-by-piece solution using {graphics} only:

1. # Data for plotting
   
2. head(dfPlt)
   
3. #       mnt  hrl              prv
   
4. # 北京 1650 16.9 <U+5317><U+4EAC>
   
5. # 上海 1620 14.0 <U+4E0A><U+6D77>
   
6. # 广东 1550 15.0 <U+5E7F><U+4E1C>
   
7. # 新疆 1520 15.2 <U+65B0><U+7586>
   
8. # 天津 1500 15.0 <U+5929><U+6D25>
   
9. # 江苏 1480 13.0 <U+6C5F><U+82CF>
   
10. 
    
11. # Function to add a fan-shape with text inside and outside in the plot
    
12. pltFan <-
    
13.   function(cnt.arc, arc, r.in, r.out, text.in = character(0), text.out = character(0), text.in.cex = 0.5, text.out.cex = text.in.cex, ...)
    
14.   {
    
15.     # Fan-shape
    
16.     polygon(c(r.in*cos(cnt.arc - arc/2), r.in*cos(cnt.arc + arc/2), r.out*cos(cnt.arc + arc/2), r.out*cos(cnt.arc - arc/2)), c(r.in*sin(cnt.arc - arc/2), r.in*sin(cnt.arc + arc/2), r.out*sin(cnt.arc + arc/2), r.out*sin(cnt.arc - arc/2)), ...)
    
17.     # Text inside
    
18.     text((r.in - par()$cxy[2]*text.in.cex)*cos(cnt.arc), (r.in - par()$cxy[2]*text.in.cex)*sin(cnt.arc), text.in, adj = c(0.5, 0.5), srt = cnt.arc/2/pi*360 - 90, cex = text.in.cex)
    
19.     # Text outside
    
20.     text((r.out + par()$cxy[2]*text.out.cex)*cos(cnt.arc), (r.out + par()$cxy[2]*text.out.cex)*sin(cnt.arc), text.out, adj = c(0.5, 0.5), srt = cnt.arc/2/pi*360 - 90, cex = text.out.cex)
    
21.   }
    
22. 
    
23. # Number of records to plot
    
24. iNbar <- nrow(dfPlt)
    
25. # Optionally, map additional font - Windows only
    
26. windowsFonts(KT = windowsFont("KaiTi"))
    
27. # Plot in .png format
    
28. png("salary.png", width = 960, height = 960)
    
29. # Maximizing the plot area
    
30. par(mar = c(0, 0, 0, 0))
    
31. # Call plot.new with a few settings
    
32. plot(c(-2, 2), c(-2, 2), type = "n", axes = F, xlab = "", ylab = "", asp = 1, xlim = c(-2, 2), ylim = c(-2, 2))
    
33. # Do the plot
    
34. sapply( 0:(iNbar - 1), function(i)
    
35.   {
    
36.   # Text inside
    
37.   pltFan(cnt.arc = pi/2 - (i*2*pi/(iNbar + 1) + 2*pi/(iNbar + 1)*(1/4 + 1/8 + 1/4)/2), arc = 2*pi/(iNbar + 1)*(1/4 + 1/8 + 1/4), r.in = 1, r.out = 2, text.in = dfPlt[i + 1, "prv"], text.out = NULL, col = NULL, border = NA, text.in.cex = 1)
    
38.   # Monthly data
    
39.   pltFan(cnt.arc = pi/2 - (i*2*pi/(iNbar + 1) + 2*pi/(iNbar + 1)/8), arc =  2*pi/(iNbar + 1)/4, r.in = 1, r.out = dfPlt[i + 1, "mnt"]/1000, text.in = NULL, text.out = as.character(dfPlt[i + 1, "mnt"]), text.out.cex = 1, col = "orange", border = NA)
    
40.   # Hourly data
    
41.   pltFan(cnt.arc = pi/2 - (i*2*pi/(iNbar + 1) + 2*pi/(iNbar + 1)*(2 + 1 + 1)/8), arc =  2*pi/(iNbar + 1)/4, r.in = 1, r.out = dfPlt[i + 1, "hrl"]/20 + 1, text.out.cex = 1, text.in = NULL, text.out = as.character(dfPlt[i + 1, "hrl"]), col = "yellow", border = NA)
    
42.   } )
    
43. # Legend and other texts in the plot
    
44. legend(2, -1.5, c("月最低工资标准", "小时最低工资标准"), col = c("orange", "yellow"), bty = "n", cex = 1.25, pt.cex = 1.2, pch = c(15, 15), xjust = 1, yjust = 1) -> lgnd
    
45. # If no font mapping is done, remove the family option
    
46. text(lgnd$rect$left, lgnd$rect$top - lgnd$rect$h, "注：数字为截至2014年\n5月1日的执行标准", adj = c(0, 1), cex = 1.25, family = "KT")
    
47. text(-2.1, 1.75, "全国各省（区、直辖市）\n2014年最低工资标准", cex = 1.75, adj = c(0, 1))
    
48. text(-2.1, 1.75 - par()$cxy[2]*4, "（单位：元）", cex = 1.25, adj = c(0, 1))
    
49. dev.off()

复制代码

Result:

不错的设计，就是图的中间空的太多，占地儿啊。希望这个图对你有用：http://excelpro.blog.sohu.com/262660652.html

有点类似circos图，这种属于简单的了，circos软件上简单些
曾经纯R画过，借助极坐标，然后画一圈的矩形就可以了

yywan0913 发表于 2014-5-25 16:14
有点类似circos图，这种属于简单的了，circos软件上简单些
曾经纯R画过，借助极坐标，然后画一圈的矩形就可 ...

可以把你以前用R画的数据和代码分享下不

zyz0329 发表于 2014-5-25 17:04
可以把你以前用R画的数据和代码分享下不

不知道扔哪去了   你提供数据的话，我还是能给你搞定的

zhangyangsmith 发表于 2014-5-27 18:54
A piece-by-piece solution using {graphics} only:Result:

这图画的真不错
dfPlt这个数据可以贴上来吗？谢谢

zhangyangsmith 发表于 2014-5-25 14:49
A piece-by-piece solution using {graphics} only:Result:

可以把全部数据贴出来吗 谢谢  我在stata区也有悬赏 你可以贴在那 论坛币给你

zyz0329 发表于 2014-5-28 19:46
可以把全部数据贴出来吗 谢谢  我在stata区也有悬赏 你可以贴在那 论坛币给你

I grabbed the data from the following link:
http://wenku.baidu.com/view/8914ed206bd97f192279e95c.html
You need to do a bit reformatting though. I simply selected the highest standard within each province as I thought the plot is just a demo...

zhangyangsmith 发表于 2014-5-30 06:27
I grabbed the data from the following link:
http://wenku.baidu.com/view/8914ed206bd97f192279e95c. ...

好的 但是我还是没明白啊 里面代码复制后有很多错误  可能是格式问题 但是我改不来  能不能发个txt或者word格式的运行命令附件  真的很感谢 谢谢