「円」の意味R

Question

月のデータセットを考えると、月が円形であることを考慮して、「平均」月を計算するにはどうすればよいですか？

months = c(1,1,1,2,3,5,7,9,11,12,12,12) 
mean(months) 
## [1] 6.333333 

このダミーの例では、平均値は1月または12月です。私は円形統計のパッケージがあることを知っていますが、ここで私のニーズに合っているかどうかはわかりません。

Answer 1

私は

months <- c(1,1,1,2,3,5,7,9,11,12,12,12) 
library("CircStats") 
conv <- 2*pi/12 ## months -> radians 

は今、数ヶ月からラジアンに変換する円弧平均を計算し、数ヶ月に戻って変換すると思います。私は結果は-0.3457である

(res1 <- circ.mean(conv*(months-1))/conv) 

... 1月/ 12時「0ラジアン」であると仮定し、ここで1を引いています。あなたが必要になることがあります。

11.65を与える

(res1 + 12) %% 12 

、すなわち途中12月を通じて（私たちは0 = 1月、11 = 12月規模で残っているので）

を私はこれが正しいと思うが、持っていないあまりにも注意深くそれをチェックした何が価値があるために

、CircStats::circ.mean機能は非常に簡単です - それは、これはあなたが必要とするすべてである場合、パッケージをロードするのオーバーヘッドの価値ではないかもしれません：

からA.Webbの巧妙な代替@取り入れ

function (x) 
{ 
    sinr <- sum(sin(x)) 
    cosr <- sum(cos(x)) 
    circmean <- atan2(sinr, cosr) 
    circmean 
} 

コメント：

m <- mean(exp(conv*(months-1)*1i)) 
12+Arg(m)/conv%%12 ## 'direction', i.e. average month 
Mod(m)    ## 'intensity'