• 한국농림기상학회지
• /
• 제25권3호
• /
• pp.129-141
• /
• 2023

작물 모형은 작물의 유전적 특성을 나타내는 품종모수를 요구하며, 품종모수는 작물의 개별 품종별로 추정되어야 한다. 품종모수의 추정에는 고품질의 많은 생육 자료가 요구되지만, 자료의 생산에 상당한 비용이 필요하다. 비교적 낮은 품질의 가용성이 높은 자료를 활용하는 대신, 대량의 랜덤 모수를 생성하고 이를 평가하여 품종모수를 추정할 수 있다. 본 연구에서는 SIMPLE 작물 모델의 불확도를 최소화하기 위해 품종모수 추정 방식을 비교하고, 두 앙상블 방식과 대한 비교를 하였다. 모수 추정을 위한 Metropolis-Hastings (MH) 알고리즘에 대한 목적함수로 로그 가능도(log-likelihood: LL)와 generic composite similarity measure (GCSM)를 사용하였다. 또한 품종모수의 평균값을 사용한 예측(E_pm)과 개별 모수들로부터 얻어진 추정값의 평균값(E_em)의 일치도를 분석하여 앙상블 방식에 따른 불확도 변화를 파악하였다. 국내에서 재배되는 사료용 벼 품종인 조우 벼와 영우 벼를 대상으로 품종모수를 추정하였다. 2013년, 2014년, 2016년에 대한 수원, 전주, 나주, 익산에 위치한 실험포장에서 얻은 수량 관측 자료를 사용하였다. 또한 2016년부터 2018년까지 수원에서 보고된 별도의 수량 관측 자료를 사용하였다. 목적함수에 따라 추정된 품종모수의 분포에 차이가 있었다. LL을 통해 얻은 품종모수는 GCSM으로 얻은 품종모수보다 좁은 범위에 분포하였다. 두 가지 앙상블 접근법은 통계적으로 유의한 차이가 나타나지 않음을 확인하였다. GCSM의 상대적으로 높은 불확도는 수용확률을 조정하여 낮출 수 있다고 사료되고, E_pm의 결과는 기존과 다른 앙상블 방식을 통해 적은 연산을 통해 불확도를 낮출 수 있음을 보인다.

• 응용통계연구
• /
• 제22권6호
• /
• pp.1355-1364
• /
• 2009

본 논문에서는 다변량 자료의 위치모수에 대한 로버스트 검정 방법으로 유사등변성과 고붕괴성을 만족하는 MVE와 MCD 추정량에 근거한 로버스트 검정방법을 제안하였다. 일반적으로 이들 추정방법은 낮은 효율성으로 인하여 통계적 추론보다는 잠재적 이상치의 발견과 같은 탐색적분석에서 사용된다. 우리는 검정력을 높이기 위하여 MVE와 MCD 추정량에 근거한 일단계 재가중절차를 사용했는데, 가중치 선정과 관련된 임계값을 조절함으로써 현실적으로 사용가능한 높은 효율성과 정확성을 갖춘 검정방법을 제시하였다. 모의실험 결과 본 연구에서 제안한 검정법은 모분포에 관계없이 모두 명목유의수준을 제대로 유지하고 검정력도 높게 나타났으며, 이상치를 포함하고 있는 사례를 이용하여 실제로 모평균에 대한 가설검정을 수행한 결과 기존 방법과는 달리 영향을 받지 않았다.

• Journal of Animal Science and Technology
• /
• 제45권5호
• /
• pp.689-694
• /
• 2003

