拉普拉胯 2022-10-18 01:15 采纳率: 100%
浏览 652
已结题

GEO平台无对应R包且下载GPL后无table信息怎么办?

GEO下载数据,在ID转换步骤卡死,根据正常方法:
1、有对应R包直接用R包——没有,放弃该法
2、直接下载GPL soft文件用R提取symbol等信息:
出现问题如下

gpl <- getGEO('GPL24676', destdir=".")
colnames(Table(gpl)) 

gpl <- getGEO('GPL24676', destdir=".")
Using locally cached version of GPL24676 found here:
./GPL24676.soft.gz
colnames(Table(gpl))
character(0)

如图所示,正常GPL应该会出现:
[1] "ID" "SPOT_ID" "CONTROL_TYPE"
[4] "REFSEQ" "GB_ACC" "GENE"
[7] "GENE_SYMBOL" "GENE_NAME" "UNIGENE_ID"
[10] "ENSEMBL_ID" "TIGR_ID" "ACCESSION_STRING"
[13] "CHROMOSOMAL_LOCATION" "CYTOBAND" "DESCRIPTION"
[16] "GO_ID" "SEQUENCE"
此类信息
但GPL24676运行出来col为0,直接卡死
3、搜芯片来源,发现其为Illumina NovaSeq 6000 (Homo sapiens),0相关贴,卡死
4、尝试过这种注释信息不全的万能法(如下图)了,也失败了,卡在提取下载的GPL数据的步骤上,它是还是0

BiocManager::install("clusterProfiler")
BiocManager::install("org.Hs.eg.db")

求问,这种数据怎么进行基因ID转换!

img

上图是GSE的矩阵,可以看到和别的数据不一样,它的start和end中间没有任何数据。

SOS!

  • 写回答

1条回答 默认 最新

  • Bioinfo Guy R语言领域新星创作者 2022-10-18 09:09
    关注

    作者都放在补充文件里面了,在geo官网下载,是SARS-CoV-2的吧,我看到有raw文件

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 10月31日
  • 已采纳回答 10月23日
  • 创建了问题 10月18日

悬赏问题

  • ¥20 wireshark抓不到vlan
  • ¥20 关于#stm32#的问题:需要指导自动酸碱滴定仪的原理图程序代码及仿真
  • ¥20 设计一款异域新娘的视频相亲软件需要哪些技术支持
  • ¥15 stata安慰剂检验作图但是真实值不出现在图上
  • ¥15 c程序不知道为什么得不到结果
  • ¥40 复杂的限制性的商函数处理
  • ¥15 程序不包含适用于入口点的静态Main方法
  • ¥15 素材场景中光线烘焙后灯光失效
  • ¥15 请教一下各位,为什么我这个没有实现模拟点击
  • ¥15 执行 virtuoso 命令后,界面没有,cadence 启动不起来