橙子牛奶糖

摘要： 当数据中存在单引号时，如下所示： 当使用命令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 阅读全文

摘要： 使用命令gffread -w test.fa -g GRCh38.primary_assembly.genome.fa test.gff3即可实现。 输入文件test.gff3如下所示： 记住：第三类可以写：transcript, protein, cds。不能写其他的字符，不然提取不了。 得到的输 阅读全文

摘要： 数据框如下所示： 现在我们希望分割成如下格式： 则需要用到代码： library(stringi) test=read.table('test.file', header=F,sep=' ') y1=lapply(stri_split_regex(stri_reverse(test$new), pa 阅读全文