Untitled

#Gender Check

library(Seurat)

## Warning: package 'Seurat' was built under R version 4.3.3

## Loading required package: SeuratObject

## Warning: package 'SeuratObject' was built under R version 4.3.3

## Loading required package: sp

## 
## Attaching package: 'SeuratObject'

## The following objects are masked from 'package:base':
## 
##     intersect, t

library(ggplot2)

## Warning: package 'ggplot2' was built under R version 4.3.3

sickle_cell <- readRDS("C:/Users/skacar/OneDrive - Indiana University/PULMONARY POST DOC/SICKLE_PROJECT/sickle_cell_anemia_integrated_1_14_25.Rds")
mouse_sex_genes <- c(
  # Y-linked (Male-specific) genes
  "Ddx3y", "Eif2s3y", "Kdm5d", "Uty", "Zfy1", "Zfy2", "Sry", "Usp9y", 
  "Rbmy1a1", "Tspy1")
  
  mouse_sex_genes <- c(# X-linked genes (higher expression in females)
  "Xist", "Eif2s3x", "Kdm6a", "Rbm3", "Jarid1c", "Tspyl2", "Bcor", "Tmem29", 
  "Pir", "Eif1ax")
  
  mouse_sex_genes <- c(  # Autosomal genes with sex-biased expression
  "Foxl2", "Sox9", "Esr1", "Ar", "Dmrt1"
)


DotPlot(sickle_cell, mouse_sex_genes, group.by = "group", scale.min = 0) + coord_flip() + theme(axis.text.x = element_text(angle = 45, hjust = 1)) + ggtitle(" Autosomal genes with sex-biased expression")

## Warning: The following requested variables were not found: Foxl2

## Warning: Scaling data with a low number of groups may produce misleading
## results

exists("sickle_cell")

## [1] TRUE