摘要：
运行以下命令时出现报错： python ldsc.py --h2 file.sumstats.gz --ref-ld-chr baseline.,eQTL. --w-ld-chr weights.hm3_noMHC. --overlap-annot --print-coefficients --fr 阅读全文
摘要：
原始数据如下所示： 现在想对第一列和第二列进行排序，得到如下结果： 则可以使用代码：sort=ori[order(as.numeric(sub("\\chr+", "", ori$V1)),ori$V2),] 阅读全文
摘要：
当数据中存在单引号时，如下所示： 当使用命令class <- read.table("file",header=T,sep="\t")读取时，就会出现如下所示情况： 所以，正确的读取方式是：class <- read.table("file",header=T,sep="\t",quote = "" 阅读全文
摘要：
运行perl good.pm出现的报错，但问题在bin下面已经有ha/haha.pm文件了。说明perl运行的时候没有找到ha/haha.pm的正确路径。 因此只需要在good.pm中加入haha.pm所在的路径即可，具体如下所示： use lib 'bin'; 阅读全文
摘要：
现有字符串如下所示： merged_file "HAH_chr14_111111_222222_-_HE_111116_chr14_OK_for" 现在我们想把HAH_chr14_111111_222222_-提取出来，那么可以用到命令： library(stringr) split_file <- 阅读全文
摘要：
运行命令plink --file test1 --het --out test_het出现的报错："Error: Line 135 of .ped file has fewer tokens than expected." 出现这种报错有两种原因： 1、map文件的SNP数量与ped的数量对应不上； 阅读全文
摘要：
现有数据框如下所示： 我们需要去除 col1 和 col2 列重复的数据，变成如下所示： 则可以通过命令： library("dplyr") x <- c(0:2,0:2,13:22) y <- c(10:12,10:12,1:10) z <- c("A","B","C","D","A","B"," 阅读全文
摘要：
使用 Rstudio 安装 R 包install.packages("ggsignif")时弹出的错误。 解决办法：改成在 R 下安装就可解决，注意不是在 Rstudio 界面。 阅读全文
摘要：
分别添加完 CluePedia 和 ClueGO 插件后，输入基因名，按照如下顺序即可显示通路中的基因名。 阅读全文
摘要：
解决方案： 1、换个IP地址； 2、换个URL地址； 比如原先是https://www.google.com，改成https://www.google.com.hk/?gws_rd=ssl 阅读全文