In irradiation process, instead of traverse on the targeted cells, there is side
effect happens to non-targeted cells. The targeted cells that had been irradiated with
ionizing radiation emits damaging signal molecules to the surrounding and then, dam-
age the bystander cells. The type of damage considered in this work is the number of
double-strand breaks (DSBs) of deoxyribonucleic acid (DNA) in cell’s nucleus. By us-
ing mathematical approach, a mechanistic model that can describe this phenomenon is
developed based on a structured population approach. Then, the accuracy of the model
is validated by its ability to match the experimental data. The Particle Swarm (PS)
optimization is employed for the data fitting procedure. PS optimization searches the
parameter value that minimize the errors between the model simulation data and exper-
imental data. It is obtained that the mathematical modelling proposed in this paper is
strongly in line with the experimental data.
Symmetric methods such as the implicit midpoint rule (IMR), implicit trapezoidal
rule (ITR) and 2-stage Gauss method are beneficial in solving Hamiltonian problems
since they are also symplectic. Symplectic methods have advantages over non-symplectic
methods in the long term integration of Hamiltonian problems. The study is to show
the efficiency of IMR, ITR and the 2-stage Gauss method in solving simple harmonic
oscillators (SHO). This study is done theoretically and numerically on the simple harmonic
oscillator problem. The theoretical analysis and numerical results on SHO problem
showed that the magnitude of the global error for a symmetric or symplectic method
with stepsize h is linearly dependent on time t. This gives the linear error growth when
a symmetric or symplectic method is applied to the simple harmonic oscillator problem.
Passive and active extrapolations have been implemented to improve the accuracy of the
numerical solutions. Passive extrapolation is observed to show quadratic error growth
after a very short period of time. On the other hand, active extrapolation is observed to
show linear error growth for a much longer period of time.
Homotopy continuation methods (HCMs) can be used to find the solutions
of polynomial equations. The advantages of HCMs over classical methods such as the
Newton and bisection methods are that HCMs are able to resolve divergence and starting
value problems. In this paper, we develop Super Ostrowski-HCM as a technique to
overcome the starting value problem. We compare the performance of this proposed
method with Ostrowski-HCM. The results provide evidence of the superiority of Super
Ostrowski-HCM.
A mechanistic model has been used to explain the effect of radiation. The
model consists of parameters which represent the biological process following ionizing
radiation. The parameters in the model are estimated using local and global optimiza-
tion algorithms. The aim of this study is to compare the efficiency between local and
global optimization method, which is Pattern Search and Genetic Algorithm respectively.
Experimental data from the cell survival of irradiated HeLa cell line is used to find the
minimum value of the sum of squared error (SSE) between experimental data and sim-
ulation data from the model. The performance of both methods are compared based on
the computational time and the value of the objective function, SSE. The optimization
process is carried out by using the built-in function in MATLAB software. The parameter
estimation results show that genetic algorithm is more superior than pattern search for
this problem.
Abstract Demographers and actuaries are very much conscious of the trend of mortality in their own country or in the world in general. This is because mortality is the basis for longevity risk evaluation. Mortality is showing a declining trend and it is expected to further decline in the future. This will lead to continuous increase in life expectancy. Several stochastic models have been developed throughout the years to capture mortality and its variability. This includes Lee Carter (LC) model which has been extended by various researchers. This paper will be focusing on comparing LC model and another mortality model proposed by Cairns, Blake and Dowd (CBD). The LC uses the log of central rate of mortality and CBD uses logit of the mortality odds as dependent variable. Analysis of comparison is done using a few techniques including Akaike information criteria (AIC) and Bayesian information criterion (BIC). From the overall results, there is no model better than the other in every aspect tested. We illustrate this via visual inspection and in sample and outof sample analysis using Malaysian mortality data from 1980 to 2017.
In this paper, we consider the system of Volterra-Fredholm integral equations
of the second kind (SVFI-2). We proposed fixed point method (FPM) to solve
SVFI-2 and improved fixed point method (IFPM) for solving the problem. In addition,
a few theorems and two new algorithms are introduced. They are supported by
numerical examples and simulations using Matlab. The results are reasonably good
when compared with the exact solutions.
In this paper, we look at the propagation of internal solitary waves over three
different types of slowly varying region, i.e. a slowly increasing slope, a smooth bump and
a parabolic mound in a two-layer fluid flow. The appropriate mathematical model for this
problem is the variable-coefficient extended Korteweg-de Vries equation. The governing
equation is then solved numerically using the method of lines. Our numerical simulations
show that the internal solitary waves deforms adiabatically on the slowly increasing slope.
At the same time, a trailing shelf is generated as the internal solitary wave propagates
over the slope, which would then decompose into secondary solitary waves or a wavetrain.
On the other hand, when internal solitary waves propagate over a smooth bump or a
parabolic mound, a trailing shelf of negative polarity would be generated as the results of
the interaction of the internal solitary wave with the decreasing slope of the bump or the
parabolic mound. The secondary solitary waves is observed to be climbing the negative
trailing shelf.
In this work, a non-abelian metabelian group is represented by G while represents conjugacy class graph. Conjugacy class graph of a group is that graph associated with the conjugacy classes of the group. Its vertices are the non-central conjugacy classes of the group, and two distinct vertices are joined by an edge if their cardinalities are not coprime. A group is referred to as metabelian if there exits an abelian normal subgroup in which the factor group is also abelian. It has been proven earlier that 25 non-abelian metabelian groups which have order less than 24, which are considered in this work, exist. In this article, the conjugacy class graphs of non-abelian metabelian groups of order less than 24 are determined as well as examples of some finite groups associated to other graphs are given.
The selection criteria play an important role in the portfolio optimization
using any ratio model. In this paper, the authors have considered the mean return as
profit and variance of return as risk on the asset return as selection criteria, as the first
stage to optimize the selected portfolio. Furthermore, the sharp ratio (SR) has been
considered to be the optimization ratio model. In this regard, the historical data taken
from Shanghai Stock Exchange (SSE) has been considered. A metaheuristic technique
has been developed, with financial tool box available in MATLAB and the particle swarm
optimization (PSO) algorithm. Hence, called as the hybrid particle swarm optimization
(HPSO) or can also be called as financial tool box particle swarm optimization (FTBPSO).
In this model, the budgets as constraint, where as two different models i.e. with
and without short sale, have been considered. The obtained results have been compared
with the existing literature and the proposed technique is found to be optimum and better
in terms of profit.
In this paper, Maxwell fluid over a flat plate for convective boundary layer
flow with pressure gradient parameter is considered. The aim of this study is to compare
and analyze the effects of the presence and absence of λ (relaxation time), and also the
effects of m (pressure gradient parameter) and Pr (Prandtl number)on the momentum
and thermal boundary layer thicknesses. An approximation technique namely Homotopy
Perturbation Method (HPM) has been used with an implementation of Adam and Gear
Method’s algorithms. The obtained results have been compared for zero relaxation time
and also pressure gradient parameter with the published work of Fathizadeh and Rashidi.
The current outcomes are found to be in good agreement with the published results.
Physical interpretations have been given for the effects of the m, Pr and β (Deborah
number) with λ. This study will play an important role in industrial and engineering
applications.
In this paper we consider a harvesting model of predator-prey fishery in which
the prey is directly infected by some external toxic substances. The toxic infection is
indirectly transmitted to the predator during the feeding process. The model is a modified
version from the classic Lotka-Volterra predator-prey model. The stability and bifurcation
analyses are addressed. Numerical simulations of the model are performed and bifurcation
diagrams are studied to investigate the dynamical behaviours between the predator and
the prey. The effects of toxicity and harvesting on the stability of steady states found in
the model are discussed.
Optimization is central to any problem involving decision making. The area
of optimization has received enormous attention for over 30 years and it is still popular
in research field to this day. In this paper, a global optimization method called Improved
Homotopy with 2-Step Predictor-corrector Method will be introduced. The method in-
troduced is able to identify all local solutions by converting non-convex optimization
problems into piece-wise convex optimization problems. A mechanism which only consid-
ers the convex part where minimizers existed on a function is applied. This mechanism
allows the method to filter out concave parts and some unrelated parts automatically.
The identified convex parts are called trusted intervals. The descent property and the
global convergence of the method was shown in this paper. 15 test problems have been
used to show the ability of the algorithm proposed in locating global minimizer.
Prediction analysis has drawn significant interest in numerous field. Taguchi’s T-Method is a prediction tool that developed practically but not limited to small sample analysis. It was developed explicitly for multidimensional system prediction by relying on historical data as the baseline model and adapting the signal to noise ratio (SNR) as well as zero proportional concepts in strengthening its robustness. Orthogonal array (OA) in T-Method is a variable selection optimization technique in improving the prediction accuracy as well as help in eliminating variables that may deteriorate the overall performance. However, the limitation of OA in dealing with higher multidimensionality restraint the optimization accuracy. Binary particle swarm optimization used in this study helps to cater to the limitation of OA as well as optimizing the variable selection process to better prediction accuracy. The results show that if the historical data consist of samples with higher correlation of determination (R2) value for the model creation, the optimization process in reducing the number of variables would be much reliable and accurate. Comparing between T-Method+OA and T-Method+BPSO in four different case study, it shows that T-Method+BPSO performing better with greater R2 and means relative error (MRE) value compared to T-Method+OA.
The effect of oil shock on the global economy is evident through many studies. However, the effect is heterogeneous over time. One of the reasons that lead to such different impacts is due to the oil source that is either the oil shock is demand or supply- driven. Applying the structural vector autoregressive (SVAR) model to generate the three oil shocks based on the three oil sources (oil supply, oil demand and oil specific- demand), we extended the examination on the effect of oil shock on the global economy using the threshold regression. Our results reveal the threshold effects of oil directly and indirectly on the global economy. The impacts of oil shocks differ across sectors, implying oil intensity, as well as oil sources, are the factors that determine the impact of oil shocks on the global economy. Overall, the oil specific-demand shock is more influential among the three oil shocks. Hence, the global economy is oil demand-driven. Besides that, the impact of oil is relatively large in the energy sector when compared to the non-energy sector and precious metals industry. Despite that, the impact of oil shocks is small if compared to the non-oil shocks such as exchange rate changes and global consumer price inflation shock. Consequently, non-oil shocks are the main determinants of the global economic fluctuation. The study leads to a better understanding of the transmission of oil shock and its sources, the interaction between oil and economic indicators and the policy implication due to oil dependency/ intensity.
1Malaysia
2 (UKM)43600 Bangi, Selangor, Malaysia
∗Corresponding author:
Numerous studies have linked biodiversity with zoonotic disease control. However, researchers have warned against simply believing that the increase in biodiversity can reduce the infection disease in the community. They proposed that amplification effect (increase in biodiversity accompanied by an increase in disease prevalence) might sometimes occur. Thus, we formulated a deterministic model to consider the impact of an amplification or dilution agent on the SNV transmission in the deer mouse population. Bifurcation analysis was carried out to examine the combined influences of the environmental carrying capacity, the interspecific competition strength and the impact of amplification or dilution agent on the deer mouse population. Our results showed that the system with amplification agent required a higher carrying capacity or stronger interspecific strength to compensate for its amplification effect in suppressing the SNV prevalence; this situation explains the lack of reduction in SNV prevalence despite the presence of high biodiversity in some empirical studies. In this study, we highlight the importance of investigating the roles of the additional species in an assemblage to better understand their relationship with the SNV prevalence in deer mouse population.
A new method to construct the distinct Hamiltonian circuits in complete
graphs is called Half Butterfly Method. The Half Butterfly Method used the concept
of isomorphism in developing the distinct Hamiltonian circuits. Thus some theoretical
works are presented throughout developing this method.
Medical diagnosis is the extrapolation of the future course and outcome of a disease and a sign of the likelihood of recovery from that disease. Diagnosis is important because it is used to guide the type and intensity of the medication to be administered to patients. A hybrid intelligent system that combines the fuzzy logic qualitative approach and Adaptive Neural Networks (ANNs) with the capabilities of getting a better performance is required. In this paper, a method for modeling the survival of diabetes patient by utilizing the application of the Adaptive NeuroFuzzy Inference System (ANFIS) is introduced with the aim of turning data into knowledge that can be understood by people. The ANFIS approach implements the hybrid learning algorithm that combines the gradient descent algorithm and a recursive least square error algorithm to update the antecedent and consequent parameters. The combination of fuzzy inference that will represent knowledge in an interpretable manner and the learning ability of neural network that can adjust the membership functions of the parameters and linguistic rules from data will be considered. The proposed framework can be applied to estimate the risk and survival curve between different diagnostic factors and survival time with the explanation capabilities.
The three dimensional free convection boundary layer flow near a stagnation point region is embedded in viscous nanofluid with the effect of g-jitter is studied in this paper. Copper (Cu) and aluminium oxide (Al2O3) types of water base nanofluid are cho- sen with the constant Prandtl number, Pr=6.2. Based on Tiwari-Das nanofluid model, the boundary layer equation used is converted into a non-dimensional form by adopting non- dimensional variables and is solved numerically by engaging an implicit finite-difference scheme known as Keller-box method. Behaviors of fluid flow such as skin friction and Nusset number are studied by the controlled parameters including oscillation frequency, amplitude of gravity modulation and nanoparticles volume fraction. The reduced skin friction and Nusset number are presented graphically and discussed for different values of principal curvatures ratio at the nodal point. The numerical results shows that, in- crement occurs in the values of Nusset number with the presence of solid nanoparticles together with the values of the skin friction. It is worth mentioning that for the plane stagnation point there is an absence of reduced skin friction along the y-direction where as for axisymmetric stagnation point, the reduced skin friction for both directions are the same. As nanoparticles volume fraction increased, the skin friction increased as well as the Nusset number. The results, indicated that skin frictions of copper are found higher than aluminium oxide.
The constraint of two ordered extreme minima random variables when one
variable is consider to be stochastically smaller than the other one has been carried
out in this article. The quantile functions of the probability distribution have been
used to establish partial ordering between the two variables. Some extensions and
generalizations are given for the stochastic ordering using the important of sign of the
shape parameter.
The box plot has been used for a very long time since 70s in checking the existence
of outliers and the asymmetrical shape of data. The existing box plot is constructed
using five values of statistics calculated from either the discrete or continous data. Many
improvement of box plots have deviated from the elegant and simplier approach of exploratory
data analysis by incorporating many other statistic values resulting the turning
back of the noble philosophy behind the creation of box plot. The modification using
range value with the minimum and maximum values are being incorporated to suit the
need of selected discrete distribution when outliers is not an important criteria anymore.
The new modification of box plot is not based on the asymmetrical shape of distribution
but more on the spreading and partitioning data into range measure. The new propose
name for the box plot with only three values of statistics is called range-box plot.