Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Adding Figure 4 - Leave One Image Out Analysis #40

Merged
merged 7 commits into from
Oct 26, 2023

Conversation

gwaybio
Copy link
Member

@gwaybio gwaybio commented Oct 3, 2023

main_figure_4_loio

Note: This figure will change. I am retraining the models (performing a new grid search for each individual image left out) and adding shuffled models. The shuffled models will likely belong in a supplement, not the main figure

@review-notebook-app
Copy link

Check out this pull request on  ReviewNB

See visual diffs & provide feedback on Jupyter Notebooks.


Powered by ReviewNB

Copy link
Member

@jenna-tomkinson jenna-tomkinson left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Very interesting figure! I mainly added comments and questions to the figure itself and not the R code as it might be easier for me to describe what I see and not what I am reading haha.

I am excited to get to hear back on my questions!

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I have a few questions and comments here:

  1. For panel A, what is "rank of prediction"? Does 1 mean the best prediction group or something? Also, why is rank 1 have high prediction probability for both correctly predicted phenotype and not? Wouldn't we expect that yes(correct prediction) should be the highest all the across?
  2. For panel B, what is a good output here? Also, what is the blue dotted line represent? Is each point a feature from each data feature space? Do we want to see more points on the left side of the blue line?
  3. For panel C, this seems the most clear that with the LOIO that there are more correct prediction than not but I don't know what "strict threshold" means. Can you explain what that is?
  4. For panel D, the x-axis title looks wrong to me... If the legend is yes or no if the phenotypic class was predicted correctly, then why would the x-axis represent "cell proportions of incorrect phenotype predictions". This seems contradictory to me.

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Sorry for my delay in answering these questions. Thanks for asking them. I will do my best to address below. We will definitely need to clarify these in the future paper.

  1. For each cell, we receive 15 phenotype predictions (p values per phenotype). The x axis is simply this prediction, rank-ordered by probability. The color means if the prediction was correct or not. For example, a cell could be labeled by apoptosis and we predict interphase with 0.8, apoptosis with 0.15 and all the rest add up to 0.5. This cell is in the plot 15 different times, yellow in each case except it is blue in x-axis rank 2.
  2. Yep, correct. A good output is everything to the left of the dotted line in the Shuffled: FALSE models.
  3. I added a comment describing what "high_threshold" was in the plot. This should help with interpretation, but we will definitely also add clarity in the manuscript!
  4. Ah! Thanks for bringing this up. The important distinction here is "Categories". These are the higher-order categories that mitocheck assigned (so groups of individual phenotypes). Your comment here made me change the legend to be more clear! Thanks


output_fig_loio <- file.path("figures", "main_figure_4_loio.png")

# Set constants
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I recommend adding more to this comment, since there is only one constant which is the threshold. What does this threshold represent and why that number.

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Good point, will do in next commits!

@gwaybio
Copy link
Member Author

gwaybio commented Oct 19, 2023

The new LOIO should be done today - I'll address your comments and update the visualization later! (sorry for the delay!)

@gwaybio
Copy link
Member Author

gwaybio commented Oct 26, 2023

New LOIO analysis completed! I will get this merged in and we can discuss at our meeting today

@gwaybio gwaybio merged commit 0c085a3 into WayScience:main Oct 26, 2023
@gwaybio gwaybio deleted the fig4 branch October 26, 2023 19:18
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants