Introduction
현재 통계청(KOSTAT)의 2018년 축산물 생산비 조사결과에 따르면 한우 비육우의 축산물 생산비(100 kg)는 전년도 보다 88,000원(8.7%) 상승한 1,106,000원이며, 경영비(100 kg)는 76,000원 증가한 955,000원으로 조사되었다. 한우 두당 사육비는 2017년 7,672,000원에서 2018년 8,406,000원으로 상승하였다. 한우 비육우 순수익은 전년도 보다 189,000원 하락한 -57,000원으로 전년대비 사육비가 총 수입보다 높게(2.6%) 상승했다. 축산물 등급제 시행 이후 한우 도매가격과 수입량은 역대 최고 수준을 기록하고 있으며 한우 자급률 동향과 향후 전망 검토가 필요한 시점이다(KREI, 2016). 현재 국내 한우 자급률은 2015년도 46.1%, 2016년도 39.0%, 2017년도 41.0%, 2018년도 36.4%로 감소의 추세를 보이고 있다(MAFRA, 2018). 국내 한우사업의 발전을 위해 한우 도체형질 개량에 환경요인이 미치는 영향에 대한 연구(Lee et al., 2018), 한우의 경락 가격에 대한 요인별 기여도 분석(Sun et al., 2012; Eum et al., 2016; Kang et al., 2019), 한우 도체형질과 경락 가격에 대한 도축월령 및 환경요인 효과분석(Hwang et al., 2016), 한우의 도체 형질이 경제성에 미치는 영향(Lee et al., 2011), 경남지방 한우의 경락 가격에 관여하는 요인별 기여도 분석(Park et al., 2015), 평창·영월·정선지역 출하 우 경락 가격에 영향을 미치는 요인 분석(Lee et al., 2013) 등 한우사업 관련 연구가 진행되었다. 하지만 각 지역별 한우에 관한 연구는 많이 부족하여 지역 한우 브랜드 가치를 높여 소비 증감을 얻기 위한 노력이 필요하다. 그 중 경기지역은 전국에서 가장 많은 소비량과 유입량이 이루어지는 수도권 지역에 가장 빠르고 신선한 유통이 가능하다(KAPE, 2015). 소비량이 가장 많은 지역이며 유통에 유리한 지역인 경기지역 한우의 개량을 위하여 도체형질(등지방두께, 등심단면적, 도체중, 근내지방도)이 경락 가격과 총 가격에 기여하는 정도와 도체형질에 미치는 환경요인의 효과를 분석하였고, 경기 외 지역과 비교 분석을 통하여 경기지역 한우의 농가 발전과 소득 증진으로 인하여 자급률을 높이고, 경기지역 한우만의 특색과 강점을 찾기 위하여 연구를 진행하였다.
Materials and Methods
공시재료
본 연구는 2012년부터 2018년도 경기지방에서 도축된 한우 총 11,610두의 자료에서 빈도수가 상대적으로 부족한 2012년에 도축된 개체의 데이터를 제거하고, 전체 년도의 수소 데이터를 제거한 후 필요한 요인들의 데이터에 결측값과 표준편차의 3배 이상되는 값을 이상치로 간주하고 제거하여 7,145두를 데이터로 선정하였다. 본 연구의 효과별 빈도는 Table 1과 같다. 또한, 경기 외 지방을 비교 데이터로 선정하고 같은 기간 데이터를 기준으로 결측값이 있는 개체를 제거하고 표준편차의 삼배수를 벗어나는 자료까지 제거하여 최종 106,779두의 자료를 이용하였다.
Table 1. Number of records of Hanwoo by sex, season at slaughter and year at slaughter.
|
|
NDF, neutral detergent fiber; ADF, acid detergent fiber. |
분석방법
환경요인의 효과
본 연구에서 조사한 도체중, 등심단명적, 등지방두께 및 근내지방도에 영향을 미치는 성, 도축계절 및 도축년도의 효과를 추정하기 위해 다음과 같은 선형혼합모형을 이용하여 최소제곱법으로 분산분석을 실시하였다.
Yijkl = μ + Sexi + SSeasonj + SYeark + eijkl (1)
여기서,
Yijkl: i번째 성의 j번째 도축 년도의 k 번째 도축 계절에 대한 측정치
μ: 전체평균
Sexi: i번째 성의 효과(i = 1, 2)
SSeasonj: j번째 도축계절의 효과(j = 1, 2, 3, 4)
SYeark: k번째 도축년도의 효과(k = 1, 2, 3, 4, 5, 6)
eijkl: 임의오차의 효과
본 연구에서 설정한 Linear Model은 PC용 SAS Package (version 9.4, SAS Institute, Cary, NC, USA)를 이용하였고, GLM (generalized linear model)분석결과 제공되는 4가지 제곱합 중에서 본 논문에서 이용되는 요인들과 같은 불균형된 자료에 적합한 TYPE Ⅲ 제곱합을 이용하여 분산분석 하였으며, 최소제곱평균치 간의 유의성 검정을 위하여 다음과 같은 귀무가설을 설정하고 유의수준 5%로 각각 검정하였다.
H0; LSM(i) = LSM(j) (2)
여기서, LSM(i(j)): i(j)번째 효과의 최소제곱평균치(i ≠ j)
요인별 기여도 분석
한우의 가격은 도체형질들에 의해 결정된다. 그러므로 경기지방과 경기 외 지방의 도체 형질들이 얼마나 영향을 미치는지 요인별로 기여도를 알아보기 위하여 경락 가격을 종속 병수로 하고 도체형질들을 독립변수로 하여 아래와 같은 선형모형을 이용하여 다중회귀분석(multiple regression analysis)을 실시하였다.
Y = a + b1BF + b2EMA + b3CW + b4MS + e (3)
여기서,
BF: Backfat thickness (등지방두께)
EMA: Eye muscle area (등심단면적)
CW: Carcass weight (도체중)
MS: Marbling score (근내지방도)
여기서, Y는 종속변수로 이용된 경락 가격과 총 가격이며, a는 Y절편을 의미한다. 등지방두께, 등심단면적, 도체중, 근내지방도 네 가지 도체형질을 독립변수로 설정하였고, b1, b2, b3 및 b4는 각 형질들의 회귀계수이며 e는 임의 오차이다. 또한 본 연구에서는 다중회귀분석에 부분상관보다 준부분상관이 중요하게 다루어지기 때문에, 각 독립변수의 기여도를 따지는데 많이 활용되고 있는 준부분상관자승(squared semi-partial)을 분석 이용하였다.
Results and Discussion
일반성적 및 분산분석 결과
일반성적
본 연구에서 조사된 경기도 지방과 경기 외 지방의 한우 가격 및 도체형질에 대한 일반성적은 Table 2와 Table 3과 같다. 본 연구에서 조사된 결과를 살펴보면 경기지방 한우의 등지방두께, 등심단면적, 도체중, 근내지방도, 경락 가격 및 총 가격의 평균과 표준편차는 각각 13.24 ± 5.15 mm, 85.44 ± 11.81 cm2, 372.69 ± 63 kg, 4.42 ± 2.01점, 16,230 ± 3,181원 및 6,137,488 ± 1,878,406원으로 나타났다. 같은 기간 도축된 경기 외 지방 한우의 등지방두께, 등심단면적, 도체중, 근내지방도, 경락 가격, 총 가격의 평균과 표준편차는 각각 13.97 ± 5.23 mm, 85.19 ± 12.03 cm2, 365.41 ± 58.9 kg, 4.23 ± 1.96점, 15,886 ± 3,162원 및 5,881,521 ± 1,767,011원으로 나타났다. 경기지방은 타 지역에 비해 등지방두께가 얇은 강점을 나타냈고, 경락 가격과 총 가격이 높은 것으로 조사되었다.
Table 2. Simple statistics for carcass traits and price of Hanwoo in Gyeonggi-do province.
|
|
SD, standard deviation; Max, maximum; Min, minimum. z Total price = auction price × carcass weight. |
분산분석
Table 4는 각 요인에 대한 분산분석 결과인데, 성에 따른 효과에서는 조사된 모든 형질에서 고도의 유의성이 인정되었고, 도축계절의 효과에서도 근내지방도에 대한 도축계절의 효과를 제외하고는 모든 형질에서 고도의 유의성을 보였으며, 도축 년도의 효과에는 모든 형질에서 고도의 유의성을 보였다.
Table 4. Mean squares and test of significance of carcass traits in Hanwoo.
|
|
d.f, degree of freedom. **p < 0.01; NS, not significant at 0.05 level of significance. |
환경요인의 효과
성의 효과
Table 5는 한우의 도체형질에 대한 성의 효과를 나타낸 것이다. 모든 형질에서 거세우가 유의적으로 높게 나타났으며, 암소는 보다 낮은 성적을 나타냈다.
도축계절의 효과
한우의 도체형질에 대한 도축계절의 효과는 Table 6과 같다. 도축계절에 따른 효과를 살펴보면, 등지방두께는 여름에 가장 얇게 나타나고, 등심단면적은 봄에 가장 높은 성적을 보였다. 도체중은 봄과 겨울에 가장 높았다. 하지만 근내지방도는 사계절 내에서 유의적인 차이를 보이지 않았다.
도축년도의 효과
한우의 도체형질에 대한 도축년도의 효과는 Table 7과 같다. 도축년도의 효과에 대하여 살펴보면, 등지방두께, 등심단면적, 도체중은 시간이 지남에 따라 성적이 유의적으로 높아졌지만 2018년에서 근내지방도는 조금 감소하였다. 도체형질의 성적은 대체적으로 시간이 지남에 따라 성적이 개선되었지만, 등지방두께도 같이 높아져 줄이기 위한 개선이 필요하며 근내지방도도 개선이 요구된다.
표현형 상관
도체 형질들 간의 표현형상관을 Table 8에 표시하였다. 모든 분석 형질들과 경락 단가 및 총 가격 간의 상관이 존재하는 것으로 나타났다. 경락 가격에서는 등지방두께의 상관에 대해 -0.02로서 부의상관을 보였으며, 등지방두께를 제외한 나머지 형질들은 정의상관을 나타냈는데, 그 중 근내지방도가 0.71로 가장 높은 정의상관을 보였다. 총 가격에서는 모든 형질들이 정의상관을 보였고 그 중 도체중이 0.82로 가장 높은 정의상관을 나타냈다.
Table 8. Phenotype correlation coefficients among carcass traits in Hanwoo.
|
|
z Total price = auction price carcass weight. ** p < 0.01. |
요인별 기여도 분석
Table 9는 요인별 기여도를 나타낸 것이다. 독립변수의 측정 단위가 다른 경우 종속변수에 미치는 기여도 계산은 준부분상관제곱(squared semi-partial correlation) 값을 이용하는데 본 연구에서도 이 값을 이용하였으며, 준부분상관제곱 분석에서 제시하는 TypeⅡ의 값을 이용하였다. 먼저 회귀계수 추정 값을 살펴보면 경락단가에서는 등지방두께, 등심단면적, 도체중 및 근내지방도가 각각 -126.16, 25.39, 4.87 및 1,014.65원으로 나타났고, 절편 값은 9,304.29원으로 나타났다. 총 가격에서는 등지방두께, 등심단면적, 도체중, 근내지방도가 각각 -50,509, 9,584.44, 17,453 및 373,245원으로 나타났고 절편 값은 -2,259,859원으로 나타났다. 경락단가에서는 회귀식의 R-Square가 0.57로 나타나 57%를 설명하고 있음을 확인할 수 있었고, 등지방두께, 등심단면적, 도체중 및 근내지방도의 기여도가 각각 10.42, 1.40, 0.92 및 86.70%로 나타나 경락단가에서 근내지방도의 기여도가 절대적인 것으로 나타났다. 총 가격에서는 회귀식이 R-Square가 0.83로 나타나 83%를 설명하고 있음을 확인할 수 있었고, 등지방두께, 등심단면적, 도체중 및 근내지방도의 기여도가 각각 6.49, 0.77, 45.95 및 45.60%로 나타나 도체중과 근내지방도의 영향을 많이 받은 것으로 나타났다.
Table 10은 경기 외 지역에서 도축된 한우의 경락 가격과 총 가격에 도체형질별 기여도를 나타낸 표이다. 먼저 회귀계수 추정 값을 살펴보면 경락단가에서는 등지방두께, 등심단면적, 도체중 및 근내지방도가 각각 -119.87, 24.32, 5.43 및 1,055.15원으로 나타났고, 절편 값은 8,985.29원으로 나타났다. 총 가격에서는 등지방두께, 등심단면적, 도체중, 근내지방도가 각각 -46,634, 8,616.53, 17,323원 및 382,168원으로 나타났고 절편 값은 -2,203,486원으로 나타났다. 경락단가에서는 회귀식의 R-Square가 0.57로 나타나 57%를 설명하고 있음을 확인할 수 있었고, 등지방두께, 등심단면적, 도체중 및 근내지방도의 기여도가 각각 9.06, 1.24, 0.99 및 88.53%로 나타나 경락단가에서 근내지방도의 기여도가 절대적인 것으로 나타났다. 총 가격에서는 회귀식이 R-Square가 0.81로 나타나 81%를 설명하고 있음을 확인할 수 있었고, 등지방두께, 등심단면적, 도체중 및 근내지방도의 기여도가 각각 5.86, 0.66, 43 및 49.64%로 나타나 도체중과 근내지방도의 영향을 많이 받은 것으로 나타났다.
경락 가격에 대해서는 경기지역과 경기 외 지역 모두 근내지방도의 기여도가 가장 높았으며, 총 가격은 경기지역에서는 도체중, 근내지방도 순으로 높았으며, 경기 외 지역에서는 근내지방도, 도체중 순으로 가장 높은 영향력을 보여 경기지역의 총 가격 결정 시 근내지방도의 영향력이 경기 외 지역보다는 낮다는 것을 알 수 있었다. 또한 경기지역의 등지방두께에 대한 기여도가 경락 가격, 총 가격 모두에서 높은 기여도를 보인 특징을 나타냈다. 이는 경기지역 한우만의 특징을 나타낸다.
Conclusion
경락 가격에 대해서는 경기지역과 경기 외 지역 모두 근내지방도의 기여도가 가장 높았으며, 총 가격은 경기지역에서는 도체중, 근내지방도 순으로 높았으며, 경기 외 지역에서는 근내지방도, 도체중 순으로 가장 높은 영향력을 보여 경기지역의 총 가격 결정 시 근내지방도의 영향력이 경기 외 지역보다는 낮다는 것을 알 수 있었다. 또한 경기지방에서 등지방두께는 경락 단가, 총 가격 모두 경기 외 지방보다 높은 기여도를 나타냈다. 이상의 연구결과를 통해 경기지방 한우의 경락 단가, 총 가격에 대한 요인별 기여도를 산출하였고, 경기 외 지방과 비교 분석을 할 수 있었다. 이상의 분석결과를 통해 경기지방 한우는 타 지방 보다 등지방두께가 얇음을 보인 강점을 나타냈고, 경락 단가, 총 가격에 대한 기여도 또한 다른 지방보다 높게 나타내는 특징을 이용하여, 경기 지방 한우의 가치를 증진시키고, 상관 분석을 통하여 도체중이 증가 시 등지방두께도 같이 증가함을 고려하여 도체성적 개량에 이용되어야 할 것이다. 경락 단가와 총 가격의 기여도를 고려하여 가장 높은 수익을 얻도록 전략적으로 출하시기를 선택하는데 활용하고자 하며, 각 형질들의 환경요인에 따른 분산분석 결과를 통해 한우농가 수익증대와 한우 개량 방향제시를 위한 기초자료로 활용하고자 한다.
Authors Information
Jewoong Yun, https://orcid.org/0000-0002-9129-6059
Yoseph Kim, https://orcid.org/0000-0002-7987-8726
Jieun Lee, https://orcid.org/0000-0003-1271-9090
Tea Hun Kang, https://orcid.org/0000-0002-9344-8386
Myunghoo Kim, https://orcid.org/0000-0002-8444-6952
Jakyeom Seo, https://orcid.org/0000-0002-9176-5206
Seong-Keun Cho, https://orcid.org/0000-0001-9309-2327
Byeong-Woo Kim, https://orcid.org/0000-0003-3300-0173