출산기술연구소 대관령지소에서 출생한 한우 암소로부터 시간적인 간격을 두고 조사된 체중측정 기록에 대해 비선형의 성장곡선 모형을 적용하여 한우 암소의 개체별 성장곡선 모수를 추정하고 개체별로 추정된 모수를 이용하여 한우 암소의 정확한 성장형태를 추정하기 위해 실시하였다. Gompertz 모형으로 추정한 개체별 성장 곡선 모수 A, b 및 k의 평균치는 각각 383.42${\pm}$97.29kg, 2.374${\pm}$0.340 및 0.0037${\pm}$0.0012였으며, 개체별 변곡점의 평균은 255.63${\pm}$109.09일, 변곡점에서 체중의 평균은 141.05${\pm}$35.79kg, 변곡점에서의 일당증체량의 평균은 0.500${\pm}$0.123 kg이었고, 성장곡선 모수 A, b 및 k의 변이계수는 각각 25.3, 14.3 및 32.4%이었다. von Bertalanffy 모형으로 추정한 개체별 성장곡선 모수 A, b 및 k의 평균치는 각각 410.47${\pm}$117.98kg, 0.575${\pm}$0.057 및 0.003${\pm}$0.001이었고, 개체별 변곡점의 평균은 211.02${\pm}$105.53일, 변곡점에서 체중의 평균은 121.62${\pm}$34.94kg, 최대 증체율의 평균은 0.504${\pm}$0.124kg이었으며, 성장곡선 모수 A, b 및 k의 변이계수는 각각 28.7, 9.9 및 33.3%로 추정되었다. Logistic 모형으로 추정한 개체별 성장곡선 모수 A, b 및 k의 평균치는 각각 347.64${\pm}$97.29kg, 6.73${\pm}$0.34 및 0.006${\pm}$0.0018이었고, 개체별 변곡점의 평균은 324.47${\pm}$3.87일, 변곡점에서 체중의 평균은 173.82${\pm}$1.13kg, 최대 증체율의 평균은 0.508${\pm}$0.003kg이었으며, 성장곡선 모수 A, b 및 k의 변이계수는 각각 27.9, 5.0 및 30.0%이었다. Gompertz 모형, von Bertalanffy 모형 및 Logistic 모형에 의해 추정된 생시체중의 평균치들은 각각 36.40${\pm}$0.33, 31.59${\pm}$0.34 및 49.09${\pm}$0.88kg으로 추정되어 실제체중 23.73${\pm}$0.10kg 보다 컸으며 36개월 체중의 경우 세 모형이 각각 349.73${\pm}$2.01, 356.41${\pm}$2.07 및 336.18${\pm}$1.91kg으로 실제체중 363.67${\pm}$2.08kg보다 작았다. 그러나 세 모형으로 추정된 변곡점이 한우 암소의 실질적인 변곡점인지에 대해서는 좀더 검토가 필요한 부분인데, 최(2001)는 한우 암소의 경우 3개월령에서 4개월령 사이의 일당증체량이 0.590으로서 이 시기에 최대의 성장이 이루어지며 3개월령 체중과 4개월령 체중은 각각 79.31 및 98.91kg이라고 보고하고 있으며, 본 연구에서 Table 4에 제시된 결과를 봐도 3개월령과 6개월령 사이 90일 동안에 50.89kg의 증체가 이루어져 일당증체량으로 환산하면 약 0.56kg이 되어 다른 어느 시기보다도 가장 증체속도가 빠른 시기인 것으로 나타났다. 변곡점이 일당증체량이 최대인 시점으로 해석한다면 각 모형에서 고정되어 있는 변곡점의 위치들인 성숙체중의 38.6%(Gompertz 모형), 29.6% (von Bertalanffy 모형) 및 50%(Logistic 모형)는 한우 암소의 실질적인 변곡점보다 늦는 것으로 판단할 수 있는데, 실제로 세 모형 중 변곡점의 위치가 제일 빠른 von Bertalanffy 모형의 적합도가 제일 좋은 것도 이러한 이유 때문일 수 있다. 따라서 한우 암소의 정확한 성장특성을 파악하기 위해서는 변곡점의 위치가 고정되어 있지 않은 기존의 모형들을 이용하는 방안과 아울러 본 연구에 이용된 모형들을 변형시켜 활용하는 방안에 대한 연구가 필요한 것으로 판단된다.

• 한국산림과학회지
• /
• 제83권2호
• /
• pp.205-210
• /
• 1994

