나는 출현 횟수가 달 단위 인 데이터 세트가 있으며 분포/평균 월의 중심을 계산하려고합니다. 가능하다면 신뢰 구간도 갖고 싶습니다.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 월이어야한다.
부트 스트랩 제안에 감사드립니다. 코드가 오류없이 실행되지만 두 번째 더미 데이터 집합에 대한 대답이 올바르지 않다는 것에 동의합니다. 내가 어디로 잘못 가고 있는지에 대한 통찰력? –
나는 왜 당신이 정확한 결과를 얻지 못했는지 완전히 확신하지 못한다. 날짜/연도 변환이 오류를 발생시키기 쉽기 때문에'obsRadians2'와'mean2'가 적절히 정의되었는지 확인하는 것이 좋습니다. –