--- a
+++ b/R/SeekGtexOrgan.R
@@ -0,0 +1,39 @@
+#' Load and Process GTEX Phenotype Data to Retrieve Primary Site Counts
+#'
+#' This function reads the GTEX phenotype data from a specified path, renames its columns for better readability,
+#' and then returns a table of primary site counts.
+#'
+#' @param path The path to the GTEX phenotype data file. Default is "./download_data/GTEX_phenotype".
+#'
+#' @return A table representing the count of samples per primary site.
+#' @importFrom utils read.table
+#' @examples
+#' # Get the file path to the example data in the package
+#' path <- system.file("extdata", "GTEX_phenotype_test", package = "TransProR")
+#' # Call the `seek_gtex_organ` function with the path and print the result
+#' SeekGtexOrgan <- seek_gtex_organ(path = path)
+#'
+#' @export
+
+seek_gtex_organ <- function(path = "./download_data/GTEX_phenotype") {
+  # Read GTEX phenotype data
+  # gtex.phe <- data.table::fread(path, header = TRUE, sep = '\t', data.table = FALSE)
+
+  gtex.phe <- utils::read.table(path,
+                        header = TRUE,
+                        sep = '\t',
+                        stringsAsFactors = FALSE,
+                        check.names = FALSE)
+
+  rownames(gtex.phe) <- gtex.phe$Sample
+
+  # Rename columns
+  colnames(gtex.phe) <- c("Sample", "body_site_detail (SMTSD)", "primary_site", "gender", "patient", "cohort")
+
+  # Create table of primary sites
+  primary_site_counts <- table(gtex.phe$primary_site)
+
+  return(primary_site_counts)
+}
+
+