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