• 응용통계연구
• /
• 제30권5호
• /
• pp.775-791
• /
• 2017

이진 데이터는 일상 생활에서 자주 접할 수 있는 데이터이다. 이진 데이터를 회귀 분석하는 방법으로 로지스틱(Logistic), 프로빗(Probit), Cauchit, Complementary log-log 모형이 주로 쓰이는데, 이 방법 이외에도 Liu(2004)가 제시한 t 분포를 이용한 로빗(Robit) 모형, Kim 등 (2008)에서 제시한 일반화 t-link 모형을 이용한 방법 등이 있다. 유연한 분포를 이용하면 유연한 회귀 모형이 가능해지는 점에 착안하여, 이 논문에서는 Theodossiou(1998)에서 제시된 기운 일반화 t 분포 (Skewed Generalized t Distribution)의 이용하여 우도 함수를 최대로 하는 이진 데이터 회귀 모형을 소개한다. 기운 일반화 t 분포를 R glm 함수, R sgt 패키지를 연결하여 이 논문에서 제시한 방법을 R로 분석할 수 있는 방법을 소개하고, 피마 인디언(Pima Indian) 데이터를 분석한다.

• Communications for Statistical Applications and Methods
• /
• 제26권3호
• /
• pp.315-323
• /
• 2019

Panel data sets have recently been developed in various areas, and many recent studies have analyzed panel, or longitudinal data sets. Often a dichotomous dependent variable occur in survival analysis, biomedical and epidemiological studies that is analyzed by a generalized linear mixed effects model (GLMM). The most common estimation method for the binary panel data may be the maximum likelihood (ML). Many statistical packages provide ML estimates; however, the estimates are computed from numerically approximated likelihood function. For instance, R packages, pglm (Croissant, 2017) approximate the likelihood function by the Gauss-Hermite quadratures, while Rchoice (Sarrias, Journal of Statistical Software, 74, 1-31, 2016) use a Monte Carlo integration method for the approximation. As a result, it can be observed that different packages give different results because of different numerical computation methods. In this note, we discuss the pros and cons of numerical methods compared with the exact computation method.

• 응용통계연구
• /
• 제33권1호
• /
• pp.47-59
• /
• 2020

본 논문에서는 마코프 이항 회귀 모형의 시차가 알려져 있거나 그렇지 않은 경우일 때, t-링크 함수를 갖는 종단적 마코프 이항 회귀 모형을 제시한다. 일반적으로, 이항 회귀 모형에서는 로직 모형이나 프로빗 모형이 주로 사용된다. t-링크 함수는 t 분포가 자유도가 커질수록 정규분포로 근사하기 때문에 프로빗 모형을 대신 더 많은 유연성을 위해 사용될 수 있다. 게다가 마코프 회귀모형은 종단 자료에 대해 사용될 수 있다. 우리는 마코프 회귀 모형의 시차를 결정하기 위해 베이지안 방법을 제시하고자 한다. 특히, 각 모델의 차수에 대해 알고 있는 경우에는 DIC를 기준으로 모델 비교를 실시하였다. 모델의 차수에 대해 모르는 경우에는 가능한 모델들의 사후 확률을 이용하였다. 복잡한 베이지안 계산을 해결하기 위하여 Albert와 Chib (1993), Kuo와 Mallick (1998)과 Erkanli 등 (2001)의 방법을 이용하여 모델을 재설정하였다. 제안하는 방법은 시뮬레이션 데이터와 Somer 등 (1984)에 의해 조사된 인도네시아 어린이 종단 데이터에 적용했다. 마코프 이항 회귀모형의 순서에 대해서 아는 경우와 모르는 경우를 각각 가정하여 최적의 모델을 알아보기 위해 MCMC 방법을 사용하였다. 또한, 매트로폴리스 해스팅 알고리즘의 수렴성을 점검하기 위해 Gelman과 Rubin의 진단을 이용했다.

