omicsgene
安生水
使用DNA序列做主成分分析（PCA）——R语言adegenet包
提到主成分分析，一般我们都是使用Plink，GCTA等软件基于SNP数据来操作，那么如何用DNA序列做主成分分析呢？ 思路是先比对，之后使用R语言的adegenet包把比对的数据转换成snp数据，用到的函数是fasta2genlight()，再进行PCA分析及绘图。
- 0
- 0
- 星莓
- 发布于 2022-10-31 11:33
- 阅读 ( 4888 )
星莓
python操作大文件删除某行
删掉文件其中一行，非常简单，不过如果一个文件很大以至于无法读到内存里，又该怎么操作呢 其实我们可以使用 open() 方法把需要修改的文件打开为两个文件，然后逐行读入内存，找到需要删除的行...
- 0
- 0
- xun
- 发布于 2022-10-24 16:10
- 阅读 ( 1528 )
xun
Perl方法对文档全文进行字符串对应替换
举个例子：从数据库中下载基因组数据，基因组.gff文件中染色体ID较为复杂（第一列） 从全基因组序列.genome.fa文件中找到染色体ID对应简写 想将gff文档中染色体ID全部替换成第二列LC*命名...
- 0
- 0
- rzx
- 发布于 2022-10-24 14:38
- 阅读 ( 1431 )
rzx
sort 按照染色体排序后按位置坐标排序
经常会遇到对bed文件或其他相关数据的排序，记录一下如何先按照染色体号排序，然后按照坐标位置排序。 sort -k1,1 -k2n file1 > file2 排序完成后如下：
- 0
- 0
- 星莓
- 发布于 2022-10-24 09:55
- 阅读 ( 2041 )
生信阿姨
GeMoMa-1.9.jar 报错mbol: g_settings_backend_get_type
mbol: g_settings_backend_get_type
- 0
- 0
- omicsgene
- 发布于 2022-10-19 15:19
- 阅读 ( 1282 )
linux删除某用户下所有的进程命令
1. ps列出用户的pid，然后依次kill掉，比较繁琐 kill -9 `ps -uf username |awk '{ print $2 }'|grep -v PID` 或 ps -ef | grep username | awk '{ print $2 }' | xargs kill -9 2. pkill方...
- 0
- 0
- 星莓
- 发布于 2022-10-14 10:21
- 阅读 ( 1408 )
awk统计文件某列中某关键词出现次数
有时我们需要统计文件的某列中某些词出现的次数，可以在Linux下使用sort | uniq -c，那用awk怎么操作呢?
- 0
- 0
- 星莓
- 发布于 2022-10-12 18:29
- 阅读 ( 1659 )
AI保存PNG/PDF/JPG时“出现了未知错误”的解决办法
AI保存PNG/PDF/JPG时“出现了未知错误”的解决办法
- 0
- 0
- 生信阿姨
- 发布于 2022-10-12 15:08
- 阅读 ( 17071 )
关于CodeML分析中四种常见模型假设
codeml是PAML软件包下的一个程序，在估算蛋白编码序列同义替换和非同义替换速率以及检测序列是否经受正选择分析中受到广泛应用。 对于一棵系统发育进化树，在考虑其不同序列间（site）和不同支...
- 0
- 0
- rzx
- 发布于 2022-10-10 15:41
- 阅读 ( 3712 )
linux将一个文件复制到多个文件夹
cp命令可以简单的将多个路径的文件复制到一个文件夹内，但反过来就会麻烦很多 不过和其他命令搭配也可以实现 echo 路径1 路径2|xargs -n 1 cp -v 需要复制的文件 当然，路径可以一直加，也可...
- 0
- 0
- xun
- 发布于 2022-10-09 09:33
- 阅读 ( 1765 )
