Tutorial: create_analytic_cohort
Source:vignettes/create_analytic_cohort_vignette.Rmd
create_analytic_cohort_vignette.RmdIntroduction
The create_analytic_cohort() function allows the user to
create a cohort from the GENIE BPC data based on cancer diagnosis
information such as cancer cohort, treating institution, histology, and
stage at diagnosis, as well as cancer-directed regimen information
including regimen name and regimen order. This function returns a list
of data frames containing each of the clinical and genomic data files
for patients that met criteria for the analytic cohort.
This vignette will walk a user through the
create_analytic_cohort() function.
Setup
Before going through the tutorial, install and load the {genieBPC}
library and log into Synapse using the
set_synapse_credentials() function. For more information on
set_synapse_credentials(), refer to the
Tutorial: pull_data_synapse() vignette.
Modifying Function Arguments
The create_analytic_cohort() function includes several
input parameters:
| Argument | Description |
|---|---|
data_synapse |
The item from the nested list returned from pull_data_synapse() corresponding to the cancer cohort of interest. |
index_ca_seq |
Index cancer sequence. Default is 1, indicating the patient's first index cancer. The index cancer is also referred to as the BPC Project cancer in the GENIE BPC Analytic Data Guide; this is the cancer that met the eligibility criteria for the project and was selected at random for PRISSMM phenomic data curation. Specifying multiple index cancer sequences, e.g. index_ca_seq = c(1, 2), will return index cancers to patients with 1 index cancer and will return the first AND second index cancers to patients with multiple. |
institution |
GENIE BPC participating institution. Must be one of 'DFCI', 'MSK', 'UHN', or 'VICC' for NSCLC, BLADDER, Prostate, and PANC cohorts; must be one of 'DFCI', 'MSK', 'VICC' for CRC and BrCa. Default selection is all institutions. This parameter corresponds to the variable `institution` in the Analytic Data Guide. |
stage_dx |
Stage at diagnosis. Must be one of 'Stage I', 'Stage II', 'Stage III', 'Stage I-III NOS', 'Stage IV'. Default selection is all stages. This parameter corresponds to the variable `stage_dx` in the Analytic Data Guide. |
histology |
Cancer histology. For all cancer cohorts except for BrCa (breast cancer), this parameter corresponds to the variable `ca_hist_adeno_squamous` and must be one of 'Adenocarcinoma', 'Squamous cell', 'Sarcoma', 'Small cell carcinoma', 'Carcinoma', 'Other histologies/mixed tumor'. For BrCa, this parameter corresponds to the variable `ca_hist_brca` and must be one of 'Invasive lobular carcinoma', 'Invasive ductal carcinoma', 'Other histology'. Default selection is all histologies. |
regimen_drugs |
Vector with names of drugs in cancer-directed regimen, separated by a comma. For example, to specify a regimen consisting of Carboplatin and Pemetrexed, specify regimen_drugs = 'Carboplatin, Pemetrexed'. Acceptable values are found in the `drug_regimen_list` dataset provided with this package. This parameter corresponds to the variable `regimen_drugs` in the Analytic Data Guide. |
regimen_type |
Indicates whether the regimen(s) specified in `regimen_drugs` indicates the exact regimen to return, or if regimens containing the drugs listed in `regimen_drugs` should be returned. Must be one of 'Exact' or 'Containing'. The default is 'Exact'. |
regimen_order |
Order of cancer-directed regimen. If multiple drugs are specified, `regimen_order` indicates the regimen order for all drugs; different values of `regimen_order` cannot be specified for different drug regimen. If multiple values are specified, e.g. c(1, 2), then drug regimens that met either order criteria are returned. |
regimen_order_type |
Specifies whether the `regimen_order` parameter refers to the order of receipt of the drug regimen within the cancer diagnosis (across all other drug regimens; 'within cancer') or the order of receipt of the drug regimen within the times that that drug regimen was administered (e.g. the first time carboplatin pemetrexed was received, out of all times that the patient received carboplatin pemetrexed; 'within regimen'). Acceptable values are 'within cancer' and 'within regimen'. |
return_summary |
Specifies whether a summary table for the cohort is returned. Default is FALSE. The `gtsummary` package is required to return a summary table. |
This function returns a list of data frames containing clinical and next generation sequencing information for patients that met the specified criteria. Optionally, if return_summary = TRUE, the list also includes summary tables for the number of records per dataset (‘tbl_overall_summary’) as well as tables of key cancer diagnosis (‘tbl_cohort’), cancer-directed regimen (‘tbl_drugs’) and next generation sequencing (‘tbl_ngs’) variables.
Example using package test data
Create analytic cohort using package test data
create_analytic_cohort(data_synapse = genieBPC::nsclc_test_data)Examples using data from pull_data_synapse()
Examples 2-4 create analytic datasets based on the GENIE BPC NSCLC
cohort. To begin, the user must first pull the NSCLC data from Synapse
and save it to the local environment, as shown below using the
pull_data_synapse function:
nsclc_2_0 <- pull_data_synapse("NSCLC", version = "v2.0-public")Example 2
Create a cohort of all patients with stage IV NSCLC adenocarcinoma and also return all of their corresponding cancer-directed drugs
nsclc_stg_iv_adeno <- create_analytic_cohort(data_synapse = nsclc_2_0$NSCLC_v2.0,
stage_dx = "Stage IV",
histology = "Adenocarcinoma")
ls(nsclc_stg_iv_adeno)
#> [1] "cohort_ca_drugs" "cohort_ca_dx"
#> [3] "cohort_ca_dx_non_index" "cohort_cna"
#> [5] "cohort_fusions" "cohort_mutations_extended"
#> [7] "cohort_ngs" "cohort_prissmm_imaging"
#> [9] "cohort_prissmm_md" "cohort_prissmm_pathology"
#> [11] "cohort_pt_char"Example 3
Create a cohort of all NSCLC patients who received Cisplatin, Pemetrexed Disodium or Cisplatin, Etoposide as their first drug regimen for their first index NSCLC
create_analytic_cohort(data_synapse = nsclc_2_0$NSCLC_v2.0,
regimen_drugs = c("Cisplatin, Pemetrexed Disodium","Cisplatin, Etoposide"),
regimen_order = 1,
regimen_order_type = "within cancer")Example 4
Create a cohort of all NSCLC patients who received Cisplatin, Pemetrexed Disodium at any time throughout the course of treatment for their cancer diagnosis, but in the event that the patient received the drug multiple times, only select the first time.
nsclc_cisplatin_pem_any = create_analytic_cohort(data_synapse = nsclc_2_0$NSCLC_v2.0,
regimen_drugs = c("Cisplatin, Pemetrexed Disodium"),
regimen_order = 1,
regimen_order_type = "within regimen",
return_summary = TRUE)
nsclc_cisplatin_pem_any$tbl_overall_summary| Characteristic | N = 182 patients1 |
|---|---|
| Number of diagnoses per patient in cohort_ca_dx data frame | |
| 1 | 182 (100%) |
| Number of regimens per patient in cohort_ca_drugs data frame | |
| 1 | 182 (100%) |
| Number of CPTs per patient in cohort_ngs data frame | |
| 1 | 171 (94%) |
| 2 | 10 (5.5%) |
| 4 | 1 (0.5%) |
| 1 n (%) | |
nsclc_cisplatin_pem_any$tbl_cohort| Characteristic | N = 182 Diagnoses1 |
|---|---|
| Cohort (cohort) | |
| NSCLC | 182 (100%) |
| Institution (institution) | |
| DFCI | 71 (39%) |
| MSK | 92 (51%) |
| UHN | 7 (3.8%) |
| VICC | 12 (6.6%) |
| Stage at diagnosis (stage_dx) | |
| Stage I | 19 (10%) |
| Stage I-III NOS | 1 (0.5%) |
| Stage II | 42 (23%) |
| Stage III | 71 (39%) |
| Stage IV | 49 (27%) |
| Histology (ca_hist_adeno_squamous) | |
| Adenocarcinoma | 142 (89%) |
| Carcinoma | 14 (8.8%) |
| Other histologies/mixed tumor | 2 (1.3%) |
| Squamous cell | 1 (0.6%) |
| Unknown | 23 |
| 1 n (%) | |
nsclc_cisplatin_pem_any$tbl_drugs| Characteristic | N = 182 Regimens1 |
|---|---|
| Cohort (cohort) | |
| NSCLC | 182 (100%) |
| Institution (institution) | |
| DFCI | 71 (39%) |
| MSK | 92 (51%) |
| UHN | 7 (3.8%) |
| VICC | 12 (6.6%) |
| Drugs in regimen (regimen_drugs) | |
| Cisplatin, Pemetrexed Disodium | 182 (100%) |
| 1 n (%) | |
nsclc_cisplatin_pem_any$tbl_ngs| Characteristic | N = 195 Cancer Panel Tests1 |
|---|---|
| Cohort (cohort) | |
| NSCLC | 195 (100%) |
| Institution (institution) | |
| DFCI | 73 (37%) |
| MSK | 103 (53%) |
| UHN | 7 (3.6%) |
| VICC | 12 (6.2%) |
| OncoTree code (cpt_oncotree_code) | |
| LCLC | 1 (0.5%) |
| LUAD | 187 (96%) |
| LUPC | 1 (0.5%) |
| LUSC | 3 (1.5%) |
| NSCLCPD | 3 (1.5%) |
| Sequence assay ID (cpt_seq_assay_id) | |
| DFCI-ONCOPANEL-1 | 1 (0.5%) |
| DFCI-ONCOPANEL-2 | 40 (21%) |
| DFCI-ONCOPANEL-3 | 32 (16%) |
| MSK-IMPACT341 | 1 (0.5%) |
| MSK-IMPACT410 | 51 (26%) |
| MSK-IMPACT468 | 51 (26%) |
| UHN-48-V1 | 5 (2.6%) |
| UHN-50-V2 | 2 (1.0%) |
| VICC-01-SOLIDTUMOR | 4 (2.1%) |
| VICC-01-T7 | 8 (4.1%) |
| 1 n (%) | |