Commit c33a42c4 authored by Robert Phillips's avatar Robert Phillips
Browse files

Update comment

parent b9e59a5e
......@@ -360,7 +360,7 @@ for(i in names(Cluster_Lists)){
#Create a data.frame where the rownames are the cells within the clusters
group.info <- data.frame(row.names = c(WhichCells(object = VTA_subset,idents = i),
WhichCells(object = VTA_subset,idents = as.character(unique(Idents(VTA_subset))[which(unique(Idents(VTA_subset))!=i)]))))
#Create a group column where cocaine is Group1 and saline is group2. The way to do that is to search the rows names for Cocaine and Saline cell identities
#Create a group column where cell type of interest is Group1 and all other cells are group2.
group.info[WhichCells(object = VTA_subset,idents = i), "group"] <- "Group1"
group.info[WhichCells(object = VTA_subset,idents = as.character(unique(Idents(VTA_subset))[which(unique(Idents(VTA_subset))!=i)])), "group"] <- "Group2"
#Make the group column a factor so it can be tested with the wilcox
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment