找到约 15 条结果

文章 免疫 侵润 分析 ESTIMATE分析

...A_download/TCGA-STAD_gene_info.csv",  stringsAsFactors=FALSE, row.names=1,check.names = F) sample.info<-read.csv("../01.TCGA_download/TCGA-STAD_sample_info.csv",  stringsAsFactors=FALSE, row.names=1,check.names = F) geneExpSet.fpkm<-read.csv("../01.TCGA_download/TCGA-STAD_gene_expressi...

文章 解决基因组当中蛋白质序列ID和gff中ID不一致的问题

...form.gff3#提取cds序列/share/work/biosoft/TransDecoder/latest/util/gff3_file_to_proteins.pl  --gff3  $species.longest_isoform.gff3 --fasta $genome  --seqType CDS >$species.cds.fa#提取pep序列/share/work/biosoft/TransDecoder/latest/util/gff3_file_to_proteins.pl  --gff3  $species.longest...

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

... 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 <-as.vector(...

文章 qiime2自己建立分类器 greengene

...,可直接在仓库中下载( http://kronos.pharmacology.dal.ca/public_files/taxa_classifiers/qiime2-2020.2_classifiers/ ),没有的话则需要自己手动建立分类器: •16S V4/V5 region (classifier_silva_132_99_16S_V4.V5_515F_926R.qza) 16S V3/V4 region (classifier_silva_132_99_16S...

文章 快速高效KM生存曲线绘制-脚本使用

...etadata input the path to the clinical information file containing survival data[required] --surv_cut Whether to group by surv_cutpoint function[optional,default:False] --median_cut Whether to group by me...

问题 老师您好,我在用MEGA11建系统发育树时出现了错误的提示

MEGA failed to parse the input file:Invalid base found:1(in line 4)

文章 绘制HIC热图

...ENOME_FRAGMENT,bed文件路径(绝对路径) config-hicpro.txt # Please change the variable settings below if necessary########################################################################### Paths and Settings  - Do not edit !###########################################################...

文章 R 安装包:configure: error: proj_api.h not found in standard or given locations.

hecking for strings.h... yeschecking for sys/stat.h... yeschecking for sys/types.h... yeschecking for unistd.h... yeschecking for gdal.h... yeschecking GDAL: linking with --libs only... yeschecking GDAL: /usr/share/gdal/pcs.csv readable... nochecking GDAL: checking whether PROJ is available for link...

文章 安装picrust2 保存:pkg_resources.DistributionNotFound: The 'cython' distribution was not found and is required by biom-format

...error of the above failed command: Traceback (most recent call last):   File "/share/work/biosoft/conda/Miniconda3.8/envs/picrust2/bin/picrust2_pipeline.py", line 4, in <module>     __import__('pkg_resources').require('PICRUSt2==2.3.0b0')   File "/share/work/biosoft/conda/Miniconda3.8/...

文章 multi_cox.r 基因表达量做多因素cox分析

...     show this help message and exit  -i data, --data data  input data file path[required]  -t time, --time time  set suvival time column name [required]  -e event, --event event                        set event column name must 0 or 1 code format                      ...

文章 单细胞转录组-inferCNV拷贝数变异分析介绍

...s/infercnv.r -i ../BCC_GSE123813.qs  \   -r "T_cells"   --annotations_file  ../cellanno_selected.tsv  --gene_location ../hg38_gencode_v27.txt \      --cpu 20 --hmm  --denoise  --analysis_mode samples    如果group_by_cluster = TRUE ,左边的树状图是所有观察细胞的,...

文章 利用vcftools分析玉米和大刍草基因组区域或基因区域是否受到选择

...1000 --out maize_1000.pi 3. R软件可视化该区域 zemaz = read.table(file='~/maize_hapmap/maize_1000.pi',header = T) dachucao = read.table(file='~/maize_hapmap/dachucao_1000.pi',header = T) zemaz$species ='maize' dachucao$species ='dachucao' df_plot = rbind(zemaz,dachucao) library(ggplot2...

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

... set event column name [required] -i immu, --immu immu input data file path[required] -m method, --method method a character string indicating which correlation coefficient (or covariance) is to be computed. One of "p...

文章 split_sample.r根据生存数据对样本随机分组

...et,default:0.6 -o outdir, --outdir outdir output file directory [default cwd] -n name, --name name output file name prefix [default demo] 参数说明: -i    输入矩阵文件 -e    文件中生存事件列名,默认为EVENT -t    文件中生存时间...

问题 这个脚本如果直接从ensamble下载cds和pep还需要去提取吗,是不是可以直接转换格式

...   #提取cds序列   /share/work/biosoft/TransDecoder/latest/util/gff3_file_to_proteins.pl  --gff3  $species.longest_isoform.gff3 --fasta $genome  --seqType CDS >$species.cds.fa   #提取pep序列   /share/work/biosoft/TransDecoder/latest/util/gff3_file_to_proteins.pl  --gff3  $speci...