drawing Venn graph
0
0
Entering edit mode
2.0 years ago
Andy ▴ 120

Good morning everyone,

I am intend to draw a Venn graph. But I have a huge trouble to find the elements to draw it. In my dataset, there are two conditions of the samples. I want to check the different expression of gene markers in these two condition patients, but I failed to find the elements.

I hope you guys could give me some suggestions on how to find the same elements in my dataset.

Best

Andy

Venn R • 1.1k views
ADD COMMENT
0
Entering edit mode

Well, what did you try?

ADD REPLY
0
Entering edit mode

I tried to use findallmarker() function, but I got an empty .csv file.

>fibroblast_class<-subset_data
 DefaultAssay(fibroblast_class)<-"RNA" 
 fibroblast_markers<-FindAllMarkers(fibroblast_class, only.pos = T, test.use = "MAST", min.pct = 0.25, logfc.threshold = 0.25, ident=c("organ")) #need to install MAST package 
write.csv(fibroblast_markers, "fibrb.markers_0.25_0.251116.csv") 
top10 <- fibroblast_markers %>% group_by(cluster) %>% top_n(n = 10, wt = avg_log2FC) 
fibroblast_class<-ScaleData(fibroblast_class) 
``

I also tried the way like we find the sub cluster, but I did not know why it says cannot find the idents, but I use table(data@active.ident) to check, and my ident is there.

ADD REPLY
0
Entering edit mode
table(subset_data@active.ident)

pancreatic     kidney      liver       lung      skine 
      1017        461        137       1116       9583 

fibroblast_venn<-subset(subset_data, idents="organ")
Error in WhichCells.Seurat(object = x, cells = cells, idents = idents,  : 
  Cannot find the following identities in the object: organ
`

This is the code I used for find sub cluster.

ADD REPLY
1
Entering edit mode

Well, you don't seem to have a category called "organ", so what exactly are you expecting that subset command to do?

ADD REPLY
0
Entering edit mode

Since I am thinking my organ ident include pancreatic, kidney, liver, lung and skine, and I should set idents="organs".

I think you meant I should set idents="pancreatic"? like this to find the subset?

ADD REPLY
0
Entering edit mode

I got your point!

I got the data finally!

😭

ADD REPLY

Login before adding your answer.

Traffic: 2140 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6