2017-02-05 11 views
1

나는 출현 횟수가 달 단위 인 데이터 세트가 있으며 분포/평균 월의 중심을 계산하려고합니다. 가능하다면 신뢰 구간도 갖고 싶습니다.R 원형 - 분포 중심 - 평균

나는 circular과 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 월이어야한다.

답변

0

코드에 문제가없는 것처럼 보입니다. 당신은 "신뢰 구간"을 물었지만, 통계 테스트 또는 부트 스트랩 절차를 수행 한 것처럼 보이지 않습니다. 나는이 결과가 없으면 어떻게 신뢰 구간을 계산할 것인지 잘 모르겠습니다. 옳은 일을했다면 수학/통계 관련 StackExchange 커뮤니티에서이 질문을 제기 할 수 있습니다.

+0

부트 스트랩 제안에 감사드립니다. 코드가 오류없이 실행되지만 두 번째 더미 데이터 집합에 대한 대답이 올바르지 않다는 것에 동의합니다. 내가 어디로 잘못 가고 있는지에 대한 통찰력? –

+0

나는 왜 당신이 정확한 결과를 얻지 못했는지 완전히 확신하지 못한다. 날짜/연도 변환이 오류를 발생시키기 쉽기 때문에'obsRadians2'와'mean2'가 적절히 정의되었는지 확인하는 것이 좋습니다. –