小编use*_*664的帖子

我有一个包含100个样本的ExpressionSet对象:

> length(sampleNames(eset1))
100

我还有一个75个样本名称的向量(不是数据本身):

> length(vecOf75)
75

如何eset1根据75个样本名称进行子集化(并保存)？也就是说,我想忽略那些eset1名字未列出的样本vecOf75.请记住,对应于75个样本名称的一些样本可能不在eset1.从而,

> length(sampleNames(eset1))

现在应该给出<75.

我想在单个绘图窗口中创建多个绘图，其中每个绘图的面板宽度与每个绘图的 xlim 成正比。

目前我使用：

layout(matrix(c(1:8,10,9), 5, 2, byrow = FALSE), widths=2)       
layout.show(10)

基本上，我想width单独应用于每个图而不是列中的所有图。这样做的最佳方法是什么？

1. 我有一段图表的代码: barplot(as.vector(t(mat1[1,3:ncol(mat1)])),las=2) 我想改变,以便x轴被线y = 2 取代; 有效地将x轴向上移动2个单位,如下图所示.

   我需要从2开始的条形图,以便:

   • 值为3的条形从y = 2线开始,并在y = 3处上升到结束.
   • 值为0的条形从y = 2行开始,并在y = 0处下降到结束

2. 我怎样才能使列名mat1我的X轴类别？

我有两个相关的问题.

问题1:我目前正在使用下面的代码生成叠加密度图的直方图:

hist(x,prob=T,col="gray")
axis(side=1, at=seq(0,100, 20), labels=seq(0,100,20))
lines(density(x))

我贴的数据(即x上图)在这里.

我对代码有两个问题:

• 最后一个刻度线和x轴的标签(100)没有出现在直方图/图上.我怎么能把它们放在上面？
• 我希望y轴是计数或频率而不是密度,但我想将密度图保留为直方图上的叠加图.我怎样才能做到这一点？

问题2:对问题1使用类似的解决方案,我现在想要覆盖三个密度图(不是直方图),再次使用y轴上的频率而不是密度.这三个数据集位于:

1. http://pastebin.com/z5X7yTLS

2. http://pastebin.com/Qg8mHg6D

3. http://pastebin.com/aqfC42fL