专栏名称: 生信宝典
生物信息分析入门、晋级和经验分享。Linux、R、Python学习教程;高通量测序数据分析学习教程;生信软件安装教程。所有内容均为原创分享,致力于从基础学习到提高整个过程。
目录
相关文章推荐
BioArt  ·  Cancer Cell丨靶向IL-15R信号 ... ·  3 小时前  
生物学霸  ·  Nature ... ·  19 小时前  
BioArt  ·  Neuron ... ·  2 天前  
BioArt  ·  EMBO ... ·  3 天前  
51好读  ›  专栏  ›  生信宝典

R语言学习 - 线图一步法

生信宝典  · 公众号  · 生物  · 2017-07-10 06:30

正文

请到「今天看啥」查看全文


  • 外层引号与内层引号不能相同

  • 凡参数值中包括了 空格 括号 逗号 等都用引号括起来作为一个整体。

完整参数列表如下:

ct@ehbio:~ $sp_lines.sh

***CREATED BY Chen Tong (chentong_biology@163.com)***

Usage:

/MPATHB/self/s-plot/sp_lines.sh options

Function:

This script is used to draw a line or multiple lines using ggplot2.
You can specify whether or not smooth your line or lines.

Two types of input files are supported, normal matrix or melted matrix format. Column separator for both types of input files is **tab**.

Here is an example of normal matrix format. The first column will be treated as X-axis variables and other columns represents each type of lines. The number of columns is unlimited and names of columns is unlimited.

**Set** column is not needed. If given, (multiple plots in one page) could be displayed.

------------------------------------------------------------
Pos    H3K27ac    CTCF    Enhancer    H3K4me3    polII
-5000    8.71298    10.69130    11.7359    10.02510    8.26866
-4000    8.43246    10.76680    11.8442    9.76927    7.78358
-3000    8.25497    10.54410    12.2470    9.40346    6.96859
-2000    7.16265    10.86350    12.6889    8.35070    4.84365
-1000    3.55341    8.45751    12.8372    4.84680    1.26110
0    3.55030    8.50316    13.4152    5.17401    1.50022
1000    7.07502    10.91430    12.3588    8.13909    4.88096
2000    8.24328    10.70220    12.3888    9.47255    7.67968
3000    8.43869    10.41010    11.9760    9.80665    7.94148
4000    8.48877    10.57570    11.6562    9.71986    8.17849
------------------------------------------------------

------------With SET------------------------------------------
Pos    H3K27ac    CTCF    Enhancer    H3K4me3    polII    Set
-5000    8.71298    10.69130    11.7359    10.02510    8.26866    1
-4000    8.43246    10.76680    11.8442    9.76927    7.78358    1
-3000    8.25497    10.54410    12.2470    9.40346    6.96859    1
-2000    7.16265    10.86350    12.6889    8.35070    4.84365    1
-1000    3.55341    8.45751    12.8372    4.84680    1.26110    1
0    3.55030    8.50316    13.4152    5.17401    1.50022    1
1000    7.07502    10.91430    12.3588    8.13909






请到「今天看啥」查看全文