library(readr)library(dplyr)library(ggplot2)library(tidyverse)library(DESeq2)mycounts<-read.csv("C:/Users/~/Desktop/datacounts.csv")metadata<-read.csv("C:/Users/~/Desktop/datametadata.csv")countData<-mycounts[!(duplicated(mycounts[[1]])|duplicated(mycounts[[1]],fromLast=TRUE)),]rownames(mycounts)=make.names(mycounts$gene,unique=TRUE)mycounts<-as.data.frame(mycounts)...[阅读全文]
longteng
随笔记录, 工作日记.
博文
(2024-03-23 15:45:30)
library(readr)library(dplyr)library(ggplot2)library(tidyverse)library(DESeq2)mycounts<-read.csv("C:/Users/~/Desktop/datacounts.csv")metadata<-read.csv("C:/Users/~/Desktop/datametadata.csv")countData<-mycounts[!(duplicated(mycounts[[1]])|duplicated(mycounts[[1]],fromLast=TRUE)),]rownames(mycounts)=make.names(mycounts$gene,unique=TRUE)mycounts<-as.data.frame(mycounts)...[阅读全文]
(2024-03-17 13:12:15)
library(dplyr)library(Seurat)library(patchwork)library(ggplot2)dat1<-Read10X("C:/Users/~/Desktop/data")#data(matrix.mtx.gz,barcode.tsv.gz,features.tsv.gz)dat1<-CreateSeuratObject(counts=dat1,project="NT",min.cells=3,min.features=200)dat1dat1[["percent.mt"]]<-PercentageFeatureSet(dat1,pattern="^MT-")VlnPlot(dat1,features=c("nFeature_RNA",&q...[阅读全文]
(2024-03-02 08:43:31)
(2024-02-17 08:02:23)
(2024-02-17 08:01:48)
(2024-02-17 08:00:54)
(2024-02-17 07:59:24)
(2024-02-09 08:41:56)
library(Seurat)library(tibble)library(magrittr)library(patchwork)#Readindatadata<-data.table::fread("~/Downloads/~.txt.gz",data.table=F)#Genenamesareinthefirstcolumnsoweneedtomovethemtorownamesdata<-data%>%column_to_rownames("GENE")#CreateSeuratNT1<-CreateSeuratObject(counts=data,min.cells=3,min.features=350,project="NT1")data2<-data.table::fread("...[阅读全文]
