| Title: | gDRtestData - R data package with testing dose response data |
|---|---|
| Description: | R package with internal dose-response test data. Package provides functions to generate input testing data that can be used as the input for gDR pipeline. It also contains qs2 files with MAE data processed by gDR. |
| Authors: | Aleksander Chlebowski [aut], Bartosz Czech [aut] (ORCID: <https://orcid.org/0000-0002-9908-3007>), Arkadiusz Gladki [cre, aut] (ORCID: <https://orcid.org/0000-0002-7059-6378>), Marc Hafner [aut] (ORCID: <https://orcid.org/0000-0003-1337-7598>), Marcin Kamianowski [aut], Sergiu Mocanu [aut], Natalia Potocka [aut], Dariusz Scigocki [aut], Janina Smola [aut] (ORCID: <https://orcid.org/0009-0007-4347-7748>), Allison Voung [aut] |
| Maintainer: | Arkadiusz Gladki <[email protected]> |
| License: | Artistic-2.0 |
| Version: | 1.11.10 |
| Built: | 2026-10-10 10:43:52 UTC |
| Source: | https://github.com/bioc/gDRtestData |
Add concentrations
add_concentration(df_layout, concentrations = 10^(seq(-3, 1, 0.5)))add_concentration(df_layout, concentrations = 10^(seq(-3, 1, 0.5)))
df_layout |
data.table that should contains the cell line, drug, concentration, and replicate columns along with the annotations that needs to be propagated |
concentrations |
vector of numeric concentrations that will be added to df_layout |
data.table with concentrations
cell_lines <- create_synthetic_cell_lines() add_concentration(cell_lines)cell_lines <- create_synthetic_cell_lines() add_concentration(cell_lines)
Add data replicates
add_data_replicates(df_layout)add_data_replicates(df_layout)
df_layout |
data.table that should contains the cell line, drug, concentration, and replicate columns along with the annotations that needs to be propagated |
data.table with replicates
cell_lines <- create_synthetic_cell_lines() add_data_replicates(cell_lines)cell_lines <- create_synthetic_cell_lines() add_data_replicates(cell_lines)
Add data with day 0
add_day0_data(df_merged, noise_level = 0.05)add_day0_data(df_merged, noise_level = 0.05)
df_merged |
data.table with merged data |
noise_level |
numeric scalar with the level of noise added to the data |
data.table with day0 data
cell_lines <- create_synthetic_cell_lines() drugs <- create_synthetic_drugs() df_merged <- prepareData(cell_lines[seq_len(2), ], drugs[seq_len(4), ]) df_merged$Duration <- 72 df_merged$ReadoutValue <- 0 add_day0_data(df_merged)cell_lines <- create_synthetic_cell_lines() drugs <- create_synthetic_drugs() df_merged <- prepareData(cell_lines[seq_len(2), ], drugs[seq_len(4), ]) df_merged$Duration <- 72 df_merged$ReadoutValue <- 0 add_day0_data(df_merged)
Create data.table with synthetic cell lines
create_synthetic_cell_lines()create_synthetic_cell_lines()
data.table with synthetic cell lines
create_synthetic_cell_lines()create_synthetic_cell_lines()
Create data.table with synthetic drugs
create_synthetic_drugs()create_synthetic_drugs()
data.table with synthetic drugs
create_synthetic_drugs()create_synthetic_drugs()
Genome-wide CRISPR/Cas9 knockout dependency scores.
Matrix with cell lines as rows and genes as columns
data.table::fread(system.file("depmap_data/CRISPRGeneEffect.csv.gz",
package = "gDRtestData"))
Rows: Cell line identifiers
Columns: NCBI gene IDs (Entrez format)
Values: Dependency scores (-1 to +1); lower = more essential
NA indicates insufficient screen coverage
Downloaded May 26, 2026 from DepMap Portal (version 24Q4).
Citation: 24Q4 DepMap Release, including CRISPR Screens, PRISM Drug Screens, Copy Number, Mutation, Expression, and Fusions DepMap, Broad (2024). DepMap 24Q4 Public. Figshare+. Dataset. https://doi.org/10.25452/figshare.plus.27993248.v1
Cell line information from DepMap Public 24Q4.
Table
data.table::fread(system.file("depmap_data/Model.csv.gz", package = "gDRtestData"))
Rows: Cell line identifiers
Columns: Metadata columns (see below)
Values: Cell line annotations and donor demographics
Columns:
Unique cell line identifier (ACH-XXXXX format)
Cell line name from Cancer Cell Line Encyclopedia (CCLE)
Common cell line name (e.g., "A375", "HCT-116")
Abbreviated ID for model type. For cancer models, from Oncotree; for other disease types, generated by DepMap
Primary cancer lineage classification (from Oncotree for cancer models, generated by DepMap for others)
Primary disease of model (from Oncotree for cancer models, generated by DepMap for others)
Subtype classification of model (from Oncotree for cancer models, generated by DepMap for others)
Oncotree code for cancer models; left blank if no corresponding code exists
Site of tumor origin (Primary/Metastatic/Recurrence/Other/Unknown)
Age at time of tissue sampling
Age category at time of sampling (Adult/Pediatric/Fetus/Unknown)
Sex at time of sampling (Female/Male/Unknown)
Patient/clinical indicated race (not derived)
Downloaded May 26, 2026 from DepMap Portal (version 24Q4).
Citation: 24Q4 DepMap Release, including CRISPR Screens, PRISM Drug Screens, Copy Number, Mutation, Expression, and Fusions DepMap, Broad (2024). DepMap 24Q4 Public. Figshare+. Dataset. https://doi.org/10.25452/figshare.plus.27993248.v1
Gene-level copy number estimates from SNP microarray or WES.
Matrix with cell lines as rows and genes as columns
data.table::fread(system.file("depmap_data/OmicsCNGene.csv.gz", package = "gDRtestData"))
Rows: Cell line identifiers
Columns: NCBI gene IDs
Values: Log2 ratio relative to diploid reference
Typical range: -2 (deletion) to +3 (amplification)
0 = diploid (2 copies)
Downloaded May 26, 2026 from DepMap Portal (version 24Q4).
Citation: 24Q4 DepMap Release, including CRISPR Screens, PRISM Drug Screens, Copy Number, Mutation, Expression, and Fusions DepMap, Broad (2024). DepMap 24Q4 Public. Figshare+. Dataset. https://doi.org/10.25452/figshare.plus.27993248.v1
RNA-seq based gene expression for protein-coding genes.
Matrix with cell lines as rows and genes as columns
data.table::fread(
system.file("depmap_data/OmicsExpressionProteinCodingGenesTPMLogp1.csv.gz",
package = "gDRtestData"))
Rows: Cell line identifiers
Columns: NCBI gene IDs (Entrez format)
Values: Log2(TPM + 1) transformed expression
Only protein-coding genes included
Downloaded May 26, 2026 from DepMap Portal (version 24Q4).
Citation: 24Q4 DepMap Release, including CRISPR Screens, PRISM Drug Screens, Copy Number, Mutation, Expression, and Fusions DepMap, Broad (2024). DepMap 24Q4 Public. Figshare+. Dataset. https://doi.org/10.25452/figshare.plus.27993248.v1
Binary matrix of damaging mutations (frame-shift, stop-gain, splice-site).
Matrix with cell lines as rows and genes as columns
data.table::fread(system.file("depmap_data/OmicsSomaticMutationsMatrixDamaging.csv.gz",
package = "gDRtestData"))
Rows: Cell line identifiers
Columns: NCBI gene IDs
Values: Binary (0 = no mutation, 1 = damaging mutation)
High-confidence loss-of-function mutations
Downloaded May 26, 2026 from DepMap Portal (version 24Q4).
Citation: 24Q4 DepMap Release, including CRISPR Screens, PRISM Drug Screens, Copy Number, Mutation, Expression, and Fusions DepMap, Broad (2024). DepMap 24Q4 Public. Figshare+. Dataset. https://doi.org/10.25452/figshare.plus.27993248.v1
Binary matrix of known hotspot mutations from whole exome sequencing.
Matrix with cell lines as rows and genes as columns
data.table::fread(system.file("depmap_data/OmicsSomaticMutationsMatrixHotspot.csv.gz",
package = "gDRtestData"))
Rows: Cell line identifiers
Columns: NCBI gene IDs
Values: Binary (0 = no mutation, 1 = hotspot mutation)
Hotspots defined by COSMIC and OncoKB
Downloaded May 26, 2026 from DepMap Portal (version 24Q4).
Citation: 24Q4 DepMap Release, including CRISPR Screens, PRISM Drug Screens, Copy Number, Mutation, Expression, and Fusions DepMap, Broad (2024). DepMap 24Q4 Public. Figshare+. Dataset. https://doi.org/10.25452/figshare.plus.27993248.v1
Calculate E inf metric
generate_e_inf(drugs, cell_lines)generate_e_inf(drugs, cell_lines)
drugs |
data.table with drugs |
cell_lines |
data.table with cell lines |
matrix with random E inf
generate_e_inf(create_synthetic_drugs(), create_synthetic_cell_lines())generate_e_inf(create_synthetic_drugs(), create_synthetic_cell_lines())
Calculate EC50 metric
generate_ec50(drugs, cell_lines)generate_ec50(drugs, cell_lines)
drugs |
data.table with drugs |
cell_lines |
data.table with cell lines |
matrix with random EC50
generate_ec50(create_synthetic_drugs(), create_synthetic_cell_lines())generate_ec50(create_synthetic_drugs(), create_synthetic_cell_lines())
Generate hill coefficient
generate_hill_coef(drugs, cell_lines)generate_hill_coef(drugs, cell_lines)
drugs |
data.table with drugs |
cell_lines |
data.table with cell lines |
matrix with random hill coefficient
generate_hill_coef(create_synthetic_drugs(), create_synthetic_cell_lines())generate_hill_coef(create_synthetic_drugs(), create_synthetic_cell_lines())
Generate response data
generate_response_data(df_layout, noise_level = 0.1)generate_response_data(df_layout, noise_level = 0.1)
df_layout |
data.table that should contains the cell line, drug, concentration, and replicate columns along with the annotations that needs to be propagated |
noise_level |
numeric scalar with the level of noise added to the data |
data.table with response data
cell_lines <- create_synthetic_cell_lines() drugs <- create_synthetic_drugs() df_layout <- prepareData(cell_lines[seq_len(2), ], drugs[seq_len(4), ]) generate_response_data(df_layout)cell_lines <- create_synthetic_cell_lines() drugs <- create_synthetic_drugs() df_layout <- prepareData(cell_lines[seq_len(2), ], drugs[seq_len(4), ]) generate_response_data(df_layout)
Returns named vector of absolute paths to test datasets.
get_test_dataset_paths(datasets_dir = NULL, pattern = "finalMAE_")get_test_dataset_paths(datasets_dir = NULL, pattern = "finalMAE_")
datasets_dir |
path to directory with datasets (default |
pattern |
used to: (1) filter to qs2 files from the dataset_dir path and (2) prettify the labels of the files |
named vector of absolute paths
Kamil FoltyĆski [email protected]
get_test_dataset_paths() path <- system.file("testdata", package = "gDRtestData", mustWork = TRUE) get_test_dataset_paths(path)get_test_dataset_paths() path <- system.file("testdata", package = "gDRtestData", mustWork = TRUE) get_test_dataset_paths(path)
Create data.table with input co-dilution data containing noise for testing purposes
prepareCodilutionData( cell_lines, drugs, drugsIdx2 = 1, conc = 10^(seq(-3, 1, 0.5)), noise = 0.1 )prepareCodilutionData( cell_lines, drugs, drugsIdx2 = 1, conc = 10^(seq(-3, 1, 0.5)), noise = 0.1 )
cell_lines |
data.table with cell line info |
drugs |
data.table with drug info |
drugsIdx2 |
numeric vector of ids for secondary drug (in |
conc |
vector of doses |
noise |
number indicating level of noise |
data.table with input data for testing
prepareCodilutionData(create_synthetic_cell_lines()[seq_len(2), ], create_synthetic_drugs()[seq_len(4), ])prepareCodilutionData(create_synthetic_cell_lines()[seq_len(2), ], create_synthetic_drugs()[seq_len(4), ])
Create data.table with input combination data containing noise for testing purposes
prepareComboMergedData( cell_lines, drugs, drugsIdx1 = 2:4, drugsIdx2 = c(26, 26, 26), concentration = c(0, 0.2, 1), noise = 0.1, modifyDf2 = FALSE )prepareComboMergedData( cell_lines, drugs, drugsIdx1 = 2:4, drugsIdx2 = c(26, 26, 26), concentration = c(0, 0.2, 1), noise = 0.1, modifyDf2 = FALSE )
cell_lines |
data.table with cell line info |
drugs |
data.table with drug info |
drugsIdx1 |
numeric vector of ids for primary drug |
drugsIdx2 |
numeric vector of ids for secondary drug |
concentration |
numeric vector of doses |
noise |
number indicating level of noise |
modifyDf2 |
Boolean indicating if the table should me modified to keep reverse single agent data |
data.table with input data for testing
prepareComboMergedData(create_synthetic_cell_lines(), create_synthetic_drugs())prepareComboMergedData(create_synthetic_cell_lines(), create_synthetic_drugs())
Create data.table with input data for testing purposes
prepareData(cell_lines, drugs, conc = 10^(seq(-3, 1, 0.5)))prepareData(cell_lines, drugs, conc = 10^(seq(-3, 1, 0.5)))
cell_lines |
data.table with cell line info |
drugs |
data.table with drug info |
conc |
vector of doses |
data.table with input data for testing
prepareData(create_synthetic_cell_lines(), create_synthetic_drugs())prepareData(create_synthetic_cell_lines(), create_synthetic_drugs())
Create data.table with input data containing noise for testing purposes
prepareMergedData(cell_lines, drugs, noise = 0.1)prepareMergedData(cell_lines, drugs, noise = 0.1)
cell_lines |
data.table with cell line info |
drugs |
data.table with drug info |
noise |
number indicating level of noise |
data.table with input data for testing
prepareMergedData(create_synthetic_cell_lines(), create_synthetic_drugs())prepareMergedData(create_synthetic_cell_lines(), create_synthetic_drugs())
The single-agent screening data from the PRISM platform - with one drug and 774 cell lines.
A MultiAssayExperiment object
containing a single-agent experiment with 5 assays.
gDRutils::get_synthetic_data("prism")
Experimental data matrices tracking cell line viability across processing stages:
Raw intensity values for treated cell line pools.
Raw intensity values for DMSO and killing controls.
Viability scores relative to controls.
Replicate-averaged viability scores.
Dose-response curve parameters (e.g., IC50, AUC).
Feature metadata describing the cell lines (models) across assays:
Unique Broad Institute cell line identifier.
Publicly recognized cancer cell line name.
Primary tissue of origin.
Identifier for the parental cell line.
Histological or molecular subtype.
Estimated doubling time in hours.
Feature metadata describing the treatment compounds across assays:
Genentech compound identifier.
Public or commercial name of the tested compound.
Mechanism of action / biological target of the drug.
Exposure time of the cell lines to the compound (e.g., 72h).
Hagenbeek, T.J., Zbieg, J., Smith, R. et al. Covalent pan-TEAD inhibitors block YAP activity and demonstrate brain penetrance in a Hippo-dependent cancer model. Nat Commun (2026). https://doi.org/10.1038/s41467-026-74722-5