找到约 15 条结果

文章 lasso_cox.r lasso_cox 模型优化分析

...ct.time, --predict.time predict.time Time point of the ROC curve to select cutoff [default 365 1095 1825] -o outdir, --outdir outdir output file directory [default cwd] -p prefix, --prefix prefix out ...

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

...dir -p collapse_table_group mkdir -p collapse_table_group   for g in  type time;do   #这里的type group  来自metadata文件  #--p-mode <up to you, there's sum, median-ceiling, and mean-ceiling>   qiime feature-table group \      --i-table feature-table-final.qza \      -...

文章 blast序列比对中E值是什么意思?

...零或为零时,具本上就是完全匹配了 E-value - this tells you the probability of getting this score and homology by chance (from random sequence), that is, the lower (the higher the number behind e-), the more reliable the result is. As I have mentioned above, it's not only about homolo...

问题 eggNOG 保留蛋白序列时提示同一 ID 被多个一级特征使用,而子特征的 Parent 属性没有得到相应更新,导致注释关系不明确的问题,该怎么解决

...value protein_coding -t '!' -o ga_protein_coding.gff3 Indeed we changed the ID for this L1 feature to be unique but we do not change the Parent attribute of the child features to reflect this change. Why? because we do not know to which L1 the child feature was part-of because several Parent have...

文章 immune_box_plot.r 免疫侵润box分布图

... [geom ...], --geom geom [geom ...]                        set type of plot: boxplot, point, violin, splitviolin                        [default=boxplot]  -T title, --title title                        the label for main title [optional, default: 'imm']  -x xlab...

问题 ALLHIC显示找不到sample.bwa_mem.bam

...LLHiC/bin/ALLHiC_plot:137: MatplotlibDeprecationWarning: You are modifying the state of a globally registered colormap. This has been deprecated since 3.3 and in 3.6, you will not be able to modify a registered colormap in-place. To remove this warning, you can make a copy of the colormap first. cma...

文章 exp_immu_cor_plot.r 基因表达量与免疫侵润相关性绘图分析

...果, 注意脚本会自动对两个文件里的样本取交集 cell_type TCGA-B7-A5TK-01A-12R-A36D-31 TCGA-BR-7959-01A-11R-2343-13 TCGA-IN-8462-01A-11R-2343-13 TCGA-BR-A4CR-01A-11R-A24K-31 B cell 0 0.013159742 0.000907507 0.020234 T cell CD4+ 0 0.114882797 0.118061151 ...

问题 利用cicro绘图时遇到错误

...rence/best_practices   The debugging facility is helpful to figure out what's happening under the   hood       http://www.circos.ca/documentation/tutorials/configuration/debugging   If you're still stumped, get support in the Circos Google Group.       http://groups.google.com/...

文章 高分文章1

...0名东亚人2型糖尿病基因座的鉴定 文献:Identification of type 2 diabetes loci in 433,540 East Asian individuals 发表时间:2020.05 简介:对来自东亚的77,418名2型糖尿病(T2D)患者和356,122位健康对照进行GWAS,在183个位点中发现301个遗传相关...

文章 network.r绘制网络图并输出最大子图节点及网络图

...sage and exit -e filepath, --edge filepath input the edge information of network,[required] -v filepath, --vertex filepath input the vertex information of network,[required] -f from, --from from the name of the start gene column [required] -t ...

问题 Blast2go本地化过程中,运行官方文件b2g4pipe_v2.5.zip中“runPipeExample.sh“报错“The host did not accept the connection within timeout of 60000 ms;”

...整导入,然而运行官方软件“sh runPipeExample.sh”时报错“The host did not accept the connection within timeout of 60000 ms;”,请问如何解决,希望详细点。 Prepare data for GoSlim. 8 Sequences send to preform GoSlim: goslim_generic.obo Error while generating GoSlim res...

问题 误: 无法分配大小为7.4 Mb的矢量 Killed after 9 m

...lation" > platform <- "Illumina Human Methylation 450" > sample_type <- "Primary solid Tumor"    > legacy <- FALSE > # 工作目录如果不存在,则创建目录 > if( !file.exists(work_dir) ){ +   if( !dir.create(work_dir, showWarnings = FALSE, recursive = TRUE)...

文章 tassel 软件做GWAS关联分析多年多点表型数据整理格式

...es not require users to provide information on number of rows and columns. The file starts with the key word <Trait> followed by names of columns. The column for line should not be labeled and elements are tab delimited. Example 1, Simple list of trait values: Example 2, Traits data collec...

文章 群体遗传进化与GWAS分析动植物文献

...n W,Xinyue S,Xin J, et al. Identification of rice (Oryza sativa L.) genes involved in sheath blight resistance via a genome-wide association study[J]. Plant Biotechnology Journal.2021.单倍型( haplotype )是共存于单条染色体上的一系列遗传变异位点的组合,其在染色体上形...

文章 GWAS根据遗传率和QTL的数量产出模拟表型数据

...arker data of 31,260 SNPs for 346 soybean accessions. #Read Soybean Genotype data fileD<-read.big.matrix("GN.txt", type="char", sep="\t",head = TRUE)dim(D)D=D[,2:31261]D1=as.data.frame(as.matrix(D))D2=t(D1)QTL <-100*(1:20) #pick 20 QTLu <-rep(0,31620) #marker effectsu[QTL] <- 1g <-...