1

我需要绘制微阵列时程数据,以便在图表上显示前 100 个基因的基因表达谱,并用复制探针标记曲线(信息如下)。

我能够使用正确的标题和轴标签绘制图形框架,但曲线上的标签没有正确的名称,使用以下代码:

  for(j in 1:12){
  pch_value = as.character(targets$grp[6*j])
  points(c(0, 6, 24, 48, 72,80), exprs.row[(6*j-1):(6*j)],
  type='b', pch=pch_value) 

对比:

mc = makeContrasts("hr6-hr0","hr24-hr0","hr48-hr0","hr72-hr0","hr80-hr0","hr72-hr6", levels = X)



colnames(X) = c('hr0', 'hr24', 'hr48', 'hr6', 'hr72', 'hr80')

复制信息:

Replicates (hr0 - 6 reps, hr 6 - 6reps, hr24-4reps, hr48-4reps,hr72-6reps, hr80-2reps) are in the following order :
hr0
hr0
hr6
hr6
hr24
hr24
hr48
hr48
hr72
hr72
hr0
hr0
hr6
hr6
hr24
hr24
hr48
hr48
hr72
hr72
hr80
hr80
hr0
hr0
hr6
hr6
hr72
hr72

我尝试过 maSigPro 和 timecourse;我无法使用这两个程序设置正确的对比度。

帮助将不胜感激。

4

0 回答 0