2

我有以弧度为单位的风向圆形数据。在这里,我们在 10 米处有 11 月的部分方向:

 nov<-c(5.4977871, 5.4977871, 5.4977871 ,5.4977871, 5.4977871 ,5.4977871, 5.4977871, 5.1050881, 5.4977871, 5.4977871, 5.8904862, 5.8904862 ,6.2831853,5.8904862, 5.8904862, 5.8904862, 6.2831853, 5.8904862, 5.8904862, 5.8904862, 5.8904862, 5.8904862, 5.8904862, 5.8904862, 5.8904862, 5.8904862)

我已经为风向绘制了时间序列和玫瑰图。这是代码:

library(circular)
par(mfrow = c(2, 1))  
ts.plot(nov, xlab="Hour", ylab="radians",main="Time series of wind direction") 
rose.diag(nov[-c(100:713)],bins=18, main="Rose Diagram",) 

在此处输入图像描述

如您所见,底部有问题,在 3pi/2 中。“2”未完全显示。我已经完成了基于 ggplot2 的玫瑰图,但我需要保存这种图。CirStats 包有类似的情节,但以度为单位。有解决办法吗?

4

2 回答 2

1

玫瑰图上的剪裁很大程度上取决于我的图形窗口有多大。如果我缩小我的默认大小的图形窗口,我可以让我的情节像你的一样。如果我使用更大的窗口,我可以让它看起来不错。

这个情节的最终目的地在哪里?如果是 PDF,那么您在该表单中可能没有问题。

使用参数更改字符大小cex可能会有所帮助:

par(mfrow = c(2, 1))  
ts.plot(nov, xlab="Hour", ylab="radians",main="Time series of wind direction") 
rose.diag(nov[-c(100:713)],bins=18, main="Rose Diagram",cex=0.5) 

给我更适合的更小的文本标签。但它仍然有点难看,文本可能与圆圈和刻度线发生冲突。

于 2014-08-29T16:46:41.333 回答
0

这可能是由于裁剪到绘图区域。2 的底部进入边缘区域,因此默认情况下被剪裁。由于绘图上没有笛卡尔坐标轴,因此允许文本延伸到页边距不会有什么坏处。在创建绘图之前尝试运行以下命令:

par(xpd=TRUE)

这将更改剪辑(对于该设备上的所有未来绘图)。有关更多详细信息,请参阅?par并向下滚动到有关部分xpd

于 2014-08-29T16:57:26.577 回答