2006년부터 2009년까지 축산물품질평가원에서 등급 판정된 도체성적이 있는 개체와 한국종축개량협회에 등록이 확인된 231,382 두의 자료를 이용하여 각 도체형질의 자료구조 분석 등을 살펴보고, 각 요인별로 도체형질에 어떠한 영향을 주는 지에 대한 환경효과 등을 분석하고, 변수선택법 등을 이용하여 향후 유전능력평가에서 사용 가능한 통계모형 찾아내고자 실시하였다. 한우 암, 수, 거세를 포함한 전체집단에서 출하체중, 도체율, 등 지방두께, 근내지방도에서 정규분포를 보이고 있었고, 다른 형질은 정규분포에 근접함을 확인할 수 있었다. 각 요인에 대한 분산분석에서는 한우 도체형질에 대하여 성별, 출생년도-계절, 도축년도-계절, 출생지역, 도축일령의 모두 고도의 유의성(p<0.01)이 인정되었다. 변수선택법을 이용한 모형적합도 검정에서는 출하체중, 도체중, 도체율, 배최장근단면적은 다섯 가지 환경효과를 선택하는 것이 최적이었고, 등지방두께와 근내지방도는 네 가지 환경효과를 고려하는 것이 최적이었지만, 다형질 분석으로 한 번에 분석을 하기 위해서는 다섯 가지 환경효과를 함께 고려하는 것이 적합하다고 사료된다.
2006년부터 2009년까지 축산물품질평가원에서 등급 판정된 도체성적이 있는 개체와 한국종축개량협회에 등록이 확인된 231,382 두의 자료를 이용하여 각 도체형질의 자료구조 분석 등을 살펴보고, 각 요인별로 도체형질에 어떠한 영향을 주는 지에 대한 환경효과 등을 분석하고, 변수선택법 등을 이용하여 향후 유전능력평가에서 사용 가능한 통계모형 찾아내고자 실시하였다. 한우 암, 수, 거세를 포함한 전체집단에서 출하체중, 도체율, 등 지방두께, 근내지방도에서 정규분포를 보이고 있었고, 다른 형질은 정규분포에 근접함을 확인할 수 있었다. 각 요인에 대한 분산분석에서는 한우 도체형질에 대하여 성별, 출생년도-계절, 도축년도-계절, 출생지역, 도축일령의 모두 고도의 유의성(p<0.01)이 인정되었다. 변수선택법을 이용한 모형적합도 검정에서는 출하체중, 도체중, 도체율, 배최장근단면적은 다섯 가지 환경효과를 선택하는 것이 최적이었고, 등지방두께와 근내지방도는 네 가지 환경효과를 고려하는 것이 최적이었지만, 다형질 분석으로 한 번에 분석을 하기 위해서는 다섯 가지 환경효과를 함께 고려하는 것이 적합하다고 사료된다.
This study was conducted to study the environment effects on live weight, carcass weight, dressing percentage, eye muscle area, backfat thickness, and marbling score, which are the carcass traits of Hanwoo, based on the estimates and all the possible regression for the selection of variable and sign...
This study was conducted to study the environment effects on live weight, carcass weight, dressing percentage, eye muscle area, backfat thickness, and marbling score, which are the carcass traits of Hanwoo, based on the estimates and all the possible regression for the selection of variable and significance test for 231,382 heads that underwent the carcass measurements. The average and standard deviation for the live weight, carcass weight, dressing percentage, eye muscle area, backfat thickness, and marbling score were 654.79${\pm}$91.61 kg, 362.30${\pm}$67.15 kg, 59.52${\pm}$0.03%, 81.79${\pm}$12.21 $cm^2$, 11.39${\pm}$5.40 mm, 4.38${\pm}$2.29, respectively. The live weight, carcass weight, dressing percentage, eye muscle area, backfat thickness, and marbling score for cow were 532.79${\pm}$78.38 kg, 313.40${\pm}$44.90 kg, 56.50${\pm}$0.03%, 75.24${\pm}$10.69 $cm^2$, 11.82${\pm}$5.10 mm, 4.30${\pm}$2.06, respectively, while for bull were 619.74${\pm}$93.27 kg, 376.89${\pm}$48.62 kg, 58.61${\pm}$0.02%, 85.61${\pm}$10.46 $cm^2$, 5.64${\pm}$2.71 mm, 1.41${\pm}$0.83, respectively, and for steer were 681.78${\pm}$70.72 kg, 415.23${\pm}$49.43 kg, 60.19${\pm}$0.02%, 88.29${\pm}$10.27 $cm^2$, 12.71${\pm}$5.23 mm, 5.42${\pm}$1.99, respectively. In the environmental variables selection based on the variables selection method, the examination by carcass traits suggested that the most appropriate model could be determined when five variables were selected for the live weight, carcass weight, dressing percentage, eye muscle area, and four variables for backfat thickness, and marbling score. When they were considered at a time altogether based on multiple traits, it was deemed to be desirable to insert all five variables into the variables for analysis. In addition, high significance was found by carcass traits.
This study was conducted to study the environment effects on live weight, carcass weight, dressing percentage, eye muscle area, backfat thickness, and marbling score, which are the carcass traits of Hanwoo, based on the estimates and all the possible regression for the selection of variable and significance test for 231,382 heads that underwent the carcass measurements. The average and standard deviation for the live weight, carcass weight, dressing percentage, eye muscle area, backfat thickness, and marbling score were 654.79${\pm}$91.61 kg, 362.30${\pm}$67.15 kg, 59.52${\pm}$0.03%, 81.79${\pm}$12.21 $cm^2$, 11.39${\pm}$5.40 mm, 4.38${\pm}$2.29, respectively. The live weight, carcass weight, dressing percentage, eye muscle area, backfat thickness, and marbling score for cow were 532.79${\pm}$78.38 kg, 313.40${\pm}$44.90 kg, 56.50${\pm}$0.03%, 75.24${\pm}$10.69 $cm^2$, 11.82${\pm}$5.10 mm, 4.30${\pm}$2.06, respectively, while for bull were 619.74${\pm}$93.27 kg, 376.89${\pm}$48.62 kg, 58.61${\pm}$0.02%, 85.61${\pm}$10.46 $cm^2$, 5.64${\pm}$2.71 mm, 1.41${\pm}$0.83, respectively, and for steer were 681.78${\pm}$70.72 kg, 415.23${\pm}$49.43 kg, 60.19${\pm}$0.02%, 88.29${\pm}$10.27 $cm^2$, 12.71${\pm}$5.23 mm, 5.42${\pm}$1.99, respectively. In the environmental variables selection based on the variables selection method, the examination by carcass traits suggested that the most appropriate model could be determined when five variables were selected for the live weight, carcass weight, dressing percentage, eye muscle area, and four variables for backfat thickness, and marbling score. When they were considered at a time altogether based on multiple traits, it was deemed to be desirable to insert all five variables into the variables for analysis. In addition, high significance was found by carcass traits.
* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.
문제 정의
2006년부터 2009년까지 축산물품질평가원에서 등급 판정된 도체성적이 있는 개체와 한국종축개량협회에 등록이 확인된 231,382두의 자료를 이용하여 각 도체형질의 자료구조 분석 등을 살펴보고, 각 요인별로 도체형질에 어떠한 영향을 주는 지에 대한 환경효과 등을 분석하고, 변수선택법 등을 이용하여 향후 유전능력평가에서 사용 가능한 통계모형 찾아내고자 실시하였다.
따라서 본 연구는 축산물품질평가원에서 2006～2009년까지 등급판정한 개체 중 한국종축개량협회에 등록된 개체 231,382두에 대하여 각 도체형질 자료를 구조 분석하고 각 요인별로 도체형질에 어떠한 영향을 주는 지에 대한 환경효과를 분석하고, 변수선택법 등을 이용하여 향후 유전능력평가에서 사용 가능한 통계모형 찾아내고자 실시하였다. 그리고 더 나아가 도체정보와 개체식별번호 정보를 수집하고 이를 한국종축개량협회의 혈통등록 정보와 결합하여 농가보유 암소에 대한 유전능력평가를 위한 기초자료로 제공하고자 이번 연구를 실시하였다.
하지만 쇠고기 생산이력제가 도입되면서 등급판정결과와 개체번호를 결합하고 이 정보를 다시 한국종축개량협회의 혈통등록 정보와 결합시킬 경우 암소 자기 자신 또는 암소의 자손(거세, 비거세)에 대한 도축성적을 얻을 수 있는 제도적 장치가 마련되었다. 따라서 본 연구는 축산물품질평가원에서 2006～2009년까지 등급판정한 개체 중 한국종축개량협회에 등록된 개체 231,382두에 대하여 각 도체형질 자료를 구조 분석하고 각 요인별로 도체형질에 어떠한 영향을 주는 지에 대한 환경효과를 분석하고, 변수선택법 등을 이용하여 향후 유전능력평가에서 사용 가능한 통계모형 찾아내고자 실시하였다. 그리고 더 나아가 도체정보와 개체식별번호 정보를 수집하고 이를 한국종축개량협회의 혈통등록 정보와 결합하여 농가보유 암소에 대한 유전능력평가를 위한 기초자료로 제공하고자 이번 연구를 실시하였다.
제안 방법
도체형질의 일반분석에서는 각 형질별 일반성적을 분석하여 암, 수, 거세, 그리고 성별을 통합한 전체자료에 대한 기술통계량을 살펴보고, 자료의 정규성을 파악하기 위하여 각 형질별 분포를 나타내어 정규분포에 근사되는지를 파악하였다.
변수선택에 포함된 환경효과는 앞에서 분석한 성별, 출생년도, 출생계절, 출생지역과 추가적으로 도축년도와 도축계절 그리고 도축월령을 삽입하여 모형적합도 검정을 실시하였다. 또한 분석을 위하여 출생년도와 출생계절은 출생년도－계절, 도축년도와 도축계절은 도축년도－계절로 묶어서 분석을 실시하였다.
변수선택에 포함된 환경효과는 앞에서 분석한 성별, 출생년도, 출생계절, 출생지역과 추가적으로 도축년도와 도축계절 그리고 도축월령을 삽입하여 모형적합도 검정을 실시하였다. 또한 분석을 위하여 출생년도와 출생계절은 출생년도－계절, 도축년도와 도축계절은 도축년도－계절로 묶어서 분석을 실시하였다.
각 모형에 이용된 자료 등에 대한 일반 정보는 Table 1에 표시한 바와 같다. 분석한 형질은 출하체중(Live weight (LW), kg), 도체중(Carcass Weight (CW), kg), 도체율(Dressing Percentage (DP), %), 배최장근단면적(Eye Muscle Area (EMA), cm2), 등지방두께(Backfat Thickness(BF), mm), 근내지방도(Marbling Score(MS), score)이다. 2006년부터 2009년까지의 등급 판정 성적을 이용한 것은 쇠고기생산이력제가 2004년 10월부터 일부 브랜드 경영체와 시ㆍ군을 중심으로 추진되었고, 2006년에 전국 20여개 브랜드 경영체와 3개 시, 군이 참여하여 사업을 확대 추진하였기 때문에 개체 확인 및 등록여부가 사업 실시 이전의 도체성적 자료보다 일치되는 개체가 많았기 때문이다.
본 연구에 이용된 재료는 2006년부터 2009년까지 축산물품질평가원에서 등급 판정된 개체 중 한국종축개량협회에서 등록된 개체 확인 후 이상치를 제외한 231,382두의 기록을 이용하였다. 평가를 위하여 성별, 출생년도, 출생계절, 도축년도, 도축계절, 출생지역으로 나누어 분석하였다. 각 모형에 이용된 자료 등에 대한 일반 정보는 Table 1에 표시한 바와 같다.
환경효과 분석에서 중요한 것 중의 하나가 환경효과의 수가 많은 경우에 고려하고 있는 모든 환경효과를 모형에 포함시킬 것인가 하는 문제가 발생하여 환경효과를 모두 회귀모형에 포함하였으며, 환경요인들 중 일부만을 포함하는 회귀모형보다 결정계수의 값이 크고 설명력을 최대화 되도록 분석하였다. 하지만 가능한 모든 환경효과들을 회귀모형에 포함시키면 환경효과들 간의 상관관계가 커지는 단점이 있어 모형의 안정성과 신뢰성에 의문이 생길 수 있기 때문에 R2와 MSE (Mean squared error) 그리고 Cp 통계량으로 모든 가능한 회귀(All possible regression) 방법을 통하여 모든 가능한 변수들의 조합을 회귀 분석한 다음, Cp 통계량을 기준으로 하고 R2와 MSE를 참조하여 적정 모형을 선택하였다.
환경효과 분석에서 중요한 것 중의 하나가 환경효과의 수가 많은 경우에 고려하고 있는 모든 환경효과를 모형에 포함시킬 것인가 하는 문제가 발생하여 환경효과를 모두 회귀모형에 포함하였으며, 환경요인들 중 일부만을 포함하는 회귀모형보다 결정계수의 값이 크고 설명력을 최대화 되도록 분석하였다. 하지만 가능한 모든 환경효과들을 회귀모형에 포함시키면 환경효과들 간의 상관관계가 커지는 단점이 있어 모형의 안정성과 신뢰성에 의문이 생길 수 있기 때문에 R2와 MSE (Mean squared error) 그리고 Cp 통계량으로 모든 가능한 회귀(All possible regression) 방법을 통하여 모든 가능한 변수들의 조합을 회귀 분석한 다음, Cp 통계량을 기준으로 하고 R2와 MSE를 참조하여 적정 모형을 선택하였다.
대상 데이터
본 연구에 이용된 재료는 2006년부터 2009년까지 축산물품질평가원에서 등급 판정된 개체 중 한국종축개량협회에서 등록된 개체 확인 후 이상치를 제외한 231,382두의 기록을 이용하였다. 평가를 위하여 성별, 출생년도, 출생계절, 도축년도, 도축계절, 출생지역으로 나누어 분석하였다.
데이터처리
일반 성적을 분석하기 위하여 SAS MEANS Procedure를 이용하여 두수, 최소, 최대, 평균, 변이계수를 추정하였다. 정규분석은 SAS UNIVARIATE Procedure의 Checking Variables for Normality를 이용하여 Kolmogorov-Smirnov 통계량(D)으로 검정을 하였고 Kolmogorov-Smirnov 통계량(D) 추정을 위한 공식은 다음과 같다.
일반 성적을 분석하기 위하여 SAS MEANS Procedure를 이용하여 두수, 최소, 최대, 평균, 변이계수를 추정하였다. 정규분석은 SAS UNIVARIATE Procedure의 Checking Variables for Normality를 이용하여 Kolmogorov-Smirnov 통계량(D)으로 검정을 하였고 Kolmogorov-Smirnov 통계량(D) 추정을 위한 공식은 다음과 같다.
이론/모형
분석을 위해서 SAS REG Procedure의 Selection option을 활용하였다.
성능/효과
각 요인에 대한 분산분석에서는 한우 도체형질에 대하여 성별, 출생년도－계절, 도축년도－계절, 출생지역, 도축일령의 모두 고도의 유의성(p<0.01)이 인정되었다.
각 환경효과에 대한 분산분석 결과를 살펴보면 출하체중(LW), 도체중(CW), 도체율(DP), 배최장근단면적(EMA), 등지방두께(BF), 근내지방도(MS)는 성별, 출생년도, 출생계절, 출생지역, 도축일령에서 모두 고도의 유의성(p<0.01)이 인정되었다.
63으로 모형이 가장 적합하게 나타났다. 그리고 등지방두께와 근내지방도는 네 가지 환경효과(성별, 출생년도－계절, 출생지역, 출하월령)를 적용하였을 경우 각각 R2=0.06, Cp=4.01, MSE=27.18과 R2=0.14, Cp=5.65, MSE = 4.14로 이때 가장 적합한 모형을 보이는 것으로 판단되었다. 출하체중, 도체중, 도체율, 배최장근단면적은 다섯 가지 환경효과를 선택하는 것이 최적이었고, 등지방두께와 근내지방도는 네 가지 환경효과를 고려하는 것이 최적이었지만, 다형질 분석으로 한 번에 분석을 하기 위해서는 다섯 가지 환경효과를 함께 고려하는 것이 적합하다고 사료된다.
07로 최적의 모형을 설정할 수 있다. 도체중과 도체율, 배최장근단면적에서도 다섯 가지 환경효과를 적용하였을 경우 각각 R2 = 0.42, Cp = 6.00, MSE =2,371.22과 R2 = 0.25, Cp = 6.00, MSE = 5.68 그리고 R2 = 0.18, Cp = 6.00, MSE = 112.63으로 모형이 가장 적합하게 나타났다. 그리고 등지방두께와 근내지방도는 네 가지 환경효과(성별, 출생년도－계절, 출생지역, 출하월령)를 적용하였을 경우 각각 R2=0.
변수선택법을 이용한 모형적합도 검정에서는 출하체중, 도체중, 도체율, 배최장근단면적은 다섯 가지 환경효과를 선택하는 것이 최적이었고, 등지방두께와 근내지방도는 네 가지 환경효과를 고려하는 것이 최적이었지만, 다형질 분석으로 한 번에 분석을 하기 위해서는 다섯 가지 환경효과를 함께 고려하는 것이 적합하다고 사료된다.(주제어: 한우, 도체형질, 변수선택, 통계모형)
전체자료를 살펴보면 출하체중, 도체중, 도체율, 배최장근단면적, 등지방두께, 근내지방도는 각각 654.79±91.61 kg, 362.30±67.15 kg, 59.52±0.03%, 81.79±12.21 cm2 , 11.39±5.40 mm, 4.38±2.29로 추정되었다.
한우 암소집단에서는 도체율과, 등지방두께 그리고 근내지방도에서 정규성을 나타내고 있으며, 수소 집단에서는 배최장근단면적과 등지방두께 그리고 근내지방도에서 정규성이 나타났고, 거세우 집단에서는 등지방두께와 근내지방도에서 정규성이 나타났으며, 도체율은 정규성에 근접하였다. 전체적으로는 출하체중, 도체율, 등지방두께, 근내지방도에서 정규성이 나타났으며, 다른 형질에서도 정규성에 근접함을 확인할 수 있었다.
14로 이때 가장 적합한 모형을 보이는 것으로 판단되었다. 출하체중, 도체중, 도체율, 배최장근단면적은 다섯 가지 환경효과를 선택하는 것이 최적이었고, 등지방두께와 근내지방도는 네 가지 환경효과를 고려하는 것이 최적이었지만, 다형질 분석으로 한 번에 분석을 하기 위해서는 다섯 가지 환경효과를 함께 고려하는 것이 적합하다고 사료된다. 하지만, 자료별 구조 및 특성이 다르기 때문에 모든 한우 도체형질 자료를 상기 모형에 적용하는 것에는 지속적인 분석과 해석이 이루어져야 한다고 판단되었다.
한우 암, 수, 거세를 포함한 전체집단에서 출하체중, 도체율, 등지방두께, 근내지방도에서 정규분포를 보이고 있었고, 다른 형질은 정규분포에 근접함을 확인할 수 있었다. 각 요인에 대한 분산분석에서는 한우 도체형질에 대하여 성별, 출생년도－계절, 도축년도－계절, 출생지역, 도축일령의 모두 고도의 유의성(p<0.
). 한우 암소집단에서는 도체율과, 등지방두께 그리고 근내지방도에서 정규성을 나타내고 있으며, 수소 집단에서는 배최장근단면적과 등지방두께 그리고 근내지방도에서 정규성이 나타났고, 거세우 집단에서는 등지방두께와 근내지방도에서 정규성이 나타났으며, 도체율은 정규성에 근접하였다. 전체적으로는 출하체중, 도체율, 등지방두께, 근내지방도에서 정규성이 나타났으며, 다른 형질에서도 정규성에 근접함을 확인할 수 있었다.
후속연구
하지만, 자료별 구조 및 특성이 다르기 때문에 모든 한우 도체형질 자료를 상기 모형에 적용하는 것에는 지속적인 분석과 해석이 이루어져야 한다고 판단되었다. 또한 향후 한우 도체형질의 유전능력평가에 필요한 유전모수와 암소의 개량 정도를 고찰하고 암소검정사업 등에서 농가 보유 암소 개체별 유전능력평가를 할 수 있는 유전능력평가 컴퓨터 프로그램 개발이 추가적으로 이루어져야 할 것으로 사료된다.
질의응답
핵심어
질문
논문에서 추출한 답변
암소에 대한 농가차원의 능력검정을 위한 최선의 정보는 무엇인가?
현재 농가에서 사육하고 있는 암소에 대한 농가차원의 능력검정은 거의 이루어지지 않고 있고, 농림수산식품부의 한우개량사업에 참여하고 있는 농가 보유 암소에 대한 혈통이 현재까지 암소사육농가에서 얻을 수 있는 최선의 정보이다. 이러한 상황에서 농가의 암소에 대한 유전능력평가를 실시한다는 것은 모래위에 성을 쌓는 것과 같다.
도체형질의 자료구조 분석 등을 살펴보고 각 요인별로 도체형질에 어떤 영향을 주는지 분석하여 변수선택법 등을 이용하여 유전능력평가에서 사용 가능한 통계모형을 찾아내고자 실시한 연구의 결과는 어떻게 나타났는가?
한우 암, 수, 거세를 포함한 전체집단에서 출하체중, 도체율, 등지방두께, 근내지방도에서 정규분포를 보이고 있었고, 다른 형질은 정규분포에 근접함을 확인할 수 있었다. 각 요인에 대한 분산분석에서는 한우 도체형질에 대하여 성별, 출생년도－계절, 도축년도－계절, 출생지역, 도축일령의 모두 고도의 유의성(p<0.01)이 인정되었다.
변수선택법을 이용한 모형적합도 검정에서는 출하체중, 도체중, 도체율, 배최장근단면적은 다섯 가지 환경효과를 선택하는 것이 최적이었고, 등지방두께와 근내지방도는 네 가지 환경효과를 고려하는 것이 최적이었지만, 다형질 분석으로 한 번에 분석을 하기 위해서는 다섯 가지 환경효과를 함께 고려하는 것이 적합하다고 사료된다.(주제어: 한우, 도체형질, 변수선택, 통계모형)
한우 도체형질에 영향을 미치는 요인에는 어떤 것들이 있는가?
한우 도체형질에 영향을 미치는 성별, 출생년도, 출생계절, 출생지역, 도축일령의 환경효과를 추정하기 위하여 다음과 같은 선형 모형을 사용하였다.
참고문헌 (10)
Choi, T. J., Kim, S. D., Agapita J. Salces and Baik, D. H. 2006. Genetic Parameter Estimation on the Growth and Carcass Traits in Hanwoo (Korean Cattle). J. Anim. Sci. & Technol. (Kor.) 48(6):759-766.
Crews, D. H., Jr., Pollak, E. J. and Quaas, R. L. 2004. Evaluation of Simmental carcass EPD estimated using live and carcass data. J. Anim. Sci. 82:661-667.
Eriksson, S., Nasholm, A., Johansson, K. and Philipsson, J. 2004. Genetic relationships between calving and carcass traits for Charolais and Hereford cattle in Sweden. J. Anim. Sci. 82: 2269-2276.
Jeong, Mi Hwang, Kim, Si dong, Choy, Yun Ho, Yoon, Ho Baek and Park, Cheol Jin. 2008. Genetic Parameter Estimation of Carcass Traits of Hanwoo Steers. J. Anim. Sci. & Technol. (Kor.) 50(5):613-620.
Kolmogorov, A. N. 1933. "On the Empirical Determination of a Distribution Function", vol.4, pp.83-91.
Koo, Yang Mo, Kim, Jung Il, Song, Chi Eun, Shin, Jae Young, Lee, Jong Heon, Cho, Byoung Dai, Kim, Byeong Woo and Lee, Jung Gyu. 2008. A Study on Genetic Parameters of Carcass Weight and Body Type Measurements in Hanwoo Steer. J. Anim. Sci. & Technol. (Kor.) 50(2):157-166.
Roh, S. H., Kim, B. W., Kim, H. S., Min, H. S., Yoon, H. B., Lee, D. H., Jeon, J. T. and Lee, J. G. 2004. Comparison between REML and Bayesian via Gibbs Sampling Algorithm with a Mixed Animal Model to Estimate Genetic Parameters for Carcass Traits in Hanwoo (Korean Native Cattle). J. Anim. Sci. & Technol. (Kor.) 46(5):719-728.
SAS. Inst. Inc. 2004. SAS/STAT 9.1 user's Guide. SAS Inst. Inc., Cary, NC, USA.
Smirnov, N. V. 1939. "On the Estimation of the Discrepency between Empirical Curves of Distribution for Two Independent Samples", Bulletin of Moscow, vol.2, pp.3-16.
Yoon, H. B., Kim., S. D., Na, S. H., Chang, U. M., Lee, H. K., Jeon, G. J. and Lee, D. H. 2002. Estimation of Genetic Parameters for Carcass Traits in Hanwoo Steer. J. Anim. Sci. & Technol. (Kor.) 44(4):383-390.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.