Участок угловых частот с использованием графика роз, R

Question

Я пытаюсь составить диаграмму / график углов / частот для роз, но не знаю точно, как этого добиться, используя R. В моих данных также можно изобразить углы в виде круговой гистограммы, но они на самом деле не представляют частоты. У меня есть другой столбец для частот, которые я не знаю, как использовать. Также я заметил, что метки не совпадают с центром столбцов на графике.

Вот образец моих данных:

library(circular)

data<-read.table(text="phase  freq
                       315  3
                       0    0
                       315  4
                       0    0
                       0    13
                       45   7
                       45   4
                       45   0
                       45   0
                       90   0
                       90   0
                       90   2
                       135  2
                       135  71
                       135  59
                       135  0
                       180  1
                       180  11
                       180  4
                       180  77
                       225  0
                       0    1
                       0    1
                       45   5
                       45   13
                       45   1
                       45   0
                       90   0
                       90   0
                       90   0
                       135  7
                       135  85
                       135  4
                       135  13
                       180  0
                       180  0
                       180  0
                       180  22
                       225  0
                       225  15
                       225  0
                       90   2
                       90   0
                       90   0
                       135  0
                       135  0
                       135  0
                       135  0
                       180  0
                       180  0
                       180  0
                       180  0
                       225  104
                       225  7
                       225  17
                       225  14
                       270  0
                       270  0
                       270  1
                       315  0
                       315  0
                       315  8
                       315  1
                       0    11
                       0    0
                       0    0
                       0    0
                       45   0
                       45   0
                       45   0
                       45   0
                       90   0
                       90   0
                       90   0
                       135  12
                       135  1
                       135  7
                       135  18
                       180  48
                       180  0
                       180  0
                       180  73
                       225  13
                       225  26
                       225  0
                       225  12
                       270  0
                       270  0
                       270  57
                       315  4
                       315  0
                       315  0
                       315  8
                       0    0
                       0    0
                       0    0
                       45   0
                       45   0
                       45   0
                       ",header=T)

head(data)

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

x<-circular(data$phase,type="angles",units="degrees",rotation="clock",zero=0)

rose.diag(x, bins = 16,shrink=1.5,prop=2,col="gray70",axes=F,zero=c(rad(90)),
              ticks=TRUE)
points(mean(x), col = "red",cex=2)
text(0,1.2,"Full", cex=1); text(0,-1.2,"New moon",cex=1); 
text(1.5,0,"Third quater",cex=1); text(-1.5,0,"First quater",cex=1)

Любые предложения будут оценены.