비운영 참조

평균으로 당기기

계층적 경험 베이즈 축소

한 줄 결론 — 적은 표본의 극단값을 공통 평균 쪽으로 완화하는 후보다. 채택되지 못한 채 2026-08-19 결정으로 비운영 참조가 됐다.

아래 해설의 실측 수치는 연구 구간 기록(M1–M3 · 2026-08-15 동결) 기준이다. 오늘 수치는 위 표가 매일 갱신한다.

Q. 이 모델은 지금 무엇을 계산하나?

관측이 적은 노선·정류장·시간 집단에서도 출발 뒤 만석 확률 추정이 한두 사건에 과도하게 출렁이지 않게 하는 문제를 다룬다.

Q. 어떻게 계산하나?

이 방식의 강점

• 표본이 적은 집단의 관측 비율이 한두 사건 때문에 극단으로 치우치는 현상을 줄인다.
• 모든 집단을 하나로 합치지 않으면서 공통 자료를 나눠 쓴다.
• 확률이 평균 쪽으로 이동한 이유를 집단 표본과 사전분포로 설명할 수 있다.

Q. 만석을 한 번도 못 본 셀은 위험 0%인가?

아니다. 코퍼스 학습분의 노선×정류장×30분 셀은 833개다. 이 중 281개(33.73%)는 표본이 10건 이하다. 719개(86.31%)는 만석을 한 번도 못 봤다. 학습 만석률 자체가 705/12,065=5.84%라서, 표본이 몇 건뿐이면 0회 관측이 흔하다.

위치와 시간은 버리기 아까운 신호다. 아무것도 안 본 기준선 AUC는 .5000, 앞차 한 대만 본 기준선은 .7187, 정류장·30분 셀 기준선은 .9613이었다. 문제는 셀을 쓰느냐가 아니다. 얇은 셀을 얼마나 믿느냐다.

무축소 셀 비율은 만석 0건인 셀에 확률 0을 냈다. 검증 log loss는 .1318이었다. 노선 평균에 해당하는 의사표본 1개를 섞자 .1033으로 줄었다. 이 진단 때문에 축소추정이 본 후보로 올라왔다.

비유하면 학생 한 명만 시험을 본 반의 평균을 그 한 점수로 확정하지 않는 것과 같다. 학생이 적은 반은 학교 평균을 조금 섞는다. 학생이 많아지면 그 반의 점수를 그대로 믿는다.

Q. 어디로, 얼마나 당기나?

Q. 한 칸은 무엇인가?

한 셀은 노선×station_id×로컬 30분대다. 셀 안에서 출발 후 잔여석이 0인 사건 수 m과 전체 사건 수 n을 센다. 도착 시점 0석은 쓰지 않는다. 하차 뒤 자리가 생길 수 있어서다. remainSeatCnt=-1도 결측이라 뺀다.

Q. 당김의 기준은 무엇인가?

먼저 전체 학습분에서 전역 만석률을 구한다. 만석과 비만석에 .5건씩 더하는 Jeffreys 보정을 쓴다. 노선 만석률은 전역 만석률 쪽으로 의사표본 4개만큼 당긴다. 마지막으로 각 셀을 그 노선 만석률 쪽으로 당긴다.

셀 예측식은 하나다.

α는 노선 평균이 몇 건의 표본처럼 작용하는지 정한다. 실제 후보는 .5, 1, 2, 4, 8, 16, 32, 64를 날짜 순방향 교차검증으로 겨뤘다. 선택값은 .5였다. 노선 평균을 반 건만 섞는 약한 축소다.

Q. 왜 ‘경험적’ 베이즈인가?

사전평균 p_route를 사람이 정하지 않고 같은 학습 데이터의 노선·전역 빈도에서 구하기 때문이다. beta-binomial 방식으로 노선 평균을 의사표본처럼 보고 셀의 만석/비만석 횟수와 합친다. 사후분포 전체가 아니라 사후평균 하나만 내보낸다.

Q. 처음 보는 셀은 어떻게 하나?

정확한 셀이 없으면 노선×정류장 → 노선×30분대 → 노선 → 전역 순서로 내려간다. 이미 본 셀은 표본이 1건이어도 이 폴백을 쓰지 않는다. 그 셀의 비율과 노선 평균만 섞는다.

출력은 만석 확률 p_full이다. 사용자에게 보이는 탈 확률은 별도 모델이 아니라 1-p_full이다.

Q. 실제 숫자는 어떻게 바뀌나?

1650번 태영아파트(station_id 209000243), 18:00~18:30 셀을 보자. 전부 학습 데이터에서 계산한 값이다.

1. 전체 학습분은 12,065건, 만석 705건이다. Jeffreys 보정을 거친 전역 만석률은 5.847%다.
2. 학습분의 1650번은 4,331건 중 298건이 만석이다. 전역 평균을 의사표본 4개만큼 섞은 노선 만석률은 6.880%다.
3. 이 셀은 5건 중 만석 0건이다. 무축소 값은 0%다.
4. 선택된 α=.5를 넣으면 (0 + .5×.068797) / (5+.5) = .006254다. 만석 확률 0.625%, 탈 확률 99.375%다.

