8.10论文研究

第一步--类似教程学习

用到的包有:affyPLM,
affy,
RColorBrewer,
immpute,
limma,
pheatmap,
ggplot2

ReadAffy函数:

读取当前工作目录下所有cel格式文件,并生成表达矩阵

rma算法:

一种Affmetrix芯片数据预处理算法参考链接

exprs函数:提取表达矩阵信息
wirte.table函数:

注意各部分的意义参考链接

read.table函数:注意各部分的意义参考链接

match函数:

左边的值在右边的位置

示例.png

user_id为1,2,3,4,每一个数字在w中的位置依次是4,1,2,

is.na函数:

检测是否为Not Available,即是不是缺失值


示例.png

image.png
!(感叹号):

表示“与、或、非”中的“非”的意思

which函数:

找到满足条件的下标


image.png
as.matrix函数:

一般是将整个数据框(各列数据类型必须相同,否则进行强制转换)转换为维数相同的矩阵,目的是为了利用那些矩阵运算的函数


image.png
image.png
方括号"[]":

表示从元素中索引


a与上图同.png
factor(因子):

名义型变量和有序类别变量在R中称为因子,函数以一个整数向量的方式储存类别值。
https://www.cnblogs.com/larryzzZ/p/6650715.html

apply函数:

把一个function作用到array或者matrix的margins(可以理解为数组的每一行或者每一列)中,返回值时vector、array、list。
简单的说,apply函数经常用来计算矩阵中行或列的均值或其他函数值。
https://www.cnblogs.com/nanhao/p/6674063.html

image.png

rnorm函数:

(n, mean = 0, sd = 1)

tapply函数:

通常有三个参数,第一个指定输入数据,第二个是指定输入数据如何分组,第三个参数指定在每一个分组内用什么函数,整体功能概括为把数据按照某种依据分组,在每个组内进行某个运算。
https://www.zhihu.com/question/39843392

自定义函数function:

https://blog.csdn.net/xxzhangx/article/details/52924851

limma函数:

一个神奇的在线工具

aggregate函数:

一种整合函数,可用于多个探针对应同一个基因的情况时对数据进行分类整合

函数参数解析

一个综合使用的例子(sjz自上传数据集-空气质量年月日)
另一个综合使用的例子(mtcars数据集)

colMeans函数:

计算矩阵中列的平均值
综合示例

dplyr包之排列: arrange()
help文档解析

熟悉的mtcars数据集使用范例
split函数

split(x, f, drop = FALSE, sep = ".", lex.order = FALSE, ...)
-x: 一个待分组的向量或者data frame
-f: 函数,一个factor或者list(如果list中元素交互作用于分组中),以此为规则将x分组
-drop: 逻辑值,如果f中的某一个level没有用上则被弃用
-value: 一个储存向量的list,其形式类似于分组完成之后返回的那个list
函数split()可以按照分组因子,把向量,矩阵和数据框进行适当的分组。它的返回值是一个列表,代表分组变量每个水平的观测。这个列表可以使用sapply(),lappy()进行处理(apply – combine步骤),得到问题的最终结果。
只是分组,既可以对向量分组,也可以对数据框分组。

一个简单示例

impute包:

使用KNN算法处理缺失值
蝴蝶迁移性配对研究范例
详情可查阅参考手册

limma包:

1.帮助文档查询

>library(limma)
>limmaUsersGuide()
 # or alternatively
>help.start()
  #and follow the links to the limma package help.
#输出:如果什么都不发生的话,你应该自己打开‘http://127.0.0.1:25351/doc/html/index.html’

错误代码示例

#错误
write.table(probe_exprs,file="cancer.probeid.exprs.txt",sep="\t")
#正确
write.table(probe_exprs,file="cancer.probeid.exprs.txt",sep="\t",quote=F,row.names = F)
错误.png

正确.png

输出时自动调整第一列为列号,而使得原始的第一列的数据变成第二列

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,776评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,527评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,361评论 0 350
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,430评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,511评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,544评论 1 293
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,561评论 3 414
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,315评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,763评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,070评论 2 330
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,235评论 1 343
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,911评论 5 338
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,554评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,173评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,424评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,106评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,103评论 2 352