12、MA图的计算过程

为了简化问题,假设有3张芯片,每组数有9个探针:

Data:

2,4,6,7,9,10,4,7,8,3

9,5,3,2,5,7,9,10,3,12

6,4,3,2,7,8,1,2,6,9

一、给3组数取2的对数

Log2Data:

1.0,2.0,2.5849626,2.807355,3.169925,3.321928,2.0,2.807355,3.0

3.169925,2.321928,1.5849625,1.0,2.321928,2.807355,3.169925,3.321928,1.5849625

2.5849626,2.0,1.5849625,1.0,2.807355,3.0,0.0,1.0,2.5849626

二、上面的结果有9列,给每一列取中位数:

Median:

2.5849626,2.0,1.5849625,1.0,2.807355,3.0,2.0,2.807355,2.5849626

四、M的值为每一列减去该列中位数:

M = Log2Data-Median:

-1.5849625,0.0,1.0,1.8073549,0.36257008,0.32192808,0.0,0.0,0.4150375

0.5849625,0.32192808,0.0,0.0,-0.4854268,-0.19264507,1.169925,0.51457316,-1.0

0.0,0.0,0.0,0.0,0.0,0.0,-2.0,-1.8073549,0.0

五、A的值是每一列的值加上该列中位数后*0.5

A = 0.5*(data_log+medianchip)

1.7924813,2.0,2.0849626,1.9036775,2.98864,3.160964,2.0,2.807355,2.7924812

2.8774438,2.160964,1.5849625,1.0,2.5646415,2.9036775,2.5849626,3.0646415,2.0849626

2.5849626,2.0,1.5849625,1.0,2.807355,3.0,1.0,1.9036775,2.5849626

六、以上的M、A值每一行代表一个样品,以M为y轴坐标,A为x轴坐标,绘制出来的散点图就是MA图

七.每个样品的Sigma=每个样品的M的中位数

0.32192808

0

0

八、每个样品的IQR=每个样品的M的IQR

0.4150375

0.7072182

0

时间: 2024-10-09 12:45:39

12、MA图的计算过程的相关文章

09、RNA降解图的计算过程

RNA降解是影响芯片质量的一个很重要的因素,因为RNA是从5’开始降解的,所以理论5’的荧光强度要低于3’.RNA降解曲线可以表现这种趋势. 以样品GSM286756.CEL和GSM286757.CEL为例: library(affy) rawData<-ReadAffy("GSM286756.CEL","GSM286757.CEL") deg<-AffyRNAdeg(rawData) plotAffyRNAdeg(deg) 结果: RNA降解图是只用P

润乾报表教程-集算报表优化计算过程

报表作为数据统计分析软件,当它出现性能问题需要对数据源计算进行优化时,执行路径难以控制是阻碍报表优化的难题之一.这是由于数据库执行路径不透明,程序员很难甚至无法干预执行路径,也就难以提高数据库的性能.而一般报表工具不具备强计算能力,大部分计算仍然要依靠数据库进行,这就导致很多报表优化效果不理想. 不同于一般报表工具,润乾集算报表内置了专门用于数据计算的集算引擎,开发人员可以通过编写集算脚本完成报表数据源准备.与数据库执行SQL路径不可控相比,集算脚本的执行过程是可控的,开发人员可根据实际情况编写

rhgb Red Hat 图形化启动过程

rhgb:Red Hat 图形化启动过程 HEL系列的跟Ubuntu不一样,需要加参数vga=0x317. rhgb Red Hat 图形化启动过程

10、探针组检测值的计算过程

一组探针是否能被检测到,可以用检测值来表示.simpleaffy包的detection.p.val方法默认alpha1=0.05,alpha2=0.065.检测值小于alpha1为A(无,Absent),介于alpha1和alpha2之间为M(不确定,Marginal Present),大于alpha2为P(有,Present). 以样品GSM286756.CEL为例,计算它的231211_s_at探针组的检测值: library(affy) library(simpleaffy) rawDat

实战MYSQL 8.0.12 主主复制配置过程

实战MYSQL 8.0.12 主主复制配置过程 搭建环境: Server name IP mysql1 192.168.200.1 mysql2 192.168.200.2 服务器版本:CentOS Linux release 7.5.1804 (Core)MYSQL版本:8.0.12 # 采用源码安装方式, 此过程略,或者参考 http://blog.51cto.com/snowlai/2140451 由于MYSQL采用的是源码安装方式,没有生成 /etc/my.cnf 文件,需要手动创建,创

随机森林之oob的计算过程

随机森林有一个重要的优点就是,没有必要对它进行交叉验证或者用一个独立的测试集来获得误差的一个无偏估计.它可以在内部进行评估,也就是说在生成的过程中就可以对误差建立一个无偏估计. 随机森铃在生成每颗决策树时,会随机且有放回的抽取样本,每棵决策树会有大概1/3的样本未抽取到,这些样本就是每棵树的oob样本.具体计算过程如下: 根据这种特点,我们可以对其进行oob估计,步骤如下: 1.计算决策树对其对应的oob样本的分类情况(约有1/3棵树参与oob估计) 2.以投票的方式确定该样本的分类 3.计算o

ES6中reduce的计算过程

reduce计算过程如下 原文地址:https://www.cnblogs.com/malong1992/p/12128672.html

gantt甘特图的制作过程

甘特图主要是用来做项目管理的,可以清楚的看到任务间的逻辑关系,任务与时间关系和任务间并行关系. 在甘特图中,横轴方向表示时间,纵轴方向并列着活动列表.图表内可以用线条.数字.文字代号等来表示计划(实际)所需时间,计划(实际)输出,计划(实际)开工或完工时间等. 1.根据计划列出所有任务 对于每一个任务,不管这项任务是并行任务或者是线性任务,列出任务的启动时间,估算完成该任务所需要花费的时间.如果本项任务是线性任务,同时也列出这项任务所依赖的任务,最后使用表格的方式将所有的任务按先后顺序列出.下面

12.5.1 声明计算类型

在这个示例中,计算的类型(在 Haskell 术语中的单子类型)称为 ValueWrapper<T>,它只是存储泛型类型参数T 的值,并不会为类型增加任何其他功能.这就是说,计算只处理标准值,但是,我们能够用 C# 中的查询表达式, F# 中的计算表达式,写代码(来增加功能). 清单 12.17 显示了在 C# 和 F# 中的类型声明.在 C# 中,我们创建简单的类,在 F# 中,将使用一个简单的只有一种情况的差别联合. 清单 12.17 C# 和 F# 中的计算值 // C# 的类声明 cl