查看: 5748|回复: 9|关注: 0
从EXCEL读入的日期在画图时怎么样标注在横坐标上？
我从EXCEL文件读入一列日期，出来的是一个单元数组，例如 { ''& & ''& & '' ..................},我想在plot作图时候把它放在横坐标上可是怎么试都不行，不知道怎么解决呢？希望大家帮帮我非常感谢。
set(gca,'XTick',-pi:pi/2:pi)
set(gca,'XTickLabel',{'-pi','-pi/2','0','pi/2','pi'})
注意你的XTick的元素个数要与XTickLabel的个数要一致。
你这个例子我试过了不行的，因为你第一句set(gca,'XTick',-pi:pi/2:pi)& &里面的,-pi:pi/2:pi 怎么弄成我从EXCEL读过来日期呢？&&这是问题是关键&&后面的好说
因为只是画图，你的数据是确定的，而且你的
日期间隔是固定的(在图形中都是距离为1），所以你的横坐标数据不是很重要。即使你没有
横坐标数据都没关系，直接设定横坐标为
x = 1:length(y)
楼上是不是脑残啊？ 我不是为了要日期横轴我干嘛来这里发帖子！ 不会就别回答，捣什么乱啊？！
请自己先试下在评论好不？
y = rand(1,3)
plot(1:3,y)
str ={ ''& & ''& & '' };
set(gca,'XTick',1:3)
set(gca,'XTickLabel',str)
关注者： 1
回复 5# wjjyinxin 的帖子
自己看不懂就骂人。。。牛。
这个方法我早就试过了 我要的不是这个效果明白吗？ 你们一群白痴！！！
回复 8# wjjyinxin 的帖子
做图后加上datetick这个命令试试……
学习了，替楼主感谢各位回帖的朋友。力争创建一个和谐的论坛环境。
Powered by查看: 2998|回复: 9
阅读权限32威望0 级论坛币3125 个学术水平0 点热心指数0 点信用等级0 点经验18320 点帖子816精华0在线时间657 小时注册时间最后登录
积分 2303, 距离下一级还需 1297 积分
权限: 自定义头衔, 签名中使用图片, 设置帖子权限, 隐身, 设置回复可见
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡, 沉默卡下一级可获得
道具: 千斤顶
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 2 天连续签到: 1 天[LV.1]初来乍到
下面是我用R的boxplot画的一张图，从横坐标轴上可以看到，刻度上的值是隔次出现，如果我想完整显示刻度值，该怎么设置？
19:42:43 上传
载入中......
(10.38 KB)
19:41:51 上传
阅读权限22威望0 级论坛币51 个学术水平0 点热心指数0 点信用等级0 点经验673 点帖子143精华0在线时间34 小时注册时间最后登录
积分 183, 距离下一级还需 77 积分
权限: 自定义头衔
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡下一级可获得
权限: 签名中使用图片
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
阅读权限30威望0 级论坛币10 个学术水平3 点热心指数18 点信用等级0 点经验24747 点帖子963精华0在线时间412 小时注册时间最后登录
积分 1453, 距离下一级还需 772 积分
权限: 自定义头衔, 签名中使用图片, 设置帖子权限, 隐身
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡下一级可获得
权限: 设置回复可见道具: 沉默卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 271 天连续签到: 1 天[LV.8]以坛为家I
最好的医生是自己，最好的药物是时间……
阅读权限32威望0 级论坛币3125 个学术水平0 点热心指数0 点信用等级0 点经验18320 点帖子816精华0在线时间657 小时注册时间最后登录
积分 2303, 距离下一级还需 1297 积分
权限: 自定义头衔, 签名中使用图片, 设置帖子权限, 隐身, 设置回复可见
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡, 沉默卡下一级可获得
道具: 千斤顶
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 2 天连续签到: 1 天[LV.1]初来乍到
trier2006 发表于
试试rug这是我生成这个图的代码，
png(&china.png&,width=1260,height=660)
boxplot(cpi[,2:33])
请问命令该加在那里？
阅读权限30威望0 级论坛币10 个学术水平3 点热心指数18 点信用等级0 点经验24747 点帖子963精华0在线时间412 小时注册时间最后登录
积分 1453, 距离下一级还需 772 积分
权限: 自定义头衔, 签名中使用图片, 设置帖子权限, 隐身
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡下一级可获得
权限: 设置回复可见道具: 沉默卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 271 天连续签到: 1 天[LV.8]以坛为家I
peter 发表于
这是我生成这个图的代码，
png(&china.png&,width=1260,height=660)
boxplot(cpi[,2:33])boxplot后面
最好的医生是自己，最好的药物是时间……
阅读权限26威望0 级论坛币3085 个学术水平4 点热心指数5 点信用等级3 点经验3443 点帖子225精华0在线时间352 小时注册时间最后登录
积分 598, 距离下一级还需 202 积分
权限: 自定义头衔, 签名中使用图片
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 隐身
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
哼哼签到天数: 9 天连续签到: 1 天[LV.3]偶尔看看II
这个主要是因为横坐标字符太多，不大可能显示所有名称，可以尝试减少box plots的数量，分两个图来画，或者把字体大小调小。
阅读权限32威望0 级论坛币3125 个学术水平0 点热心指数0 点信用等级0 点经验18320 点帖子816精华0在线时间657 小时注册时间最后登录
积分 2303, 距离下一级还需 1297 积分
权限: 自定义头衔, 签名中使用图片, 设置帖子权限, 隐身, 设置回复可见
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡, 沉默卡下一级可获得
道具: 千斤顶
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 2 天连续签到: 1 天[LV.1]初来乍到
trier2006 发表于
boxplot后面请问是这么加：
png(&china.png&,width=1260,height=660)
boxplot(cpi[,2:33],rug)
还是这么加：
png(&china.png&,width=1260,height=660)
boxplot(cpi[,2:33]);rug
第一种情况下出来的图形还是原来的样子，没什么变化。
第二种情况下图形也没变，出来一些错误提示：
& png(&china.png&,width=1260,height=660)
& boxplot(cpi[,2:33]);rug
function (x, ticksize = 0.03, side = 1, lwd = 0.5, col = par(&fg&),
& & quiet = getOption(&warn&) & 0, ...)
& & x &- as.vector(x)
& & ok &- is.finite(x)
& & x &- x[ok]
& & if (!quiet) {
& && &&&u &- par(&usr&)
& && &&&u &- if (side%%2 == 1) {
& && && && &if (par(&xlog&))
& && && && && & 10^u[1L:2]
& && && && &else u[1L:2]
& && &&&else {
& && && && &if (par(&ylog&))
& && && && && & 10^u[3:4]
& && && && &else u[3:4]
& && &&&if (any(x & u[1L] | x & u[2L]))
& && && && &warning(&some values will be clipped&)
& & Axis(side = side, at = x, labels = FALSE, lwd = 0, lwd.ticks = lwd,
& && &&&col.ticks = col, tck = ticksize, ...)
&bytecode: 0x024fc15c&
&environment: namespace:graphics&
& dev.off()
null device
& && && & 1
阅读权限32威望0 级论坛币3125 个学术水平0 点热心指数0 点信用等级0 点经验18320 点帖子816精华0在线时间657 小时注册时间最后登录
积分 2303, 距离下一级还需 1297 积分
权限: 自定义头衔, 签名中使用图片, 设置帖子权限, 隐身, 设置回复可见
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡, 沉默卡下一级可获得
道具: 千斤顶
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 2 天连续签到: 1 天[LV.1]初来乍到
毓闻 发表于
这个主要是因为横坐标字符太多，不大可能显示所有名称，可以尝试减少box plots的数量，分两个图来画，或者把 ...请问如何把刻度上的字调小？
阅读权限26威望0 级论坛币3085 个学术水平4 点热心指数5 点信用等级3 点经验3443 点帖子225精华0在线时间352 小时注册时间最后登录
积分 598, 距离下一级还需 202 积分
权限: 自定义头衔, 签名中使用图片
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 隐身
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
哼哼签到天数: 9 天连续签到: 1 天[LV.3]偶尔看看II
peter 发表于
请问如何把刻度上的字调小？可以在boxplot函数里面添加参数cex.axis=0.8之类的，参数值小于1就会使字体缩小，但是能不能解决你的问题我还不清楚。
阅读权限32威望0 级论坛币3125 个学术水平0 点热心指数0 点信用等级0 点经验18320 点帖子816精华0在线时间657 小时注册时间最后登录
积分 2303, 距离下一级还需 1297 积分
权限: 自定义头衔, 签名中使用图片, 设置帖子权限, 隐身, 设置回复可见
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡, 沉默卡下一级可获得
道具: 千斤顶
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 2 天连续签到: 1 天[LV.1]初来乍到
毓闻 发表于
可以在boxplot函数里面添加参数cex.axis=0.8之类的，参数值小于1就会使字体缩小，但是能不能解决你的问题 ...谢谢指点！
我把cex.axis设为0.5就可以完整显示了。
只是可惜字迹太小太模糊，不过知道了一个方法，蛮开心的。谢谢！
论坛好贴推荐下次自动登录
现在的位置:
R语言绘图符号
语言画图的时候，所能利用的pch参数值如下图，可有下面的代码生成
plot(rep(1,10),ylim=c(-2,1.2),pch=1:10,cex=3,axes=F,xlab=&&,ylab=&&)
text(rep(0.6,10),as.character(1:10))
points(rep(0,10),pch=11:20,cex=3)
text(rep(-0.4,10),as.character(11:20))
points(rep(-0.8,5),pch=21:25,cex=3)
text(rep(-1.2,5),as.character(21:25))
points(6:10, rep(-0.8,5),pch=c(&*&,&?&,&X&,&x&,&&&),cex=3)
text(6:10,rep(-1.2,5),c(&*&,&?&,&X&,&x&,&&&))
【上篇】【下篇】
您可能还会对这些文章感兴趣！
2000人的生物信息学QQ群
生物信息学①群 :(已满)
生物信息学②群 :(开放中)
请先在注册，凭PLoB的用户名验证入群。
生物信息快速问答社区
百度站内搜索
生物信息培训班推荐
我的收藏夹 &&----> 基础绘图工具1
基础绘图工具1
&&&&R语言的绘图功能及应用案例&&&&&&&&深圳华大基因研究院闻博wenbo@日&&&&&&&&主要内容&&&&一、基础绘图工具1.绘图函数2.绘图参数介绍3.画图面板分割4.图形保存二、图形案例三、其它图形包&&&&&&&&一、基础绘图工具&&&&1.2.3.4.绘图函数绘图参数介绍画图面板分割图形保存&&&&&&&&一、基础绘图工具&&&&1.绘图函数&&&&在R中有两种绘图函数：&&&&1.1高级绘图函数1.2低级绘图函数创建一个新的图形在现有的图形上添加元素&&&&&&&&demo(graphics)demo(persp)&&&&&&&&一、基础绘图工具&&&&1.绘图函数&&&&1.1高级绘图函数创建一个新的图形&&&&&&&&函数名plot(x)plot(x,y)pie(x)boxplot(x)hist(x)barplot(x)&&&&pairs(x)coplot(x~y|z)matplot(x,y)qqnorm(x)image(x,y,z)heatmap(x)&&&&&&&&功能以x的元素值为纵坐标、以序号为横坐标绘图x与y的二元作图饼图盒形图（也称箱线图）x的频率直方图x的值的条形图如果x是矩阵或是数据框，作x的各列之间的二元图&&&&关于z的每个数值（或数值区间）绘制x与y的二元图二元图，其中x的第一列对应y的第一列，依次类推正态分位数-分位数图x,y,z三元图热图&&&&&&&&一、基础绘图工具&&&&1.绘图函数&&&&1.2低级绘图函数&&&&函数名points(x,y)lines(x,y)text(x,y,labels,…)mtext(text,side=3,line=0,…)segments(x0,y0,x1,y1)arrows(x0,y0,x1,y1,...)abline(a,b)abline(h=y)abline(v=x)abline(lm.obj)rect(x1,y1,x2,y2)polygon(x,y)legend(x,y,legend)title()axis(side,vect)box()rug(x)locator(n,type=n,…)&&&&&&&&在现有的图形上添加元素&&&&功能添加点添加线在(x,y)处添加用labels指定的文字在边空添加用text指定的文字从(x0,y0)各点到（x1,y1）各点画线段同上，但添加箭头绘制斜率为b和截距为a的直线在纵坐标y处画水平线在横坐标x处画垂直线画出lm.obj确定的回归线绘制长方形，(x1,y1)为左下角，(x2,y2)为右上角绘制连接各x,y坐标确定的点的多边形在点(x,y)处添加图例，说明内容由legend给定添加标题，也可添加一个副标题画坐标轴在当前的图上加边框在x轴上用短线画出x数据的位置获取鼠标在图中点击处的坐标&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?高级绘图函数共同参数选项：&&&&&&&&选项add=FALSE&&&&axes=TRUE&&&&&&&&功能如果是TRUE，叠加图形到前一个图上(如果有的话)&&&&如果是FALSE，不绘制轴与边框指定图形的类型，“p”:点，“l”:线，“b”:点连线，“o”:同上，但是线在点上，“h”:垂直线，“s”:阶梯式，垂直线顶端显示数据，“S”：同上，但是垂直线底端显示数据&&&&&&&&type=p&&&&&&&&xlim=,ylim=xlab=,ylab=main=sub=&&&&&&&&指定轴的显示范围坐标轴的标签主标题副标题&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?高级绘图函数共同参数选项：&&&&&&&&maintype=o&&&&&&&&ylim=c(20,100)&&&&&&&&ylab&&&&&&&&axesxlim=c(7,20)sub&&&&&&&&xlab&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?高级绘图函数共同参数选项：&&&&&&&&plot(cars,font.lab=2,main=主标题,sub=副标题,type=p,xlim=c(7,20),ylim=c(20,100))plot(cars,font.lab=2,main=主标题,sub=副标题,type=p,xlim=c(7,20),ylim=c(20,100),axes=F)&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?高级绘图函数共同参数选项：&&&&&&&&plot(cars,font.lab=2,type=p)plot(cars,font.lab=2,type=l)plot(cars,font.lab=2,type=b)plot(cars,font.lab=2,type=o)plot(cars,font.lab=2,type=h)plot(cars,font.lab=2,type=s)plot(cars,font.lab=2,type=S)&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?其它常用绘图参数help(par)&&&&参数adjbgbtycexcolfontlasltylwdmarmfcol,mfrowpchpsptytckxaxt,yaxt功能控制关于文字对齐方式背景色图形边框形状符号和文字大小,cex.axis,cex.lab,cex.main,cex.sub颜色,col.axis,col.lab,col.main,col.sub文字字体，font.axis,font.lab,font.main,font.sub坐标轴刻度数字标记方向连线的线性线的宽度图形边空的大小c(nr,nc)的向量，前者按列分隔绘图窗口,后者按行分隔绘图窗口绘图符号的类型文字大小绘图区域类型轴刻度长度，tck=1时，则绘制grid如果xaxt=“n”,设置x轴不显示,如果yaxt=“n”,设置y轴不显示&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?其它常用绘图参数help(par)&&&&par(mgp=c(1.6,0.6,0),mar=c(3,3,2,1))X=1:10;Y=1:10cex_type=c(0.7,0.8,0.9,1,2,3,4,5,6,7)plot(X,Y,col=blue,pch=16,cex=cex_type,#绘图符号大小，1是正常大小cex.axis=1.1,#坐标轴刻度数字大小cex.lab=1.2,#坐标轴标签大小cex.main=1.4,#标题文字大小#cex.sub=1.3,#副标题文字大小font.axis=3,#坐标轴刻度数字字体，斜体font.lab=2,#坐标轴标签字体，粗体font.main=4,#标题字体，粗斜体#font.sub=1,#副标题字体，正常main=字体,xlim=c(1,12),ylim=c(1,12))text(X,Y,adj=0,labels=paste(cex=,cex_type),col=red,cex=1.2,font=3)&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?其它常用绘图参数help(par)&&&&?R的绘图符号：pch=1:25&&&&&&&&?R的绘图符号：任意字符par(mgp=c(1.6,0.6,0),mar=c(3,3,2,1))pch_type=c(1:5,15:19)X=1:10;Y=rep(6,10)plot(X,Y,col=1,pch=pch_type,cex=2,main=pch,font.lab=2)text(X,Y,adj=-0.5,labels=paste(pch=,pch_type),srt=90)&&&&注：背景色bg参数只对符号21-25有作用&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?其它常用绘图参数help(par)&&&&?颜色名称：white,black,red,yellow,…&&&&&&&&colors()&&&&?数字代号：1,2,3,4,5,6,7,8?颜色函数：&&&&&&&&rainbow,heat.colors,terrain.colors,topo.colors,cm.colors,colors()&&&&par(mgp=c(1.6,0.6,0),mar=c(3,3,2,1))pch_type=c(1:5,15:19)X=1:10;Y=rep(6,10)plot(X,Y,col=1:10,pch=pch_type,cex=2,main=pch,font.lab=2)text(X,Y,col=1:10,adj=-0.5,labels=paste(pch=,pch_type),srt=90)&&&&&&&&一、基础绘图工具&&&&2.绘图参数介绍&&&&?其它常用绘图参数help(par)&&&&?控制连线的类型：lty1:实线，2:虚线，3:点线，4:点虚线，5:长虚线，6:双虚线）?控制线的宽度：lwd1为正常大小，1为减小宽度，1为增加宽度&&&&&&&&par(mgp=c(1.6,0.6,0),mar=c(3,3,2,1))X=1:10;Y=1:10line_type=c(1:6)line_width=seq(from=0.6,by=0.6,length=6)plot(X,Y,col=blue,pch=16,type=n,main=linetypeandlinewidth,font.lab=2,cex.lab=1.2)abline(h=3:8,lty=line_type,col=colors()[120:126],lwd=line_width)text(3:8,3:8,adj=0,labels=paste(lty=,line_type,,lwd=,line_width),col=red,cex=1.2,font=3)&&&&&&&&一、基础绘图工具&&&&3.画图面板分割&&&&?在一个面板中画多张图?par中参数mfrow和mfcol#par(mfrow=c(2,2))source(mfrow.R)?layout函数生成复杂的图形排列layout(mat,widths=rep(1,ncol(mat)),heights=rep(1,nrow(mat)),respect=FALSE)source(layout.R)?其它函数?split.screen&&&&&&&&一、基础绘图工具&&&&3.画图面板分割&&&&?在一个面板中画多张图图形边距控制：par(mar=c(5,4,4,2)+0.1)par(mai=c(bottom,left,top,right))par(oma=c(bottom,left,top,right))par(mgp=c(3,1,0))&&&&&&&&一、基础绘图工具&&&&4.图形保存&&&&?输出到屏幕?windows,X11windows()X11()?输出到文件?pdf,postscript,xfig,bitmap,pictex,cairo_pdf,svg,png,jpeg,bmp,tiffpdf(plot.pdf,width=4,height=4)png(plot.png,width=400,height=600)dev.off()#绘制完图形后关闭图形设备dev.list()通过菜单命令保存图形&&&&&&&&?&&&&&&&&二、图形案例&&&&&&&&二、图形案例&&&&1.散点图――点线混合图&&&&require(stats)head(cars)speeddist225816plot(cars[,1],cars[,2],pch=16,col=red,xlab=speed,ylab=dist,type=p,font.axis=2,font.lab=2,cex.lab=1.5)x=seq(5,25,by=0.1);y=10*sin(x)+50lines(x,y,col=“green,lwd=2)a=5:25;b=2*a+10*cos(a)points(a,b,pch=15,col=“blue)legend(“topleft”,pch=c(16,-1,15),lty=c(-1,1,-1),col=c(“red”,”green,blue),legend=c(plotwithpoint,lines,points))&&&&&&&&二、图形案例&&&&1.散点图――并列散点图&&&&a=read.table(“acgt.txt”,head=T)matplot(a,type=o,pch=15:18,lty=1,bg=yellow,lwd=2,xlab=Position(nt),ylab=Frequency,font=2,font.lab=2,cex.lab=1.5)legend(topright,pch=15:18,lty=1,merge=T,legend=names(a),col=1:4)&&&&&&&&01020
&&&&&&&&A0.......&&&&&&&&C0.......&&&&&&&&G0.......&&&&&&&&T/U0.......&&&&&&&&二、图形案例&&&&1.散点图――坐标对数化&&&&TagTCACAGTGAACCGGTCTCTTTAGAGGTAGTAGGTTGCATAGTTTCACAGTGAACCGGTCTCTTTTTGAGGTAGTAGGTTGTATGGTTTGAGGTAGTAGGTTGTATAGTTTGAGGTAGTAGGTTGTATGTAGCACCATCTGAAATCGGTTATGAGGTAGTAGATTGTATAGTTTAAGGCACGCGGTGAATGTGAGGTAGTAGGTTGTATGGTTTR1RS90&&&&&&&&a=read.table(“replicate_experiment_sRNA.txt,head=T)plot(a[,2],a[,3],log=xy,cex.lab=1.4,cex.axis=1.2,font.lab=2,font.axis=2,xlab=R1,ylab=RS,pch=.,col=red,cex=1.1,axes=F)axis(1,at=axTicks(1),labels=as.integer(axTicks(1)))axis(2,at=axTicks(2),labels=as.integer(axTicks(2)))title(main=“R1VSRS”)box()&&&&&&&&二、图形案例&&&&2.条形图――标准条形图&&&&&&&&&&&&a-read.table(lengthVSnumber.txt)#pdf(lengthVSnumber.pdf)bar=barplot(a[,2],ylim=c(0,1.2*max(a[,2])),cex.lab=1.5,names.arg=a[,1],xlab=Length(nt),font.lab=1.5,ylab=NumberofsmallRNAs,space=0.9,col=blue,main=Lengthdistribution)ratio=paste(sprintf(%.2f,a[,2]/sum(a[,2])*100),%,sep=)text(bar,a[,2],labels=ratio,adj=c(0.5,0),cex=0.8,font=2,pos=3)abline(h=axTicks(2),lty=2,col=rgb(0,0,0,0.5));box()&&&&&&&&二、图形案例&&&&2.条形图――堆积条形图&&&&80.......10.......30.......2915T/U0.......&&&&&&&&h=read.table(“acgt.txt”,head=T)barplot(t(as.matrix(h)),ylim=c(0,1),xlab=Position,ylab=Frequency,main=Homosapiens,col=rainbow(4),font=2,cex.lab=1.4)box()legend(topright,legend=c(A,G,C,T),pch=15,col=rainbow(4))&&&&&&&&二、图形案例&&&&2.条形图――邻接条形图&&&&t0t1t2t3t40...4554&&&&&&&&read.table(fGC.txt)-fpread.table(“tGC.txt)-tptable(cut(fp[,2],breaks=seq(0,1,by=0.1)))-fpftable(cut(tp[,2],breaks=seq(0,1,by=0.1)))-tpfdata.frame(tpf,fpf)-tfp.gctfp.gc-tfp.gc[,c(2,4)]dat-data.frame(tfp.gc[,1]/sum(tfp.gc[,1]),tfp.gc[,2]/sum(tfp.gc[,2]))row.names(dat)-seq(0.1,1,by=0.1)a-barplot(t(as.matrix(dat)),xlab=GCcontent,ylab=Frequency,main=Pre-piRNAsGCcontent,legend=c(Positivesample,Negativesample),yaxt=n,col=c(green,cyan),ylim=c(0,0.5),beside=T,args.legend=list(bty=n,cex=0.8));sp-spline(a[1,],dat[,1],n=100)lines(sp$x,sp$y,lwd=2,col=green)sp-spline(a[2,],dat[,2],n=100)lines(sp$x,sp$y,lwd=2,col=cyan);box()axis(2,seq(0,0.5,by=0.1),labels=c(0%,10%,20%,30%,40%,50%))&&&&&&&&二、图形案例&&&&2.条形图――ErrorBar条形图&&&&?显示均值和标准差的条形图&&&&&&&&data(ToothGrowth)head(ToothGrowth)tg.m-tapply(ToothGrowth$len,ToothGrowth$dose,mean)tg.sd-tapply(ToothGrowth$len,ToothGrowth$dose,sd)tg-cbind(tg.m,tg.sd)par(mgp=c(1.6,0.6,0))bp-barplot(tg.m,xlab=Dose,ylab=Growth,cex.lab=1.5,ylim=c(0,35),space=2,col=1)arrows(bp,tg.m,bp,tg.m+tg.sd,lwd=1.5,angle=90,length=0.1)box()&&&&&&&&二、图形案例&&&&3.文氏图&&&&&&&&library(plotrix)plot(0:10,seq(0,10,length=11),type=n,axes=F,xlab=,ylab=)draw.circle(2,5,2,col=rgb(154/255,0/255,205/255,0.6))draw.circle(4,5,2,col=rgb(21/255,3/255,252/255,0.6))text(1,5,labels=10.12%,col=white,font=2)text(5,5,labels=40.38%,col=white,font=2)text(3,5,labels=49.5%,col=white,font=2)legend(6.2,5,pch=15,xjust=0,yjust=0.5,bty=n,cex=1.3,col=c(rgb(154/255,0/255,205/255),rgb(74/255,2/255,233/255),rgb(21/255,3/255,252/255)),legend=c(sample1uniq,samplesample2,sample2uniq))text(3.5,7.5,labels=Vennchartforuniq_sRNAs,font=2,cex=1.5)&&&&&&&&二、图形案例&&&&4.饼图――二维饼图&&&&&&&&dat=read.table(“pie.txt)ratio=sprintf(%.2f,100*dat[,2]/sum(dat[,2]))ratio=paste(ratio,%,sep=)label=paste(dat[,1],ratio,sep=\n)pie(dat[,2],col=1:6,main=PiechartformiRNAsuniq,border=purple,labels=label,font=2)&&&&&&&&miR1miR2miR3miR4miR5miR6&&&&&&&&0800700&&&&&&&&二、图形案例&&&&4.饼图――三维饼图&&&&dat=read.table(dat.txt)ratio=sprintf(%.2f,100*dat[,2]/sum(dat[,2]))ratio=paste(ratio,%,sep=)label=paste(dat[,1],ratio,sep=\n)pie3D(dat[,2],col=1:6,main=Pie3DchartformiRNAsuniq,border=purple,labels=label,font=2,labelcex=1,,explode=0.1,radius=0.95)&&&&&&&&miR1miR2miR3miR4miR5miR6&&&&&&&&0800700&&&&&&&&二、图形案例&&&&5.经验分布图&&&&&&&&?经验分布图&&&&&&&&x=c(40.5,20.8,16.9,20.2,17.3,32.4,34.5,43.6,19.5,26.4,41.4)plot(ecdf(x),verticals=T,col=2,main=Empiricaldistributionfunctionx,font=2,cex.lab=1.5)lines(14:46,pnorm(14:46,mean(x),sd(x)),lwd=2,col=blue)&&&&&&&&二、图形案例&&&&6.盒形图(箱线图)&&&&?箱线图&&&&frequncypair10C:G/P16C:G/P17A:U/N9A:U/N14A:U/P4A:U/P20C:G/P16C:G/P10A:U/P7C:G/N&&&&&&&&data=read.table(“pair.txt”,head=T)boxplot(frequncy~pair,data=data,ylab=Frequncy,font.lab=2,cex.lab=1.5)&&&&&&&&二、图形案例&&&&7.频率直方图&&&&?直方图x=rnorm()x=x[x0]hist(x,xlim=range(x),xaxs=i,yaxs=i,main=HistofMFE,xlab=MFE,nclass=30,freq=F,font=2,cex.lab=1.4)lines(density(x),col=blue,lwd=2)lines(-30:0,dnorm(30:0,mean(x),sd(x)),col=red,lwd=2)legend(-30,0.15,pch=c(15,-1,-1),lty=c(1,1,1),col=c(gray,blue,red),legend=c(Histogram,Densityline,Normaldensityline))&&&&&&&&8.QQ图&&&&&&&&二、图形案例&&&&&&&&QQ图：鉴别样本的分布是否近似于某种类型的分布&&&&x=c(75.0,64.0,47.4,66.9,62.2,62.2,58.7,63.5,66.6,64.0,57.0,69.0,56.9,50.0,72.0)qqnorm(x,pch=21,col=red,bg=yellow,font=2,cex.lab=1.4);qqline(x,lwd=2,col=blue)x=rnorm(200,2,2)qqnorm(x,pch=21,col=red,bg=yellow,font=2,cex.lab=1.4);qqline(x,lwd=2,col=blue)&&&&&&&&二、图形案例&&&&9.热图&&&&?方法1&&&&x-read.table(cluster.txt,head=T)library(pheatmap)pheatmap(x[,-1],border_color=NA)&&&&&&&&二、图形案例&&&&9.热图&&&&?方法2&&&&library(gplots)x-read.table(cluster.txt,head=T)heatmap.2(as.matrix(x[1:100,]),col=redgreen(75),cexCol=0.9,key=T,symkey==none,trace=none,labRow=NA)&&&&&&&&二、图形案例&&&&9.热图&&&&?方法3&&&&x-read.table(cluster.txt,head=T)heatmap(as.matrix(x[,-1]),cexCol=0.9,labRow=NA)&&&&&&&&二、图形案例&&&&10.多组图&&&&x-read.table(cluster.txt,head=T)pairs(x[,-1])&&&&&&&&三、其它图形包&&&&gplotsggplot2lattice&&&&&&&&
分享给好友::
All Rights Reserved 图宝贝
本站声明：本站转所转载之内容，无任何商业意图，如本网站转载稿件涉及版权、著作权等问题，请您来函与本站管理员取得联系,友情链接请加QQ,要求PR 或者 BR >=2（联系方式：QQ ）