PCR仪 » 讨论区 » 经验共享 » 做完荧光定量后如何统计数据【转自 丁香园论坛】

采购询价

点击提交代表您同意 《用户服务协议》 《隐私政策》

 
需要登录并加入本群才可以回复和发新贴

标题:做完荧光定量后如何统计数据【转自 丁香园论坛】

细水长流[使用道具]
二级
Rank: 2


UID 70450
精华 0
积分 142
帖子 63
信誉分 100
可用分 1006
专家分 0
阅读权限 255
注册 2011-8-10
状态 离线
21
 
关注这个话题很久了,但是大家好像没有说出个什么明确的思路。
我最近也在处理PCR的数据,就相对定量方法说说我的一些理解和看法。
1.我们的目的:研究一种或者是多种干预对一种疾病(现象)的影响。那么相对定量就需要得出我们的目的基因在某一种(多种)样品中的表达量相对于对照样品的上调(下调)倍数。这是对研究的终极目的---统计学上的总体而言。
2.实际操作:我们研究常常是对一个样本进行的。这里就需要从我们研究的样本数据经过统计分析而实现对终极目的的推断。
3.具体的操作:上面好多的战友都说了几种方法,都是有道理的,但是还是说的很粗略。包括我在学位论文数据库中检索了国内多所所谓的“牛”校他们的数据处理方法都是五花八门。而且你很难得看的明白具体的思路!!
例如:20只大鼠随机等分为两组,干预组以“善宁”治疗,对照组不用任何处理。观察“善宁”对肝脏中目的基因白介素-1(IL-1)的mRNA表达的影响。我们的终极目的就是想得到“善宁”是上调(下调)该基因表达的疗效。而此处20只大鼠仅仅是一个样本。要得到这个终极目标,就要用这20只大鼠的样本统计量去统计推断总体参数,才能实现这个终极目标。我们选用B-actin作为内参基因。

前提:目的基因和内参基因的扩增效率和相关系数要满足不同分析方法的前提。不需多说,上面的战友说的很多了。
那么
首先,干预组和对照组中,IL-1基因在相同的样品(肝脏)中的表达需要用内参基因B-actin归一化(人为的假定为“1”)。主要是消除我们取材时候“量”的差异(不能保证相同的细胞数)。
其次,(“善宁”)干预措施的效果是相对(没有处理)而言的,就是说从基因表达角度看,应该理解为:经过干预后,目的基因的表达是上调(下调)了。反应在数据上就是相对于对照样品,干预样品的目的基因上调(下调)了多少倍(这就是“相对”定量的确切含义)。此处就是以对照样品目的基因的表达量再一次的归一化。而实际上,如果在此用对照样品目的基因IL-1的表达量进行归一化处理(人为的假定为“1”),我们就只能得到这10只“善宁”干预组大鼠相对与没有任何处理的对照组的10只大鼠的目的基因IL-1的表达的倍数关系(对样本而言),还不能就说“善宁”能上调(下调)IL-1基因表达。还需要进一步的统计分析才能这样说!因为这一结论是我们的终极目的(对总体而言)。那么此处是两组样本的比较就可以选择t检验达到统计推断的目的。(当然如果是多组的比较可以选择单因素方差分析)。
而恰恰在进行统计推断时,我们面临一个问题!是什么呢?就是我们能得到“善宁”干预组肝脏中IL-1基因相对与对照组的表达倍数的一组数值(10个数据),可是要完成t检验就还需要对照组的另外10个数据,此处我们已经把对照组的IL-1的表达量假定为“1”了,总不能以均数为“1”,标准差为“0”进行统计分析吧!怎么办呢?比较好的办法就是:在上面进行对照组的归一化处理时,我们是用对照组的目的基因IL-1表达的“均数”来假定为“1”的。因此,可以以这个“均数”再次归一化处理对照组中不同样品的IL-1的表达量。这样就可以得到10个数据(此时“均数”不是10个“1”,标注差也不是10个“0”)。以此就可以完成统计分析了。
其实我们大可不必完成后两种归一化处理。直接把“善宁”干预组和“未处理”对照组中已经经过内参归一化处理的目的基因IL-1的表达量(此处的表达量是各个组织中的目的基因“绝对”表达量)进行t检验,得出我们要的结论------“善宁”是否上调(下调)IL-1的表达。
这里,需要说明的是:因为相对定量是相对而言的,所以在进行实验操作时“一定”要把各个样品一一对应起来,这样得出的CT值才能够进行相关的数据处理。理解也很简单,就是说干预组中1号样品的目的基因就必须和1号样品的内参基因的CT值进行相关数据处理,不能把1号样品的目的基因的CT值和2号样品的内参基因CT值进行数据处理!假如这样做就没有起到“归一化”的作用!但是,在用对照组的的目的基因来归一化处理干预组的目的基因时就需要把所有的对照组的目的基因的表达考虑进去了。否则不能得出干预的效果。这也是把对照组目的基因表达的“均数”用来归一化处理干预组目的基因表达的原因。
说了这么多,不知道对大家有没有帮助。请大家指正!
顶部
不懂[使用道具]
三级
Rank: 3Rank: 3


UID 70203
精华 0
积分 208
帖子 195
信誉分 100
可用分 1661
专家分 0
阅读权限 255
注册 2011-8-6
状态 离线
22
 

我也是real-time PCR的初做者,看了各位的发言,感触颇多。现有几个疑问,请大家帮忙解答:
我用的是SYBR染料法,拿到引物后首先要看看引物的扩增效率吧?(对cDNA倍比稀释做标准曲线,R2的意义是加样准确与否吗?斜率可以反映扩增效率)如果目的基因和管家基因扩增效率不同,可以用Pfaffl方法做统计,但是这种方法还需要加做一步---看同一基因在不同处理组是否具有相同的扩增效率。那么这一步每个处理组都要做吗?还是只看一两个处理组就可以了呢?
请大家指教!万分感谢!



