小编use*_*000的帖子

我正在尝试在ggplot2中的R中创建一个简单的密度图.这是我的代码很好用.

d <-  ggplot(result, aes(x=result$baseMeanA)) 
d + geom_density(colour="darkgreen", size=2, fill="darkgreen") + 
scale_x_log10() + scale_y_continuous(limits = c(0, 0.45))

问题是我不能像我想的那样将x轴调整为负数.

scale_x_log10(limits= c(1, 10000))

效果很好,但是

scale_x_log10(limits= c(-1, 10000))

根本不起作用!它给了我这个错误:

if(zero_range(range)){:缺少值需要TRUE/FALSE时出错

请帮忙!

我真的需要你的 R 技能。已经研究这个情节好几天了。我是 R 新手，所以这可以解释它。

我有染色体的序列覆盖数据（基本上是每个染色体长度上每个位置的值，使向量的长度达到数百万）。我想为我的阅读内容制作一个很好的覆盖图。这是我到目前为止得到的：

看起来不错，但我缺少 y 标签，这样我就可以知道它是哪条染色体，而且我在修改 x 轴时也遇到了麻烦，所以它在覆盖范围结束的地方结束。此外，我自己的数据要大得多，使得这个图特别需要很长时间。这就是我尝试这个 HilbertVisplotLongVector 的原因。它有效，但我不知道如何修改它、x 轴、标签、如何记录 y 轴，以及矢量在绘图上都获得相同的长度，即使它们的长度不相等。

source("http://bioconductor.org/biocLite.R")
biocLite("HilbertVis")
library(HilbertVis)
chr1 <- abs(makeRandomTestData(len=1.3e+07)) 
chr2 <- abs(makeRandomTestData(len=1e+07)) 

par(mfcol=c(8, 1), mar=c(1, 1, 1, 1), ylog=T)

# 1st way of trying with some code I found on stackoverflow
# Chr1
plotCoverage <- function(chr1, start, end) { # Defines coverage plotting function.
  plot.new()
  plot.window(c(start, length(chr1)), c(0, 10))
  axis(1, labels=F) 
  axis(4)
  lines(start:end, log(chr1[start:end]), type="l")
}
plotCoverage(chr1, start=1, end=length(chr1)) # Plots coverage result.

# Chr2
plotCoverage <- …

我一直在努力做这个盒子图一整天.我很快就放弃了.请帮忙.我只想为我的3组A,B和C提供一个简单的箱线图.这是我简化的数据:

> mydata 
   A01 A02 A03 B01 B02 B03 C01 C02 C03
id0001    1    2    3    6    7    8    11    12    13

A01,A02和A03应该弥补的置信区间为A,等等.我怎么能有3组在x轴一个非常简单的箱线图？虽然盒子将是巨大的,我会很感激,如果有人可以帮助我让我的情节相似,这一个.最后,我想打一个包含多个ID的图形,所以它看起来像这一个,但误差棒像第一个.

如果有人可以帮助我或给我如何进行的提示,我会非常感激!

谢谢.