검증 라벨은 이 계산에 넣지 않았다. 뒤 검증에서는 같은 셀 4건 중 1건이 만석이었다. 한 셀의 4건으로 모델 전체를 판정할 수는 없다. 다만 학습의 0/5가 위험 0이라는 뜻은 아니었다.

표본이 두꺼우면 덜 움직인다. 3330번 롯데백화점.범계역 07:00~07:30 셀은 학습 23건 중 14건이 만석이었다. 무축소 60.87%가 축소 뒤 59.69%가 됐다.

재료 한눈에

쓰는 재료는 네 가지다. 출처 필드까지 아래 표가 정본이다.

쓰는 재료 4개 — 출처 필드까지

지금은 안 쓰는 재료 5개

• 차량 정원·점유율 — TBBMSVEHINFOM.SEAT_RDNG_PSN_CNT; 점유율 후보식 1-msgBody.busLocationList[].remainSeatCnt/TBBMSVEHINFOM.SEAT_RDNG_PSN_CNT20일 코퍼스 차량키가 plateNo 우선 HMAC이라 외부 정원과 조인된 차량이 0/77이다. 기존 231,048행의 capacity_exact도 모두 0이다. 관측 최대 잔여석은 정원의 하한일 뿐 고정 정원이 아니다. 라벨은 출발 후 좌석이 0인지에 정원이 필요 없다. 프로토타입의 defaultSeatCapacity=80은 분포 배열 상한이며 점유율 분모가 아니다.
• 혼잡도 — GBIS getBusLocationListv2 msgBody.busLocationList[].crowded공식 제공 노선유형은 13·15·23이고 대상 노선은 11이다. 직접 수집분에서 같은 잔여석 값에 다른 코드가 붙은 값이 26개 나와 잔여석만의 전역 함수라는 해석은 반증됐다. 다만 같은 (vehId, remainSeatCnt)에서는 충돌이 0셀이고 180초 이하 인접 관측 81,147쌍에서 crowded만 바뀐 연속쌍도 0건이라 차량별 임계화일 가능성이 남는다. 산식은 미확정이고 아직 특징으로 써 보지 않았다.
• 현재 상류·지난 버스 상태 — msgBody.busLocationList[].remainSeatCnt, msgBody.busLocationList[].stateCd, msgBody.busLocationList[].stationSeq, msgBody.busLocationList[].vehId에서 만드는 S1 upstream_*, passed_*, previous_passage_*, route_seat_*이 후보는 역사적 시공간 위험만 비교하려고 실시간 스냅샷 재료를 읽지 않는다.
• 도착 ETA — GBIS getBusArrivalListv2 msgBody.busArrivalList[].predictTime1, msgBody.busArrivalList[].predictTime2, msgBody.busArrivalList[].predictTimeSec1, msgBody.busArrivalList[].predictTimeSec2; getBusArrivalItemv2 msgBody.busArrivalItem.predictTime1, msgBody.busArrivalItem.predictTime2, msgBody.busArrivalItem.predictTimeSec1, msgBody.busArrivalItem.predictTimeSec2공급자가 만든 예측값을 다시 설명변수로 배우는 구조가 된다. 분·초 필드 관계도 실측 6쌍 중 2쌍에서 맞지 않았다. v1 후보 입력에서 뺐다.
• API queryTime 파생 시간 — GBIS getBusLocationListv2 msgHeader.queryTime; HTTP 응답 헤더 Date응답 수신 시각과 어긋나는 반증 실측이 있어 감사 원문만 보존하고 나이·속도·헤드웨이·분할에 쓰지 않는다. 실측 수치는 평가 규약에 있다.

Q. 보정은 실제로 적용됐나?

이 모델 자체에는 보정층을 얹지 않는다. 축소로 만든 확률을 그대로 낸다. 확률 보정은 별도 후보 숫자 다듬기 · Platt 확률 보정이 맡는다. 보정 눈금이 뒤 날짜로 옮겨가지 않았던 실측은 아래 ‘언제 실패하나’에 있다.

Q. 언제 실패하나?

• 학습 날짜가 적거나 만석이 특정 날짜에 몰리면 공통 경향 자체가 불안정해진다.
• 실제로 다른 집단을 같은 사전분포로 묶으면 중요한 노선·정류장·시간 차이를 지나치게 줄일 수 있다.
• 운영 수집기 실패나 검열이 특정 구간에 몰리면 관측된 공통 평균이 운영 전체를 대표하지 못한다.
• 운행 정책, 노선, 정류장 순서, 수집 방식이 바뀌면 이전 공통 경향이 새 구간에 맞지 않을 수 있다.
• 처음 보는 집단의 확률은 공통 평균에 크게 의존하므로 해당 집단의 직접 성능 근거로 쓰면 안 된다.
• 평가 구간에 만석 또는 비만석이 없으면 순위를 가르는 일부 지표는 정의되지 않는다.