• 터널과지하공간
• /
• 제31권1호
• /
• pp.52-65
• /
• 2021

본 연구에서는 멕시코시티 외곽에 있는 프로판 저장기지인 PEMIX 터미널에서 발생한 프로판 누출에 따른 증기운 폭발을 분석하였다. 누출된 4750 kg의 프로판에 대한 TNT 등가량은 9398 kg으로 평가되었다. 폭원으로부터 40~400 (m) 떨어진 지점에서의 최대과압, 양의 압력 지속시간, 충격량과 같은 폭발변수를 TNT 등가법과 다중에너지법을 적용하여 구하였다. 폭발 변수들을 이용하여 구한 프로빗 함수를 적용하여 폐 손상, 고막 파열, 머리 충격, 전신 전위 충격으로 인한 손상 확률을 평가하였다. 고려한 모든 거리에서 다중에너지법을 이용하여 구한 최대과압이 TNT 등가법을 적용하여 구한 최대과압보다 큰 것으로 나타났으나, 200 m 이후 지점부터는 큰 차이가 없는 것으로 평가되었다. 다중에너지법에 의해 구해진 최대과압을 적용하여 구조물 손상 범위를 평가한 결과 폭원으로부터 100 m 이내에 있는 구조물의 경우 완전히 붕괴될 것으로 예측되고, 400 m 떨어진 구조물의 유리창도 거의 파손될 것으로 추정되었다. 폐 손상에 의한 사망 확률은 충격파 진행방향으로 위치하고 있는 인체의 자세에 따라 달라지는 것으로 나타났으며, 인체 주변에 반사면이 있는 경우 사망 확률이 가장 큰 것으로 평가되었다. 충격파가 폐 손상, 고막 파열, 머리 충격, 전신 전위 충격에 미치는 영향을 평가한 결과 전신 전위 충격 < 폐 손상 < 고막파열 < 머리 충격 순으로 영향을 미치는 것으로 나타났다.

• 화약ㆍ발파
• /
• 제41권4호
• /
• pp.1-8
• /
• 2023

2020년 8월 4일에 베이루트항 저장 창고에 저장되었던 2750 ton의 질산암모늄이 폭발하였다. 이 폭발은 지금까지의 질산암모늄 폭발 중 가장 규모가 큰 것으로 알려졌다. TNT 등가법을 적용하여 2750 ton의 질산암모늄의 폭발 에너지에 상응하는 TNT 등가량을 구한 결과 856 ton으로 나타났다. Kingery-Bulmash 폭발 특성 계산기 툴을 활용하여 폭원으로부터 3600 m 까지의 범위에서 과압과 충격량을 산정하였다. 폭원으로부터 멀어짐에 따라 과압과 충격량은 지수적으로 감소하지만 과압이 더 크게 감소하여 과압이 충격량보다 거리에 따른 영향을 더 크게 받는 것으로 나타났다. 과압과 충격량이 구조물에 미치는 영향을 평가하기 위해서 구조물의 손상 기준을 적용한 결과 구조물의 부분적 붕괴, 심각한 손상, 가벼운 손상이 발생하는 임계거리는 폭원으로부터 각기 약 500, 800, 2200 m로 나타났다. 구조물과 인체의 손상 확률을 평가하기 위해서 프로빗 함수를 적용하였다. 구조물의 붕괴, 구조물의 심각한 손상, 구조물의 가벼운 손상, 창유리의 파손 가능성이 50% 이상이 되는 지점은 각기 약 500, 810, 2200, 3200 m가 되는 것으로 나타났다. 폭원으로부터 200 m 이내 지점에 있는 사람의 경우 폐 손상으로 인해 사망할 확률이 99% 이상인 것으로, 고막 파열이 발생할 확률이 50%인 지점은 약 300 m인 것으로 나타났다. 전신 이동에 따른 두개골 파열과 전신 충격에 의해 사망할 확률이 100%인 지점은 각기 300, 100 m인 것으로 평가되었다.

