1000论坛币求作一个图

tmdxyz

5834

收藏 2017-05-22

悬赏 1000 个论坛币已解决

前段时间看到“柳叶刀”杂志上有一篇文章（原文献见附件），其中的一个图（Figure 3）很有意思。感觉是用ggplot做的（stack column + 极坐标）。

为此我假设了一套数据（见附件），试图仿做一下（见附图）。

只可惜
（1）画不出最外头那一圈代表疾病的分类，“A,B,C...”
（2）中间那个空白的圆圈怎么画
（3）那一圈缺刻，里头的年份，不知怎么画进去
（4）那些疾病种类的文本，字号能不能缩小一些，且旋转一定的角度？因为若不然，下方有些就会显示不全。

（5）整个图形的颜色填充，能不能设定为浅色？我附图太暗了

本人初学R，实在不会，故此悬赏1000个论坛币。请您把根据您的代码所画出来的图贴上，我若觉得确实可行，就将购买您设为1000论坛币的代码。

2.jpeg

原图尺寸 311.94 KB

1.png

原图尺寸 141.21 KB

data.rar

大小:634.05 KB

马上下载

本附件包括：

data.xlsx

paper.pdf

大小:715.89 KB

马上下载

最佳答案

cheetahfly 查看完整内容

水平所限，逐渐试错拼凑出的来的。首先，数据放在变量a中，数据变形的过程就不啰嗦了，最终变成了tidy data，共四个变量：“type”，“name”，“year”，“value”，其中前两个是factor，“year”是integer。最后的效果如下：另外，对于避免字重合的问题，我首先尝试着旋转45度，发现必然会有某个位置上的字存在着重合的风险，而且影响阅读体验，因此，我的建议是将比较短的标签放在正底部附近，人为调整一下顺序，就可 ...

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

cheetahfly

2017-5-22 10:20:29

水平所限，逐渐试错拼凑出的来的。
首先，数据放在变量a中，数据变形的过程就不啰嗦了，最终变成了tidy data，共四个变量：“type”，“name”，“year”，“value”，其中前两个是factor，“year”是integer。

复制代码

最后的效果如下：

另外，对于避免字重合的问题，我首先尝试着旋转45度，发现必然会有某个位置上的字存在着重合的风险，而且影响阅读体验，因此，我的建议是将比较短的标签放在正底部附近，人为调整一下顺序，就可以解决该问题了。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zerofung

2017-5-22 16:18:21

复制代码

附件列表

TIM截图20170522161819.png

原图尺寸 48.84 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zerofung

2017-5-22 16:32:57

代码里的美元符号显示不好。。。。有些细节其实还是用PS擦掉或添加，比花心思去写代码更划算，做事毕竟要考虑时间投入产出比的。针对你的问题回答一下
1.外圈分类，PS吧也只是几个字了。
2.空白圈和缺刻是通过人为增加不要的值来实现的，像上面代码那样。
3.字体大小角度当然可以改，在theme里可以慢慢调，但是由于你给的数据只有那么短，也不知道具体怎么才合适你的，所以就提供个思路。（另外你说的重叠问题，注意ggplot作图导出的时候图片大小直接影响字体的重叠，所以你要自己再试试）
4.颜色填充在scale_fill_gradient这个函数里自己改，取色我是随意的，教你个好方法，打开QQ，用Ctrl+Alt+A开启截图模式，移动到你想要的颜色地方，小框就会显示那里的RGB值，再按住ctrl，可以看“#F77400”这种值，可以直接应用到代码里

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-22 16:57:12

zerofung 发表于 2017-5-22 16:18

谢谢！但是还有一点小问题：
（1）那一圈代表疾病的分类，“A,B,C...”，怎么添加？
（2）代表年份的那个缺刻，角度好像有些歪呀？
（3）我在执行代码到前面几行，就出现了错误提示（见下面截图），不知为何？
（4）代码中有一些“%>%”是什么意思？

我把数据整理了一下，您再帮忙看看，好吗？

附件列表

www.jpg

原图尺寸 70.55 KB

ccc.rar

大小:3.45 KB

马上下载

