今天本打算看3个视频,结果应该是算看了2个,第二个没有看完,而第一个只是挑挑拣拣的看了一部分,可能第二个视频的内容太多了,每个新的知识点都需要我去延伸,虽然这只是对以前看过的东西的复习,但是可见以前看得有多么囫囵吞枣,很多内容都没有理解透彻。
1.就拿向量化循环来说,之前来来回回看了有2遍,都没有注意到这个概念,在r中,当对向量运算时,不需要写循环,自动对向量里面的每个值进行循环预算。
2.在r中的循环函数中,apply(a,1,function(x){}), 老师说老老实实的记住就行了,其中1是对行进行循环计算,而2是对列进行循环计算,但是对于取矩阵中符合要求的行,比如:a[apply(a,1,function(x),sum(x)>floor(ncol(a)/50)),],其中里面那个循环计算出来应该是每一排有一个逻辑值,那么a[x,y],x还可以是逻辑值? 疑惑。
3.还有之前一直觉得CPM很复杂,一看到数学的计算公司就头疼,今天就像忽然醒悟一样,不就是稍微换算一下,把基因表达值除个reads数,再乘个一百万就可以了,去除了基因间的文库大小的差异。
4.看rpkm值也是,听老师讲了一下倒很懂,但是去发现网上有各种各样的定义,后来发现,只看一个定义实在是不好的,还是要多看看,综合去理解的,rpkm值既然是错误的,那么最好还是不用吧,但是我的那个数据一下载下来就是rpkm值,我也没办法。
5.至于spike-in和ERCC,也是标准化范畴里面的概念,spike-in就和housekeepinggene一样,是对基因表达数据进行标准化的,而ERCC,external rna control cortium,是spike-in里,