查看: 300  |  回复: 6

小马科研部落

新虫 (小有名气)

[交流] R语言之初体验,导个数据,画个图已有4人参与

每天都学一点点,每天都进步一点点。
----------------------------------------------------
之前以为R跟SPSS一样是一个统计软件,
由于不小心跨行到基因测序行业,发现在生物信息分析的某端,序列拼接、注释完后,竟然需要R软件来分析、作图。
由于在技术支持岗位上,万一以后有客户咨询到图是怎么做的,万一答不出来,就“惨”喽!
而公司又没有培训,只能自己自学R语言。
毕竟有点SPSS基础,所以对于R语言的统计分析还算是比较容易上手。
-------------------------------------------------------
与SPSS的鼠标操作相比,R语言需要输入代码,但是运行速度很快。
简单的一个代码/命令,就可以省略SPSS的“点点点点点”。
唯一的缺点就是要记代码。
----------------------------------------------------------
R的下载、安装就不说了。跟普通软件的安装一样。

对于SPSS来说,拿到数据后,把数据导入,之后进行运算处理。R语言也一样。
但要注意2点:
1 路径。
在【文件】--【改变工作目录】中,修改到 数据所在的目录。
假设你的数据在D:\data里面,那么工作目录就改为 D:\data。
改变好了路径后,可以用【getwd()】这个命令来看一下当前工作路径。

2 数据格式。
xls的格式R 是不认的,需要将其转为csv格式。用read.csv来带入。
------------------------------------------------------

做完上面的准备,就可以用R做普通分析和画图了。
以【数量生态学—R语言】中的数据为例,画个图。

由于数据中有3列,第一列为样本名称,一般第一列应该是样本名,所以应该用row.names=1

由于有的数据行数比较多,你想看一下有多少行,可以用这个命令:dim()
-------------------------------------------------------
画个图
其中main为图片的主题,xlab为横坐标,ylab为纵坐标。
plot(spa,main='地图',xlab='x轴',ylab='y轴')

图上的点的标志也都可以换。如pch=2
plot(spa,main='地图',xlab='x轴',ylab='y轴',pch=2,col='red')
------------------------------------------------------
具体换法如下图所示

当然,也可以在图中加上线
lines(spa,col='blue') 记住,为lines,有S。

-------------------------------------------------------
可以把每个样本的编号标上。
text(spa,row.names(spa),cex=0.8,col='black')

由于该图,是由画大框、连线、标字母,3部分组成。
为了怕出错,可以将代码写入notepad++,将3行代码,直接复制到R中运行。
--------------------------------------------------------
具体的图片,可以去微信公众号【小马科研部落】里查看。
每天学习一点点,进步一点点。

R语言之初体验,导个数据,画个图
加点的注释.png
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

更多内容,请关注微信公众号【小马科研部落】,有问题可找QQ308091253
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

deephill

铁杆木虫 (职业作家)

送红花一朵
顶个
2楼2017-09-04 23:06:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一个人不够

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主有没有包的介绍呀?我最近在学R包的muma,学不会。。。
3楼2017-09-12 20:26:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小马科研部落

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 一个人不够 at 2017-09-12 20:26:15
楼主有没有包的介绍呀?我最近在学R包的muma,学不会。。。

这个我还没有用到哈,我主要学习跟二代高通量测序相关的分析。
更多内容,请关注微信公众号【小马科研部落】,有问题可找QQ308091253
4楼2017-09-14 16:23:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪竹小木虫

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我想问一下亲爱的楼主,R语言能分析差异表达基因吗?跪求!急!

发自小木虫Android客户端
5楼2017-10-17 21:04:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

8023zx

新虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
Excel格式文件R是可以识别的吧?

发自小木虫Android客户端
6楼2017-12-01 19:58:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小马科研部落

新虫 (小有名气)

引用回帖:
6楼: Originally posted by 8023zx at 2017-12-01 19:58:33
Excel格式文件R是可以识别的吧?

可以识别的,excel格式存为csv,在R语言用read.csv进行调用

发自小木虫IOS客户端
更多内容,请关注微信公众号【小马科研部落】,有问题可找QQ308091253
7楼2017-12-01 21:51:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小马科研部落 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 读研感觉 +18 zaj123 2018-09-18 18/900 2018-09-19 14:58 by 胖胖的懒虫虫
[考博] 导师竟然没有实验室 +37 永远的208 2018-09-16 43/2150 2018-09-19 13:03 by 无视军团
[有机交流] 聚四氟乙烯盖子怎么会融化了? 5+3 北方sdy 2018-09-19 5/250 2018-09-19 11:14 by 天涯不归路
[硕博家园] 博士想退学,不知道自己读博的目标是什么? +48 凌古风 2018-09-13 65/3250 2018-09-19 10:47 by theseachem
[找工作] 关于博士毕业后的去处 +3 zhangnaeve 2018-09-18 3/150 2018-09-19 07:59 by zfuh
[论文投稿] 有偿数据处理 +4 mmmmm1010 2018-09-18 15/750 2018-09-18 22:48 by linjx2000
[硕博家园] 写给毕业找工作 想转行的你 +8 瑾瑜兰珞 2018-09-13 19/950 2018-09-18 16:34 by 星辰00544
[海外博后] 分享一个美国博后信息 +9 15165266591 2018-09-15 12/600 2018-09-18 12:25 by F5-BOY
[教师之家] 餐厅的饭菜越来越质次价高 +3 舒克8265 2018-09-15 6/300 2018-09-18 11:47 by 舒克8265
[有机交流] 钯碳氢气还原硝基把吡嗪环也还原了 20+4 asone77 2018-09-16 10/500 2018-09-18 08:43 by qq123neng
[基金申请] NSFC限5篇代表作预示:数论文数项目数帽子时代的终结! +23 路口狮子 2018-09-13 33/1650 2018-09-17 18:09 by jzhang_cs
[论文投稿] 计算机核心期刊哪些比较好中呀? +7 苏醒266 2018-09-15 7/350 2018-09-17 16:37 by char1es
[考博] 博导回复 +11 懒懒的夏日 2018-09-12 19/950 2018-09-17 15:36 by 懒懒的夏日
[硕博家园] 有一点动心 +5 小花888 2018-09-16 8/400 2018-09-16 23:09 by 小花888
[硕博家园] 计算机专业 +3 桃之妖妖呀 2018-09-16 3/150 2018-09-16 12:02 by gorgan
[硕博家园] 辅导员打电话叫我直博 +12 doraiba 2018-09-16 25/1250 2018-09-16 10:56 by 郝华伟皇甫
[硕博家园] 有木有不住学校,在外边租房子学习的研二或者研三同学啊 +9 小花888 2018-09-15 12/600 2018-09-15 13:32 by emu_yu
[考博] 十万火急~心凉凉之后冉冉升起~ +11 小蜗牛高速 2018-09-15 12/600 2018-09-15 12:09 by 纪文朋
[版块工场] 【早起签到贴】2018年9月14日(q q 48626260 进群申请:小木虫——您个人ID) +35 8475 2018-09-14 35/1750 2018-09-14 13:56 by coke
[考博] 要是没什么文章,有机会读博吗? +18 科研的杜公子 2018-09-13 22/1100 2018-09-14 12:53 by 十早韦
信息提示
请填处理意见