R中无报错使用mas5calls函数？miRNA芯片CDF包缺失求助

问题重现 #

执行代码 #

setwd("C:/A549_ALI/4_tert-Butanol (22)/")
list.celfiles()
my.affy=ReadAffy()
dim(exprs(my.affy))

正常输出 #

Show in New Window
     [1] "(46) 22-B1-1_(miRNA-4_0).CEL"   
     [2] "(47) 22-B1-2_(miRNA-4_0).CEL"   
     [3] "(48) 22-B1-3_(miRNA-4_0).CEL"   
     [4] "(49) 22-R1-1_(miRNA-4_0).CEL"   
     [5] "(50) 22-NEC 1-1_(miRNA-4_0).CEL"
     [6] "(51) 22-B2-1_(miRNA-4_0).CEL"   
     [7] "(52) 22-B2-2_(miRNA-4_0).CEL"   
     [8] "(53) 22-B2-3_(miRNA-4_0).CEL"   
     [9] "(54) 22-R2-1_(miRNA-4_0).CEL"   
    [10] "(55) 22-NEC 2-1_(miRNA-4_0).CEL"
    [11] "(56) 22-B3-1_(miRNA-4_0).CEL"   
    [12] "(57) 22-B3-2_(miRNA-4_0).CEL"   
    [13] "(58) 22-B3-3_(miRNA-4_0).CEL"   
    [14] "(59) 22-R3-1_(miRNA-4_0).CEL"   
    [15] "(60) 22-NEC 3-1_(miRNA-4_0).CEL"
    [1] 292681     15

调用mas5calls的错误信息 #

background correction: mas 
PM/MM correction : mas 
expression values: mas 
background correcting...'getOption("repos")' replaces Bioconductor standard repositories, see '?repositories' for details

replacement repositories:
    CRAN: https://cran.rstudio.com/

Error in getCdfInfo(object) : 
  Could not obtain CDF environment, problems encountered:
Specified environment does not contain miRNA-4_0
Library - package mirna40cdf not installed
Bioconductor - mirna40cdf not available

替代解决方案 #

• 方法1：手动指定CDF文件路径
  如果能获取到与miRNA-4_0芯片匹配的.cdf格式文件，读取CEL时直接指定路径：

  my.affy <- ReadAffy(cdfname = "/path/to/miRNA-4_0.cdf")
  mas5calls(my.affy)
  

  CDF文件需与芯片型号严格匹配，可从实验原始数据或Affymetrix官方资源中获取

• 方法2：使用oligo包处理
  oligo包对Affymetrix芯片兼容性更强，无需依赖特定CDF包，步骤如下：

  1. 安装并加载包

     if (!requireNamespace("BiocManager", quietly = TRUE))
         install.packages("BiocManager")
     BiocManager::install("oligo")
     library(oligo)
     

  2. 读取文件并执行MAS5调用

     setwd("C:/A549_ALI/4_tert-Butanol (22)/")
     celFiles <- list.celfiles()
     rawData <- read.celfiles(celFiles)
     mas5Calls <- mas5calls(rawData)
     

• 方法3：尝试安装旧版Bioconductor中的CDF包
  R 4.2.2对应Bioconductor 3.16，可尝试从该版本归档安装：

  BiocManager::install("mirna40cdf", version = "3.16")
  

内容的提问来源于stack exchange，提问作者BerKa