找到约 15 条结果

文章 linux将一个文件复制到多个文件夹

...用特殊符号表示 比如想给~/test/中所有文件夹都加一个~/file.txt的文件 echo ~/test/*/|xargs -n 1 cp -v ~/file.txt 需要复制的文件夹多的时候还是比较方便的

文章 loom 文件转换成Seurat对象

...0220526/ann1/3.seurat_cluster/B0507-old/cell_ann")brain.loom <- Connect(filename = "l5_all.loom", mode = "r")brain <- as.Seurat(brain.loom)saveRDS(brain, file = "Adolescent.l5_all.rds")write.table(data.frame(barcode=row.names(brain@meta.data),brain@meta.data),file="Adolescent.l5_all.metadata.t...

文章 服务器docker防止退出之后容器杀死导致任务中断

...ontainer af05fafe31f6bfb00c2599255c47e35813ecf5af9bbe6760ae8a4abffd343627: file exists ERRO[0000] Error refreshing container b4620633d99f156bb59eb327a918220d67145f8198d1c42b90d81e6cc29cbd6b: error acquiring lock 2 for container b4620633d99f156bb59eb327a918220d67145f8198d1c42b90d81e6cc29cbd6b: file e...

文章 单细胞转录组数据挖掘流程记录-BRCA乳腺癌(E-MTAB-8107)

...格按sample ID合并,注意保存txt格式 Rscript $scripts/merge_tsv_files.r -i GSE210347_meta.txt GSE210347_study_metadata.tsv -b SampleID -p metadata_all​ #筛选乳腺癌数据的metadata cat metadata_all.tsv |awk 'NR==1|| ($0~"E-MTAB-8107" && $0~/Breast/){print $0}'>E-MTAB-810...

文章 基因组版本 Primary 和 Alternate的区别

...th long stretches of phased blocks. The concept has been used by GRC. BAC-to-BAC assemblies can all be regarded as primary assemblies. Falcon-unzip is perhaps the first to produce such assemblies for whole-genome shotgun data.Alternate assembly: an incomplete assembly consisting of haplotigs in het...

文章 “特朗普”成为科研新热点,100多篇CNS都有他......

...ure在线发表了一篇研究报告,标题是——Trump’s plan would make government stupid 翻译一下:特朗普的计划会让政府变得愚蠢 这位作者可是地地道道的美国人,不得不说,他们“研究”起自己的大Boss可真是毫不留情面, 事实上,...

文章 借助perl脚本进行文件整理 | 将按第一列信息整理的文件改成按照第二列信息整理

... Cyellow  B有perl脚本如下:print "perl $0  <id_list> <out_file>\n" and die unless(@ARGV==2);open IN,"$ARGV[0]" or die "$!; can't open file $ARGV[0]\n";open OUT,">$ARGV[1]" or die "$!; can't open file $ARGV[1]\n";my%id;my%hash;while(<IN>){        chomp;        my@...

文章 R语言基础入门—函数

...泛使用的 R 函数。 # Create a sequence of numbers from 32 to 44.print(seq(32,44))# Find mean of numbers from 25 to 82.print(mean(25:82))# Find sum of numbers from 41 to 68.print(sum(41:68)) 当我们执行上面的代码,它产生以下结果 : [1] 32 33 34 ...

问题 请教老师,Gemma分析时的报错问题?

[root@iZ2e]# ./gemma-0.98.5-linux-static-AMD64 -bfile c -k output/result1.sXX.txt -lmm 1 -p p.txt -c c.txt GEMMA 0.98.5 (2021-08-25) by Xiang Zhou, Pjotr Prins and team (C) 2012-2021 ERROR: Enforce failed for file_kin output/result1.sXX.txt: open file in src/param.cpp at line 934 in CheckParam

文章 MA plot 绘图计算公式

...面) ############## Microarray data is often normalized within arrays to control for systematic biases in dye coupling and hybridization efficiencies, as well as other technical biases in the DNA probes and the print tip used to spot the array.[3] By minimizing these systematic variations, tru...

文章 食管癌 10X多样本(GSE145370)

...et -c "https://www.ncbi.nlm.nih.gov/geo/download/?acc=GSE145370&format=file" -O GSE145370_RAW.tar#解压tar xvf GSE145370_RAW.tar 准备map.txt 文件: #accession SampleID gender age pathological_stage tissue_type patient GSM4317409 S133A Male 60 IIIB T S133 GSM4317410 S134A Female 53 II...

文章 从tar压缩文件夹中仅提取某一个目标文件

...成不清楚的话,可以使用tar tvf指令进行查看: tar tvf file.tar 比如我们需要从中提取某个MD5文件,那么执行如下指令: tar xvf 1.tar ./MD5/file.md5

问题 用deseq_analysis.r分析差异表达分析时报如下的错误

[1] "abstract_factors is start" Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec,  :   line 2 did not have 3 elements Calls: abstract_factors -> read.table -> scan Execution halted

文章 遗传图绘制:

...utfile,ruler = F,lg.col = "lightblue1")## draw tickmarks at each cM from 0 to largest position of linkage groups to be drawnmaxpos <- floor(max(carrot$position))at.axis <- seq(0, maxpos)## put labels on ruler at every 10 cMaxlab <- vector()for (lab in 0:maxpos) {  if (!lab %% 10) {    ax...

文章 linux命令之cp

...果没有文件复制的权限,则系统会显示出错信息。将文件file复制到目录/usr/men/tmp下,并改名为file1cp file /usr/men/tmp/file1将目录/usr/men下的所有文件及其子目录复制到目录/usr/zh中cp -r /usr/men /usr/zh交互式地将目录/usr/men中的以m打头的...