Some questions about mixed model

[HelloWorldLTY]

Hi, I intend to use mixed model method to identify marker gene across disease/control + male/female group.

My code is like: `library(Seurat) library(fgsea) source('sex_specificity_helper.R') library(stats) library(Libra)

data = readRDS("file.rds")

data$replicate = colnames(data)

data$cell_type = data$gender data$label = data$condition

DE = run_de(data, de_family = 'mixedmodel', de_method = "negbinom", n_threads = 12)

`

My gender is like "Female", "Male". My condition is like "disesae_female". It is single-cell dataset. But my running time and memory usage are quite large. Is it normal? More than 50 GB and 5 hours. Thanks.

[jordansquair]

Why are you setting replicate to the barcode? Replicate should be the animal/sample. Likely that is the issue here.

[HelloWorldLTY]

Hi, I think referring from the paper, it is: So I treate replicates are barcodes. Thanks for your remind.