24小时热门版块排行榜     先丰纳米 买纳米材料就点这里

热门推荐

查看: 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的回帖
相关版块跳转 我要订阅楼主 小马科研部落 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 他一点儿都看不懂的文章都敢要一作 +32 偏振片 2018-01-17 61/3050 2018-01-17 19:06 by caoximin
[硕博家园] 博导要求每个月汇报至少四次 +32 笑得开心 2018-01-16 33/1650 2018-01-17 16:26 by shifukuan
[找工作] 汕头大学如何 +15 owen1394 2018-01-11 15/750 2018-01-17 10:36 by oyzc
[硕博家园] 深度学习 +3 追忆时光 2018-01-15 3/150 2018-01-16 07:30 by uerth
[有机交流] 格式试剂溶剂用乙醚和四氢呋喃有区别吗? 5+5 暮野星空 2018-01-13 17/850 2018-01-16 01:23 by hocky
[硕博家园] 你们东北真的是太冷了 +57 眼睛666 2018-01-11 67/3350 2018-01-15 20:11 by hollyword
[考研] 考研复试英语应该怎么准备 +4 卢森林 2018-01-12 4/200 2018-01-15 13:50 by lyn311
[微生物] 发酵罐入口空气过滤用什么 +4 爱菜樱花? 2018-01-11 5/250 2018-01-15 11:16 by 爱菜樱花?
[版块工场] 【早起签'到'贴】2018年1月15日(q q 48626260 进群申请:小木虫——您个人ID) +80 8475 2018-01-15 80/4000 2018-01-15 09:06 by flypanzer
[教师之家] 怎样解决迟到 (金币+3) +6 没有如果… 2018-01-13 9/450 2018-01-14 13:23 by 可爱可敬的人
[考研] 325求调剂 +11 丫丫686 2018-01-13 14/700 2018-01-14 12:40 by 2226471723
[教师之家] 非以诚相待就足矣 +7 行走的虫虫 2018-01-11 12/600 2018-01-13 19:57 by @小小芳
[金属] 炉中焊,同时出现晶粒长大和裂缝(如图)的因素有哪些? 10+3 拓海1627 2018-01-12 14/700 2018-01-13 16:02 by 拓海1627
[硕博家园] 论老师要求数据分析到2017的痛苦性 +6 学禾 2018-01-12 12/600 2018-01-13 09:53 by shindowfeng
[考研] 调剂 +5 shero_ 2018-01-11 6/300 2018-01-12 23:17 by Nitara
[硕博家园] 有用吗 (金币+2) +7 雪儿1993 2018-01-10 13/650 2018-01-12 21:48 by 764765360
[硕博家园] 学校对文章排名认可不专业丫 +4 Me诺言liu 2018-01-12 6/300 2018-01-12 14:25 by guo854639263
[硕博家园] 我认为自己做出了原创性的东西 (金币+2) +18 物理所院士 2018-01-11 20/1000 2018-01-12 12:50 by mincarrot
[基金申请] 山东省基金2018年面上项目大约多少经费呢? +7 nenu-jlu 2018-01-11 7/350 2018-01-12 12:21 by gaozheng81
[版块工场] 【早起签'到'贴】2018年1月11日(q q 48626260 进群申请:小木虫——您个人ID) +60 8475 2018-01-11 60/3000 2018-01-11 08:43 by phoo2000
信息提示
请填处理意见