Input format for `combinations` in `hashedDrops`
1
0
Entering edit mode
@mikhaelmanurung-17423
Last seen 2.5 years ago
Netherlands

Has anyone tried using the combinations argument in hashedDrops? Unfortunately, I cannot fully understand the explanation in the help docs:(

Thanks in advance!

MIkhael

DropletUtils • 833 views
ADD COMMENT
1
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 46 minutes ago
The city by the bay

Let's have a look at what ?DropletUtils says:

An integer matrix specifying valid _combinations_ of HTOs. Each row corresponds to a single sample and specifies the indices of rows in ‘x’ corresponding to the HTOs used to label that sample.

So let's say we have a matrix like:

x <- matrix(rpois(1000, 10), ncol=100) # made up
rownames(x) <- LETTERS[1:10]
x[,1:5]
##   [,1] [,2] [,3] [,4] [,5]
## A    9    9   16    8    9
## B   12   15    7   12    9
## C   10    7   14   13   16
## D   13    5   13    7    5
## E    8    7   13   10    6
## F   10   10   12    7    7
## G   12   15    9    6    8
## H    9   13   12   11   11
## I    5   10   11    6    4
## J    9    7   12   11    6

I might then have this combinations:

combinations <- rbind(
    Sample1 = c(1, 5), # Sample 1 is marked by A and E
    Sample2 = c(1, 4), # Sample 2 is marked by A and D
    Sample3 = c(2, 3)  # Sample 3 is marked by B and C
)

and so on. This particular example uses combinations of 2, but you could also imagine higher-order combinations if your experimental design was more ambitious. I do recall there being some constraints though - something like the number of HTOs used in each sample must be less than half the total number of HTOs.

ADD COMMENT

Login before adding your answer.

Traffic: 405 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