找到约 15 条结果

问题 采用RSeQC 对bam 文件进行质控分析时,用inner_distance.py报错。

...实验质量评估,采用RSeQC 对bam 文件进行质控分析时,用inner_distance.py报错。 报错信息: manager@bl8vbox[QC] inner_distance.py -i /home/manager/rnaseq/work/3.alignments/hisat2/UHR.bam    -r stringtie_merged.bed  -o UHR_inner_size Traceback (most recent call last): ...

文章 使用R语言安装包编译时gcc标准报错——“error: for loop initial declarations are only allowed in C99 mode”

使用BiocManager安装Biostrings包的时候报错了: BiocManager::install("Biostrings",force=TRUE) 报错信息提示是gcc标准的问题,应该使用C99标准,可以使用withr包指定使用标准: with_makevars(c(PKG_CFLAGS = "-std=c99"), + BiocManager::install("Biostrings",force=T...

问题 非模式物种建立注释库问题

...column types or set `show_col_types = FALSE` to quiet this message. Error in file(file, "rt") : cannot open the connection Calls: read.table -> file In addition: Warning message: In file(file, "rt") :   cannot open file 'NA/pathway2name.tsv': No such file or directory Execution halted 已...

文章 GATK外显子测序数据肿瘤体细胞突变数据分析

...rences_hg38_v0_Homo_sapiens_assembly38.dbsnp138.vcf \ --known-sites Mills_and_1000G_gold_standard.indels.hg38.vcf.gz -I N154.sort.dedup.bam -O N154.sort.dedup.bam.tablegatk BaseRecalibrator -R Homo_sapiens_assembly38.fasta \ --known-sites references_hg38_v0_Homo_sapiens_assembly38.dbsnp138.vc...

文章 MYSQL数据库常用语句命令

...    # 返回集合中最小的值1SELECT mod(5,3);    # 余数2SELECT rand();    # 返回0到1内的随机值,每次不一样SELECT rand(5);   # 提供一个参数(种子)使RAND()随机数生成器生成一个指定的值。SELECT round(1415.1415);   # 四舍五入1415SELECT round(1...

问题 基因家族解析报错

根据课程学习中,执行代码orthofinder.py  -f pep  -S diamond  -M msa  -A muscle  -T fasttree  -t 1 -a 1 后,后台nohup文件报错如下,但是查看了文件权限和文件的实际路径都有该文件,且文件大小不为0。这几天运行了好几次,都是卡在这个...

文章 FGF21突变让你吃多甜食也不发胖?

...gar Intake Is Associated with Body Shape, Lower Total Body-Fat Percentage, and Higher Blood Pressure[J]. Cell Reports, 2018, 23(2): 327-336.

问题 老师,请问做Conet 构建的时候只有两个样本,但是有多个环境因子能做图吗?我导入数据报错了“After filtering, less than 3 columns remain in the input matrix!”请问该如何解决呢?

...图? Error message: java.lang.IllegalArgumentException: After filtering, less than 3 columns remain in the input matrix! at be.ac.vub.bsb.cooccurrence.conversion.MatrixFilterer.filter(MatrixFilterer.java:376) at be.ac.vub.bsb.cooccurrence.conversion.Preprocessor.processMatrix(Preprocessor....

文章 转录组+sRNA+降解组=7分文章

...文献:Han X, Yin H, Song X, et al. Integration of smallRNAs, degradome and transcriptome sequencing in hyperaccumulatorSedum alfrediiuncovers a complex regulatory network and provides insights into cadmium phytoremediation:[J]. Plant Biotechnology Journal, 2016, 14(6):1470-1483.

问题 WGCNA分析中,导入txt的性状数据时(数据为12行样本,10列性状数据),报错Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : line 1 didot have 12 elements。但是数据如果是12行样本,8列性状,不会报错

...xt",row.names=1,header=T,comment.char = "",check.names=F)但是报错Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec,  :    line 1 did not have 12 elements。但是数据如果是12行样本,8列性状,就不会报错。谢谢

问题 在运行这段代码时出现Invalid value for '--m-metadata-file': Metadata column is of type 'numeric', but expected Categorical.

...自metadata文件  #--p-mode <up to you, there's sum, median-ceiling, and mean-ceiling>   qiime feature-table group \      --i-table feature-table-final.qza \      --p-axis sample \      --m-metadata-file $metadata \      --m-metadata-column $g \      --p-mode sum  \ ...

问题 请问QQ图绘制时报错,文件明明在,却找不到文件怎么回事啊?

> library(qqman) > results_log = read.table("logistic_results.assoc.logistic", head=TRUE) Error in file(file, "rt") : cannot open the connection In addition: Warning message: In file(file, "rt") :   cannot open file 'logistic_results.assoc.logistic': No such file or directory

文章 正选择分析之 Site Models

...作用的检测,而是用于位点间ω值是否一致的检测。 M1 and M2 以及M7 and M8是用于正选择作用的检测,作者推荐使用这两组比较进行LRT检验来验证正选择。不过Prof.Yang认为,The M1-M2 comparison 与 the M7- M8 comparison相比,更加的稳定。M7...

问题 检查数据是否正确

...=1,header=T,comment.char = "",check.names=F) > allowWGCNAThreads() Allowing multi-threading with up to 8 threads. > ALLOW_WGCNA_THREADS=4 > memory.limit(size = 20000) [1] 32457 Warning message: In memory.limit(size = 20000) : cannot decrease memory limit: ignored > # 查看部分样本...

问题 vcftools --gzvcf all.varFilter.vcf.gz --recode --recode-IN.............els |gzip - > clean.vcf.gz这一步输出的vcf.gz文件,用 bcftools query -l clean.vcf.gz | wc -l查看样本名出现报错

Failed to read from clean.vcf.gz: unknown file type 0 这是为什么呢