找到约 15 条结果

问题 执行fastqc $datadir/*.gz -o $workdir/1.fastqc报错

...据质控这一节时候，执行上述命令报错，报错为 Exception in thread "Thread-3" java.lang.OutOfMemoryError: GC overhead limit exceeded at java.lang.String.toCharArray(String.java:2899) at uk.ac.babraham.FastQC.Modules.BasicStats.processSequence(BasicStats.java:123) at uk.ac.babrah...

问题 在RStudio中安装WGCNA包后调用WGCNA包提示错误

运行 install.packages("WGCNA") 后提示安装GO.db，IRanges，S4Vectors，Biobase，BiocGenerics，DBI，AnnotationDbi等一些列包，手动安装好这些包后 运行 library(WGCNA) 载入需要的程辑包：dynamicTreeCut 载入需要的程辑包：fastcluster 载入程辑包：‘f...

文章 R语言-Cox比例风险模型

... 同样还是用到以下两个包： library("survival") library("survminer") 计算Cox模型主要用到的是coxph()函数，但需要先用Surv()函数产生一个生存对象；另外coxph()函数支持的方法有：exact，breslow以及exact，默认是exact 以lung数据集为例，看...

文章 遗传图绘制：

代码参考：https://cran.r-project.org/web/packages/LinkageMapView/vignettes/LinkageMapView.html library(LinkageMapView)#示例数据data(carrot)head(carrot)carrot=carrot[order(carrot$group),]outfile="genetic_map_ruler.pdf"lmv.linkage.plot(carrot,outfile,ruler = TRUE,lg.col = "lightblue1")ou...

问题 WGCNA中性状与模块相关性图总出现colorMatrix[, c] : incorrect number of dimensions错误

之前一切运行正常 Error in colorMatrix[, c] : incorrect number of dimensions 我就不知道哪里出问题了。求助 啊

问题 在用tassel对vcf进行sort时报的错，有哪位大佬知道怎么操作吗

...的错，有哪位大佬知道怎么操作吗 ERROR net.maizegenetics.plugindef.ThreadedPluginListener - Out of Memory: SortGenotypeFilePlugin could not complete task: Current Max Heap Size: 29809 Mb Use -Xmx option in start_tassel.pl or start_tassel.bat to increase heap size. Included with tassel ...

文章 R语言跑循环，如何跳过报错信息？

...数： for () {   err = try(想要执行的函数体)}#比如：for (i in 1:length(flies)) {   err = try(eQTpLot(GWAS.df = GWASdata, eQTL.df = eQTLdata , Genes.df = genedata, gene =max_r_gene,           gbuild = "hg19",  trait = phe, tissue = "all" , CollapseMethod = "min",LD.df = LDdata...

文章 GATK Hard-filter 过滤变异结果推荐阈值

...pipeline. These distributions are sharped by both the pipeline methodology and the underlying physical properties of the sequence data; so for a given pairing of data generation technology + analysis pipeline, you can derive filtering thresholds based on what the distributions look like for the trut...

文章 R语言基础入门—数组

... the array.print(result[3,,2])# Print the element in the 1st row and 3rd column of the 1st matrix.print(result[1,3,1])# Print the 2nd Matrix.print(result[,,2]) 运行以上的代码，输出结果如下： COL1 COL2 COL3   5    9   12[1] 11     COL1 COL2 COL...

问题 你好，我是初次接触网络分析绘图，在R中运行以下代码总是报错

...[which(rowSums(table)>=1),]  Abu<-table.generalist总是出现Error in V<-(*tmp*, value = *vtmp*) : invalid indexing这种报错，我想可能是数据格式的问题，但是具体是什么问题也不知道，能麻烦不吝赐教吗，感谢

文章 2区SCI赏析—木麻黄中MYB基因家族鉴定及其在盐胁迫中的响应

...，IF=3.497）题为“Genome-wide analysis of MYB transcription factors and their responses to salt stress in Casuarina equisetifolia.”的基因家族分析文章，看看它是如何发到2区期刊的。 文章简介 MYB是植物中最大的转录因子家族之一，是一类能与真核基...

问题 群体遗传学GWAS课程docker很多命令无法运行

[root@4d714c8732b0  14:14:49 /work/my_gwas2/04.LDdecay]# PopLDdecay -InVCF  ../00.filter/all.varFilter.vcf.gz -SubPop popid.txt -MaxDist 500 -OutStat ld.stat the Number of subPop samples[found in VCF] is 0 sub Group Population szie is too small, SubGroup sample size: 0 ##begin pair-wise R^2 cal...

文章 numpy从文件直接读入得到数组对象

...据读取与写出； 测试数据：https://archive.ics.uci.edu/ml/machine-learning-databases/iris/ import numpy as npimport osos.chdir("D://python_script//")a = np.loadtxt("iris.data", dtype = float64, delimiter=',', usecols = [0, 1, 2, 3])print(a)b = np.loadtxt("iris.data", dtype = "S", d...

问题 重测问题

...再你们服务器上运行一些重测序命令的时候，例如，for i in $(cat $workdir/data/data.txt); do echo "RUN CMD: fastp --thread 50 --qualified_quality_phred 10 \    --unqualified_percent_limit 50 \    --n_base_limit 10 \    -i $datadir/${i}_1.fastq.gz \    -I $datadir/${i}_...

问题 运行go_enrich.r报错

...OG/go_enrich.r -l Private.gene.list -d eggnog/GO_orgdb/ -o Private.go WARNING: ignoring environment value of R_HOME --> Q&A for bioinformatics, please visit the website: https://www.omicsclass.com/ --> R beginners ? I suggest your  learning  R language: https://study.omicsclass.com/in...