我有以弧度为单位的风向圆形数据。在这里,我们在 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 包有类似的情节,但以度为单位。有解决办法吗?