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)...[
阅读全文]
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...[
阅读全文]
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("...[
阅读全文]