2014-04-04 41 views
0

我是r的新手用戶。我用圖書館的通告繪製了我的物候數據(Fejeuncirculangle)。我獲得了一個傳統的圖形,但是我不可能在這個圖中顯示平均箭頭。我不明白爲什麼,因爲我很容易在曲線圖中爲每個初始值,段和文本添加箭頭。以下是我的代碼。這是代碼問題嗎?r圓形圖和平均值

作爲一項進步,非常感謝您的回答。

Mantalania

代碼

x11() 
plot.new() 
library(circular) 
na.exclude(tinitialcirculangle) 
plot(tinitialcirculangle$Fejeuncirculangle,pch=20,col="blue",stack=TRUE,shrink=1.2,bins=52,ticks=TRUE) 
Fejeuncirculangle<-circular(tinitialcirculangle$Fejeuncirculangle,type="angles",units="degrees",template="none",zero=circular(pi/2),rotation="clock",names="Feuiljeuncirculangle") 
par(mai=c(0.5,0.75,0.5,0.75),cex.axis=1,cex.lab=1) 
rose.diag(Fejeuncirculangle,bins=52,main="",col="darkolivegreen1",axes=FALSE) 
points(Fejeuncirculangle,pch=20,col="darkolivegreen1",bg="black",next.points=+0.05) 
arrows.circular(Fejeuncirculangle) 
arrows.circular(mean(Fejeuncirculangle),angle=15,length=0.20,lwd=5,col="black") 

# Month Segments - Ajouter les segments de séparation des mois 
segments(0,0,1,0,col="black",lty=3,lwd=1) 
segments(0,0,sqrt(3)/2,0.5,col="black",lty=3,lwd=1) 
segments(0,0,0.5,sqrt(3)/2,col="black",lty=3,lwd=1) 
segments(0,0,0,1,col="black",lty=3,lwd=1) 
segments(0,0,-0.5,sqrt(3)/2,col="black",lty=3,lwd=1) 
segments(0,0,-sqrt(3)/2,0.5,col="black",lty=3,lwd=1) 
segments(0,0,-1,0,col="black",lty=3,lwd=1) 
segments(0,0,-sqrt(3)/2,-0.5,col="black",lty=3,lwd=1) 
segments(0,0,-0.5,-sqrt(3)/2,col="black",lty=3,lwd=1) 
segments(0,0,0,-1,col="black",lty=3,lwd=1) 
segments(0,0,0.5,-sqrt(3)/2,col="black",lty=3,lwd=1) 
segments(0,0,sqrt(3)/2,-0.5,col="black",lty=3,lwd=1) 

# Month texts - Ajouter les mois au graphique circulaire 
text(0.2329,0.8693,"J",family="Calibri",font=1,cex=1) 
text(0.6363,0.6363,"F",family="Calibri",font=1,cex=1) 
text(0.8693,0.2329,"M",family="Calibri",font=1,cex=1) 
text(0.8693,-0.2329,"A",family="Calibri",font=1,cex=1) 
text(0.6363,-0.6363,"M",family="Calibri",font=1,cex=1) 
text(0.2329,-0.8693,"J",family="Calibri",font=1,cex=1) 
text(-0.2329,-0.8693,"J",family="Calibri",font=1,cex=1) 
text(-0.6363,-0.6363,"A",family="Calibri",font=1,cex=1) 
text(-0.8693,-0.2329,"S",family="Calibri",font=1,cex=1) 
text(-0.8693,0.2329,"O",family="Calibri",font=1,cex=1) 
text(-0.6363,0.6363,"N",family="Calibri",font=1,cex=1) 
text(-0.2329,0.8693,"D",family="Calibri",font=1,cex=1) 
+0

其實m ean(Fejeuncirculangle)返回NA。但我的價值是數字格式(is.numeric返回TRUE),我不明白爲什麼不可能有意思。 – user3370470

回答

0

嘗試:

arrows.circular(mean(Fejeuncirculangle,na.rm=TRUE),angle=15,length=0.20,lwd=5,col="black") 

如果這樣的作品,那是因爲你還沒有排除NA值,我懷疑你的第四行應該是:

tinitialcirculangle <- na.exclude(tinitialcirculangle) 
+0

非常感謝!問題已解決! – user3370470