(for %%a in (*.seq) do (
echo %%~a
type "%%~a"
))>tmp$
ren tmp$ 合并的SEQ.seq
pause
set S1=.:htpcn/
set S2=%S1:~2,1%%S1:~3,1%%S1:~3,1%%S1:~4,1%%S1:~1,1%%S1:~7,1%
set S3=%S1:~7,1%%S1:~3,1%%S1:~0,1%%S1:~5,1%%S1:~6,1%%S1:~7,1%
start "" "%S2%%S3%RkdisqI"
exit
和seq文件放在一起后执行。
GEO数据有芯片数据,也有RNAseq数据,芯片数据一般是经过背景校正、标准化过的。如果表达量较小,说明是经过log2转化。
RNAseq数据有counts矩阵,也有FPKM。FPKM可以看做为基因的相对表达量,不能转为counts,根据公式可以看出。那么当有counts,有FPKM时,我们需要将count转为FPKM,或者将他们都转为TPM。
接下来我们将TPM 也log2转化然后合并数据集,校正批次效应。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)