问题标签 [plotrix]

请考虑以下示例极坐标图：

我想添加角度为 30 和 330 以及 150 和 210（从中心到外部）的线。我尝试了 line 功能，但无法让它工作。

我有一个类似于以下格式的时间序列面板数据集，其中 V1 是字符（这里是 A、B、C...），V2 是数字（5,8,6...）。我在 V1 中有 65 个不同的项目（总共大约 50000 个奇怪的观察值）。我正在尝试找到一种方法来绘制 V2 提供的 V1=A 或 B 等等。我可以简单地做的是plot(V1[1:1065])，plot([V1:1066:2085])等等，但我一直在寻找一种强大的方法来做到这一点。尝试过这样的事情：

但由于 R 中的 'if' 不接受向量，它会给出错误消息，"the condition has length > 1 and only the first element will be used"

有人可以提出一些建议吗？

这种创建人口金字塔的方法曾经对我有用，现在不知何故冻结了我的 RStudio，所以我每次都需要重新启动它，我没有进步......

如果有人能让我知道是什么让 R 停止，将不胜感激......

如果我将以下内容添加到代码中，它会起作用：laxlab=seq(0,1e7,by=1e6),raxlab=seq(0,1e7,by=1e6)) 如这篇文章中所述： R pyramid plot very slow if fed更大的数字

我正在尝试在泰勒图中绘制多个模型，并且在代码中略有挣扎。我设法制作了图表（见图），但无法弄清楚如何减少轴，因为它们太大了，标准化轴以标记为 1,2,3,4 并在相关性上添加刻度线 - 带有刻度标记我希望每 0.1 有一次主要刻度线，每 0.05 次有次要刻度线，直到 0.9，之后我试图在 0.95 有一次主要刻度线，此时每 0.01 次有次要刻度线（如果这有意义的话）。上述任何帮助/建议都会有所帮助。我在“plotrix”包中使用了“taylor.diagram”（并阅读了它的指南——但我对 R 相对缺乏经验）并附上了我的（有点基本的）代码，但我的情节看起来相当混乱。谢谢

我最近发现了 JESS3 对肥胖和饥饿之间联系的可视化：http: //jess3.com/obesity-and-hunger-infographic/

由于我喜欢信息/空间比，我试图在 R 中绘制一个类似的图表（即一个在中心有散点的极坐标图，在中心的末端，一个堆叠的条形图从 0 开始的新轴到x.但是，我没有找到一种方法来实现plotrix或ggplot2中不同极坐标图版本的堆叠。

为了说明，在一个非常基本的版本中，两个单独的图表应该是什么样子，我用 mtcars 包创建了一个简单的例子：

是否有可能，如果是这样，我怎样才能实现新轴开始并将条形添加到散点图之外？

非常感谢您，非常感谢所有提示！

我成功地创建了一个脚本，它在 RStudio 中产生了一个漂亮的情节：

abscissa = nrow(LoopVariable) library(plotrix) plotCI(1:abscissa, LoopVariable[ ,13], uiw=LoopVariable[ ,"F_corrected_normed_error"], err="y", main=(LoopVariable[1,"Samples..Sample_ID"]), xlab="", ylab="correct Y axis", pch=19) abline(h=LoopVariable[1,"mean_fraction-modern"], lwd=2, col="red") mtext("no additional error added")

到目前为止，一切都很好。然后我写了这段代码：

png(filename=paste0("I:/C14Data/QualityControl/Johannes/code_output/", RNumberList[i], "_", date, " (R#_yyyy-mm-dd_hh-mm-ss).png"), type="cairo", units="cm", width=abscissa, height=10, pointsize=12, res=96) plotCI(1:abscissa, LoopVariable[ ,13], uiw=LoopVariable[ ,"F_corrected_normed_error"], err="y", main=(LoopVariable[1,"Samples..Sample_ID"]), xlab="", ylab="fraction modern", pch=19) dev.off() 这就像我想要的那样导出一个 png 文件，BUT：如何将以下两行添加到导出函数并让它看起来与前面的代码完全相同？： abline(h=LoopVariable[1,"mean_fraction-modern"], lwd=2, col="red") mtext("no additional error added")

我正在尝试制作一个“相反的堆积条形图”，并且pyramid.plot从plotrix包装中发现似乎可以完成这项工作。（我很感激ggplot2将成为你们中的一些人的首选解决方案，但我希望base在这个上坚持使用图形。）

不幸的是，当我尝试将限制设置为非整数值时，它似乎对 x 轴做了一件奇怪的事情。如果我让它自动定义限制，它们是整数，在我的情况下，它只会留下太多的空白。但是将它们定义为xlim=c(1.5,1.5)会产生下面的奇怪结果。

如果我从文档中正确理解，则无法传递其他图形参数以例如抑制轴并稍后添加它，或者更不用说定义刻度点等。有没有办法让它更灵活？

这是一个用于生成下图的最小工作示例。

以防万一其他人感兴趣，我不是在做人口金字塔，而是尝试使用一些值为负数的堆积条形图。上面的代码包含一个“技巧”，我使用它可以使每边有不同数量的条形集，即向矩阵添加空列，希望有人会发现这很有用 - 很抱歉，工作示例不是本来可以最小的！

我有以下数据，

SampleID Pos Dep Pvalues
sample_1 849 62 0.02755358
sample_1 859 63 0.07406833
sample_1 864 63 0.00351564
sample_1 883 60 0.02780868
sample_1 893 58 0.00451450
sample_1 895 58 0.03600795
sample_2 54 66 0.11864407
sample_2 55 67 0.01515152
sample_2 71 91 0.02712367
sample_2 78 97 0.00077325

我已经生成了 P 值的直方图，每个条形图的顶部都有频率值。下面，是代码

pval_at_site <- read.table("samples.pval")
s <- hist(pval_at_site$Pvalues, xlab="Pval",cex=0.8)
text(s$mids,s$counts,s$count,srt=90, pos = 3,offset=1,cex=0.6)

现在，我想做的是，连同 P 值频率，在每个条形顶部添加样本数。

例如，如果我在第一个间隔中有 1000 个数据点，并且这些值来自 20 个唯一样本，我希望我的绘图在第一个条形顶部显示“1000,20”。

请让我知道我应该怎么做。希望我已经说清楚了。

谢谢。

我一直在尝试使用 plotrix 在极坐标图中创建误差线。虽然散点函数没有问题，但我无法添加误差线，因为当我添加例如具有弧度模式的第二个极坐标图的线时，我无法将它们从绘图中心移开以从散点开始点...我最好的猜测是在顶部添加第二个极坐标图，散点图点的弧度+误差线，然后在其顶部添加一条白线，散点图点 - 误差线，覆盖误差线到下限...这显然不仅在代码中非常难看，而且在图表中也是如此，所以我想知道是否有更好的解决方案。到目前为止，我的代码（使用 mtcars 作为一个无意义的示例）如下所示：

我想做一个带有自定义 x 轴的堆积面积图。具体来说，我想用 cex.axis=0.8 标记 x 轴上的刻度线（使其更小以适合所有值），并且我想使用红色突出显示一些 x 值（例如，周末日期与工作日形成对比）。

plotrix 包中的 stackpoly 在创建绘图方面做得很好。但我无法控制 cex.axis 或抑制 x 轴的绘图。

这是我完成这项工作的两次尝试：