R/intermediate_feather_file_creator-cellimage_network_files.R

In CRI-iAtlas/iatlas-feather-files: iAtlas Data

build_cellimage_network_files <- function(){

  hgnc_labels <- "syn21781366" %>%
    synapse_feather_id_to_tbl %>%
    dplyr::select("value" = "Variable", "x", "y")

  hgnc_scaffold <- "syn21782167" %>%
    synapse_feather_id_to_tbl %>%
    dplyr::select("from" = "From", "to" = "To", "label" = "interaction") %>%
    dplyr::mutate("number" = 1:dplyr::n()) %>%
    tidyr::pivot_longer(-c("label", "number"))

  tcga_genes <- "syn22133677" %>%
    synapse_feather_id_to_tbl %>%
    tidyr::drop_na() %>%
    dplyr::mutate("entrez" = as.character(entrez))

  genes <- "syn22240716" %>%
    synapse_feather_id_to_tbl %>%
    tidyr::drop_na() %>%
    dplyr::mutate("entrez" = as.character(entrez))

  cell_features <- c(
    "B_cells",
    "Dendritic_cells",
    "Eosinophils",
    "Macrophage",
    "Mast_cells",
    "NK_cells",
    "Neutrophils",
    "T_cells_CD4",
    "T_cells_CD8"
  )

  other_features <- dplyr::tibble(
    "value" = "Tumor_cell",
    "new_value" = "Tumor_fraction"
  )

  hgnc_translation <- hgnc_scaffold %>%
    dplyr::select("value") %>%
    dplyr::distinct() %>%
    dplyr::inner_join(genes, by = c("value" = "hgnc")) %>%
    dplyr::rename("new_value" = "entrez")

  cell_translation <-
    hgnc_scaffold %>%
    dplyr::filter(value %in% c(cell_features)) %>%
    dplyr::select("value") %>%
    dplyr::distinct() %>%
    dplyr::mutate("new_value" = stringr::str_c(value, "_Aggregate2")) %>%
    dplyr::bind_rows(other_features)


  translation <-
    dplyr::bind_rows(
      hgnc_translation, cell_translation
    )

  entrez_scaffold <- hgnc_scaffold %>%
    dplyr::left_join(translation) %>%
    dplyr::select(-value) %>%
    dplyr::rename(value = new_value) %>%
    tidyr::pivot_wider() %>%
    tidyr::drop_na() %>%
    dplyr::select(-number)

  entrez_labels <- hgnc_labels %>%
    dplyr::inner_join(translation, by = "value") %>%
    dplyr::select("node" = "new_value", "x", "y")

  synapse_store_feather_file(
    entrez_scaffold,
    "cellimage_scaffold.feather",
    "syn23518511"
  )

  synapse_store_feather_file(
    entrez_labels,
    "cellimage_node_labels.feather",
    "syn23518511"
  )

}