我试图以黑白方式为出版物完成一些图形,并且在轴上遇到了一些问题。
我需要绘图轴(不是网格)中的数据的子标签,我通过使用 annotation_logticks(base = 2,sides =trbl)基2,因为我的数据是不相关的日志。
但似乎不适用于日期格式类型(xts)的y轴。
到目前为止,我得到的图是我所需要的,就是删除背景网格(该部分应该很容易)并在y上添加刻度轴(这对我来说是困难的部分)。
是我用来做图的代码:
ratio_plot_start_time =01:45 ratio_plot_end_time =02: (x =比率* 100,y =发病))+ geom_point(aes(y =发病))p < - ggplot(onsetratios,aes (粘贴(date,ratio_plot_start_time,sep =)),as.POSIXct(paste (表达(分钟=起始误差,ymax =开始+误差),宽度= 0.0) 17)))+主题(aspect.ratio = 2 /(1 + sqrt(5)))p <-p + theme_bw(base_size = 14,base_family =Helvetica)+ annotate(text ,x = max(onsetratios $ ratio),y = as.POSIXct(paste(date,ratio_plot_end_time,sep =)),vjust = 4,hjust = 2,label = paste(F,channelNo,' ))p< - p + annotation_logticks(base = 2,sides =trbl)用于构建图的数据为:
> onsetratios 比例发病误差 1 0.01092106 2012-05-17 01:49:03 19.05407 2 0.02092106 2012-05-17 01:49:35 14.89534 3 0.03092106 2012- 05-17 01:51:01 12.82914 4 0.04092106 2012-05-17 01:50:22 21.45099 5 0.05092106 2012-05-17 01:50:22 21.45099 6 0.06092106 2012 -05-17 01:51:13 20.09877 7 0.07092106 2012-05-17 01:51:13 20.09877 8 0.08092106 2012-05-17 01:51:13 20.09877 9 0.09092106 2012-05-17 01:51:13 20.09877 10 0.10092106 2012-05-17 01:52:34 15.53187 11 0.11092106 2012-05-17 01:53:53 14.79303 12 0.12092106 2012-05-17 01:53:53 14.79303 13 0.13092106 2012-05-17 01:53:32 17.12521 14 0.14092106 2012-05-17 01:53:32 17.12521 15 0.15092106 2012-05-17 01:53:32 17.12521 16 0.16092106 2012-05-17 01:53:32 17.12521 17 0.17092106 2012-05-17 01:46:30 36.96600 18 0.18092106 2012-05-17 01:46:30 36.96600 19 0.19092106 2012-05-17 01:46:30 36.96600解决方案
我不建议使用 annotation_logticks 。这是做这件事的预期方法。我删除了标题,因为我没有安装字体。
编辑:我借用函数此处,它可以正常工作,通过手工拆分结构提供:
onsetratios < - read.table(clipboard,head = T ) onsetratios $ onset < - as.POSIXct(onsetratios $ onset) ratio_plot_start_time =01:45 ratio_plot_end_time =02:10 channelNo = 1 start_date< - as.POSIXct(paste(date,ratio_plot_start_time,sep =)) end_date< - as.POSIXct(paste(date,ratio_plot_end_time,sep =))日期< - 2012-05-17p < - ggplot(onsetratios,aes(x = ratio * 100,y = onset))+ geom_point(aes(y = onset))+ geom_errorbar(aes(ymin =发作误差,ymax =发作+误差),宽度= 0.0)+ xlab(表达式(峰值计数的百分比))+ ylab(表达式(UT Time(2012 May 17 )))+主题(方面.ratio = 2 /(1 + sqrt(5)))+ theme_bw(base_size = 14)+ annotate(text,x = max(onsetratios $ ratio),y = as.POSIXct (paste(date,ratio_plot_end_time,sep =)),vjust = 4,hjust = 2,label = paste(F,channelNo,'))+ ylim(c(start_date,end_date) ) insert_minor< - function(major_labs,n_minor){labs < - c(sapply(major_labs,function(x)c(x,rep(,4))) ) labs [1:(length(labs)-n_minor)]} date_br1 < - seq(from = start_date,to = end_date,by =1 min) date_br5 < - seq(from = start_date,to = end_date,by =5 min)p + scale_x_continuous(breaks = 0:20,labels = insert_minor(seq(0,20,by = 5) ,4))+ scale_y_datetime(breaks = date_br1,labels = insert_minor(format(date_br5,%H:%M),4))+ theme(panel.grid.minor = element_blank ),panel.grid.major = element_blank())
I'm trying to get some graphics done for a publication in black and white, and I'm having some trouble with the axis.
I need to plot the subticks for the data in the axis (not the grid), and I found a trick around it by using the annotation_logticks(base=2, sides="trbl") with base 2, since my data is not relevant in log.
But that doesn't seem to work for the y axis, which is in date format type (xts).
The graph I got so far is the one bellow, all I need is to remove the background grid (that part should be easy) and to add the ticks on the y axis (that's the difficult part for me).
This is the code I used to do the graph:
ratio_plot_start_time = "01:45" ratio_plot_end_time = "02:10" channelNo = 1 p <- ggplot(onsetratios, aes(x=ratio*100, y=onset)) + geom_point(aes(y=onset)) p <- p + ylim(as.POSIXct(paste(date, ratio_plot_start_time, sep=" ")), as.POSIXct(paste(date, ratio_plot_end_time, sep=" "))) p <- p + geom_errorbar(aes(ymin=onset-error, ymax=onset+error), width=0.0) p <- p + xlab(expression("percent of peak counts"))+ ylab(expression("UT Time (2012 May 17)")) + theme(aspect.ratio = 2/(1+sqrt(5))) p <- p + theme_bw(base_size = 14, base_family = "Helvetica") + annotate("text",x=max(onsetratios$ratio),y=as.POSIXct(paste(date, ratio_plot_end_time, sep=" ")),vjust=4,hjust=2,label=paste("F", channelNo, "'")) p <- p + annotation_logticks(base=2, sides="trbl")And the data used for building the graph is:
> onsetratios ratio onset error 1 0.01092106 2012-05-17 01:49:03 19.05407 2 0.02092106 2012-05-17 01:49:35 14.89534 3 0.03092106 2012-05-17 01:51:01 12.82914 4 0.04092106 2012-05-17 01:50:22 21.45099 5 0.05092106 2012-05-17 01:50:22 21.45099 6 0.06092106 2012-05-17 01:51:13 20.09877 7 0.07092106 2012-05-17 01:51:13 20.09877 8 0.08092106 2012-05-17 01:51:13 20.09877 9 0.09092106 2012-05-17 01:51:13 20.09877 10 0.10092106 2012-05-17 01:52:34 15.53187 11 0.11092106 2012-05-17 01:53:53 14.79303 12 0.12092106 2012-05-17 01:53:53 14.79303 13 0.13092106 2012-05-17 01:53:32 17.12521 14 0.14092106 2012-05-17 01:53:32 17.12521 15 0.15092106 2012-05-17 01:53:32 17.12521 16 0.16092106 2012-05-17 01:53:32 17.12521 17 0.17092106 2012-05-17 01:46:30 36.96600 18 0.18092106 2012-05-17 01:46:30 36.96600 19 0.19092106 2012-05-17 01:46:30 36.96600解决方案
I do not advise to use annotation_logticks. Here's the intended way to do that. I removed the title since I do not have the font installed.
Edit:I borrowed the function from here, and it works fine, provided by manual breaks construction:
onsetratios <- read.table("clipboard", head=T) onsetratios$onset <- as.POSIXct(onsetratios$onset) ratio_plot_start_time = "01:45" ratio_plot_end_time = "02:10" channelNo = 1 start_date <- as.POSIXct(paste(date, ratio_plot_start_time, sep=" ")) end_date <- as.POSIXct(paste(date, ratio_plot_end_time, sep=" ")) date <- "2012-05-17" p <- ggplot(onsetratios, aes(x=ratio*100, y=onset)) + geom_point(aes(y=onset)) + geom_errorbar(aes(ymin=onset-error, ymax=onset+error), width=0.0) + xlab(expression("percent of peak counts"))+ ylab(expression("UT Time (2012 May 17)")) + theme(aspect.ratio = 2/(1+sqrt(5))) + theme_bw(base_size = 14) + annotate("text",x=max(onsetratios$ratio),y=as.POSIXct(paste(date, ratio_plot_end_time, sep=" ")),vjust=4,hjust=2,label=paste("F", channelNo, "'")) + ylim(c(start_date, end_date)) insert_minor <- function(major_labs, n_minor) {labs <- c( sapply( major_labs, function(x) c(x, rep("", 4) ) ) ) labs[1:(length(labs)-n_minor)]} date_br1 <- seq(from = start_date, to = end_date, by = "1 min") date_br5 <- seq(from = start_date, to = end_date, by = "5 min") p + scale_x_continuous(breaks = 0:20, labels = insert_minor(seq(0, 20, by=5), 4)) + scale_y_datetime(breaks = date_br1, labels = insert_minor(format(date_br5, "%H:%M"), 4)) + theme(panel.grid.minor=element_blank(), panel.grid.major=element_blank())更多推荐
在日期格式轴中绘制ggplot2中的次要刻度线(不是网格刻度线)
发布评论