本附件包括：

ccc.csv

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zerofung

2017-5-22 17:55:08

tmdxyz 发表于 2017-5-22 16:57
谢谢！但是还有一点小问题：
（1）那一圈代表疾病的分类，“A,B,C...”，怎么添加？
（2）代表年份的那 ...

1，ABC用PS加比研究出代码怎么加更快；2，年份角度你可以仔细看看你给的lancet图也是歪的，他把空白放在第一，即十二点方向的右边，我放在最后，即十二点的左边，仅此而已；3，是我说的这个论坛美元符号显示问题，你把a和Val，a和Year中间加个美元符号就好。另外你的数据我改过列名，类型我用的英语Type，名称我用的Name，尽量避免在R里用中文。4，%>%是dplyr里面的管道操作符，有空可以研究一下dplyr，对R的学习很有帮助

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

tmdxyz

2017-5-23 08:26:27

cheetahfly 发表于 2017-5-23 00:20
水平所限，逐渐试错拼凑出的来的。
首先，数据放在变量a中，数据变形的过程就不啰嗦了，最终变成了tidy da ...

非常感谢您。现在您这个图形几乎完美。但是我还有两个小小的问题：
（1）那个空白的放年份字符的缺刻，稍微有些宽，请问在哪个地方调节？请原谅我有一点强迫症
（2）您这个code太复杂了，我实在是蒙圈了。不过，这不是您的问题，是因为我太小白了
（3）您做这个图所依据的数据文件是哪一个？
（4）待我重复一下您的代码，再把您的设为精彩回复哈。请稍等

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 08:57:33

在美元符号之前的那一个语句
a <- df %>% gather(key = "Year", value = "Val", - c(Type, Name))
就出现了这个错误提示
Error in eval(expr, envir, enclos) : object 'Type' not found

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 09:01:47

唉，回复您的问题，还需要其他一些麻烦的手续？
您现在画的这个图，几乎完美。但是：
（1）那个空白的缺刻，其宽度怎么调节，语句在什么地方呀？
（2）总的来说，您这个code实在是让我蒙圈。我原以为ggplot+coord_polar就可以了呢
（3）您作这个图，用的是哪一个数据呢？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cheetahfly

2017-5-23 09:48:39

*****xyz 发表于 2017-5-23 09:01
唉，回复您的问题，还需要其他一些麻烦的手续？
您现在画的这个图，几乎完美。但是：
（1）那个空白的缺刻 ...

我不知道你是问的我还是楼上那位，我从自己的角度回答一下：
我的数据是你最开始的数据；
调节圆圈的缺口宽度和相对角位置，主要靠两个函数：
scale_x_discrete(expand = c(0, 1.5)) # 留出缺口
coord_polar(start = pi/360*15.5) # 保持一边竖直
这个图为什么代码这么多，因为它不是统计研究中的标准统计图表，而是将结果展示的商业个性图表，是需要花时间来细细调整，以达到自己想要的效果。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 10:24:11

最后一个问题。由于您的代码比较复杂，我一下子也学不会。我现在仿照您第一次回答我的那段代码，没有用tidyverse的包，而仅仅用ggplot2，采用以下代码
ggplot(dat)+geom_tile(aes(x=type,y=year,fill=death))
+ylim(c(2001,2014.5))
+annotate("text",x="e1",y=2004:2013,label=as.character(2004:2013))
+labs(x=NULL,y=NULL)+coord_polar()+theme_classic()
+theme(axis.line=element_blank(),axis.text.y=element_blank(),axis.ticks=element_blank())
+coord_polar(theta="x")
+scale_fill_gradient(high="#F77400",low="yellow",na.value="green")
得到的图形是歪的（我指的是年份那一圈缺刻）。请问，仅就我这一段代码而言，修改哪个地方，能使其扶正？

附件列表

new.jpeg

原图尺寸 291.38 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cheetahfly

2017-5-23 10:35:00

tmdxyz 发表于 2017-5-23 10:24
最后一个问题。由于您的代码比较复杂，我一下子也学不会。我现在仿照您第一次回答我的那段代码，没有用tidy ...

