找到约 15 条结果

文章 小尺度的种群历史分析——∂a∂i原理解析

...文章 “A genomic variation map provides insights into the genetic basis of cucumber domestication and diversity”，在这个文章中，作者对来自印度、西双版纳、欧洲和东亚四个地方的黄瓜进行了种群分化的研究，同时通过遗传多样性和选择性清除定位...

问题 比较基因组绘图时报错

我的配置文件： seqid 1,2,3,4,5A01,A02,A03,A04,A05,A06,A07,A08,A09,A10,A11,A12,A13,D01,D02,D03,D04,D05,D06,D07,D08,D09,D10,D11,D12,D13 layout # y, xstart, xend, rotation, color, label, va,  bed.8,     .2,    .75,       0,      red, Gbar, top, gbar.bed.2,     .2, ...

文章 pip报错

...d conflicting behaviour with the system package manager. It is recommended to use a virtual environment instead: https://pip.pypa.io/warnings/venv 以“root”用户运行pip可能导致权限中断，可以尝试使用虚拟环境 1，创建环境 python3 -m venv tutorial-env 如出现报错...

问题 关于图形泛基因组中的Minigraph-Cactus的问题。

...的三种文件，其中 clip graph: This is the default graph. Stretches of sequence >10kb that were not aligned to the underlying SV/minigraph are removed. "Dangling" nodes (ie that don't have an edge on each side) that aren't on the reference path are also removed, so that each chromosome only...

问题 重测序var_density.R脚本报错

... in seq.default(0, chorm.maxlen/bp, round(xticks[2])) :   invalid '(to - from)/by' Calls: CMplot -> DensityPlot -> seq -> seq.default Execution halted

问题 重测序var_density.R脚本报错

... in seq.default(0, chorm.maxlen/bp, round(xticks[2])) :   invalid '(to - from)/by' Calls: CMplot -> DensityPlot -> seq -> seq.default Execution halted

问题 R包安装

 Warning message: package ‘edgeR’ is not available for this version of R A version of this package for your version of R might be available elsewhere, see the ideas at https://cran.r-project.org/doc/manuals/r-patched/R-admin.html#Installing-packages    请问老师，出现这种情况...

问题 RNA 速率分析Pseudotime vs latent time

...-Seq analysis, and I am doing RNA velocity and pseudotime analyses as part of my analysis pipeline. Now, RNA velocity is carried out using scvelo, and this package generates a metric known as latent time. I am having a very had time seeing the difference between latent time and standard pseudotime a...

文章 R语言is.null()函数

...都没有，变量为空。     判断向量是否空值的函数：is.null() 2. R语言中，NA代表位置上的值为空。     判断向量中的元素是否缺失值：is.na() 示例： $var_name=NULL;if (is.null($var_name)){print "变量var_name为NULL";}else{print "变量var_name...

文章 eggNOG 注注释报near "WITHOUT": syntax error

... for CentOS7 users.Because CentOS7 users should install sqlite(>=3.8.2) from source, CentOS7 users need to set export PATH=/share/work/biosoft/SQLite/latest/bin:$PATH export LD_LIBRARY_PATH=/share/work/biosoft/SQLite/latest/lib/:$LD_LIBRARY_PATH

文章 R包 CandiHap 进行单倍型分析

...Phenotype.txt",sep="\t")hap <- snp2hap(pheno,ovl, hapname="Hap")## want to extract resultsHap_results <- results(hap, gene="Si9g49990")write.table(Hap_results, file = "Hap_results.txt", sep = "\t", quote = FALSE, col.names = T, row.names = F)## Plot Gene LDheatmapif (! require("LDheatmap")) in...

文章 MFP - Molecular Functional Portrait

... Python3.7import pandas as pdfrom portraits.clustering import clustering_profile_metrics, clustering_profile_metrics_plotfrom portraits.utils import read_gene_sets, ssgsea_formula, median_scale# Example script# Read signaturesgmt = read_gene_sets('signatures.gmt')  # GMT format like in MSIGdb# Read...

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

Table S1. Example of R script that was used to simulate a trait, having Heritability=60% and number of QTLs=20 (H60_Q20), using the marker 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...

文章 perl语言——正则表达式

正则表达式(regular expression)描述了一种字符串匹配的模式，可以用来检查一个串是否含有某种子串、将匹配的子串做替换或者从某个串中取出符合某个条件的子串等。 Perl语言的正则表达式功能非常强大，基本上是常用语言中最...

问题 无法链接到docker守护进程

docker 容器自动退出，再连就报错：Cannot connect to the Docker daemon at unix:///var/run/docker.sock. Is the docker daemon running? 我检查docker运行状态显示正在运行： docker.service - Docker Application Container Engine Loaded: loaded (/lib/systemd/system/docker.servic...