push script

Author: kpatel427

getOrthologs.R

@@ -0,0 +1,16 @@
+# script to get all genes orthologous to human in mouse
+
+library(tidyverse)
+library(homologene)
+
+# Retrieve all orthologs between human and mouse
+all_human_mouse_orthologs <- homologeneData2 |> 
+  filter(Taxonomy %in% c(9606, 10090))
+
+# Split human and mouse data
+human_genes <- all_human_mouse_orthologs |>  filter(Taxonomy == 9606)
+mouse_genes <- all_human_mouse_orthologs |>  filter(Taxonomy == 10090)
+
+
+# Merge human and mouse data by HID to get ortholog pairs
+human_mouse_pairs <- merge(human_genes, mouse_genes, by = "HID", suffixes = c("_human", "_mouse"))