본 연구는 혼파초지에서 기상변수를 기본으로 시비량, 파종량 및 조성연차 변수를 단계적으로 적용하여 정밀도(모형의 해석력 및 풀사료 생산 이론과 적합여부)가 높은 수량예측모형을 구축하는데 목적이 있다. 혼파초지 수량예측모형 구축과정은 풀사료 및 기상자료 수집(자료수집), 엑셀에 수치를 입력(1차가공), 풀사료의 오류 자료를 수정(2차가공), 풀사료자료와 기상자료를 병합하며 기상변수를 선정(3차가공), 이상점 및 결측치 제거를 통해 데이터세트를 구성(4차가공)하였다. 데이터세트는 상관분석, 기술통계분석, 다중회귀분석 및 일반선형분석 등을 이용하여 수량예측모형을 구축하고 모형의 정밀도를 확인하였다. 수량예측모형의 단계적 적용은 우선 변수탐색 과정으로서 1) 상관분석을 이용하여 기상변수의 기준을 설정하였다. 다음 단계로서 2) 기술통계량 및 히스토그램을 통해 반응변수(건물수량)의 정규분포 및 설명변수(기상, 시비량, 파종량 및 조성연차)의 분포를 확인하였다. 3) 혼파초지의 수량예측모형 구축 과정은 변수(기상, 시비, 파종량 및 조성연차)를 단계적으로 적용하여 모형 I∼VI을 구축하였다. 4) 혼파초지의 수량예측 검토과정은 구축된 모형 I∼VI 중 해석력을 고려하여 가장 우수한 모형을 선택한 뒤 풀사료 생산 이론 관점에서 검토하였다. 건물수량에 영향을 미치는 주요 기상변수는 최한월평균기온, 생육일수 0, 강수일수 0, 강수량 0 및 하고일수였고 건물수량이 정규분포를 이루는 것을 확인하였다. 단, 시비량 및 파종량은 분포가 편중되는 경향이 나타났다. 모형의 단계적 적용결과 모형 VI(기상, 시비량, 파종량 및 조성연차를 고려)의 해석력이 가장 우수하였으나 하고일수와 건물수량은 정(+)의 상관관계로 나타나 풀사료 생산 이론과 맞지 않았다. 이를 해결하기 위해 지역별 하고일수와 건물수량의 관계를 검토한 결과 12개 지역 중 6개 지역에서 부(-)의 상관관계, 3개 지역에서 정(+)의 상관관계로 나타났으며 3개 지역은 표본크기가 부족하여 상관관계를 알 수 없었다. 한편 각 지역별 표본크기가 부족하여 지역별 수량예측모형은 구축할 수 없었다. 따라서 기존 변수로는 문제를 해결할 수 없다고 판단하여 혼파초지 건물수량에 보다 밀접한 영향을 미친다고 생각된 적산기상변수(하고적산기온, 적산생육도일 및 적산일조시간)를 추가하여 수정(Modified) 수량예측모형을 구축하여 하고일수의 문제를 해결하였다. 최종적으로 선택된 모형 VI의 그룹 1은 건물수량=91.8×강수일수 0 + 15.3×인산 시비량 + 8.7×질소 시비량 + 2.5×적산생육도일 0~25 - 214.6×파종량 - 9.8×칼리 시비량 - 1.2×하고적산기온 - 1.0×적산강수량 0 - 3875.3이었으며, 모형 VI의 그룹 2는 건물수량 = 128.7×강수일수 0 + 24.7×인산 시비량 - 8669.2였다. 이상에서 최종적으로 선택된 수정 모형 VI는 해석력 및 풀사료 생산 이론을 고려하였을 때 정밀도가 높은 모형이었다. 수정 모형 VI의 그룹 1은 조성연차 1∼3년 및 그룹 2는 조성연차 4∼7년에서 이용할 수 있었다.
This research was conducted to construct a yield prediction model with high precision for mixed pasture in South Korea. The climatic variables were considered as basic explanatory variables, then fertilization, seeding rate and number of periods after pasture establishment were considered as explanatory in order to improve the precision of the yield prediction model for mixed pasture. The processes of constructing yield prediction model for mixed pasture included the following steps: forage and climatic data collection (data collection), inserting data into excel file (1st processing), correction of errors in forage data (2nd processing), selection of climatic variables after combination of forage and climatic data (3rd processing), final data set generation after elimination of the outliers and data points with missing values (4th processing). Yield prediction models were constructed through correlation analysis, descriptive statistics analysis, multiple regression analysis and general linear analysis. Meanwhile the precision of the models was checked The stages of model construction were as follows: 1) setting the standards for climatic variables selection through correlation analysis; 2) checking the normal distribution of response variable (dry matter yield) and the distributions of explanatory variables (climatic variables, fertilization, seeding rate and periods after pasture establishment); 3) constructing the yield prediction models (Model I-VI) for mixed pasture considering different variables (climatic variables, fertilization, seeding rate and periods after pasture establishment) via statistical methods in stages; 4) selecting the optimal model via considering the interpretabilities among the models (Model I-VI), then examined the selected model from the perspective of forage production theories. The selected climatic variables including Mean Temperature in January (MTJ), Growing Days 0 (GD 0), Number of Days with Precipitation 0 (NDP 0), Accumulated Amount of Precipitation 0 (AAP 0) and Days of Summer Depression (DSD) were considered as having significant effects on dry matter yield. The normal distribution of dry matter yield was confirmed. While the observed distribution of fertilization and seeding rate were tended to be skewed. Model VI (considering climatic variables, fertilization rate, seeding rate and periods after pasture establishment) had the best interpretability. However, DSD showed positive correlation with dry matter yield in this model, this is inconsistent with the forage production theories. To solve this problem, the correlation relations of DSD and dry matter yield were investigated by cultivated locations. The results showed that negative correlation were present in 6 locations, positive correlation relations were observed in 3 locations and correlation were unavailable in 3 locations because of too small data size. Meanwhile, the sample size of each cultivated location was too small to fulfill the statistical model to construct yield models by locations. Therefore, it could be judged that the problems existing in the variable DSD could not be solved in the data set used in this research. Three new variables related to accumulated temperature (Accumulated Summer Depression, Accumulated Growing Degree Days and Accumulated Duration of Sunshine) were generated and they were considered to have more direct relations with yield prediction. Then the modified yield prediction model for mixed pasture based on Model VI were constructed as follows: The precision of the modified Model VI was considered as the best among the model we made based on its interpretability and consistence with forage production theories. The equation of Group 1 of the modified Model VI could be utilized for yield prediction of mixed pasture in the 1 to 3 periods after pasture establishment. The equation of Group 2 of the modified Model VI could be utilized for yield prediction of mixed pasture in the 4 to 7 periods after pasture establishment.