combat_tumor: Process and Correct Batch Effects in Tumor Data

In TransProR: Analysis and Visualization of Multi-Omics Data

Process and Correct Batch Effects in Tumor Data

Description

The function first extracts histological types from the provided tumor data set. After displaying these types, the user is prompted to input specific types to retain. The data is then filtered based on this input.

Note: This example assumes that different tumor types represent different batches in a general sense. Users need to adjust the batch and group vectors based on real-life scenarios.

Usage

combat_tumor(
  tumor_data_path,
  CombatTumor_output_path,
  auto_mode = FALSE,
  default_input = "01,06"
)

Arguments

tumor_data_path	The path to the tumor data stored in an RDS file.
CombatTumor_output_path	A character string specifying the path where the output RDS file will be saved.
auto_mode	Logical. If set to TRUE, the function will not prompt the user for input and will instead use the values provided in default_input. Default is FALSE.
default_input	Character string. When auto_mode is TRUE, this parameter specifies the default tumor types to be retained. It should be provided as a comma-separated string (e.g., "01,06").

Details

This function takes a tumor data set, asks the user for specific tumor types to retain, and then corrects for batch effects using the ComBat_seq function from the 'sva' package.

The ComBat_seq function from the sva package is used to correct batch effects. The function requires the 'sva' package to be installed and loaded externally.

Value

A data.frame with corrected values after the ComBat_seq adjustment. Note that this function also saves the combat_count_df data as an RDS file at the specified output path.

See Also

ComBat_seq

Examples

tumor_file <- system.file("extdata",
                          "SKCM_Skin_TCGA_exp_tumor_test.rds",
                          package = "TransProR")
output_file <- file.path(tempdir(), "SKCM_combat_count.rds")

  SKCM_combat_count <- combat_tumor(
  tumor_data_path = tumor_file,
  CombatTumor_output_path = output_file,
  auto_mode = TRUE,
  default_input = "01,06"
)

head(SKCM_combat_count)[1:5, 1:5]

TransProR documentation built on April 4, 2025, 3:16 a.m.