박스 플롯에 (Mean +/- 2 Sigma) 라인 추가하기

Question

R에서 boxplot을 만들고 싶습니다. 간단하게 (여기서는 샘플 데이터를 사용했습니다) R 콘솔에서 다음 명령을

x <- rnorm(100); 

y <- rnorm(100); 

boxplot(x, y, col=rainbow(2)) 

지금 내가 정말 원하는 (2 (+/- 평균 3) * 시그마)를 보여주는 것입니다 내 데이터에 대한 표시는 위에 설정 boxplot 위. 나는 표준 편차를 알고 있으며, 평균은 강력한 통계는 아니지만, 나의 정보 세트는 그 정보가 유용 할 수있다. 나는. 의도 된 그림은 다음과 같아야합니다 :

여기서 점선은 해당 플롯의 95 % 또는 99.7 % 신뢰 구간 또는 (평균 +/- (2 또는 3) * 시그마)를 나타냅니다.

N.B. 불행히도 여기에 수치를 추가 할 수는 없으므로 충분한 평판이 없으며 누군가 위의 링크에서 업로드 할 수 있다면 감사 할 것입니다.

Answer 1

매우 간단하고 쉬운 방법은 인수를 boxplot() 함수 (?boxplot 참조)에 사용하는 것입니다. 이것은 평균에 대한 95 % 신뢰 구간과 동일하지는 않지만 대신 중앙값에 대한 95 % CI와 유사합니다 (견고성에 대한 우려를 나타냄). 설명되어 있습니다 here :

노치 (요청 된 경우)는 +/- 1.58 IQR/sqrt (n)까지 확장됩니다. 이것은 에서 1.57의 공식과 동일한 계산을 기반으로하는 으로 간주됩니다. Chambers et al. (1983, p. 62), McGill et al. (1978, 16 쪽). 이들은 중앙값의 점근 적 정규성을 기반으로하며, 비교되는 두 중앙값에 대해 대략 동일한 샘플 크기를 기반으로하며, 샘플의 기본 분포에 다소 민감하지 않다고합니다. 아이디어는 두 중간 값의 차이 인 의 대략 95 % 신뢰 구간을 제공하는 것으로 보입니다.