• 자원ㆍ환경경제연구
• /
• 제20권1호
• /
• pp.119-153
• /
• 2011

본 연구는 조건부가치측정법(CVM)의 지불의사유도 방법으로 가장 많이 적용되고 있는 양분선택형 응답자료를 활용하여 WTP 함수의 추정과 편익의 측정단계에 쟁점이 될 수 있는 경제이론적 근거, 확률분포의 가정, 확률모형의 표기, 함수형태 그리고 편익의 범위 등 분석자의 입장에서 중요한 제요소들에 대해 이론적으로 검토하고 산림보존을 위한 가치추정 자료를 사용하여 실증적으로 분석하였다. WTP 함수를 선형으로 표기하는 한, 확률효용모형과 확률지출모형은 쌍대적 접근으로 해석할 수 있었고, 또한 확률모형의 오차항의 분포를 정규분포로 하든지 로지스틱분포로 가정하든지 계수들의 추정치와 이들을 사용하여 측정된 표본 WTP도 차이가 없었다. 그러나 이중양분선택모형은 WTP 함수 추정과 편익측정 결과는 두 CV 제시금액에 대한 응답이 같은 WTP 분포로 왔다고 보기 어려워 단일양분선택형 모형이 갖는 유인일치성이 약화되었다. 나아가서 WTP를 비음의 영역으로 제한하는 WTP 함수형태의 경우 척도모수가 클 경우 표본 WTP의 평균과 중앙값이 크게 다른 것으로 나타났다.

• 대한가정학회지
• /
• 제38권9호
• /
• pp.143-157
• /
• 2000

The purpose of this study was to show the possibility of labor mobility among married women by analysing some variables to have effects on their employment status and willingness to move in the labo market, namely, whether they would willingly enter into/come out the labor market. The personal collecting data were used and one discirminant function and two probit models were employed to analyse the efficients to the independent variables. The resets of this study were as follows. The very discriminant variables to have significant effects on the employment status of married women were the age of the youngest and the number of children. The willingness to enter into the labor market of unemployed married women was high in contrast as the willingness to keep the job was high among employed married women. The subjective perception on their own health, children related variables, and financial variable were very significant to predict the possiblities to move in the labor market among married women. It was remarkable that the psychological variables, which were related with the attitudes toward the occupation of house wife and employment of married women, were also important. These results showed that first of all, the decision of whether they enter into or come out in the labor market among married women is largely on the basis of the household situation not on the personal capacities of married women themselves.

• 대한교통학회지
• /
• 제29권3호
• /
• pp.103-114
• /
• 2011

본 연구는 2000년 이후 대형교통사고 발생지점 112건의 자료를 이용, 다양한 교차 및 빈도분석을 통해 대형교통사고와 도로 기하구조의 관계를 규명하고, 이를 토대로 대형교통사고 심각도 모형을 구축하였으며, 주행안전성 향상을 위해 720회의 컴퓨터 모의실험으로 다음과 같은 결론을 얻을 수 있었다. 첫째, 교차 및 빈도분석의 결과 커브구간에서 43.7%, 종단경사 기타조건에서 60.7%, 곡선반경 0~24m 구간에서 57.2%, 편경사 0.1~2.0% 구간에서 83.9%, 편도2차로 도로에서 49.1%, 차종별로는 승용(33.0%), 화물(20.5%), 버스(14.3%) 순이었으며, 편경사 설치 유 무가 대형교통사고 발생에 가장 큰 영향을 주는 것으로 분석되었다. 둘째, 순서형 프로빗 모형(Ordered Probit Model)을 이용하여 다양한 도로조건에서의 피해 예측이 가능한 대형교통사고 심각도 모형을 개발하였으며, 개발된 모형을 기반으로 도로의 위험성을 사전에 예측하고, 대책 마련이 가능토록 기여 하였다. 셋째, 컴퓨터 모의실험(Simulation) 결과, 이미 대형교통사고가 발생한 장소에 편경사를 설치했을 경우 약 85% 이상의 지점들에서 대형교통사고가 발생하지 않는 개선효과가 있는 것으로 분석되었으며, 이 분석결과를 통해 도로의 구조 시설 기준에 관한 규칙(해설 및 지침)의 편경사 설치 예외규정을 더욱더 강화시킬 필요가 있다고 사료된다.

