Skip to content

Commit 2824705

Browse files
authored
Update README.md
1 parent a41eacd commit 2824705

File tree

1 file changed

+32
-0
lines changed

1 file changed

+32
-0
lines changed

README.md

Lines changed: 32 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -20,6 +20,38 @@ The immune defense function protecting body from invasive patho-gens is a key in
2020
- GSVA 1.50.1
2121
- clusterProfiler 4.10.1
2222
- msigdbr 7.5.1
23+
# 🧬 ssGSEA Immune Pathway Enrichment Analysis (C7 Gene Sets)
24+
25+
This R script performs **single-sample Gene Set Enrichment Analysis (ssGSEA)** using the **MSigDB C7 immunologic signature gene sets**. It takes a gene expression matrix as input and outputs enrichment scores for each sample across all C7 pathways.
26+
27+
---
28+
29+
## 📦 Required R Packages
30+
31+
Make sure the following packages are installed:
32+
33+
34+
> ⚠️ **Note**: This script only supports **human (Homo sapiens)** data due to the use of `org.Hs.eg.db`.
35+
36+
---
37+
38+
## 📁 Input Format
39+
40+
### ✅ Requirements
41+
- **File type**: CSV (comma-separated)
42+
- **Rows**: Genes (using official **HGNC gene symbols**, e.g., `TP53`, `CD3D`)
43+
- **Columns**: Samples
44+
- **Values**: Expression values
45+
- Use **raw counts** (e.g., from RNA-seq) → script uses `kcdf="Poisson"`
46+
- For **log-transformed data** (e.g., log2(TPM+1)), you must modify the script to use `kcdf="Gaussian"`
47+
48+
Example (`input.csv`):
49+
```csv
50+
,SAMPLE_A,SAMPLE_B,SAMPLE_C
51+
TP53,120,98,105
52+
CD3D,450,512,488
53+
FOXP3,320,290,310
54+
...
2355
2456
# Cite
2557

0 commit comments

Comments
 (0)