[ 本帖最后由 不懂 于 2011-8-12 16:09 编辑 ]
顶部
不懂[使用道具]
三级
Rank: 3Rank: 3


UID 70203
精华 0
积分 208
帖子 195
信誉分 100
可用分 1661
专家分 0
阅读权限 255
注册 2011-8-6
状态 离线
23
 
求助:
求助:
我是新手。我想请问运用双delta Ct方法进行数据处理后,
最终作图时的标准差值应该怎样计算阿?

具体说:
待测样品中的目的基因和内参,加上对照样品的目的基因和内参
共有4个值会根据PCR重复结果取平均值,那么最后的目的基因相
对表达量值标准差该怎样计算?

谢谢  
顶部
跳跳糖[使用道具]
三级
Rank: 3Rank: 3


UID 70449
精华 0
积分 157
帖子 73
信誉分 100
可用分 1066
专家分 0
阅读权限 255
注册 2011-8-10
状态 离线
24
 
说了这么多,我觉得还是没有人把如何统计最后结果说明白?
比如这个吧,文献中对照组都是用(1+-标准差)表示的,可是123三个对照组的均值根本就很难是1,除非凑巧是1,那好办了。

有人给解答一下么?
顶部
嘉年华[使用道具]
三级
Rank: 3Rank: 3


UID 70441
精华 0
积分 222
帖子 184
信誉分 100
可用分 1651
专家分 0
阅读权限 255
注册 2011-8-10
状态 离线
25
 


QUOTE:
原帖由 跳跳糖 于 2011-8-12 16:10 发表 bbcodeurl('http://bbs.antpedia.com/images/common/back.gif', '%s')
说了这么多,我觉得还是没有人把如何统计最后结果说明白?
比如这个吧,文献中对照组都是用(1+-标准差)表示的,可是123三个对照组的均值根本就很难是1,除非凑巧是1,那好办了。

有人给解答一下么? ...


用2-ddtCT统计。
顶部
嘉年华[使用道具]
三级
Rank: 3Rank: 3


UID 70441
精华 0
积分 222
帖子 184
信誉分 100
可用分 1651
专家分 0
阅读权限 255
注册 2011-8-10
状态 离线
26
 
我想请问一下,我要做不同组织间同蛋白表达的相对定量.可是我发现我选的内参GAPDH在不同组织间的表达量不同,那还可不可以做为内参啊?
顶部
饭团团[使用道具]
三级
Rank: 3Rank: 3


UID 70425
精华 1
积分 182
帖子 99
信誉分 102
可用分 1201
专家分 10
阅读权限 255
注册 2011-8-10
状态 离线
27
 
我也想知道无论用Delta Ct还是其他方法处理之后得到的比值如何进行统计分析。做方差分析还是做多重比较。如果是以2-Delta Ct为原始数据进行统计分析的话,单位是什么。
顶部
开心蔬菜帮[使用道具]
二级
Rank: 2


UID 70426
精华 0
积分 142
帖子 63
信誉分 100
可用分 1001
专家分 0
阅读权限 255
注册 2011-8-10
状态 离线
28
 
做定量,在普通上扩增得到很好的目的条带,无引物二聚体,但是到定量PCR仪上确得不到扩增曲线,换了另外一台定量PCR仪也是如此。各位大侠帮帮忙,找找原因。谢谢了
顶部
快乐的大脚[使用道具]
三级
Rank: 3Rank: 3


UID 70429
精华 0
积分 211
帖子 161
信誉分 100
可用分 1501
专家分 0
阅读权限 255
注册 2011-8-10
状态 离线
29
 
实验中的real-time PCR数据通过机器直接得到:CT平均值(mean CT)、CT值的标准差(CT-SD),我通过2 -△△Ct 公式计算出ddCT值(相对表达水平),但是不知道如何求ddCT值的SD和P值?

举个例子或许更容易说清楚:

House-keeping gene: GAPDH

Target gene: CCL6

Ctrl group:
GAPDH
3次CT值 18.21;17.89; 18.15
平均值:18.08
SD值:0.17222

CCL6
3次CT值 31.79; 32.08; 31.68
平均值 31.85
SD值:0.20903

Experiment group:
GAPDH
3次CT值 19.88; 19.63; 19.78
平均值 19.76
SD值0.12777

CCL6
3次CT值 33.35; 33.07; 33.63
平均值 33.35
SD值:0.27807

这样通过公式求出的CCL6的ddCT=1.132883885

但是如何求这个ddCT值的SD,及其P值大小?
顶部
章鱼小丸子[使用道具]
三级
Rank: 3Rank: 3


UID 70431
精华 0
积分 177
帖子 113
信誉分 100
可用分 1281
专家分 0
阅读权限 255
注册 2011-8-10
状态 离线
30
 


QUOTE:
原帖由 快乐的大脚 于 2011-8-12 16:32 发表 bbcodeurl('http://bbs.antpedia.com/images/common/back.gif', '%s')
实验中的real-time PCR数据通过机器直接得到:CT平均值(mean CT)、CT值的标准差(CT-SD),我通过2 -△△Ct 公式计算出ddCT值(相对表达水平),但是不知道如何求ddCT值的SD和P值?

举个例子或许更容易说清楚:

House-keeping gene: GA ...

仪器可以算出相对表达量及其标准偏差,但需要设置好仪器。比如基因,样品号,重复等等,随仪器不通而不同但大同小异。然后再根据这个进行统计分析。例如方差分析,t测验等。一般的统计分析应该都学了吧
顶部