Skip to content

eBayes error when running propeller.anova on serial data #12

@bapoorva

Description

@bapoorva

Hi,

Thanks for the wonderful package. I have single cell data on two groups (control and treatment) across 3 timepoints (T0, T1, T2) and I'm interested in seeing if there are proportional differences between the two groups across time ? I've been looking at your vignette with the PBMC data which looks into the effect of age and sex and I was thinking if the following design would be the appropriate way to do it

alltp= readRDS("alltp_integrated.RDS")
meta= alltp@meta.data %>% select(orig.ident,celltypes,Group,Timepoint)
prop_list <- getTransformedProps(
  clusters = alltp$celltypes,
  sample = alltp$orig.ident
)

sample_info <- meta %>%distinct(orig.ident, Group, Timepoint)
design_matrix <- model.matrix(~ 0+Group*Timepoint,data=sample_info)

res=propeller.anova(prop.list = prop_list,design = design_matrix,coef = 1:2) 

I keep getting this error

Error in eBayes(fit[, coef[-1]], robust = robust, trend = trend) : 
  promise already under evaluation: recursive default argument reference or earlier problems?

Any idea on how to fix this issue or if there is a better way to do this ?

Thanks so much

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions