Skip to content
This repository has been archived by the owner on Jun 21, 2023. It is now read-only.

V22 subtyping - chr instab (2/N) #1434

Merged
merged 46 commits into from
Jun 6, 2022
Merged

V22 subtyping - chr instab (2/N) #1434

merged 46 commits into from
Jun 6, 2022

Conversation

jharenza
Copy link
Collaborator

@jharenza jharenza commented May 26, 2022

Purpose/implementation Section

What scientific question is your analysis addressing?

Run chromosomal instability

What was your approach?

This module also needed to be divided into script to run for subtyping and plot scripts to not run. I added that code to the bash script.

01-localization-of-breakpoints.Rmd was not set up to intake either histologies file, so I added that. However, when I tried to run the OPENPBTA_BASE_SUBTYPING=1 bash ../analyses/chromosomal-instability/run_breakpoint_analysis.sh, I get an error:

Quitting from lines 150-154 (01-localization-of-breakpoints.Rmd) 
Error: '../../data/pbta-histologies.tsv' does not exist in current working directory ('/home/rstudio/kitematic/analyses/chromosomal-instability').

So I am not sure where the logic is going wrong.

In the meantime, I changed the param to 1 inside of the Rmd and ran it, to get the analysis file, in order to continue the PRs.

What GitHub issue does your pull request address?

#1426 partial

Directions for reviewers. Tell potential reviewers what kind of feedback you are soliciting.

Which areas should receive a particularly close look?

logic for choosing histology file

Is there anything that you want to discuss further?

Also note, this 01 script requires short_histology

Is the analysis in a mature enough form that the resulting figure(s) and/or table(s) are ready for review?

Results

What types of results are included (e.g., table, figure)?

What is your summary of the results?

Reproducibility Checklist

  • The dependencies required to run the code in this pull request have been added to the project Dockerfile.
  • This analysis has been added to continuous integration.

Documentation Checklist

  • This analysis module has a README and it is up to date.
  • This analysis is recorded in the table in analyses/README.md and the entry is up to date.
  • The analytical code is documented and contains comments.

@jharenza jharenza changed the base branch from master to v22-subtyping May 26, 2022 03:11
@jharenza jharenza requested a review from jaclyn-taroni May 26, 2022 03:12
breaks_file <- file.path(output_dir, "breaks_lists.RDS")
uncallable_file <- file.path(
"..", "copy_number_consensus_call", "ref",
"cnv_excluded_regions.bed"
)

# set up histologies file
if ( params$base_run == 0 ){
Copy link
Member

Choose a reason for hiding this comment

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

Can you confirm that this runs as expected if you pass 0 as the param? I can think of other cases where params did not exactly behave as expected.

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Do you mean within a v21 branch or within v22 branch which should error out?

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

while in this branch (v22 data download), running:

root@334c2bfa1312:/home/rstudio/kitematic/scripts# OPENPBTA_BASE_SUBTYPING=0 bash ../analyses/chromosomal-instability/run_breakpoint_analysis.sh

Error: 
 Could not find needed file(s):
../../data/pbta-histologies.tsv
Check your options and set up.
Execution halted

This is expected. I will open a new PR with a v21 run from master.

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Seems the param works if 1 is default, but not 0.

Copy link
Member

Choose a reason for hiding this comment

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

I meant did you do functional testing by passing params = list(base_run = 0). I just did locally.

Copy link
Member

@jaclyn-taroni jaclyn-taroni left a comment

Choose a reason for hiding this comment

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

👍🏻 tested the logic locally

@jharenza jharenza merged commit d51f501 into v22-subtyping Jun 6, 2022
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants