IntegrateData error: number of items to replace is not a multiple of replacement length

[xinmiaoyan]

I used below code to integrate data by rPCA, but error occurred, I've tried many ways but still couldn't solve this problem. appreciate it if you could help.

##### Perform integration #######
sc.anchors <- FindIntegrationAnchors(object.list = sc.list, anchor.features = features, reduction = "rpca")
# this command creates an 'integrated' data assay
sc.combined <- IntegrateData(anchorset = sc.anchors)

Merging dataset 2 into 4 Extracting anchors for merged samples Finding integration vectors Finding integration vector weights Error in idx[i, ] <- res[[i]][[1]] : number of items to replace is not a multiple of replacement length

This error couldn't solved by adjust parameters, such has l2.norm and k.filter.... Appreciate it if you could help.

[Gesmira]

Hi, what are the sizes of your datasets in sc.list? k.weight (the number of neighbors to consider when weighting anchors) in IntegrateData() defaults to 100 so if one of your datasets is smaller than 100 it will cause this error to occur. Have you tried adjusting the k.weight parameter to the size of your smallest dataset? Otherwise, you may want to remove very small datasets.

[xinmiaoyan]

@Gesmira, yes, I tried to adjust k.weight. the smallest sc in my sc.list has 130 cells, when I adjusted k.weight to 50, this error was still here and didn't be resolved.

[Gesmira]

Have you tried also setting k.filter inFindIntegrationAnchors() before you run IntegrateData()? For example:

sc.anchors <- FindIntegrationAnchors(object.list = sc.list, anchor.features = features, reduction = "rpca", k.filter = 100)
sc.combined <- IntegrateData(anchorset = sc.anchors, k.weight = 100)

[xinmiaoyan]

yes, I tried k.filter = 30, it didn't work for this error

[Gesmira]

Just to see if the issue is due to the small dataset, are you able to run the integration through without it?

[xinmiaoyan]

I removed the smallest 130 sc, and it works this time. and the smallest one is 205, I'm wondering what's the cutoff for the smallest obj that can be used to integrate? thanks

[xinmiaoyan]

yes, I did checked the size of sc.list, the smallest one is 130, which is higher than 100. And I also adjusted k.weight to 50, also the same error.

On Fri, Aug 26, 2022 at 2:08 PM gesmira @.***> wrote:

Hi, what are the sizes of your datasets in sc.list? k.weight (the number of neighbors to consider when weighting anchors) in IntegrateData() defaults to 100 so if one of your datasets is smaller than 100 it will cause this error to occur. Have you tried adjusting the k.weight parameter to the size of your smallest dataset? Otherwise, you may want to remove very small datasets.

— Reply to this email directly, view it on GitHub https://github.com/satijalab/seurat/issues/6341#issuecomment-1228826037, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATZU2PDNHFI6M625ZHZHSWLV3EI3PANCNFSM57QEIF4Q . You are receiving this because you authored the thread.Message ID: @.***>

[decodebiology]

I got the lowest cell number "3770" for a smaller object. It is same error. tried playing with different parameters k.filter, k.anchor, k.score, no luck so far.