Skip to content

Model does not converge  #6

Open

Description

@innaaver

I tried twice, checked the input was correct but this is what I get:
image
Should I be running on all samples and genes or filter first?

Activity

  1. AlexanderAivazidis commented on Nov 22, 2022

    @AlexanderAivazidis
    Collaborator

    Hi! This usually happens when there are some unrealistically high expressing genes in the dataset. But from I have seen the optimization still converges for the remaining genes. It is just that the total or average loss is nan. So you can either filter very high expressing genes out before, or just leave it as it is.

  2. innaaver commented on Nov 23, 2022

    @innaaver
    Author

    Thank you for the reply! Is the correction calculated by the model different for every gene in one sample? I thought it just finds a per-sample correction

  3. AlexanderAivazidis commented on Nov 24, 2022

    @AlexanderAivazidis
    Collaborator

    yes it is different for every gene

  4. innaaver commented on Nov 24, 2022

    @innaaver
    Author

    Is there a way to see for which genes the model did not converge?

  5. AlexanderAivazidis commented on Nov 25, 2022

    @AlexanderAivazidis
    Collaborator

    No and I think the optimization never fully converges. However, the prior distributions and the starting point for the optimization are already sensible background binding values for all genes (specifically the mean of negative probes in each sample). So the values should be ok even after short optimisations.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions