r - 突出情节中的几个点

Question

我尝试用几个要突出显示的位置(点)绘制以下数据

plot(b$pos,b$log_p,col==ifelse(b$pos==c(14824849,13920386,14837470),90,100), pch=19, xlab='Chromosome 21 position', ylab='-log10(p)')

生成的图，仅显示一个点以红色突出显示，并带有以下警告消息:

In b$pos == c(14824849, 13920386,14837470) : longer object length is not a multiple of shorter object length

Answer 1

好的，问题很可能是你在 ifelse 中的条件。如果您在 ifelse() 之外尝试条件 (b$pos==c(14824849,13920386,14837470))，您将收到如下错误消息:

longer object length is not a multiple of shorter object length

如果将条件更改为:

b$pos %in% c(14824849,13920386,14837470)

您将获得一个 TRUE/FALSE 值向量，确定 b$pos 中的每个条目是否存在于向量 (14824849,13920386,14837470) 中，而不是 b$pos 中的条目是否等于 c(14824849,13920386) ,14837470).

x = c(49, 7, 66, 51, 43, 70, 35, 53, 6, 29)
y = c(10, 98, 44, 31, 37, 14, 64, 84,  4, 34)
x %in% c(6, 7)

[1] FALSE TRUE FALSE FALSE FALSE FALSE FALSE FALSE TRUE FALSE

plot(x, y, col=ifelse(x %in% c(6, 7), 'red', 'blue'))

现在这个数据集有 10 个 x 值，如果你这样写:

plot(x, y, col=ifelse(x == c(1, 7), 'red', 'blue'))

这可以正常工作，x 值将交替与 1 和 7 进行比较，例如:

49 == 1 ?
7 == 7 ?
66 == 1?
51 == 7? .... etc etc.

错误信息是说你的向量长度 3 没有完全进入 b$pos 的长度。