inst/scripts/make-data_20Q1.R
In UCLouvain-CBIO/depmap: Cancer Dependency Map Data Package
### =========================================================================
### Make 20Q1 depmap data
### -------------------------------------------------------------------------
###
## This scripts documents how to download and generate the data files.
## Note 1: This scipt assumes (!!!) that it is run in ./depmap/inst/scripts/ and
## saves the resulting .rda files in ./depmap/inst/data
## (e.g. *setwd(./depmap/inst/scripts/)*)
## Note 2: the Broad Institute may change the download urls for these datasets.
## If the link to one of these datasets is broken, please contact the package
## maintainer. All Depmap data can be downloaded from the API at
## https://depmap.org/portal/download/ even if the specific links may change
library("readr")
library("dplyr")
library("tidyr")
library("ExperimentHub")
##########################################
## EH3294 depmap `metadata_20Q1` dataset
##########################################
## data cleaning of `metadata` dataset
### loading data (downloading .csv file from online source)
url_28 <- "https://ndownloader.figshare.com/files/21522000"
sample_info <- read_csv(url_28)
### DepMap-2020q1-celllines.csv data renamed to `metadata`
metadata_20Q1 <- sample_info
# names(metadata_20Q1)
# [1] "DepMap_ID" "stripped_cell_line_name" "CCLE_Name"
# [4] "alias" "COSMIC_ID" "lineage"
# [7] "lineage_subtype" "lineage_sub_subtype" "lineage_molecular_subtype"
# [10] "sex" "source" "Achilles_n_replicates"
# [13] "cell_line_NNMD" "culture_type" "culture_medium"
# [16] "cas9_activity" "RRID" "sample_collection_site"
# [19] "primary_or_metastasis" "disease" "disease_subtype"
# [22] "age" "Sanger_model_ID" "additional_info"
#### Rename `metadata` columns to contain underscores and be in snake case
## note: "metadata_20Q1" has different columns than "metadata_20Q1"
names(metadata_20Q1)[1:24] <-c("depmap_id", "stripped_cell_line_name",
"cell_line", "aliases", "cosmic_id", "lineage",
"lineage_subtype", "lineage_sub_subtype",
"lineage_molecular_subtype", "sex",
"source", "Achilles_n_replicates",
"cell_line_NNMD", "culture_type",
"culture_medium", "cas9_activity", "RRID",
"sample_collection_site", "primary_or_metastasis",
"primary_disease", "subtype_disease", "age",
"sanger_id", "additional_info")
### visual check
View(metadata_20Q1)
### saving cleaned and converted `metadata` data as .rda file
save(metadata_20Q1, file = "../eh_data/metadata_20Q1.rda",
compress = "xz", compression_level = 9)
## access the data on ExperimentHub
# hub <- ExperimentHub()
# x <- query(hub, "depmap")
##########################################
## `depmap_id_to_name_20Q1 map `depmap_id` to `cell_line`
##########################################
## generation of `metadata` subset `depmap_id_to_name_20Q1`
## The subset of `metadata`, `depmap_id_to_name_Q2` is used to map `cell_line`
## and depmap_id via left_join in other depmap datasets that do not contain
## both variables. If you are generating the depmap data from scratch, you will
## need to run the following code to generate the data correctly.
### `depmap_id_to_name` to add `depmap_id` or `cell_line` to other datasets
depmap_id_to_name_20Q1 <- metadata_20Q1 %>% dplyr::select(depmap_id, cell_line)
### visual check
# head(depmap_id_to_name_20Q1)
##########################################
## EH3293 depmap `mutationCalls_20Q1` dataset
##########################################
## data cleaning of `mutationCalls` dataset
### loading data (downloading .csv file from online source)
url_29 <- "https://ndownloader.figshare.com/files/21521967"
CCLE_mutations_20Q1 <- read_csv(url_29)
### depmap_20Q1_mutation_calls data renamed to `mutationCalls`
mutationCalls_20Q1 <- CCLE_mutations_20Q1
# head(mutationCalls_20Q1)
## remove unnecessary X1 variable
mutationCalls_20Q1 <- mutationCalls_20Q1 %>% dplyr::select(-X1)
#### rename last column to depmap_id
# names(mutationCalls_20Q1)
# [1] "Hugo_Symbol" "Entrez_Gene_Id" "NCBI_Build"
# [4] "Chromosome" "Start_position" "End_position"
# [7] "Strand" "Variant_Classification" "Variant_Type"
# [10] "Reference_Allele" "Tumor_Seq_Allele1" "dbSNP_RS"
# [13] "dbSNP_Val_Status" "Genome_Change" "Annotation_Transcript"
# [16] "Tumor_Sample_Barcode" "cDNA_Change" "Codon_Change"
# [19] "Protein_Change" "isDeleterious" "isTCGAhotspot"
# [22] "TCGAhsCnt" "isCOSMIChotspot" "COSMIChsCnt"
# [25] "ExAC_AF" "CGA_WES_AC" "SangerWES_AC"
# [28] "SangerRecalibWES_AC" "RNAseq_AC" "HC_AC"
# [31] "RD_AC" "WGS_AC" "Variant_annotation"
# [34] "DepMap_ID"
## note: "mutationCalls_20Q1" has different columns than "mutationCalls_20Q1"
## the variable "VA_WES_AC" is no longer present in this dataset, unlike
## previous releases (e.g. 19Q1)!
names(mutationCalls_20Q1)[1:34] <- c("gene_name", "entrez_id", "ncbi_build",
"chromosome", "start_pos", "end_pos", "strand",
"var_class","var_type", "ref_allele",
"tumor_seq_allele1", "dbSNP_RS",
"dbSNP_val_status", "genome_change",
"annotation_transcript", "tumor_sample_barcode",
"cDNA_change", "codon_change", "protein_change",
"is_deleterious", "is_tcga_hotspot",
"tcga_hsCnt", "is_cosmic_hotspot",
"cosmic_hsCnt", "ExAC_AF", "CGA_WES_AC",
"sanger_WES_AC", "sanger_recalib_WES_AC",
"RNAseq_AC", "HC_AC", "RD_AC", "WGS_AC",
"var_annotation","depmap_id")
### rearrange columns into same column format as other datasets
mutationCalls_20Q1 <- mutationCalls_20Q1 %>%
dplyr::select(depmap_id, gene_name, entrez_id, ncbi_build, chromosome,
start_pos, end_pos, strand, var_class, var_type, ref_allele,
tumor_seq_allele1, dbSNP_RS, dbSNP_val_status, genome_change,
annotation_transcript, tumor_sample_barcode, cDNA_change,
codon_change, protein_change, is_deleterious, is_tcga_hotspot,
tcga_hsCnt, is_cosmic_hotspot, cosmic_hsCnt, ExAC_AF,
CGA_WES_AC, sanger_WES_AC, sanger_recalib_WES_AC, RNAseq_AC,
HC_AC, RD_AC, WGS_AC, var_annotation)
### visual check
# head(mutationCalls_20Q1)
### saving cleaned and converted `mutationCalls` data as .rda file
save(mutationCalls_20Q1, file = "../eh_data/mutationCalls_20Q1.rda",
compress = "xz", compression_level = 9)
## access the data on ExperimentHub
# hub <- ExperimentHub()
# x <- query(hub, "depmap")
##########################################
## EH3291 depmap `copyNumber_20Q1` dataset
##########################################
## data cleaning of `copyNumber` dataset
### loading data (downloading .csv file from online source)
url_30 <- "https://ndownloader.figshare.com/files/21521964"
CCLE_gene_cn_20Q1 <- read_csv(url_30)
### public_20Q1_gene_cn.csv data renamed to `copyNumber`
copyNumber_20Q1 <- CCLE_gene_cn_20Q1
### rename column first column to "depmap_id"
names(copyNumber_20Q1)[1] <- "depmap_id"
### gather into long form on columns: `depmap_id`, `gene`, `logCopyNumber`
copyNumber_20Q1_long <- gather(copyNumber_20Q1, gene, log_copy_number,
-depmap_id)
### mutate gene column into `gene_name` and `entrez_id`
copyNumber_20Q1_long <- copyNumber_20Q1_long %>%
mutate(entrez_id = gsub("&", ";", sub("\\)", "", sub("^.+ \\(", "", gene))),
gene_name = gsub("&", ";", sub(" \\(.+\\)$", "", gene)))
### left_join `copyNumber` & `depmap_id_to_name_20Q1` on `depmap_id`,
## `cell_line`
copyNumber_20Q1 <- copyNumber_20Q1_long %>%
left_join(depmap_id_to_name_20Q1, by = c("depmap_id" = "depmap_id"))
### rearrange columns into same column format as other datasets
copyNumber_20Q1 <- copyNumber_20Q1 %>%
dplyr::select(depmap_id, gene, log_copy_number,
entrez_id, gene_name, cell_line) %>%
type_convert(cols(entrez_id = "i"))
### visual check
head(copyNumber_20Q1)
### saving cleaned and converted `copyNumber` data as .rda file
save(copyNumber_20Q1, file = "../eh_data/copyNumber_20Q1.rda",
compress = "xz", compression_level = 9)
## access the data on ExperimentHub
# hub <- ExperimentHub()
# x <- query(hub, "depmap")
##########################################
## EH3290 depmap `crispr_20Q1` dataset
##########################################
## data cleaning of `crispr` dataset`
### loading data (downloading .csv file from online source)
url_31 <- "https://ndownloader.figshare.com/files/21521910"
Achilles_gene_effect_20Q1 <- read_csv(url_31)
### gene_effect_corrected.csv data renamed to `crispr`
crispr_20Q1 <- Achilles_gene_effect_20Q1
### visual check
head(crispr_20Q1)
### rename column first column to "depmap_id"
names(crispr_20Q1)[1] <-"depmap_id"
### gather cripsr into long form with columns: `depmap_id`, `gene`, `dependency`
crispr_20Q1_long <- gather(crispr_20Q1, gene, dependency, -depmap_id)
### mutate gene into `gene_name` and `entrez_id`
crispr_20Q1_long <- crispr_20Q1_long %>%
mutate(entrez_id = gsub("&", ";", sub("\\)", "", sub("^.+ \\(", "", gene))),
gene_name = gsub("&", ";", sub(" \\(.+\\)$", "", gene)))
### left_join `crispr_long` and `depmap_id_to_name` to add `cell_line` column
crispr_20Q1 <- crispr_20Q1_long %>% left_join(depmap_id_to_name_20Q1,
by = c("depmap_id" = "depmap_id"))
### rearrange columns into same column format as other datasets
crispr_20Q1 <- crispr_20Q1 %>% dplyr::select(depmap_id, gene,
dependency, entrez_id,
gene_name, cell_line) %>%
type_convert(cols(entrez_id = "i"))
### visual check
head(crispr_20Q1)
### saving cleaned and converted `crispr` data as .rda file
save(crispr_20Q1, file = "../eh_data/crispr_20Q1.rda",
compress = "xz", compression_level = 9)
## access the data on ExperimentHub
# hub <- ExperimentHub()
# x <- query(hub, "depmap")
##########################################
## EH3292 depmap `TPM_20Q1` dataset
##########################################
## data cleaning of `TPM` dataset
### loading data (downloading .csv file from online source)
url_32 <- "https://ndownloader.figshare.com/files/21521940"
CCLE_expression_full_20Q1 <- read_csv(url_32)
### CCLE_depMap_20Q1_TPM.csv data renamed to `TPM`
TPM_20Q1 <- CCLE_expression_full_20Q1
names(TPM_20Q1)
### rename column first column to "depmap_id"
names(TPM_20Q1)[1] <-"depmap_id"
### gather `TPM` into long form on columns: `depmap_id`, `gene`, `expression`
TPM_20Q1_long <- gather(TPM_20Q1, gene, expression, -depmap_id)
### mutate gene into gene_name and ensembl_id
TPM_20Q1_long <- TPM_20Q1_long %>%
mutate(ensembl_id = gsub("&", ";", sub("\\)", "", sub("^.+ \\(", "",gene))),
gene_name = gsub("&", ";", sub(" \\(.+\\)$", "", gene)))
### left_join join `TPM` and `depmap_id_to_name_20Q1` to add `cell_line` column
TPM_20Q1_long %>%
left_join(depmap_id_to_name_20Q1, by = c("depmap_id" = "depmap_id")
) -> TPM_20Q1
### rearrange columns into same column format as other datasets
TPM_20Q1 %>%
dplyr::select(depmap_id, gene, expression, entrez_id, gene_name,
cell_line) %>%
type_convert(cols(entrez_id = "i")) -> TPM_20Q1
### visual check
head(TPM_20Q1)
### saving cleaned and converted `TPM` data as .rda file
save(TPM_20Q1, file = "../eh_data/TPM_20Q1.rda", compress="xz",
compression_level=9)
UCLouvain-CBIO/depmap documentation built on March 24, 2024, 2 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
### =========================================================================
### Make 20Q1 depmap data
### -------------------------------------------------------------------------
###
## This scripts documents how to download and generate the data files.
## Note 1: This scipt assumes (!!!) that it is run in ./depmap/inst/scripts/ and
## saves the resulting .rda files in ./depmap/inst/data
## (e.g. *setwd(./depmap/inst/scripts/)*)
## Note 2: the Broad Institute may change the download urls for these datasets.
## If the link to one of these datasets is broken, please contact the package
## maintainer. All Depmap data can be downloaded from the API at
## https://depmap.org/portal/download/ even if the specific links may change
library("readr")
library("dplyr")
library("tidyr")
library("ExperimentHub")
##########################################
## EH3294 depmap `metadata_20Q1` dataset
##########################################
## data cleaning of `metadata` dataset
### loading data (downloading .csv file from online source)
url_28 <- "https://ndownloader.figshare.com/files/21522000"
sample_info <- read_csv(url_28)
### DepMap-2020q1-celllines.csv data renamed to `metadata`
metadata_20Q1 <- sample_info
# names(metadata_20Q1)
# [1] "DepMap_ID" "stripped_cell_line_name" "CCLE_Name"
# [4] "alias" "COSMIC_ID" "lineage"
# [7] "lineage_subtype" "lineage_sub_subtype" "lineage_molecular_subtype"
# [10] "sex" "source" "Achilles_n_replicates"
# [13] "cell_line_NNMD" "culture_type" "culture_medium"
# [16] "cas9_activity" "RRID" "sample_collection_site"
# [19] "primary_or_metastasis" "disease" "disease_subtype"
# [22] "age" "Sanger_model_ID" "additional_info"
#### Rename `metadata` columns to contain underscores and be in snake case
## note: "metadata_20Q1" has different columns than "metadata_20Q1"
names(metadata_20Q1)[1:24] <-c("depmap_id", "stripped_cell_line_name",
"cell_line", "aliases", "cosmic_id", "lineage",
"lineage_subtype", "lineage_sub_subtype",
"lineage_molecular_subtype", "sex",
"source", "Achilles_n_replicates",
"cell_line_NNMD", "culture_type",
"culture_medium", "cas9_activity", "RRID",
"sample_collection_site", "primary_or_metastasis",
"primary_disease", "subtype_disease", "age",
"sanger_id", "additional_info")
### visual check
View(metadata_20Q1)
### saving cleaned and converted `metadata` data as .rda file
save(metadata_20Q1, file = "../eh_data/metadata_20Q1.rda",
compress = "xz", compression_level = 9)
## access the data on ExperimentHub
# hub <- ExperimentHub()
# x <- query(hub, "depmap")
##########################################
## `depmap_id_to_name_20Q1 map `depmap_id` to `cell_line`
##########################################
## generation of `metadata` subset `depmap_id_to_name_20Q1`
## The subset of `metadata`, `depmap_id_to_name_Q2` is used to map `cell_line`
## and depmap_id via left_join in other depmap datasets that do not contain
## both variables. If you are generating the depmap data from scratch, you will
## need to run the following code to generate the data correctly.
### `depmap_id_to_name` to add `depmap_id` or `cell_line` to other datasets
depmap_id_to_name_20Q1 <- metadata_20Q1 %>% dplyr::select(depmap_id, cell_line)
### visual check
# head(depmap_id_to_name_20Q1)
##########################################
## EH3293 depmap `mutationCalls_20Q1` dataset
##########################################
## data cleaning of `mutationCalls` dataset
### loading data (downloading .csv file from online source)
url_29 <- "https://ndownloader.figshare.com/files/21521967"
CCLE_mutations_20Q1 <- read_csv(url_29)
### depmap_20Q1_mutation_calls data renamed to `mutationCalls`
mutationCalls_20Q1 <- CCLE_mutations_20Q1
# head(mutationCalls_20Q1)
## remove unnecessary X1 variable
mutationCalls_20Q1 <- mutationCalls_20Q1 %>% dplyr::select(-X1)
#### rename last column to depmap_id
# names(mutationCalls_20Q1)
# [1] "Hugo_Symbol" "Entrez_Gene_Id" "NCBI_Build"
# [4] "Chromosome" "Start_position" "End_position"
# [7] "Strand" "Variant_Classification" "Variant_Type"
# [10] "Reference_Allele" "Tumor_Seq_Allele1" "dbSNP_RS"
# [13] "dbSNP_Val_Status" "Genome_Change" "Annotation_Transcript"
# [16] "Tumor_Sample_Barcode" "cDNA_Change" "Codon_Change"
# [19] "Protein_Change" "isDeleterious" "isTCGAhotspot"
# [22] "TCGAhsCnt" "isCOSMIChotspot" "COSMIChsCnt"
# [25] "ExAC_AF" "CGA_WES_AC" "SangerWES_AC"
# [28] "SangerRecalibWES_AC" "RNAseq_AC" "HC_AC"
# [31] "RD_AC" "WGS_AC" "Variant_annotation"
# [34] "DepMap_ID"
## note: "mutationCalls_20Q1" has different columns than "mutationCalls_20Q1"
## the variable "VA_WES_AC" is no longer present in this dataset, unlike
## previous releases (e.g. 19Q1)!
names(mutationCalls_20Q1)[1:34] <- c("gene_name", "entrez_id", "ncbi_build",
"chromosome", "start_pos", "end_pos", "strand",
"var_class","var_type", "ref_allele",
"tumor_seq_allele1", "dbSNP_RS",
"dbSNP_val_status", "genome_change",
"annotation_transcript", "tumor_sample_barcode",
"cDNA_change", "codon_change", "protein_change",
"is_deleterious", "is_tcga_hotspot",
"tcga_hsCnt", "is_cosmic_hotspot",
"cosmic_hsCnt", "ExAC_AF", "CGA_WES_AC",
"sanger_WES_AC", "sanger_recalib_WES_AC",
"RNAseq_AC", "HC_AC", "RD_AC", "WGS_AC",
"var_annotation","depmap_id")
### rearrange columns into same column format as other datasets
mutationCalls_20Q1 <- mutationCalls_20Q1 %>%
dplyr::select(depmap_id, gene_name, entrez_id, ncbi_build, chromosome,
start_pos, end_pos, strand, var_class, var_type, ref_allele,
tumor_seq_allele1, dbSNP_RS, dbSNP_val_status, genome_change,
annotation_transcript, tumor_sample_barcode, cDNA_change,
codon_change, protein_change, is_deleterious, is_tcga_hotspot,
tcga_hsCnt, is_cosmic_hotspot, cosmic_hsCnt, ExAC_AF,
CGA_WES_AC, sanger_WES_AC, sanger_recalib_WES_AC, RNAseq_AC,
HC_AC, RD_AC, WGS_AC, var_annotation)
### visual check
# head(mutationCalls_20Q1)
### saving cleaned and converted `mutationCalls` data as .rda file
save(mutationCalls_20Q1, file = "../eh_data/mutationCalls_20Q1.rda",
compress = "xz", compression_level = 9)
## access the data on ExperimentHub
# hub <- ExperimentHub()
# x <- query(hub, "depmap")
##########################################
## EH3291 depmap `copyNumber_20Q1` dataset
##########################################
## data cleaning of `copyNumber` dataset
### loading data (downloading .csv file from online source)
url_30 <- "https://ndownloader.figshare.com/files/21521964"
CCLE_gene_cn_20Q1 <- read_csv(url_30)
### public_20Q1_gene_cn.csv data renamed to `copyNumber`
copyNumber_20Q1 <- CCLE_gene_cn_20Q1
### rename column first column to "depmap_id"
names(copyNumber_20Q1)[1] <- "depmap_id"
### gather into long form on columns: `depmap_id`, `gene`, `logCopyNumber`
copyNumber_20Q1_long <- gather(copyNumber_20Q1, gene, log_copy_number,
-depmap_id)
### mutate gene column into `gene_name` and `entrez_id`
copyNumber_20Q1_long <- copyNumber_20Q1_long %>%
mutate(entrez_id = gsub("&", ";", sub("\\)", "", sub("^.+ \\(", "", gene))),
gene_name = gsub("&", ";", sub(" \\(.+\\)$", "", gene)))
### left_join `copyNumber` & `depmap_id_to_name_20Q1` on `depmap_id`,
## `cell_line`
copyNumber_20Q1 <- copyNumber_20Q1_long %>%
left_join(depmap_id_to_name_20Q1, by = c("depmap_id" = "depmap_id"))
### rearrange columns into same column format as other datasets
copyNumber_20Q1 <- copyNumber_20Q1 %>%
dplyr::select(depmap_id, gene, log_copy_number,
entrez_id, gene_name, cell_line) %>%
type_convert(cols(entrez_id = "i"))
### visual check
head(copyNumber_20Q1)
### saving cleaned and converted `copyNumber` data as .rda file
save(copyNumber_20Q1, file = "../eh_data/copyNumber_20Q1.rda",
compress = "xz", compression_level = 9)
## access the data on ExperimentHub
# hub <- ExperimentHub()
# x <- query(hub, "depmap")
##########################################
## EH3290 depmap `crispr_20Q1` dataset
##########################################
## data cleaning of `crispr` dataset`
### loading data (downloading .csv file from online source)
url_31 <- "https://ndownloader.figshare.com/files/21521910"
Achilles_gene_effect_20Q1 <- read_csv(url_31)
### gene_effect_corrected.csv data renamed to `crispr`
crispr_20Q1 <- Achilles_gene_effect_20Q1
### visual check
head(crispr_20Q1)
### rename column first column to "depmap_id"
names(crispr_20Q1)[1] <-"depmap_id"
### gather cripsr into long form with columns: `depmap_id`, `gene`, `dependency`
crispr_20Q1_long <- gather(crispr_20Q1, gene, dependency, -depmap_id)
### mutate gene into `gene_name` and `entrez_id`
crispr_20Q1_long <- crispr_20Q1_long %>%
mutate(entrez_id = gsub("&", ";", sub("\\)", "", sub("^.+ \\(", "", gene))),
gene_name = gsub("&", ";", sub(" \\(.+\\)$", "", gene)))
### left_join `crispr_long` and `depmap_id_to_name` to add `cell_line` column
crispr_20Q1 <- crispr_20Q1_long %>% left_join(depmap_id_to_name_20Q1,
by = c("depmap_id" = "depmap_id"))
### rearrange columns into same column format as other datasets
crispr_20Q1 <- crispr_20Q1 %>% dplyr::select(depmap_id, gene,
dependency, entrez_id,
gene_name, cell_line) %>%
type_convert(cols(entrez_id = "i"))
### visual check
head(crispr_20Q1)
### saving cleaned and converted `crispr` data as .rda file
save(crispr_20Q1, file = "../eh_data/crispr_20Q1.rda",
compress = "xz", compression_level = 9)
## access the data on ExperimentHub
# hub <- ExperimentHub()
# x <- query(hub, "depmap")
##########################################
## EH3292 depmap `TPM_20Q1` dataset
##########################################
## data cleaning of `TPM` dataset
### loading data (downloading .csv file from online source)
url_32 <- "https://ndownloader.figshare.com/files/21521940"
CCLE_expression_full_20Q1 <- read_csv(url_32)
### CCLE_depMap_20Q1_TPM.csv data renamed to `TPM`
TPM_20Q1 <- CCLE_expression_full_20Q1
names(TPM_20Q1)
### rename column first column to "depmap_id"
names(TPM_20Q1)[1] <-"depmap_id"
### gather `TPM` into long form on columns: `depmap_id`, `gene`, `expression`
TPM_20Q1_long <- gather(TPM_20Q1, gene, expression, -depmap_id)
### mutate gene into gene_name and ensembl_id
TPM_20Q1_long <- TPM_20Q1_long %>%
mutate(ensembl_id = gsub("&", ";", sub("\\)", "", sub("^.+ \\(", "",gene))),
gene_name = gsub("&", ";", sub(" \\(.+\\)$", "", gene)))
### left_join join `TPM` and `depmap_id_to_name_20Q1` to add `cell_line` column
TPM_20Q1_long %>%
left_join(depmap_id_to_name_20Q1, by = c("depmap_id" = "depmap_id")
) -> TPM_20Q1
### rearrange columns into same column format as other datasets
TPM_20Q1 %>%
dplyr::select(depmap_id, gene, expression, entrez_id, gene_name,
cell_line) %>%
type_convert(cols(entrez_id = "i")) -> TPM_20Q1
### visual check
head(TPM_20Q1)
### saving cleaned and converted `TPM` data as .rda file
save(TPM_20Q1, file = "../eh_data/TPM_20Q1.rda", compress="xz",
compression_level=9)
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.