График R не отображается полностью с использованием кругового пакета R

Question

У меня есть круговые данные направления ветра в радианах. Здесь у нас часть ноябрьских направлений на 10 метров:

 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 имеет похожий сюжет, но в градусах. Есть ли решение?

Answer 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) 

дает мне меньшие текстовые метки, которые подходят лучше. Но все еще немного некрасиво, и текст может конфликтовать с кружком и отметками.

Answer 2

Вероятно, это связано с отсечкой по региону. Нижняя часть 2 идет в область полей и поэтому по умолчанию обрезается. Поскольку у вас нет декартовых осей на графике, не помешает тексту выйти за пределы поля. Попробуйте выполнить следующую команду перед созданием графика:

par(xpd=TRUE)

Это изменит отсечение (для всех будущих графиков на этом устройстве). Увидеть ?par и прокрутите вниз до раздела xpd Больше подробностей.