본 연구는 채종원에서 개화결실이 부진한 낙엽송 (Larix leptolepis S. et Z. Gordon)의 삽목증식법을 개발하기 위하여, 삽수의 모수령, 채취부위, 클론 및 발근촉진제가 삽목 발근에 미치는 영향을 구명하고저 실시하였다. 2, 8, 16, 30년생 모수에서 7월 하순경 당년에 생장한 삽수를 채취한 후, 발근촉진제인 IBA(Indole Butyric Acid)를 1000, 2000. 5000ppm의 세가지 농도로 5초간 처리하였다. 삽수는 peatmoss, vermiculite, perlite를 1 : 1 : 1 로 배합하여 조제된 온실내 삽목상에 삽목한 후 3개월이 경과된 10월 하순에 발근여부를 조사하였다. 발근촉진제 무처리의 경우 발근율은 2년생에서 52%, 8년생에서 48%, 16년생에서 36%, 30년생에서 20%인 반면, IBA 1000ppm 처리시 가장 높은 발근율을 얻었는데 2년생 80%, 8년생 71%, 16년생 52%, 30년생 25%로 발근율이 향상되었으며, 부정근의 수와 주근의 길이가 모수령이 낮을수록 증가하였으며, IBA 처리로도 증가하였다. 시기별 생존율은 삽목후 15-45일 사이에 현저히 감소한 반면, 그 후에는 생존율이 거의 감소하지 않았으며, IBA 처리로 생존율이 향상되었다. 모수의 클론별 발근율 조사는 8년생 모수를 사용하였는데, IBA 1000ppm 처리시 강원 49호에서 40%로 저조한 반면, 충남 6, 7호와 전북 1, 9호는 67%로 가장 커서 클론간에 발근율이 차이가 있었다. 삽수의 채취부위에 따른 위치효과는 수관 상부에서 채취한 삽수가 IBA 처리시 발근율이 평균 53%로 가장 낮았으며, 수관 하부가 56%, 수관 중부는 60%의 순으로 발근율이 증가되었다.

• 한국해안·해양공학회논문집
• /
• 제24권6호
• /
• pp.390-399
• /
• 2012

본 연구에서는 지속적인 담수유입이 존재하고 S자의 수로형태를 가진 경기만 석모수로에서 소조기 홍수시 13시간 동안 7개의 정선에서 단면유속 및 염분을 관측하였다. 각 단면의 최강 창조와 낙조시의 유속 크기와 방향을 파악하였고, 단면 내의 유속 분포 및 염분구조를 분석하였다. 또한 정선 별로 나타나는 횡 방향 흐름이 어떠한 운동량에 의해 지배되는지 파악하고자 횡 방향의 운동량 분석을 수행하였다. 운동량 분석에서는 석모수로의 S자 형태의 영향을 고려하기 위해서 원심력을 고려하였다. 분석 결과 소조기 홍수시 석모수로는 횡 방향 압력 경도력과 수직적 마찰력이 가장 우세하기 때문에 염분 분포와 성층에 의한 흐름이 주로 나타났다. 수로의 특성은 크게 북단과 남단으로 나누어 볼 수 있는데 상대적으로 조간대가 넓게 형성되고 담수의 영향이 큰 석모수로 북단의 4개 정선중에서 수심이 깊은 정선에서는 횡 방향 압력 경도력이 우세하지만 수직적으로 크기가 다르며, 수심이 낮은 정선에서는 수직 마찰항이 우세하였다. 이와는 달리 수심이 깊고 수로의 굴곡이 심한 석모수로 남단에서는 낙조시 지형학적 원인과 담수의 영향에 따라 이류 가속항과 원심력이 강해지게 된다. 이와 같은 결과를 종합할 때, 석모수로는 위치와 수심, 수로의 굴곡 등에 따라 운동량 분포가 각기 달리 나타나며 이러한 영향으로 인해 횡 방향 흐름 특성이 발생했음을 알 수 있다.

• 한국농림기상학회지
• /
• 제20권4호
• /
• pp.277-283
• /
• 2018

