-
Notifications
You must be signed in to change notification settings - Fork 0
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
Add ability to select columns for the last sample replicates #39
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
A couple questions here but I think we're moving along!! 🎉
#' @description A short description... | ||
#' @description This bar graph first uses the specified `filter_zerocount_target_col` columns to flag pgRNA constructs that have a raw count of 0 in any one of those columns/samples of interest. | ||
#' Then, it looks at the specified columns for the final day/sample replicates (`filter_replicates_target_col`) to see for pgRNAs that were flagged by the filter, how many of those replicate samples had raw counts of zeros. And it produces a bar plot reporting on this. | ||
#' Note, if you select samples/columns to check with the filter that don't have the replicate samples, this graph won't be informative. So you want there to be overlap between the columns for the two target_col parameters to have an informative graph |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Great notes. This is excellent documentation!
R/plots-qc.R
Outdated
@@ -113,18 +120,32 @@ qc_variance_hist <- function(gimap_dataset, wide_ar = 0.75){ | |||
#' @examples \dontrun{ | |||
#' gimap_dataset <- get_example_data("gimap") | |||
#' qc_constructs_countzero_bar(gimap_dataset) | |||
#' | |||
#' #or if you want to select a specific column(s) for looking at where/which samples zero counts are present for | |||
#' qc_constructs_countzero_bar(gimap_dataset, filter_zerocount_target_col = c(3:5)) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
You shouldn't need to have the c()
here?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Yeah I think you're right. I just hadn't gone through and removed them from all the PRs following the code review last week.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
pushed a change that simplifies the ones I found
Co-authored-by: Candace Savonen <cansav09@gmail.com>
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
LGTM!
This parameter for selecting columns doesn't directly impact a filter, but does impact two graphs. So there are only changes to the
R/plots-qc.R
file and theinst/rmd/gimapQCTemplate.Rmd
file. Specifically impacting theqc_variance_hist()
andqc_constructs_countzero_bar()
functions. The changes to these functions include adding to the headers (param and description information), replacing anywhere numbers are hardcoded with the parameter, and sets the parameters to a default if they're null. Within the QC Template Rmd, changes include referencing the parameters in the function calls.I can confirm that I knit the
getting-started.Rmd
successfully locally. The only oddity is that the heatmap is still showing up in the vignette rather than the QC report (as mentioned in PR #33)Also need to setup something to check the identity of the target columns if they are not null as described in PR #38