• The Journal of the Korea Contents Association
• /
• v.15 no.4
• /
• pp.350-356
• /
• 2015

The aim of this study is to perform simulations to design the detector module with crystal array by Monte Carlo simulation. For this purpose, a small animal PET scanner, employing module with 1~8 crystal array discrimination scheme, was designed. The proposed scanner has an inner diameter of 100 mm with detector modules in crystal array. Each module is composed of a 5.0 mm LSO crystal with a $2.0{\times}2.0mm^2$ sensitive area with a pitch 2.1 mm and 10.0 mm thickness. The LSO crystals are attached to the SiPM which has a dimension of $2.0{\times}2.0mm^2$. The detector module with crystal array of the designed PET detector was simulated using the Monte Carlo code GATE(Geant4 Application for Tomographic Emission). The detector is enough compensation for the loss of data in sinogram due to gaps between modules. The results showed that the high sensitivity and effectively reduced the problem about the missing data were greatly improved by using the detector module with 1 crystal array.

• Progress in Medical Physics
• /
• v.14 no.4
• /
• pp.240-248
• /
• 2003

The Monte Carlo calculation is the most accurate means of predicting radiation dose, but its accuracy is accompanied by an increase in the amount of time required to produce a statistically meaningful dose distribution. In this study, the effects on calculation time by introducing variance reduction techniques and increasing computing power, respectively, in the Monte Carlo dose calculation for a 6 MV photon beam from the Varian 600 C/D were estimated when maintaining accuracy of the Monte Carlo calculation results. The EGSnrc­based BEAMnrc code was used to simulate the beam and the EGSnrc­based DOSXYZnrc code to calculate dose distributions. Variance reduction techniques in the codes were used to describe reduced­physics, and a computer cluster consisting of ten PCs was built to execute parallel computing. As a result, time was more reduced by the use of variance reduction techniques than that by the increase of computing power. Because the use of the Monte Carlo dose calculation in clinical practice is yet limited by reducing the computational time only through improvements in computing power, introduction of reduced­physics into the Monte Carlo calculation is inevitable at this point. Therefore, a more active investigation of existing or new reduced­physics approaches is required.

• Proceedings of the Korean Statistical Society Conference
• /
• 2000.11a
• /
• pp.135-140
• /
• 2000

베이지안 신경망 모형(Bayesian Neural Networks Models)에서 주어진 입력값(input)은 블랙 박스(Black-Box)와 같은 신경망 구조의 각 층(layer)을 거쳐서 출력값(output)으로 계산된다. 새로운 입력 데이터에 대한 예측값은 사후분포(posterior distribution)의 기대값(mean)에 의해 계산된다. 주어진 사전분포(prior distribution)와 학습데이터에 의한 가능도함수(likelihood functions)를 통해 계산되어진 사후분포는 매우 복잡한 구조를 갖게 됨으로서 기대값의 적분계산에 대한 어려움이 발생한다. 이때 확률적 추정에 의한 근사 방법인 몬테칼로 적분을 이용한다. 이러한 방법으로서 Hybrid Monte Carlo 알고리즘은 우수한 결과를 제공하여준다(Neal 1996). 본 논문에서는 Hybrid Monte Carlo 알고리즘과 기존에 많이 사용되고 있는 Gibbs sampling, Metropolis algorithm, 그리고 Slice Sampling등의 몬테칼로 방법들을 비교한다.

• The Korean Journal of Applied Statistics
• /
• v.13 no.2
• /
• pp.505-514
• /
• 2000

A Markov Chain Monte Carlo method with data augmentation is developed to compute the features of the posterior distribution. For each observed failure epoch, we introduced mixture failure model of Rayleigh and Erlang(2) pattern. This data augmentation approach facilitates specification of the transitional measure in the Markov Chain. Gibbs steps are proposed to perform the Bayesian inference of such models. For model determination, we explored sum of relative error criterion that selects the best model. A numerical example with simulated data set is given.

• Journal of Korea Society of Industrial Information Systems
• /
• v.23 no.3
• /
• pp.49-57
• /
• 2018

In a previous article, for analyzing a stochastic activity network, Cho proposed a method for estimating the moments (mean, variance, skewness, kurtosis) of the project completion time under the assumption that the durations of activities are independently and normally distributed. Developed in the present article is a method for estimating those moments for stochastic activity networks which allow any type of distributions for activity durations. The proposed method uses the moment matching approach to discretize the distribution function of activity duration, and then a discrete inverse-transform method to determine activity durations to be used for calculating the project completion time. The proposed method can be easily applied to large-sized activity networks, and computationally more efficient than Monte Carlo simulation, and its accuracy is comparable to that of Monte Carlo simulation.

• The Korean Journal of Applied Statistics
• /
• v.4 no.1
• /
• pp.13-24
• /
• 1991

This paper presents a Monte Carlo evaluation of estimators for nonlinear consored simultaneous equations models. We examine the performance of the maximum likelihood estimator (MLE), the two-step quasi maximum likelihood estimator (2QMLE) proposed by Lee and Hurd (1989), and another quasi MLe using least squares at the first step (LSAE) under varying degrees of freedom and underlying distributions, Although QMLE's are not necessarily consistent, the Monte Carlo results show that the 2QMLE may be used as an alternative to MLE when MLE is not applicable in practice.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2001.05a
• /
• pp.436.2-436
• /
• 2001

• Proceedings of the Korean Nuclear Society Conference
• /
• 2001.10a
• /
• pp.299.2-299
• /
• 2001

• Proceedings of the Korean Nuclear Society Conference
• /
• 2003.10a
• /
• pp.447-447
• /
• 2003

• Progress in Medical Physics
• /
• v.21 no.3
• /
• pp.253-260
• /
• 2010

Most brachytherapy treatment planning systems employ a dosimetry formalism based on the AAPM TG-43 report which does not appropriately consider tissue heterogeneity. In this study we aimed to set up a simple Monte Carlo-based intracavitary high-dose-rate brachytherapy (IC-HDRB) plan verification platform, focusing particularly on the robustness of the direct Monte Carlo dose calculation using material and density information derived from CT images. CT images of slab phantoms and a uterine cervical cancer patient were used for brachytherapy plans based on the Plato (Nucletron, Netherlands) brachytherapy planning system. Monte Carlo simulations were implemented using the parameters from the Plato system and compared with the EBT film dosimetry and conventional dose computations. EGSnrc based DOSXYZnrc code was used for Monte Carlo simulations. Each $^{192}Ir$ source of the afterloader was approximately modeled as a parallel-piped shape inside the converted CT data set whose voxel size was $2{\times}2{\times}2\;mm^3$. Bracytherapy dose calculations based on the TG-43 showed good agreement with the Monte Carlo results in a homogeneous media whose density was close to water, but there were significant errors in high-density materials. For a patient case, A and B point dose differences were less than 3%, while the mean dose discrepancy was as much as 5%. Conventional dose computation methods might underdose the targets by not accounting for the effects of high-density materials. The proposed platform was shown to be feasible and to have good dose calculation accuracy. One should be careful when confirming the plan using a conventional brachytherapy dose computation method, and moreover, an independent dose verification system as developed in this study might be helpful.