24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huaidan0331

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 176.9
帖子: 81
在线: 26.6小时
虫号: 4003214
注册: 2015-08-03
性别: MM
专业: 环境工程

[交流] 论文画图求助信息！已有9人参与

请教大位大神，下面这个图是用什么软件画的？看着很好看，我用origin没有画出来呢

发自小木虫IOS客户端

回复此楼

» 猜你喜欢

售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有4人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有6人回复
课题组招2027级博士上海工程技术大学激光智能制造方向已经有6人回复
现代”学阀”该如何界定已经有14人回复
求合成方法已经有8人回复
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？已经有4人回复
有机合成以后会不会被AI改变？做科研的虫友怎么看已经有8人回复
我的奶奶已经有3人回复

1楼 2018-01-10 14:48:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lvchengrui

专家顾问 (正式写手)

专家经验: +12
应助: 4 (幼儿园)
金币: 929.8
散金: 255
红花: 4
帖子: 327
在线: 284.6小时
虫号: 2448875
注册: 2013-05-05
性别: GG
专业: 信息系统与管理
管辖: 程序语言

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
paperhunter: 金币+3, 鼓励交流 2018-01-11 09:49:52

ggplot2可以应付，替换下面的TS数据

> require(ggplot2)
> figuredata <- read.csv("c:/figuredata.csv",
+                      header = TRUE)
> figuredata
YearMonth    Type TS
1 2009Nov Influent  6
2 2009Nov    STW1 16
3 2009Nov    STW2 17
4 2009Nov    STW3 15
5 2010Mar Influent  5
6 2010Mar    STW1 13
7 2010Mar    STW2 15
8 2010Mar    STW3 14
9 2010Jul Influent  5
10 2010Jul    STW1 45
11 2010Jul    STW2 46
12 2010Jul    STW3 48
13 2010Nov Influent  5
14 2010Nov    STW1 18
15 2010Nov    STW2 19
16 2010Nov    STW3 17
17 2011Mar Influent  5
18 2011Mar    STW1 18
19 2011Mar    STW2 19
20 2011Mar    STW3 17
21 2011Apr Influent NA
22 2011Apr    STW1 19
23 2011Apr    STW2 19
24 2011Apr    STW3 18
25 2011May Influent NA
26 2011May    STW1 21
27 2011May    STW2 20
28 2011May    STW3 19
29  2011June Influent NA
30  2011June    STW1 21
31  2011June    STW2 21
32  2011June    STW3 21
> levels(figuredata$YearMonth) <-
+ figuredata$YearMonth[seq(1, 32, by = 4)]
>
>
> ggplot(figuredata, aes(x = YearMonth, y = TS)) +
+ geom_col(aes(color = Type, fill = Type),
+          width = 0.1,
+          position = position_dodge(width = 0.6)) +
+ geom_point(aes(color = Type, shape = Type),
+             position = position_dodge(width = 0.6)) +
+ scale_y_continuous(limits = c(0,100),
+                   breaks = seq(0,100,by = 20)) +
+ geom_vline(xintercept = 5.5) +
+ annotate("text", x = c(3, 7), y = c(90, 90),
+          label = c("Feedind period", "Resting period"

) +
+ geom_hline(yintercept = seq(0,100,by = 20), linetype = 2) +
+ theme(panel.background = element_rect(fill = "white",
+ colour = "black"

,
+       legend.background = element_blank(),
+       legend.position = c(0.1,0.8),
+       legend.title = element_blank(),
+       legend.box.background = element_rect(fill = "white",
+                                           colour = "black"

,
+ legend.key = element_rect(colour = "black"

,
+ legend.key.width = unit(0.5, units = "cm"

,
+ legend.key.height = unit(0.5, units = "cm"

,
+       axis.text.x = element_text(colour = c("darkred",
+                                              rep("darkcyan", 3),
+                                              rep("darkblue", 4)))) +
+ labs(x = NULL, y = "TS(%)"

+
+ scale_color_brewer(palette = "Spectral", type = "seq"

+
+ scale_fill_brewer(palette = "Spectral", type = "seq"