Description Usage Arguments Details See Also Examples
A list of of GatingSet
objects. This class exists for method dispatching.
use GatingSetList
constructor to create a GatingSetList from a list of GatingSet
1 | GatingSetList(x, samples = NULL)
|
x |
a |
samples |
|
Objects store a collection of GatingSets,which usually has the same gating trees and markers. Most GatingSets methods can be applied to GatingSetList.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 | ## Not run:
#load several GatingSets from disk
gs_list<-lapply(list.files("../gs_toMerge",full=T) ,function(this_folder){
load_gs(this_folder)
})
#gs_list is a list
gs_groups <- merge(gs_list)
#returns a list of GatingSetList objects
gslist2 <- gs_groups[[2]]
#gslist2 is a GatingSetList that contains multiple GatingSets and they share the same gating and data structure
gslist2
class(gslist2)
sampleNames(gslist2)
#reference a GatingSet by numeric index
gslist2[[1]]
#reference a GatingSet by character index
gslist2[["30104.fcs"]]
#loop through all GatingSets within GatingSetList
lapply(gslist2,sampleNames)
#subset a GatingSetList by [
sampleNames(gslist2[c(4,1)])
sampleNames(gslist2[c(1,4)])
gslist2[c("30104.fcs")]
#get flow data from it
gs_pop_get_data(gslist2)
#get gated flow data from a particular popoulation
gs_pop_get_data(gslist2, "3+")
#extract the gates associated with one popoulation
gs_pop_get_gate(gslist2,"3+")
gs_pop_get_gate(gslist2,5)
#extract the pheno data
pData(gslist2[3:1])
#modify the pheno data
pd <- pData(gslist2)
pd$id <- 1:nrow(pd)
pData(gslist2) <- pd
pData(gslist2[3:2])
#plot the gate
autoplot(gslist2[1:2],5)
#remove cerntain gates by loop through GatingSets
gs_get_pop_paths(gslist2[[1]])
lapply(gslist2,function(gs)gs_pop_remove("Excl",gs = gs))
#extract the stats
gs_pop_get_count_fast(gslist2)
#extract statistics by using getQAStats defined in QUALIFIER package
res<-getQAStats(gslist2[c(4,2)],isMFI=F,isSpike=F,nslaves=1)
#archive the GatingSetList
save_gslist(gslist2, path ="~/rglab/workspace/flowIncubator/output/gslist",overwrite=T)
gslist2 <- load_gslist(path ="~/rglab/workspace/flowIncubator/output/gslist")
#convert GatingSetList into one GatingSet by merge_list_to_gs
gs_merged2 <- merge_list_to_gs(gslist2)
gs_merged2
## End(Not run)
## Not run:
samleNames(gsA) # return A1, A2
samleNames(gsB) # return B1, B2
gs.list <- list(gsA, gsB)
gslist<- GatingSetList(gs.list)
sampleNames(gslist) #return A1,A2,B1,B2
#set different order when create the GatingSetList
gslist<- GatingSetList(gs.list, samples = c("A1","B1", "A2", "B2"))
sampleNames(gslist) #return A1,B1,A2,B2
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.