好像你现在还没有意识到，前后有两个人回答了你的问题

anyway......
我现在的网络看不到图，不知道你具体所指，如果是圆圈的角度要调整的话，主要调：
coord_polar(theta="x") ——> coord_polar(theta = "x", start = pi / 360 * 15.5)
最后一个数字不断微调，可以把图摆正。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 14:43:02

cheetahfly 发表于 2017-5-23 09:48
我不知道你是问的我还是楼上那位，我从自己的角度回答一下：
我的数据是你最开始的数据；
调节圆圈的缺 ...

在美元符号之前的那一个语句
a <- df %>% gather(key = "Year", value = "Val", - c(Type, Name))
就出现了这个错误提示
Error in eval(expr, envir, enclos) : object 'Type' not found

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 14:43:34

cheetahfly 发表于 2017-5-23 10:35
好像你现在还没有意识到，前后有两个人回答了你的问题
anyway......
我现在的网络看不到图，不知 ...

在美元符号之前的那一个语句
a <- df %>% gather(key = "Year", value = "Val", - c(Type, Name))
就出现了这个错误提示
Error in eval(expr, envir, enclos) : object 'Type' not found

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cheetahfly

2017-5-23 15:00:25

tmdxyz 发表于 2017-5-23 14:43
在美元符号之前的那一个语句
a % gather(key = "Year", value = "Val", - c(Type, Name))
就出现了这个 ...

这是提示你没有“Type”这个变量名，你检查一下，是不是大小写错误？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 15:21:21

cheetahfly 发表于 2017-5-23 15:00
这是提示你没有“Type”这个变量名，你检查一下，是不是大小写错误？

这个是来自于你自己的语句呀。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 15:26:30

cheetahfly 发表于 2017-5-22 10:20
水平所限，逐渐试错拼凑出的来的。
首先，数据放在变量a中，数据变形的过程就不啰嗦了，最终变成了tidy da ...

“数据放在变量a中”，这句话对应的语句在哪里呢？为什么我在运行你的代码是，会提示找不到“a”呢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 15:26:38

cheetahfly 发表于 2017-5-22 10:20
水平所限，逐渐试错拼凑出的来的。
首先，数据放在变量a中，数据变形的过程就不啰嗦了，最终变成了tidy da ...

“数据放在变量a中”，这句话对应的语句在哪里呢？为什么我在运行你的代码是，会提示找不到“a”呢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 15:39:51

zerofung 发表于 2017-5-22 16:18

虽然没有给你最佳回复，但是你也了我许多帮助。请你随便设置一个附件及其论坛币价格，我来购买吧。以示感谢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cheetahfly

2017-5-23 15:57:39

tmdxyz 发表于 2017-5-23 15:26
“数据放在变量a中”，这句话对应的语句在哪里呢？为什么我在运行你的代码是，会提示找不到“a”呢

我先将Excel中的数据变成下面的形式，注意我的变量名都是小写，我一般不主动用大小写混合。

然后，我用Crl+C将数据区域的单元格全部复制。
然后打开RStudio，输入代码：

复制代码

a就是这样来的了。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 16:56:15

cheetahfly 发表于 2017-5-23 15:57
我先将Excel中的数据变成下面的形式，注意我的变量名都是小写，我一般不主动用大小写混合。

然后，我用 ...

谢谢。这会儿按照您告诉我的，我也大体做出来了，尽管我还是不懂代码里所含的意思。

现在还有一个问题，就是我做出来的图，legend离图形过于近了（见附图），该怎么让它稍微远离一些？

附件列表

ppp.jpg

原图尺寸 284.15 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cheetahfly

2017-5-23 17:08:18

tmdxyz 发表于 2017-5-23 16:56
谢谢。这会儿按照您告诉我的，我也大体做出来了，尽管我还是不懂代码里所含的意思。

现在还有一个问题 ...

代码中的这一句：
legend.justification = c(-1, 1.5),
是微调legend的位置的，你可以自己微调一下。
另外，将图片另存为PDF格式后再截图，比直接存图或者直接截图的质量要高很多。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 17:14:05

