Создание полуполярного графика (диаграмма-роза) с круговым пакетом

Я строю распределение углов с помощью функции rose.diag из библиотеки circular. Входными данными являются радианты. Мой код MWE

    library(circular);
    dat<-read.csv(file.choose(),header=F);
    data=unlist(dat);
    rose.diag(data, bins=24)

и я получаю этот график: введите здесь описание изображения

Мне интересно показать только часть данных, от -pi/2 до pi/2, и чтобы длина самого большого тика равнялась длине радиуса круга, как показано здесь: введите здесь описание изображения

Любая помощь будет оценена по достоинству!

ИЗМЕНИТЬ

Как предложил #lawyeR, вот код с образцом данных:

    library(circular);
    data<- c(-0.188,-0.742,-0.953,-0.948,-0.953,-1.187,-0.9327200,-0.855,-  0.024,1.303,-1.041,-1.068,-1.066,1.442,1.150,0.965,0.665,0.649,0.984,-1.379,-0.584,-0.573,-0.357,-0.237,-0.287,-0.486,-0.783,-0.298,0.849,1.088,-1.003,-0.952,-0.776,-0.811,-0.880);
    rose.diag(data, bins=24);

r polar-coordinates
person michltm    schedule 16.10.2015    source источник
comment
Почему бы вам не отправить свои данные, чтобы люди могли помочь вам лучше?   -  person lawyeR    schedule 16.10.2015
comment
Спасибо за ваш ответ. Я отредактировал сообщение, чтобы добавить образец данных.   -  person michltm    schedule 16.10.2015
comment
Хорошо. Я отредактировал ваш вопрос и добавил еще один тег, надеюсь, чтобы получить ответ или два.   -  person lawyeR    schedule 16.10.2015

Ответы (1)


arrow_upward
4
arrow_downward

Возможно, вы можете сначала нарисовать полукруг без распределения, используя функции графика по умолчанию. Затем заполните распределение без его круга:

library(circular)

data<- c(-0.188,-0.742,-0.953,-0.948,-0.953,-1.187,-0.9327200,-0.855,-  0.024,1.303,-1.041,-1.068,-1.066,1.442,1.150,0.965,0.665,0.649,0.984,-1.379,-0.584,-0.573,-0.357,-0.237,-0.287,-0.486,-0.783,-0.298,0.849,1.088,-1.003,-0.952,-0.776,-0.811,-0.880)

freq <- diff(colSums(outer( data %% (2*pi), (1:24)*pi/12,"<"))) / length(data)
r.max <- sqrt(max(freq))

#-----------------------------------------------------------------
# Plot the half circle:

lab.width  <- 0.15*r.max
lab.height <- 0.15*r.max

plot( c(-r.max,r.max), c(0,0),
      axes=FALSE,
      ylim=c(0,r.max+lab.height),
      xlim=c(-r.max-lab.width,r.max+lab.width),
      xlab="", ylab="", type="l")

for ( i in 0:(5*12-1) )
{
  psi <- i*pi/(5*12)
  x1 <- r.max*cos(psi)
  y1 <- r.max*sin(psi)
  x2 <- r.max*cos(psi+pi/(5*12))
  y2 <- r.max*sin(psi+pi/(5*12))

  lines( c(x1,x2), c(y1,y2), type="l")

  if (i %% 5 == 0) { lines( x1*c(1,0.95), y1*c(1,0.95), type="l" ) }
}

par(cex=2.0)

text( x = c(-r.max,0,r.max),
      y = c(0,r.max,0),
      labels = c("-pi/2","0","pi/2"),
      pos = c(2,3,4))

#------------------------------------------------------------
# Plot the distribution, but without the circle:

rose.diag(data,
          bins = 24,
          rotation = "clock",
          tcl.text = NA,
          ticks = FALSE,
          zero = pi/2,
          control.circle = circle.control( col="white" ),
          add = TRUE )

. введите здесь описание изображения

person mra68    schedule 16.10.2015
comment
Это отличная идея! Спасибо за вашу помощь. У меня был бы последний вопрос: как лучше всего добавить -pi/2, 0, pi/2 во внешний круг? - person michltm; 19.10.2015
comment
Я отредактировал свой ответ. Вы можете использовать text для добавления ярлыков. - person mra68; 20.10.2015
comment
Большое спасибо за вашу драгоценную помощь, все будет хорошо! - person michltm; 21.10.2015