관측이 적은 구간은 공통 평균 쪽으로 더 많이 당겨진다. 그 구간의 현장 차이가 작다는 증거는 아니다. 관측이 많은 구간은 자체 비율이 더 반영되지만, 그 자체가 뒤 날짜 검증을 통과했다는 뜻은 아니다.

그 밖의 약점도 남는다.

• 정확한 셀이 한 번이라도 관측되면 정류장 평균이나 시간 평균을 더 섞지 않는다. 표본 1건 셀도 노선 평균과만 섞는다.
• 이웃 정류장과 이웃 시간의 연속성을 쓰지 않는다. 정류장 순번을 ±1 옮겨도 AUC .9100/.9113이 남았지만 이 정보는 버린다.
• 모든 노선·셀에 같은 α 하나를 쓴다. 셀마다 다른 흔들림을 따로 배우지 않는다.
• 새 셀 성능은 아직 재지 못했다. 코퍼스 validation의 새 셀은 2행·만석 0건뿐이다.
• 직접 수집분에서는 AUC 점추정이 7종 중 가장 낮았다. train 1,081건·만석 40건, validation 443건·만석 14건에서 ECE .0328, AUC .8527이고 네 후보에 유의하게 밀렸다.

Q. 축소했는데 왜 보정은 나빠졌나?

S1 층화 기준선은 의사표본 1개를 섞었다. 경험적 베이즈가 고른 α=.5는 그보다 약한 축소다. 날짜 순방향 2-fold에서 ECE·Brier·log loss의 합이 .188043으로, α=1의 .188211보다 아주 조금 낮아 선택됐다. ECE만 고른 결과가 아니다.

ECE·Brier·log loss는 낮을수록 좋고 AUC는 높을수록 좋다.

log loss는 .0009 좋아졌지만 ECE는 .0007 나빠졌다. 0% 같은 극단값을 피한 이득과 전체 확률 눈금의 작은 편향이 서로 다른 지표에 잡혔다.

후처리는 더 분명히 실패했다. OOF는 각 날짜를 그보다 앞선 날짜로만 예측해 만든 값이다. 학습 OOF 10,236건·만석 572건에서는 ECE가 raw .0128에서 Platt .0064, isotonic .0016으로 좋아졌다. 하지만 뒤 검증 10,040건·만석 583건에서는 raw .0070이 Platt .0134, isotonic .0131로 나빠졌다. 평균 예측도 raw 5.4829%에서 4.4684%·4.4963%로 내려가 관측 5.8068%를 크게 밑돌았다. 앞 구간에서 배운 눈금이 뒤 날짜로 옮겨가지 않은 것이다. 날짜별 분포 이동, 보정 표본 구성, isotonic의 60개 계단이 각각 얼마나 영향을 줬는지 분리한 실측은 아직 없다.

Q. 프로토타입과 무엇이 다른가?

• 프로토타입의 화면 탑승률은 정류장×평일/주말×30분에서 마지막 스냅샷의 비만석 비율 1-z/n이었다. 이 모델의 한 행은 버스 한 대의 정류장 통과 사건이고, 정답은 출발 후 0석이다.
• 프로토타입은 비율을 그대로 쓰고 n<10이면 추천을 멈췄다. 이 모델은 얇은 셀도 노선 평균과 섞고, 새 셀은 정류장·시간·노선·전역 순으로 폴백한다.
• 프로토타입의 화면 탑승률과 차량별 좌석 예보는 서로 다른 계산이었다. 이 모델은 만석 확률 하나만 내고 탈 확률은 1-p_full로 바꾼다. 남을 좌석 수·대기 인원·이동 추천은 계산하지 않는다.
• 프로토타입 좌석 예보는 현재 잔여석과 정류장별 수요분포를 흘렸다. 이 모델은 현재 좌석·앞차·대기열을 보지 않는 정류장·시간대 빈도 모형이다.
• 프로토타입 좌석 예보의 공통 상한 80석도 여기 없다. 이 모델은 정원·점유율·crowded를 모두 입력에서 뺐다.

Q. 개인정보와 미래 정답은 어떻게 막나?

• 공개 평가 파일에는 집계 표본, 공식 노선·정류장, 모델 성능만 남긴다.
• 차량 결합 키와 번호판 값은 공개 응답에 포함하지 않는다.
• 개별 차량의 시각·장소 이동 사례와 원본 응답은 공개하지 않는다.
• 공통 경향과 모델 설정은 평가일보다 앞선 자료로만 정한다. 평가 구간 라벨을 사전분포, 집단 정의, 하이퍼파라미터 선택에 되먹이지 않는다.

관문 기준·구간 정의·지표 정의·학습/평가 분할·라벨 규칙·검열 처리 규칙은 평가 규약에 있다.