본 실험은 경기도 소재 한우 암소 사육농가에서 3곳에서 수행되었으며, 각 실험농장에서 한우 3산 이상의 다산암소를 선발하였다. 암소 전용 배합사료는 비육전기와 비육후기로 구분하여 실험농가에 공급하고 실험 대상축에 일일 두당 \( 8 \sim10 \mathrm{kg} \) 수준에서 자유채식하도록 급여하였다. 각 실험 농장에서 실험대상 가축은 우방을 따로 관리하였다. 실험용 배합사료를 급여하기 전에 생체중을 측정하고 출하시 체중과 비교하여 실험 사료 급여 기간 동안의 일당 증체량을 계산하였다. 쇠고기 시료는 갈비 12번째 마디에서 13번째 마디의 등심 (M. Longissimus dorst) 부위의 냉장육 \( 2.5 \mathrm{kg} \)을 이용하였으며, 농장 3곳의 처리구와 대조구간 평균 성적 차이를 이용하여 사료급여효과를 분석할 때 렌덤효과모형을 적용하여 통합적인 효과를 산출하였다.
<h1>재료 및 방법</h1> <h2>실험농장 및 실험 두수</h2> <p>본 실험은 경기도 소재 3개의 한우 암소 사육농가에서 수행되었으며, 각 실험농장에서 한우 3산 이상의 다산암소를 선발하였다. 농장에서는 대조구에 3두, 처리구에 29두를 공시하였으며, 농장B에서는 대조구에 3두, 처리구에 8두를 공시하였다. 농장C에서는 대조구에 4두 그리고 처리구에 11두를 공시하였다. 결과적으로 총 3개의 농장에서 다산 한우 암소 48두(평균체중, \( 456 \pm 75 \mathrm{~kg} \))를 처리구 그리고 10두(평균체중,\( 513 \pm 21 \mathrm{~kg}) \)를 대조구로 공시하여 실험에 사용하였다.</p> <h2>실험사료 및 사양관리</h2> <p>암소 전용 배합사료는 비육전기와 비육후기로 구분하여 아래의 표(Table 1과 Table 2)와 같이 생산하여 실험농가에 공급하고 실험 대상축에 일일 두당 \( 8 \sim10 \mathrm{~kg} \) 수준에서 자유채식하도록 급여하였다. 대조구의 사료는 일반 사료회사에서 판매하는 거세우 비육후기 사료를 이용하였다. 실험기간은 농장여건에 따라서 다소간의 차이는 있었으나 평균 \( 211 \pm 30 \) 일간 급여하였고, 실험개시 후 4 개월 동안은 다산암소비육 전기사료를 급여하였고, 이후 출하시까지의 기간 동안에는 다산암소비육 후기사료를 급여하여 효과를 검정하였다. 급여사료를 제외한 나머지 사양관리 조건은 일반적으로 한우농장에서 사용하는 방법에 준하여 수행되었다.</p> <h2>증체 효과</h2> <p>각 실험 농장에서 실험대상 가축은 우방을 따로 관리하고 실험용 배합사료를 급여하기 전에 생체중을 측정하여 출하시 체중과 비교하여 실험 사료 급여 기간 동안의 일당 증체량을 계산하였다.</p> <h2>도체 성적 분석</h2> <p>쇠고기 시료는 갈비 12번째 마디에서 13번째 마디의 등심 (M. Longissimus dorst) 부위의 냉장육 \( 2.5 \mathrm{~kg} \)을 이용하였다. 농협 부천 공판장에 출하하여 도축된 공시축의 등급판정은 축산물품질평가원 소속 등급사가 측정하였다. 공시축의 육질 등급 요인(근내지방도, 등지방두께, 등심단면적, 육량등급, 육색, 지방색, 조직감, 성숙도)은 20-22시간 냉장된 반도체의 13번째 마디 갈비부위를 절개하여 측정하였다. 결과로 얻어진 육량등급A, B 및 C 등급을 각각 1, 2 및 3 점으로 나타내어 효과 검정에 사용하였다. 근내지방도는 1에서 9까지 등급을 구분하여 평가하였고, 육색은 1에서 7등급으로 나누었으며, 1등급은 선홍색을 나타내며, 7등은 검붉은색을 나타낸다. 조직감은 1에서 3등급으로, 성숙도는 1에서 9등급으로 구분하여 평가하였다. 지방색은 1에서 7등급으로 나누었으며, 1등급은 흰색 그리고 7등급은 노란색으로 구분하였다. 육질등급은 \( 1^{++} \) 등급은 1 점, \( 1^{+} \)등급은 2점, 1 등급은 3점 그리고 2등급과 3등급은 각각 4점과 5점으로 구분하여 평가하였다.</p> <h2>통계분석</h2> <p>본 연구는 3개의 서로 다른 농장에서 일반적인 사양관리 방법에 준하여 수행된 것이며 농장마다 서로 다른 기간 동안에 그리고 서로 다른 관리인 및 관리프로그램 아래서 수행되었다. 따라서 농장간 연구결과에는 이질성이 존재할 수 있다. 이에 본 연구에서는 메타분석방법을 이용하여 통합적인 처리 효과를 산출하였다. 3개 농장에서의 처리구와 대조구간 평균 성적 차이를 이용하여 사료급여효과를 분석함에 있어 렌덤효과모형을 적용하였고, 통합적인 효과를 산출하였다. 렌덤효과모형의 산출방법 및 통합효과 분석방법은 Kook 등의 방법에 준하여 실시하였다. 대략적인 방법은 다음과 같다. 다산암소 전문사료의 급여와 일반비육우사료의 급여를 통하여 얻어진 각 성적들의 평균차 \( \left(d_{i}\right) \)를 산출하였고, 렌덤효과모형을 산출하기 위하여, 각 평균차에 대한 근사적 역분산 \( \left(w_{i}\right) \)을 산출하였다. 최종적으로 통합적 효과(summary effect, \( \bar{d}) \)를 아래와 같은 수식을 이용하여 산출하였다.</p> <p>\( \bar{d}=\frac{\sum_{i}^{k} d_{i} w_{i}}{\sum_{i}^{k} w_{i}} \)</p> <p>위의 통합적 효과에서 각 실험농장에서 발생하는 편의를 보정하기 위하여 \( \hat{\tau}^{2} \) (tau-squared)를 산출하였고, 다시 렌덤효과모형의 조정된 근사적 역분산 \( \left(w_{i}^{*}\right) \)을 계산하였고 최종적으로 편의가 조정된 렌덤효과모형의 통합적 효과를 산출하였다.</p> <p>\( \hat{\tau}^{2}=\operatorname{Max}\left(0, \frac{Q-(k-1)}{\sum_{i}^{k} w_{i}-\left(\sum_{i}^{k} w_{i}^{2}\right) / \sum_{i}^{k} w_{i}}\right) \)</p> <p>\( w_{i}^{*}=\frac{1}{\left(d_{i}+\hat{\tau}^{2}\right)} \)</p> <p>\( \bar{d}^{*}=\frac{\sum_{i}^{k} d_{i} w_{i}^{*}}{\sum_{i}^{k} w_{i}^{*}} \)</p> <p>\( \operatorname{Var}\left(\bar{d}^{*}\right)=\frac{1}{\sum_{i}^{k} w_{i}} \)</p> <p>최종적으로 표준화된 효과와 효과의 표준편차를 이용하여 \( Z \)-value를 산출하였고, 정규누적분포의 양측검정을 이용하여 확률값을 산출하였다. 렌덤효과모형으로 얻어진 결과들의 편의 검증은 Egger's regression을 이용하여 검정하였다.</p>