저메탄 사료를 급여한 다음 주에 메탄 농도가 낮아졌다면 사료가 원인이라고 말할 수 있을까요? 시간상 먼저 일어난 사건과 뒤에 일어난 결과가 연결돼 보이지만, 이것만으로 인과관계가 성립하지는 않습니다. 같은 기간에 기온이 낮아졌거나 환기 설정이 바뀌고, 사육두수가 줄거나 사료 섭취량과 생산성이 달라졌을 수 있습니다. 센서 교정이나 측정 위치가 바뀐 경우도 있습니다.
인과관계를 확인한다는 것은 “사료를 적용하지 않았더라면 같은 동물 또는 같은 농장에서 메탄이 어떻게 되었을까”라는 반사실을 설계로 가깝게 추정하는 일입니다. 실제로 같은 시간의 한 개체를 처리와 비처리 상태에 동시에 둘 수 없으므로 비교군, 무작위 배정, 교차설계, 반복 측정과 통계 모형을 사용합니다.
전후 비교가 자주 틀리는 이유
가장 단순한 설계는 적용 전 평균과 적용 후 평균을 비교하는 것입니다. 현장 운영에는 편리하지만 계절 추세와 사료 효과가 겹칩니다. 예를 들어 기준기간은 한여름, 적용기간은 초가을이라면 열 스트레스, 섭취량, 환기량이 함께 달라집니다. 메탄이 10% 낮아져도 그중 얼마가 사료 때문인지 분리하기 어렵습니다.
평균 회귀도 주의해야 합니다. 우연히 메탄 값이 높은 주를 기준선으로 잡으면 이후 값은 아무 조치가 없어도 평소 수준으로 돌아갈 수 있습니다. 반대로 센서 드리프트가 적용일과 겹치면 가짜 효과가 만들어집니다. “적용 날짜와 그래프 하락이 일치한다”는 시각적 설명은 가설을 세우는 근거이지 효과 크기를 확정하는 근거가 아닙니다.
또 하나의 혼란은 지표 선택입니다. g CH₄/일, g CH₄/kg 건물섭취량, g CH₄/kg 우유, 축사 평균 ppm은 서로 다른 질문에 답합니다. 사료가 총 메탄을 낮추면서 섭취량이나 생산량도 낮춘다면 환경 성과와 생산 성과의 해석이 달라질 수 있습니다. FAO LEAP의 사료첨가제 지침은 첨가제 생산과 사용을 포함한 공급망 환경성과를 일관된 생애주기평가 방식으로 비교하도록 안내합니다. 장내발효 한 지표의 감소만으로 전체 환경성과가 자동 확정되는 것은 아닙니다.
현장에서 사용할 수 있는 인과 설계
가장 이해하기 쉬운 방식은 동시 비교군입니다. 비슷한 개체를 처리군과 대조군에 무작위로 배정하고 같은 기간, 같은 축사 운영 아래 관측합니다. 무작위 배정은 알려진 요인뿐 아니라 미처 생각하지 못한 개체차가 한쪽에 몰릴 가능성을 낮춥니다. 개체 수가 적다면 체중, 산차, 비유일수, 생산량 같은 주요 특성으로 블록을 만든 뒤 블록 안에서 배정할 수 있습니다.
개체가 자기 대조군이 되는 교차설계도 있습니다. 한 집단은 대조 사료 후 시험 사료를, 다른 집단은 반대 순서로 급여합니다. 개체 간 차이를 줄일 수 있지만 이전 처리의 효과가 다음 기간까지 남는 이월효과를 고려해야 합니다. 충분한 적응기간과 필요시 세척기간을 두고, 기간과 순서 효과를 분석에 포함합니다. 사료가 반추위 환경을 바꾸는 데 시간이 걸릴 수 있으므로 적용 직후 며칠만 골라 결론 내리면 안 됩니다.
농장 단위로 처리하는 경우에는 여러 농장을 처리군과 비교군으로 배정하는 군집설계를 검토합니다. 한 농장의 많은 센서 관측을 서로 독립인 수백 개 표본처럼 취급하면 표본 수를 부풀리는 의사반복 문제가 생깁니다. 처리가 농장에 적용됐다면 실험 단위는 농장일 수 있고, 개체와 시간 측정은 그 안에 중첩됩니다.
무작위화가 불가능한 상업 현장에서는 차이의 차이, 중단 시계열 같은 준실험 설계를 사용할 수 있습니다. 그러나 이 방법도 비교 농장이 같은 외부 추세를 따랐다는 가정이나 다른 동시 개입이 없었다는 가정이 필요합니다. 사전에 그래프와 운영 기록으로 가정을 점검하고 한계를 공개해야 합니다.
측정 계획이 실험 설계와 맞아야 한다
ICAR 지침은 필요한 결과와 현장 조건에 따라 호흡챔버, SF₆, GreenFeed, sniffer, 휴대형 축적 챔버와 대리지표 등 측정법의 적합성이 다르다고 설명합니다. 절대 배출량을 구하려는 시험과 개체 순위를 매기려는 시험에 같은 장비가 항상 최선은 아닙니다. 장비를 먼저 정한 뒤 질문을 끼워 맞추지 말고, 1차 평가변수와 필요한 정확도를 먼저 정해야 합니다.
측정은 하루의 특정 급여 직후에만 몰리지 않게 해야 합니다. 반추동물의 메탄 배출은 급여와 일주기 패턴을 따릅니다. ICAR GreenFeed 절차는 방문 측정을 하루에 분산하고 배경농도와 정량적 유량을 확보하는 운영 원칙을 제시합니다. 어떤 방법이든 처리군과 대조군의 측정 기회, 장치 접근성, 결측 규칙을 같게 유지해야 합니다.
사료 노출도 측정해야 합니다. 배합표에 들어갔다는 사실과 개체가 목표 용량을 실제 섭취했다는 사실은 다릅니다. 배치번호, 배합농도, 급여량, 잔량, 개체별 접근 기록, 중단일과 이상반응을 남깁니다. 계획대로 섭취한 집단만 사후에 골라 분석하면 순응도가 낮은 개체를 제거하는 선택 편향이 생길 수 있으므로 주 분석대상을 사전에 정합니다.
분석 전에 합의할 기준
1차 결과변수, 비교 기간, 제외 기준, 최소 데이터 가용률, 통계모형과 유의수준을 데이터 확인 전에 문서화합니다. 가능한 경우 표본 수는 기대 효과뿐 아니라 개체 간·일간 변동과 측정오차를 반영한 검정력 계산으로 정합니다. 작은 표본에서 유의하지 않았다는 결과는 “효과가 없다”와 같지 않고, 유의한 결과도 효과 크기와 신뢰구간이 넓다면 정밀하지 않을 수 있습니다.
분석에는 처리 외에도 기준값, 산차, 비유단계, 건물섭취량, 생산량, 기상과 환기처럼 사전에 정한 공변량을 포함할 수 있습니다. 다만 처리의 결과로 변한 변수를 무조건 통제하면 실제 효과 일부를 제거할 수 있습니다. 어떤 변수는 교란요인이고 어떤 변수는 매개요인인지 인과 경로를 먼저 그려 보는 것이 좋습니다.
한 지표와 한 기간에서 좋은 결과가 나왔다고 하위집단과 기간을 계속 바꿔 유리한 결과만 고르면 다중검정 문제가 커집니다. 탐색 분석은 탐색이라고 표시하고 확인 시험과 구분합니다. 효과의 재현성은 다른 계절, 다른 농장, 다른 사료 배치와 독립 분석에서도 확인해야 합니다.
실행 체크리스트
“사료가 없었을 때”를 대신할 동시 비교군 또는 타당한 준실험 비교를 둡니다.
실험 단위가 개체인지 축사인지 농장인지 명시하고 의사반복을 피합니다.
무작위 배정, 블록화, 교차 순서와 적응·세척기간을 사전에 정합니다.
1차 결과변수와 단위, 측정법, 표본 수, 제외·결측 규칙을 등록합니다.
실제 섭취량과 배치, 급여 중단, 생산성·건강 지표를 함께 기록합니다.
처리군과 대조군에 동일한 측정 일정과 품질 기준을 적용합니다.
효과 크기, 신뢰구간, 데이터 가용률과 민감도 분석을 함께 보고합니다.
장내발효 감소와 제품 단위 환경성과, 안전성·생산성을 구분합니다.
결론
저메탄 사료의 인과 효과는 적용 전후 그래프 한 장에서 나오지 않습니다. 비교 가능한 대조 상태를 만들고 개체차와 시간 변화를 분리하며, 실제 섭취와 측정 품질을 확인하고, 분석 규칙을 결과보다 먼저 정해야 합니다. 현장 시험은 연구실보다 복잡하지만 그 복잡성을 기록에서 지우는 것이 아니라 설계와 분석에 포함하는 것이 신뢰의 출발점입니다. 이 절차를 통과한 결과라야 “함께 변했다”를 넘어 “사료가 변화에 기여했다”는 주장에 가까워질 수 있습니다.