perform anova on deconvoluted single cell output

[complexgenome]

Dear music developers,

I got deconvoluted output from MuSic R library. I am interested in performing differential gene expression of the deconvoluted result.

I have following outputs upon deconvolution: Est.prop.weighted, Est.prop.allgene, Weight.gene, r.squared.full, Var.prop

I am interested in two things:

1. I have per subject proportion within in each cell type. I have weight/count per gene within each person. However, I do not know the expression of a gene from a certain cell type in any individual. How can I extract that information from MuSic output?

2. I am interested in using anova on the music output. https://xuranw.github.io/MuSiC/reference/Anova_info.html

How do I proceed for both these tasks?

Thanks.

[complexgenome]

@xuranw can you assist?

[sam0per]

@sariya how did you solve it?