小编jay*_*mee的帖子

我正在尝试生成一个条形图,但是y轴刻度太短.这是我的代码:

barplot(as.matrix(dat), log="y", ylim=c(10000,100000000), beside=TRUE,
        ylab = "Number of reads", col = c("gray","black","white"))

它离开了轴的空间(根据ylim),但没有填充实际的轴.我经历过?barplot并尝试了一些事情(从谷歌搜索我认为xpd = F, yaxs = c(10000,10000000,5)应该工作,但它没有).

我知道这是一件小事,但这正是我长期坚持的问题,而不是实际工作,所以任何帮助都会非常感激!

编辑:为输入人员干杯!

我最初没有使用ylim进行绘图,但最终会出现一个更奇怪的轴(同样的问题); 我实际上选择了我的ylim值来给它一个更好的间隔轴.

这是数据:

dat <- read.table(text="D2,D3n,D3m,D4n,D4m
21234722,34262282,31920464,25486357,20712943
35343,64403,22537,39934,46547
126646,312286,101105,81537,76944", header=TRUE, sep=",")

编辑2:@DWin没错 - 我更新了我的R,现在情节很好 - 谢谢大家!

我很欣赏,考虑到我的知识水平(或缺乏),这项任务可能有点雄心勃勃,但仍然如此.

我有一个包含16个字符串的列表,大约3000个项目,其中每个字符表示另一个数字列表.不确定我是否说得那么清楚; 它实际上是16个氨基酸长肽的列表,其中每个氨基酸(20个中的1个)可由5个数字表示.

我想遍历那个列表(肽段),然后为每个字符(氨基酸)添加相关的5个数字(Atchley因子,如果你有兴趣)到一个数组,制作一个三维数组,我的轴是:该肽中的肽(3000)x氨基酸的实例(16)x因子(5).

我不可思议地超出了我的深度,所以我不确定我所拥有的是否有用是有用的,但这里是(使用numpy):

array = np.empty(shape=(len(peptides),16,5)

for i in peptides:

    for j in str(i):

(此时我尝试了一堆不同的东西,因为我在论坛上搜索,结尾的东西有点像这样,但我确信我甚至错过了我的目标)

    if j == 'A':    
            L16Afctrs = np.append([-0.59145974, -1.30209266, -0.7330651, 1.5703918, -0.14550842], axis=1)
    elif j == 'C':
            L16Afctrs = np.append([-1.34267179, 0.46542300, -0.8620345, -1.0200786, -0.25516894], axis=1)
    ...
    elif j == 'Y':
            L16Afctrs = np.append([0.25999617, 0.82992312, 3.0973596, -0.8380164, 1.51150958], axis=1)

就像我说的那样,我老老实实地在努力,任何帮助都会非常感激.

编辑:澄清(希望如此)

我有一个大约3000个不同的16个字符串的列表,其中这些字符串中的每个字符表示另外5个数字.

我想生成一个三维数组或结构,通过查看三维数组中的给定平面(我设想的维度是原始的),我可以(最终)绘制所有3000个字符串中给定位置的5个数字.字符串x 16个字符x 5个因子).

我目前正在制作不同角色的字典,与@Winston的帖子相关,然后尝试将其折叠成3d数组.

编辑2:成功!

温斯顿的修复效果非常好!