기후변화 영향평가를 위해 작물생육모형을 폭넓게 사용하고 있지만 모형을 구동하기 위해서는 품종모수를 결정하는 것이 필수적이다. 그러나 품종모수 결정을 위한 실험은 장시간의 노력이 필요하여 대부분 작황자료 또는 지역적응 시험 자료를 많이 사용하고 있다. 그러나 밭작물의 경우 작황자료 또는 지역적응 시험을 사용하는 경우에는 포장의 관개량과 시기에 대한 자료가 없고 또한 별도의 기상관측 없이 최인접지역의 기상자료를 사용하기 때문에 문제가 발생할 수 있다. CERES-MAIZE를 이용하여 밭작물인 옥수수에 대해서 이 문제점들을 검토하였다. 출사기와 관련된 품종모수는 최대 27km 내에 기상관측 지점이 있는 경우에도 신뢰성있는 품종모수가 얻어졌다. 온도의 경우에는 지형에 따라 유효한 거리가 달라질 수 있지만 본 연구의 대상 지역에서는 품종 모수 추정에 문제가 크지는 않을 것으로 보인다. 또한 온도 이외의 요소인 강수 또는 관개량은 생물계절관련 품종모수와는 큰 영향이 없기 때문에 비교적 정확도가 높은 결과가 나온 것으로 보인다. 그러나 수량과 관련된 품종 모수 요소에서는 그렇지 못하였다. 이는 밭작물의 경우 강수량에 따라 스트레스 정도가 결정되기 때문에 관개 및 강수량 정보가 중요한데, 관개량에 대한 정보를 작황 또는 지적시험 보고서에서는 얻을 수 없기 때문이다. 더구나 강수량의 경우에 온도보다 더 가까운 위치에 기상관측소가 존재해야만 유의미한 정보를 제공할 수 있다. 그러나 시험포장에 따라 기상관측소와의 거리가 충분히 가깝지 않은 경우가 대부분이었다. 따라서, 작황 또는 지역적응 시험자료를 이용하여 옥수수의 품종모수를 결정할 때는 기상관측 지점이 최소한 20km 이내의 인접지역에서 생물계절과 관련된 모수에 대해서만 결정하는 것이 타당할 것으로 생각된다. 그 반면에 수량과 관련된 요소의 결정은 적절하지 않을 것으로 생각된다. 수량과 관련된 요소를 결정하기 위해서는 가급적 직접 기상관측망을 설치하여 해당 포장에서 관개시기와 관개량을 모두 확보한 실험한 결과를 바탕으로 얻는 것이 적절할 것으로 생각된다.

• 정보과학회 컴퓨팅의 실제 논문지
• /
• 제22권11호
• /
• pp.567-575
• /
• 2016

스마트폰의 대중화로 이용자의 위치정보를 이용한 모바일 앱이 늘어나면서 위치정보에 대한 관심이 증가하고 있다. 본 논문에서는 위치정보를 이용한 사용자의 이동속도를 나타내는 간단한 분포함수를 찾고자 한다. 스마트폰에서 제공하는 위치정보는 경우에 따라 큰 오차가 발생할 수 있기 때문에 이를 제거하는 과정이 필요하다. 또한 이동속도 분포의 경우 교통수단에 따라 여러 가지 다른 분포들의 혼합분포로 표현할 수 있다. 속도의 분포함수를 찾기 위하여 혼합분포를 이용하여 꼬리부분에 해당하는 분포를 찾고 이부분을 설명할 수 있는 모수 분포모형을 찾는다.

• Journal of the Korean Data Analysis Society
• /
• 제20권6호
• /
• pp.2805-2812
• /
• 2018

