【怪毛匠子整理】
一、使用SOFT文件
1. 確定GEO平台,比如GPL25318, 通過NCBI找到這個平台,也可以通過下面的鏈接,把GPL25318 替換成你想要的平台ID。
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL25318
2. 下載SOFT文件

點開紅框:
下載文件:

解壓文件:

R語言處理# 讀入文件
df2=read.table("GPL25318_family.soft",sep = "\t")
# 使用merge合並
【怪毛匠子整理】
二、使用GEOquery包
library(GEOquery)
GPL=getGEO("GPL25318",destdir = ".")
g1=Table(GPL)
head(exprSet$ID)
[1] "200000_s_at" "200001_at" "200002_at" "200003_s_at" "200004_at" "200005_at"
new_datasets=merge(exprSet,g1,by = "ID" ,all.x = T)
colnames(new_datasets)
[1] "ID" "GSM******" "GSM******"
.....
[247] "GSM******" "GSM******" "GSM******"
[250] "Ensembl Version" "Ensembl Genome Version" "Orientation"
[253] "Probeset mapping position" "Ensembl Gene ID" "Gene Symbol"
[256] "Chromosomal location" "Strand" "Gene Description"
[259] "Entrez Gene" "SPOT_ID"
這樣就搞定了!
【怪毛匠子整理】