cheetahfly 发表于 2017-5-23 17:08
代码中的这一句：
legend.justification = c(-1, 1.5),
是微调legend的位置的，你可以自己微调一下。

再次感谢。又有新问题：
我把代表疾病的name换成真正的名字以后（见截图10.jpg），然后全部复制，再利用语句
a <- tbl_df(read.table("clipboard", header = TRUE))
发现又有新的出错提示（见截图11.jpg）
“Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, :
line 4 did not have 13 elements”

请问我该怎么办？

附件列表

11.jpg

原图尺寸 20.1 KB

10.jpg

原图尺寸 82.44 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-23 17:20:46

zerofung 发表于 2017-5-22 16:18

您好！我利用您的代码，若果把这些疾病名称（a1,a2...d11）换成真正的名称（Leptospirosis,ECM,NT,Rabies,AD），那么画出来的图形，这些名称会自动按照字母顺序排列（AD,ECM,Leptospirosis,NT,Rabies）。可是我不想这样。该怎么办？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cheetahfly

2017-5-23 17:26:57

*****xyz 发表于 2017-5-23 17:14
再次感谢。又有新问题：
我把代表疾病的name换成真正的名字以后（见截图10.jpg），然后全部复制，再利用 ...

这是由于你的单元内容中含有空格造成的。
我们私聊吧，不要在这里“盖楼”了。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zhangyangsmith

2017-5-24 15:55:45

See the original post:
请问这张图是用R哪个包画出来的？
https://bbs.pinggu.org/forum.php? ... &from^^uid=563172

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-24 17:20:29

zhangyangsmith 发表于 2017-5-24 15:55
See the original post:
请问这张图是用R哪个包画出来的？
https://bbs.pinggu.org/forum.php?mod=viewthr ...

运行到此行时，提示出错。您能否看看？
Error in seq(from = lAngRngbgn, to = lAngRngbgn, to = lAngRngend, length.out = 101) :
  object 'lAngRngbgn' not found

# Fill in data

lapply(
  grep( "year", names(dfPlt) ),
  function(i)
  {
ivY <-
   as.integer(
gsub( "year", "", grep( "year", names(dfPlt), value = TRUE ) )
      )

iRad <- ivY + iYearMap

lapply(
   seq.int(iNclass - 1),
   function(j)
   {
lAngRng <- fIndToAng(j)

dvAng <- seq(from = lAngRngbgn,to=lAngRngbgn, to = lAngRngend, length.out = 101)

polygon(
   x = c( iRad*sin(dvAng), rev( (iRad + 1)*sin(dvAng) ) ),
   y = c( iRad*cos(dvAng), rev( (iRad + 1)*cos(dvAng) ) ),
   border = NA, col = hsv( h = iHue, s = dfPlt[j, i], v = 1 )
         )
   }
   )
  }
   )
最后的图形是这样

附件列表

zaq.jpg

原图尺寸 226.54 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zerofung

2017-5-24 22:02:22

复制代码

附件列表

TIM截图20170524220236.png

原图尺寸 64.98 KB

ccc.xlsx

大小:21.11 KB

马上下载

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zerofung

2017-5-24 22:04:46

尽力了，根据楼上各位的提示做的，我也学到很多。也许要改的只有颜色了吧，按我之前说的自己挑选配色（真的花好多时间在细节微调上）

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2017-5-24 22:52:17

zerofung 发表于 2017-5-24 22:04
尽力了，根据楼上各位的提示做的，我也学到很多。也许要改的只有颜色了吧，按我之前说的自己挑选配色（真的 ...

非常感谢您！我在RStudion上运行，会出错，提示的结果见
Error: unexpected input in:
" labs(x = NULL, y = NULL) + annotate("text", x = c(3, 11, 27, 40), y = 2014, label = LETTERS[1:4]) +
coord_polar(theta = "x"?

然后我把语句适当分成一块块，在R上运行，也是运行到这一步：coord_polar(theta = "x"，start = pi/360*9.5) 就卡住了，只能出现一个堆叠的column图。

附件列表

22.jpeg

原图尺寸 197.23 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

最佳答案

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群