找到约 15 条结果

文章 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...

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

...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 已...

问题 基因家族解析报错

根据课程学习中，执行代码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.

文章 转录组+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.

问题 老师，请问做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....

问题 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  \ ...

文章 IF=4.8 | 联合转录组与代谢组分析揭示了硒对百合生长和品质的影响

...tiers in plant science杂志发表了一篇题为“Combined transcriptome and metabolome analyses reveal the effects of selenium on the growth and quality of Lilium lancifolium”的转录组代谢组联合分析文章。” 01 — 研究背景 卷丹百合（学名：Lilium lancifolium Thunb...

文章 正选择分析之 Site Models

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

问题 请问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

问题 检查数据是否正确

...=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 这是为什么呢

问题 perl求助，基因ID的对应替换perl脚本怎么写

...的Bna基因ID对应替换上去，perl脚本该怎么写。。 #!/usr/bin/perl ;#usage：perl repalce.pl list.file fa.file out.fileopen LIST,$ARGV[0];open FAFILE,$ARGV[1];open OUT,">$ARGV[2]";undef $/;#$in=<FAFILE>;#study $in;$/="\n";#while(<LIST>){    chomp;    my($new,$old)=...

问题 perl求助，基因ID的对应替换perl脚本怎么写

...的Bna基因ID对应替换上去，perl脚本该怎么写。。 #!/usr/bin/perl ;#usage：perl repalce.pl list.file fa.file out.fileopen LIST,$ARGV[0];open FAFILE,$ARGV[1];open OUT,">$ARGV[2]";undef $/;#$in=<FAFILE>;#study $in;$/="\n";#while(<LIST>){    chomp;    my($new,$old)=...