Построение подмножеств объекта AffyRNAdeg {affy} с помощью plot AffyRNAdeg?

library(affy)    
microarrays <- ReadAffy()         # 98 CEL files are read into the same object
RNAdeg <- AffyRNAdeg(microarrays)

Теперь я хочу построить подмножества RNAdeg

plotAffyRNAdeg(RNAdeg[.......?])   # What can I do?

Я безуспешно пробовал различные циклы for.

Но если указаны цвета линии графика, plotAffyRNAdeg отображает подмножество 1:(количество указанных цветов), но я не придумал, как это эффективно использовать. Например, ниже представлены графики с первого по шестой набор данных микрочипа AffyRNAdeg (с первого по шестой файл .CEL, прочитанный ReadAffy() )

plotAffyRNAdeg(RNAdeg,col=c(2,2,2,3,3,3))

plot subset bioconductor
person matt    schedule 07.06.2013    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Итак, один из способов был найден путем запуска AffyRNAdeg() на подмножествах объекта, в котором находятся файлы CEL, и помещения полученных данных в список списков, организованных экспериментально, а затем построения списка элементов. Может быть, есть более простой способ, но это сработало (я новичок в R).

library(affy)    
library(RColorBrewer)

> sampleNames(ARTHwoundMA[,11:14])
[1] "GSE18960_05_GSM469416_trt_rep2.CEL" "GSE18960_06_GSM469418_trt_rep3.CEL"
[3] "GSE5525_GSM128715_ctrl12h.CEL"      "GSE5525_GSM128716_ctrl24h.CEL

# RNA DEG

# Indices to subset by experiment

cel_names <- substr(sampleNames(ARTHwoundMA),1,7)
unique_exp <- unique(substr(sampleNames(ARTHwoundMA),1,7))
exp_ind <- list()
for (i in 1:length(unique_exp))
  { 
  tempvec <- vector()
  for (j in 1:length(cel_names))
  {
    if (cel_names[j]==unique_exp[i])
    {
      tempvec <- append(tempvec,j)
    }
  }
  exp_ind[[(length(exp_ind)+1)]] <- tempvec
}


# Calculating

RNAdeg_exp <- list()
for(i in 1:length(exp_ind))
{
  RNAdeg_exp[[i]] <- AffyRNAdeg(ARTHwoundMA[,exp_ind[[i]]]) 
}


# Plotting

colors <- colorRampPalette(rev(brewer.pal(9, "Reds")))(length(exp_ind[[i]])
pdf(file="C:\\R working directory\\TEST\\RNAdeg_plots.pdf")
for(i in 1:length(exp_ind))
{
  par(bg="gray")
  colors <- colorRampPalette(rev(brewer.pal(9, "Reds")))(length(exp_ind[[i]]))
  plotAffyRNAdeg(RNAdeg_exp[[i]], col=colors)
  plot.new()
  legend("topleft", lty=1, lwd=2,col=colors,
     legend=paste(sampleNames(ARTHwoundMA[,exp_ind[[i]]])))
}

dev.off()
person matt    schedule 08.06.2013