• 한국산림과학회지
• /
• 제88권1호
• /
• pp.117-123
• /
• 1999

잉여가치는 일정한 시간범위에 대한 가치로 표현된다. 특히, 가상적시장에서 시간범위의 선택문제는 매우 중요하다. 본 연구에서는 상이한 기간에 따른 지불의사액의 차이원인을 할인율을 통해 살펴보고자 한다. 이를 위해 속리산 국립공원을 방문객을 대상으로 1년과 5년 두 종류의 시간범위에 대한 방문객의 휴양가치를 이용하였다. 지불의사액에 대한 질문은 양분선택형을 택하였으며 분석에서는 프로빗모형을 사용하였다. 상한치 절단에 의한 편의의 문제점을 회피하기 위하여 휴양가치에 대한 기대치는 중앙값으로 평가하였다. 분석 결과, 방문객의 1년과 5년의 휴양가치는 각각 16,569원과 27,111원으로 산출되었다. 이러한 상이한 시간범위에 따라 평가된 가치의 차이를 일치시키는 할인율은 년(年) 153%로 산출되었다. 이러한 높은 할인율은 (1)미래방문에 대한 불확실로 인한 현재가치에 대한 높은 선호 (2)시간범위와 비례하는 대체가능성 (3)시간범위에 대한 응답자의 인지부족 (4)재화의 가격비율과 효용함수에 의한 영향 (5)지불수단에 대한 기본적인 프레미엄으로 인해 1년에 대한 휴양가치가 과대 추정 등, 복합적인 원인으로 해석 가능하다.

• 한국산학기술학회논문지
• /
• 제22권6호
• /
• pp.197-206
• /
• 2021

우수한 연구성과를 상업화하여 산업발전에 기여하는 '기업가적 행위자(entrepreneurial actor)'로서의 공공 연구기관의 기능과 역할에 대한 관심이 높아지고 있다. 그러나 이들의 상용화 단계 성과는 미흡하다는 지적이 제기되고 있다. 공공기술의 기업이 활용하기에 기술성숙도 또는 재현성(repeatability)이 미흡하여 공공 연구기관으로부터 개발된 기술의 사업화가 성공적이지 못하다는 주장이 있다. 본 연구에서는 공공기술을 도입한 기업을 대상으로 조사한 자료를 활용하여 기술사업화 성공에 대한 기술성숙도(TRL) 및 사업화 활동의 영향을 실증적으로 분석한다. 이를 위해 대학 및 출연(연)으로부터 기술을 이전받은 369개 기업의 설문응답 결과를 활용하여 프로빗 모형 분석을 실시하였으며, 특히 도입시점의 TRL 수준과 현재시점 TRL 수준의 차이(TRL Gap)를 측정하여 이의 조절효과를 검증하고자 하였다. 연구 결과, 기업이 도입한 공공기술의 도입시점 TRL이 높을수록 기술사업화 성공 가능성 또한 높은 것으로 나타났다. 또한, 도입시점의 TRL 수준과 사업화 시점의 TRL 수준의 차이(TRL Gap)는 그 자체로 기술사업화 성공에 영향을 미치지 않으며, 도입시점 TRL과 결합하여 추가적인 효과를 발생시키는 것으로 나타났다. 마지막으로, 한계효과 추정을 통해 TRL이 4~6 수준인 기술을 도입하여 TRL 9 수준으로 향상시키면 기술사업화 성공 가능성이 가장 높게 나타나는 것으로 관측되었다.