我正在尝试创建一个玫瑰图,显示每个细胞子集的平均轨迹角度和距离。我希望玫瑰图周围的角度是轨迹角度,图中条形的长度是总位移。
这是每组平均角度和位移的测试数据集。
testsum<-data.frame(Group=c(1,2,3),
angle=c(0.78,1.04,2.094),
displacement=c(1.5,2,1))
当我尝试用循环方法绘制它时,我的图表看起来非常错误。
p1<-ggplot(testsum, aes(x=angle,y=displacement))+
coord_polar(theta="x",start=0)+
geom_bar(stat="identity",aes(color=Group,fill=Group),width=.01)+
scale_x_continuous(breaks=seq(0,360,60))
它为我提供了输出图。
根据数据所说的内容,它应该看起来更像这样(绘制预期的输出)。
似乎错误地放置了角度?知道我做错了吗?
答案 0 :(得分:2)
尽管MLavoie在20分钟内“击败了我”,但我认为可以使用NISTunits
包添加一些可读性:
library(ggplot2)
library(NISTunits)
testsum <- data.frame(
Group = c(1, 2, 3),
angle = c(0.78, 1.04, 2.094),
displacement = c(1.5, 2, 1)
)
testsum$angle = NISTradianTOdeg(testsum$angle)
ggplot(testsum, aes(x = angle, y = displacement)) +
coord_polar(theta = "x", start = NISTdegTOradian(-90), direction = 1) +
geom_bar(stat = "identity",
aes(color = Group, fill = Group),
width = 1) +
scale_x_continuous(breaks = seq(0, 360, 10), limits = c(0, 360))
结果:
要剪切下半部分,请检查this answer.
答案 1 :(得分:1)
也许你可以在露天包中使用windRose功能
library(openair)
testsum<-data.frame(Group=c(1,2,3),angle=c(0.78,1.04,2.094),displacement=c(1.5,2,1))
testsum$angle=180*testsum$angle/pi
windRose(testsum,ws="displacement",wd="angle",breaks=c(0.5,1,1.5,2),paddle=F,key.header = "My Rose",angle=10,
statistic="prop.mean",key.footer = "Displacement")
答案 2 :(得分:1)
也许你可以试试这个:
testsum$angle_b=180*testsum$angle/pi
#
ggplot(testsum, aes(x=angle_b,y=displacement))+
geom_bar(stat="identity",aes(color=Group,fill=Group),width=1) +
scale_x_continuous(breaks=seq(0,360,10), limits=c(0,360)) + coord_polar(direction=1)