Return the concordance between two assays (i.e. single cell and hundred cell).
The "average" of singleCellRef
(after adjusting for the number of cells) and
singleCellComp
are taken per gene, per groups
.
A data.frame
with one row per gene-groups
is returned with some additional columns.
getConcordance(singleCellRef, singleCellcomp, groups = NULL, fun.natural = expavg, fun.cycle = logmean) getwss(concord, nexp) getss(concord) getrc(concord)
concord
concordance between two assays
getwss
: getrc the sum of squares, weighted by nexp
getss
: return the sum of squares
getrc
: Return Lin's (1989) concordance correlation coefficient
data(vbetaFA) sca1 <- subset(vbetaFA, ncells==1) sca100 <- subset(vbetaFA, ncells==100) concord <- getConcordance(sca1, sca100)#>#>getss(concord)#> [1] 3.621746getrc(concord)#> [1] 0.8374755