R円 - 分布の中心 - 平均

Question

私は月ごとの出現数を持つデータセットを持っており、分布/平均月の中心を計算したいと思います。可能であれば、私も信頼区間を持っていたいと思います。

私は円形とCircStatsのマニュアルを読んでおり、同様の質問とhereを見てきました。

私はあるケースでは一見賢明な結果を得ることができましたが、他のケースでは妥当な結果を得ることができず、信頼区間を計算する方法をまだ考えていませんでした。ここでは、私のポイントを説明するために

は、いくつかのダミーデータがある：観測の最初のセットについて

library(CircStats) 

# The number of observations by month (Jan-Dec): 
obsMonths1 <- c(12,15,1,2,3,1,1,4,1,2,7,1) 
obsMonths2 <- c(1,1,1,1,2,10,11,2,1,1,2,1) 

# Convert data to radians: 
obsRadians1 <- (obsMonths1/12*2)*pi 
obsRadians2 <-(obsMonths2/12*2)*pi 

# Calculate circular mean: 
mean1 <- circ.mean(obsRadians1-1)#assume January is 0 
mean2 <- circ.mean(obsRadians2-1)#assume January is 0 

# Convert radians to months: 
mean1*12/(2*pi)+12 
mean2*12/(2*pi)+12 

答えが賢明なようだが、観測値の第二セットのためには、7〜8月でなければなりません。

Answer 1

あなたのコードに問題はないようです。あなたは "信頼区間"と尋ねましたが、統計テストまたはブートストラップ手順のようには見えません。私は最初にこれらの結果がなくても信頼区間をどのように計算するかはわかりません。あなたが正しいことをしたことが確かであれば、数学/統計特有のStackExchangeコミュニティでこの質問を提起したいかもしれません。