|
...
|
...
|
@@ -29,13 +29,17 @@ add_meta <- function(ds){
|
|
29
|
29
|
}
|
|
30
|
30
|
fc <- lapply(cg,g=g, FUN=function(x,g) g %in% x)
|
|
31
|
31
|
names(fc) <- c("Mt","coding","ribosomal")
|
|
32
|
|
- ds <- addQCPerCell(ds, subsets=fc, percent_top=c(20,50,100,200))
|
|
|
32
|
+ if(exists("addQCPerCell")){
|
|
|
33
|
+ ds <- addQCPerCell(ds, subsets=fc, percent_top=c(20,50,100,200))
|
|
|
34
|
+ }else{
|
|
|
35
|
+ ds <- addPerCellQC(ds, subsets=fc, percent_top=c(20,50,100,200))
|
|
|
36
|
+ }
|
|
33
|
37
|
ds$total_features <- ds$detected
|
|
34
|
38
|
ds$log10_total_features <- log10(ds$detected)
|
|
35
|
39
|
ds$total_counts <- ds$sum
|
|
36
|
40
|
ds$log10_total_counts <- log10(ds$sum+1)
|
|
37
|
41
|
ds$featcount_ratio <- ds$log10_total_counts/ds$log10_total_features
|
|
38
|
|
- ds$featcount_dist <- .getFeatCountDist(ds)
|
|
|
42
|
+ ds$featcount_dist <- getFeatCountDist(ds)
|
|
39
|
43
|
ds$pct_counts_top_50_features <- ds$percent_top_50
|
|
40
|
44
|
for(f in names(fc))
|
|
41
|
45
|
ds[[paste0("pct_",f)]] <- ds[[paste0("subsets_",f,"_percent")]]
|
