我知道这似乎微不足道,但由于某种原因,我无法弄清楚如何做到这一点.
我已经制作了这个情节
使用这些命令:
data1 <- read.table("data_for_question.txt", sep="\t", header=TRUE)
plot(data1, type="l", lwd=2)
Run Code Online (Sandbox Code Playgroud)
从数据我将在下面添加.
这是一个很好的情节,但刻度线和数字相距太远.我想在每个第4个数据点都有一个刻度标记,它将以20,000个时间点的增量,而不是如此处所示的50,00.
我为长长的数据列表道歉,但绘制的数据点数似乎与我在刻度标记中的问题有关.
谢谢!
我知道这似乎是重复的,但我已经检查了几个小时,其他类似的问题,并且由于某些原因,没有一个解决方案适合我.
数据:
Time Counts
0 29
5000 68
10000 55
15000 34
20000 61
25000 47
30000 42
35000 39
40000 25
45000 3
50000 30
55000 2
60000 36
65000 29
70000 13
75000 13
80000 17
85000 11
90000 10
95000 29
100000 20
105000 8
110000 0
115000 1
120000 1
125000 1
130000 1
135000 0
140000 6
145000 0
150000 …Run Code Online (Sandbox Code Playgroud) 我想做一个相对简单的情节(让人想起如下的时间线:http://www.ats.ucla.edu/stat/sas/code/timeline.gif),但不是时间在x轴上,它将成为基因组中的基础位置."时间跨度"将是DNA序列支架的覆盖距离,显示它们落入基因组的范围,它们重叠的位置和没有覆盖的位置.这是我正在寻找的粗略模型,显示rRNA的重叠群覆盖,(我遗漏了,但需要,x轴显示开始和停止的位置,以及重叠群的标记(彩色线)): http://i.imgur.com/MDABx.png,具有以下坐标:
Contig# Start1 Stop1 Start2 Stop2 Start3 Stop3 Start4 Stop4
1 1 90 90 100 120 150 200 400
2 1 100 120 150 200 400 NA NA
3 1 30 90 100 120 135 200 400
4 1 100 120 140 200 400 NA NA
5 -35 80 90 100 130 150 200 400
6 1 100 200 300 360 400 NA NA
Run Code Online (Sandbox Code Playgroud)
我很确定这可以在R中完成,可能使用ggplot2,但由于某种原因我无法弄明白.