t*g
2 楼
准备买个白色牛排32g wifi送人
apple store税后要655
amazon上有个609 这是最便宜的了吗 谢谢
apple store税后要655
amazon上有个609 这是最便宜的了吗 谢谢
S*l
3 楼
最近拿到collaborator的一些RNAseq数据。他们找了个bioinformatics公司处理。公司
给找出来一些treatment相对control significant的gene.我看这些gene都是在control
里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
达差不多的。。。
现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
多谢指点!
给找出来一些treatment相对control significant的gene.我看这些gene都是在control
里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
达差不多的。。。
现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
多谢指点!
y*g
4 楼
自己养两条狗,再养三只猫,盖住原来的味道.
C*e
7 楼
马上A6 + USB 外设了,等等吧
m*c
8 楼
I knew someone is good at this. send me a note if you need help.
H*N
12 楼
借风问一下,对mouse genome 来说,一个sample最少需要多少matched reads?
j*p
14 楼
理想条件,假设library里面有n个rna,每个rna expression 为ni,假设fragments,
sequencing等都是random,那么理论上讲,每个rna被sequence到几率应该和rna的
expression level成正比。问题是有两个不同的rna pool,一个control,一个treat,
这两个sample中rna level 的distribution 不同,total number of rna copies 也不
同,假设两个sample拿去被sequencing的总量一样,并且都被sequence出来10M reads
,那么有可能,在control里面只sequence到expression level最大的1k个rna,其他的
rna reads=0,在treat里面,却sequence到了2k个rna。假设你说的rna在control里面
排第1001,有10个copy,在treat里面也排1001,有10个copy,用sequence,一个没表
达,一个有表达。如果用pcr,两个都有表达。
我认为,不应该是map的问题,the question is how to call significant genes in
RNA-seq. RPKM, which might the way the bioinformatics company used, is
definitely not a decent method.
control
【在 S**********l 的大作中提到】![](/moin_static193/solenoid/img/up.png)
: 最近拿到collaborator的一些RNAseq数据。他们找了个bioinformatics公司处理。公司
: 给找出来一些treatment相对control significant的gene.我看这些gene都是在control
: 里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
: 达差不多的。。。
: 现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
: 多谢指点!
sequencing等都是random,那么理论上讲,每个rna被sequence到几率应该和rna的
expression level成正比。问题是有两个不同的rna pool,一个control,一个treat,
这两个sample中rna level 的distribution 不同,total number of rna copies 也不
同,假设两个sample拿去被sequencing的总量一样,并且都被sequence出来10M reads
,那么有可能,在control里面只sequence到expression level最大的1k个rna,其他的
rna reads=0,在treat里面,却sequence到了2k个rna。假设你说的rna在control里面
排第1001,有10个copy,在treat里面也排1001,有10个copy,用sequence,一个没表
达,一个有表达。如果用pcr,两个都有表达。
我认为,不应该是map的问题,the question is how to call significant genes in
RNA-seq. RPKM, which might the way the bioinformatics company used, is
definitely not a decent method.
control
【在 S**********l 的大作中提到】
![](/moin_static193/solenoid/img/up.png)
: 最近拿到collaborator的一些RNAseq数据。他们找了个bioinformatics公司处理。公司
: 给找出来一些treatment相对control significant的gene.我看这些gene都是在control
: 里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
: 达差不多的。。。
: 现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
: 多谢指点!
t*d
16 楼
难道 RNA-seq 还用这种方法判断有无表达?
microarray 里,Affy 的 P and A 内行是从来不相信的。
一般来说找差异表达的基因应该先把这种在所有样本中都低于或接近可检测下限的的基
因放弃了。
RNA-seq 的检测下限和测序深度有关,lz 的这个例子也许是因为深度不够,该基因的
数据都是噪音而已?
reads
in
【在 j*p 的大作中提到】![](/moin_static193/solenoid/img/up.png)
: 理想条件,假设library里面有n个rna,每个rna expression 为ni,假设fragments,
: sequencing等都是random,那么理论上讲,每个rna被sequence到几率应该和rna的
: expression level成正比。问题是有两个不同的rna pool,一个control,一个treat,
: 这两个sample中rna level 的distribution 不同,total number of rna copies 也不
: 同,假设两个sample拿去被sequencing的总量一样,并且都被sequence出来10M reads
: ,那么有可能,在control里面只sequence到expression level最大的1k个rna,其他的
: rna reads=0,在treat里面,却sequence到了2k个rna。假设你说的rna在control里面
: 排第1001,有10个copy,在treat里面也排1001,有10个copy,用sequence,一个没表
: 达,一个有表达。如果用pcr,两个都有表达。
: 我认为,不应该是map的问题,the question is how to call significant genes in
microarray 里,Affy 的 P and A 内行是从来不相信的。
一般来说找差异表达的基因应该先把这种在所有样本中都低于或接近可检测下限的的基
因放弃了。
RNA-seq 的检测下限和测序深度有关,lz 的这个例子也许是因为深度不够,该基因的
数据都是噪音而已?
reads
in
【在 j*p 的大作中提到】
![](/moin_static193/solenoid/img/up.png)
: 理想条件,假设library里面有n个rna,每个rna expression 为ni,假设fragments,
: sequencing等都是random,那么理论上讲,每个rna被sequence到几率应该和rna的
: expression level成正比。问题是有两个不同的rna pool,一个control,一个treat,
: 这两个sample中rna level 的distribution 不同,total number of rna copies 也不
: 同,假设两个sample拿去被sequencing的总量一样,并且都被sequence出来10M reads
: ,那么有可能,在control里面只sequence到expression level最大的1k个rna,其他的
: rna reads=0,在treat里面,却sequence到了2k个rna。假设你说的rna在control里面
: 排第1001,有10个copy,在treat里面也排1001,有10个copy,用sequence,一个没表
: 达,一个有表达。如果用pcr,两个都有表达。
: 我认为,不应该是map的问题,the question is how to call significant genes in
S*l
17 楼
最近拿到collaborator的一些RNAseq数据。他们找了个bioinformatics公司处理。公司
给找出来一些treatment相对control significant的gene.我看这些gene都是在control
里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
达差不多的。。。
现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
多谢指点!
给找出来一些treatment相对control significant的gene.我看这些gene都是在control
里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
达差不多的。。。
现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
多谢指点!
m*c
19 楼
I knew someone is good at this. send me a note if you need help.
H*N
21 楼
借风问一下,对mouse genome 来说,一个sample最少需要多少matched reads?
j*p
22 楼
理想条件,假设library里面有n个rna,每个rna expression 为ni,假设fragments,
sequencing等都是random,那么理论上讲,每个rna被sequence到几率应该和rna的
expression level成正比。问题是有两个不同的rna pool,一个control,一个treat,
这两个sample中rna level 的distribution 不同,total number of rna copies 也不
同,假设两个sample拿去被sequencing的总量一样,并且都被sequence出来10M reads
,那么有可能,在control里面只sequence到expression level最大的1k个rna,其他的
rna reads=0,在treat里面,却sequence到了2k个rna。假设你说的rna在control里面
排第1001,有10个copy,在treat里面也排1001,有10个copy,用sequence,一个没表
达,一个有表达。如果用pcr,两个都有表达。
我认为,不应该是map的问题,the question is how to call significant genes in
RNA-seq. RPKM, which might the way the bioinformatics company used, is
definitely not a decent method.
control
【在 S**********l 的大作中提到】![](/moin_static193/solenoid/img/up.png)
: 最近拿到collaborator的一些RNAseq数据。他们找了个bioinformatics公司处理。公司
: 给找出来一些treatment相对control significant的gene.我看这些gene都是在control
: 里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
: 达差不多的。。。
: 现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
: 多谢指点!
sequencing等都是random,那么理论上讲,每个rna被sequence到几率应该和rna的
expression level成正比。问题是有两个不同的rna pool,一个control,一个treat,
这两个sample中rna level 的distribution 不同,total number of rna copies 也不
同,假设两个sample拿去被sequencing的总量一样,并且都被sequence出来10M reads
,那么有可能,在control里面只sequence到expression level最大的1k个rna,其他的
rna reads=0,在treat里面,却sequence到了2k个rna。假设你说的rna在control里面
排第1001,有10个copy,在treat里面也排1001,有10个copy,用sequence,一个没表
达,一个有表达。如果用pcr,两个都有表达。
我认为,不应该是map的问题,the question is how to call significant genes in
RNA-seq. RPKM, which might the way the bioinformatics company used, is
definitely not a decent method.
control
【在 S**********l 的大作中提到】
![](/moin_static193/solenoid/img/up.png)
: 最近拿到collaborator的一些RNAseq数据。他们找了个bioinformatics公司处理。公司
: 给找出来一些treatment相对control significant的gene.我看这些gene都是在control
: 里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
: 达差不多的。。。
: 现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
: 多谢指点!
t*d
23 楼
难道 RNA-seq 还用这种方法判断有无表达?
microarray 里,Affy 的 P and A 内行是从来不相信的。
一般来说找差异表达的基因应该先把这种在所有样本中都低于或接近可检测下限的的基
因放弃了。
RNA-seq 的检测下限和测序深度有关,lz 的这个例子也许是因为深度不够,该基因的
数据都是噪音而已?
reads
in
【在 j*p 的大作中提到】![](/moin_static193/solenoid/img/up.png)
: 理想条件,假设library里面有n个rna,每个rna expression 为ni,假设fragments,
: sequencing等都是random,那么理论上讲,每个rna被sequence到几率应该和rna的
: expression level成正比。问题是有两个不同的rna pool,一个control,一个treat,
: 这两个sample中rna level 的distribution 不同,total number of rna copies 也不
: 同,假设两个sample拿去被sequencing的总量一样,并且都被sequence出来10M reads
: ,那么有可能,在control里面只sequence到expression level最大的1k个rna,其他的
: rna reads=0,在treat里面,却sequence到了2k个rna。假设你说的rna在control里面
: 排第1001,有10个copy,在treat里面也排1001,有10个copy,用sequence,一个没表
: 达,一个有表达。如果用pcr,两个都有表达。
: 我认为,不应该是map的问题,the question is how to call significant genes in
microarray 里,Affy 的 P and A 内行是从来不相信的。
一般来说找差异表达的基因应该先把这种在所有样本中都低于或接近可检测下限的的基
因放弃了。
RNA-seq 的检测下限和测序深度有关,lz 的这个例子也许是因为深度不够,该基因的
数据都是噪音而已?
reads
in
【在 j*p 的大作中提到】
![](/moin_static193/solenoid/img/up.png)
: 理想条件,假设library里面有n个rna,每个rna expression 为ni,假设fragments,
: sequencing等都是random,那么理论上讲,每个rna被sequence到几率应该和rna的
: expression level成正比。问题是有两个不同的rna pool,一个control,一个treat,
: 这两个sample中rna level 的distribution 不同,total number of rna copies 也不
: 同,假设两个sample拿去被sequencing的总量一样,并且都被sequence出来10M reads
: ,那么有可能,在control里面只sequence到expression level最大的1k个rna,其他的
: rna reads=0,在treat里面,却sequence到了2k个rna。假设你说的rna在control里面
: 排第1001,有10个copy,在treat里面也排1001,有10个copy,用sequence,一个没表
: 达,一个有表达。如果用pcr,两个都有表达。
: 我认为,不应该是map的问题,the question is how to call significant genes in
l*1
24 楼
try test this kind of Cloud mapping platform to NGS 2.0 or 3.0 version database:
//www.ncbi.nlm.nih.gov/pubmed/21645377
//sourceforge.net/projects/cloudaligner/
//cloudaligner.sourceforge.net/
control
【在 S**********l 的大作中提到】![](/moin_static193/solenoid/img/up.png)
: 最近拿到collaborator的一些RNAseq数据。他们找了个bioinformatics公司处理。公司
: 给找出来一些treatment相对control significant的gene.我看这些gene都是在control
: 里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
: 达差不多的。。。
: 现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
: 多谢指点!
//www.ncbi.nlm.nih.gov/pubmed/21645377
//sourceforge.net/projects/cloudaligner/
//cloudaligner.sourceforge.net/
control
【在 S**********l 的大作中提到】
![](/moin_static193/solenoid/img/up.png)
: 最近拿到collaborator的一些RNAseq数据。他们找了个bioinformatics公司处理。公司
: 给找出来一些treatment相对control significant的gene.我看这些gene都是在control
: 里面根本没有map一个read的。回去做pcr,果然不对,这些gene都是在两个sample里表
: 达差不多的。。。
: 现在的问题是。对那些map reads极小的。怎么知道是map的部队还是没有表达?
: 多谢指点!
a*n
25 楼
别去管那些copy数太低的, 就这么简单
RNA-seq给的信息那么多, 不会跟你hypothesis相符的基因都是这些copy数很低的吧
RNA-seq给的信息那么多, 不会跟你hypothesis相符的基因都是这些copy数很低的吧
相关阅读
张三不转行,今天喝完新年酒习皇拟取消两院生命科学部和医学部青千明年大概率取消PI现在都要傻眼了LOL (转载)加微信txm149,进入全球华人高端私密交友平台。请弟兄们帮忙下一篇论文的Supplementary Table 1。 谢谢如果不是贺建奎,而是早年的Elizabeth Holmes做出转基因人paper help贺建奎证明了科学和科学家的虚伪性科技出版都有什么诀窍?如何办好一本生物医学期刊并被SCI收录科学界的特大事故,负面影响将无法估量Re: 生物老千有研究这个的吗,叔的疑问是”为啥听音乐能产生欣(转载)千老转行花样繁多转行3d动画制作如何?寻求生信生统领域期刊审稿机会生物化学在药厂也pay的很好嘛一手拿美国钱,一手替中国搞科研,“千人计划”再惹争议hejiankui 后续?椰子油是超级食品还是心脏病杀手严肃地请教一下癌症的诊断标准: (转载)