본 연구는 생명정보학(bio-informatics) 분야 중, 특정 병에 관련된 유전자 위치를 찾고자 DNA 시퀀싱(DNA sequencing) 방법을 이용한 메틸화(methylation) 데이터의 분석에 관한 것이다. 반복적인 시퀀싱 과정을 통해 도출되는 메틸화 여부 자료를 비율로 표현한 메틸화 점수는 0과 1사이의 값을 가지게 된다. 이러한 데이터에 집단별 메틸화 점수의 차이를 검토하기 위해 t-검정을 단순히 적용하는 것은 정규분포의 가정에 위배된다. 또한 메틸화 점수 생성과정에서 시퀀싱의 반복수에 따라 결과가 달라 질 수 있으므로 이러한 오차를 고려해서 분석할 수 있는 방법도 필요하다. 이에 본 논문에서는 메틸화 데이터를 하나의 숫자 데이터가 아닌 불확실성을 포함하는 구간형(interval) 데이터로 변환하여 분석하는 심볼릭 데이터 분석(symbolic data analysis) 및 구간형 K-S 검정법을 적용하였다. 또한 구간형 데이터로 변환하는 과정에서 정규분포를 이용하지 않고 베타분포를 이용하여 메틸화 점수의 특성을 반영하여 분석할 수 있게 하였다. 자료분석을 위하여 174명의 실제 암환자 및 정상인들의 DNA 시퀀싱 데이터를 이용하여 제안한 방법의 성질을 살펴보았다. t-검정은 위치모수에 관한 검정만 가능한 반면, 구간형 K-S 통계량은 구간자료에 대해 위치모수뿐만 아니라 분포함수의 이질성에 검정할 수 있으므로 t-검정이 놓칠 수 있는 유의미한 유전자 위치를 찾아낼 수 있음을 확인하였다.

• Communications for Statistical Applications and Methods
• /
• 제16권3호
• /
• pp.397-408
• /
• 2009

마이크로어레이 기술은 수만 개 유전자의 발현 패턴을 동시에 관찰하는 것을 가능하게 하였고, 이들을 하나씩 검정하여 찾아낸 특이발현 현상을 보이는 유전자를 중심으로 질병의 진단, 치료법 정립과 신약 개발을 위한 기본 정보를 확립하였다. 그러나 개별 유전자분석의 여러 문제점이 발견되면서 유전자들을 생물학적 대사경로나 염색체 위치가 같은 것끼리 묶은 집단을 분석하여 질병의 발생이나 생존에 영향을 미치는 집단을 찾는 방법이 제시되었다. 이러한 유전자 집단의 유의성에 대한 연구는 2002년에 MIT에서 비롯되어 GSEA, SAM-GS와 중심극한 정리의 개념을 이용한 모수적 방법인 PAGE 등이 사용되고 있다. 본 논문에서는 이들 통계량의 구조적 한계를 극복하고 계산이 간단한 새로운 모수적 방법을 제안하고 자료 분석을 통하여 효율성을 보였다.

• Journal of the Korean Data and Information Science Society
• /
• 제17권3호
• /
• pp.707-716
• /
• 2006

분산함수는 회귀함수와 더불어 회귀모형의 연구에 매우 중요한 함수이며 이 함수가 불연속일 때의 연구는 Delgado and Hidalgo (2000)와 Perron (2001)은 시계열모형에서는 비모수적 추정법에 의해 분산함수의 추정을 연구하였으며 Kang and Huh (2006)은 Perron의 추정법을 회귀모형에 적용하여 분산함수의 불연속점의 추정에 대하여 연구하였고, Huh (2005)는 Kang and Huh의 잔차제곱들을 이용한 분산함수의 불연속점의 추정 대신 이차적률함수를 이용하여 분산함수의 불연속점을 추정하였다. 이는 Kang and Huh의 연구에서 잔차제곱들을 구하기 위하여 회귀함수의 추정이 우선되어야 하기에 전체적인 계산량이 늘어나게 되고, 늘어난 만큼 불연속점 추정의 정도가 떨어지게 됨으로 반응변수의 표본의 제곱을 이용하여 이차적률함수의 추정으로 불연속점을 추정하는 것이 더 용이하기 때문이다. 이러한 연구를 바탕으로 본 연구에서는 Huh의 점프의 크기 추정량의 점근분포를 이용하여 불연속점의 존재 유무에 대한 가설검정법을 제안하였다. 즉, 점프의 크기 추정량의 귀무가설 하의 점근분포가 가지고 있는 장애모수인 불연속점의 위치에서 확률밀도함수와 4차적률함수를 비모수적 방법으로 추정하는 방법을 제안하고 이들의 균일 일치성을 보여 가설검정법을 제안하였다. 불연속점의 추정에 앞서 불연속점의 존재 여부의 가설검정이 우선되어야 하기에 다른 통계적 함수에 대한 불연속점의 연구에서도 이러한 본 논문에서 연구한 방법으로 불연속점의 존재 유무에 대한 가설검정법을 제안 할 수 있을 것이다.