我需要绘制微阵列时程数据,以便在图表上显示前 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;我无法使用这两个程序设置正确的对比度。
帮助将不胜感激。