主成分分析PCA烟花状连线图

rm(list=ls())
# 生物信息学PCA分析：烟花状连线图（适配自定义数据+Lancet配色）
# 数据文件要求：
# - fireworks_group.txt：两列（Sample, Group），制表符分隔
# - fireworks_data.txt：行=基因，列=样本（第一列是基因名，制表符分隔）

# 安装依赖（若未安装）
if (!requireNamespace("ggsci", quietly = TRUE)) {
  install.packages("ggsci")
}
library(ggplot2)
library(ggsci)  # 提供Lancet配色


### 1. 数据读取与整理 ###
# 读取分组信息（Sample → Group 映射）
group_df <- read.table(
  "fireworks_group.txt", 
  header = TRUE, 
  sep = "\t", 
  stringsAsFactors = FALSE
)

# 读取基因表达数据（行=基因，列=样本；第一列作为行名）
data_df <- read.table(
  "firew