搜索
查看: 1910|回复: 1

[mRNA-seq] jimmy大神,5行代码搞定GEO学习总结版,代码的一点疑问

[复制链接]

2

主题

5

帖子

526

积分

高级会员

Rank: 4

积分
526
发表于 2017-9-26 16:09:15 | 显示全部楼层 |阅读模式
本帖最后由 quan1912 于 2017-9-27 20:10 编辑

这是jimmy大神的原帖,5行代码搞定GEO学习总结版,全程无尿点
http://www.biotrainee.com/forum. ... 12&fromuid=2929
(出处: 生信技能树)


在学习大神的代码,对于代码有一些问题请教


pdata=pData(eSet[[1]])
#这一行代码的大概意思是在title里面用空格将文字分开,然后选择了第一组
group_list=factor(unlist(lapply( pdata$title ,function(x) strsplit(as.character(x)," ")[[1]][1])))
#这是在第一组里面选择了normal作为对照组
group_list=relevel(group_list,'Normal')

我的问题是如果我有3组数据,只想计算其中的两组如何实现呢?自己想了好久,希望大神能帮助解惑,不胜感激!



上一篇:多重序列比对及系统发生树的构建
下一篇:jimmy大神,5行代码搞定GEO学习总结版,代码的一点疑问
回复

使用道具 举报

2

主题

41

帖子

385

积分

中级会员

Rank: 3Rank: 3

积分
385
发表于 2017-10-9 00:18:41 | 显示全部楼层
删除你不想比较的那一组数据就可以实现了
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树 ( 粤ICP备15016384号  

GMT+8, 2019-7-21 04:58 , Processed in 0.036605 second(s), 25 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.