【发布时间】:2012-04-30 15:49:08
【问题描述】:
我已经创建了一个绘图函数:
plotfun = function(dat, pos, tetrad.name=NULL, snp.cols=c("blue", "red", "lightgray")){
plot(1,1, type="n", xlim=c(min(pos), max(pos)), ylim=c(0.75,4.25), yaxt="n", xaxt="n", main=paste("CCT6 -",tetrad.name), xlab="", ylab="", bty="n", cex=.75)
abline(h=1:4, lty="dotted")
tmpnull = sapply(1:4, function(n) points(pos, rep((4:1)[n],length(pos)), col=snp.cols[dat[n,]], pch=16))
mtext(c("a", "b", "c", "d"), side=2, line=0.5, at=4:1, las=1, cex=.5)
}
我正在尝试使用此绘图函数来绘制我从更大列表中索引的矩阵子集(在此示例中为x)。我正在使用sapply 分别绘制所有矩阵:
sapply(1:69, function(n) plotfun(dat=matrix.transform(x[[n]]), pos=cct6.pos, tetrad.name= n ))
我的问题是我不知道如何将绘图标题读取为附加到每个矩阵的列表编号。相反,它将它们标记为1:69,从而丢失了索引输入中的信息。例如,如果第一个输入是cct6[[35]],我希望情节标题改为"CCT6-35" 而不是"CCT6-1"
抱歉,这是一个基本问题。我只是卡在这个...
【问题讨论】: