Selected Papers from the International Conference on Information, Communication, and Engineering 2013
View this Special IssueResearch Article  Open Access
Analyzing Taiwan IC Assembly Industry by GreyMarkov Forecasting Model
Abstract
This study utilizes the black swan theorem to discuss how to face the lack of historical data and outliers. They may cause huge influences which make it impossible for people to predict the economy from their knowledge or experiences. Meanwhile, they cause the general dilemma of which prediction tool to be used which is also considered in this study. For the reason above, this study uses 2009 Q1 to 2010 Q4 quarterly revenue trend of Taiwan’s semiconductor packaging and testing industry under the global financial turmoil as basis and the grey prediction method to deal with nonlinear problems and small data. Under the lack of information and economic drastic changes, this study applies Markov model to predict the industry revenues of GM(1,1) and DGM(1,1) results. The results show that the accuracy of 2010 Q1–Q3 is 88.37%, 90.27%, sand 91.13%, respectively. Besides, they are better than the results of GM(1,1) and DGM(1,1) which are 86.51%, 77.35%, 75.46% and 73.77%, 74.25%, 59.72%. The results show that the prediction ability of the grey prediction with Markov model is better than traditional GM(1,1) and DGM(1,1) sfacing the changes of financial crisis. The results also prove that the greyMarkov chain prediction can be the perfect criterion for decisionmakers judgment even when the environment has undergone drastic changes which bring the impact of unpredictable conditions.
1. Introduction
Financial turmoil has huge influences on varieties of industries, like the black swan theorem which produces a great deviation and causes adverse reactions [1]. In Taiwan, semiconductor production capacity accounts for more than half of the global semiconductor manufacturing market, and the IC packaging and testing accounts for the global foundry capacity of over 60%. The speed and environment of new products replacing older market create the competitiveness of related industries in Taiwan. The black swan theorem will produce outliers. This study found that the way to remove the deviation would be more difficult to predict regression or other methods of forecasting future trends accurately.
The grey method proposed by professor Long can provide very good prediction under shorttime historical data [2]. Among the method, GM(1,1) provides suitable reference tools to decisionmakers in the changing and competitive environment. In the previous researches, grey theory applied to the economic, financial, engineering, and other fields reached ideal results [3–5]. For example, Ken and Li [6] used a different grey prediction model which included traditional GM(1,1), RGM(1,1), and Tan’s GM(1,1), and he used chained convenience stores sales patterns as backgrounds for model selection. Hung et al. [7] also used GM(1,1) to predict under incomplete information. Compared with two prediction models of Egli and Walfischbertoni [8], it had proved that the GM(1,1) rolling model produced the lowest error. Lin and Yang [9] used GM(1,1) to predict Taiwan semiconductor optoelectronics industry, and the average residual is less than 10%, which proves the high accuracy of grey prediction to semiconductor industry. But the financial crisis causes excessive fluctuations in random sequence and it easily leads to the condition of inadequate accuracy. If it is applied to the Taiwan semiconductor packaging and testing industry revenue data, it will cause the lack of precision because of deficiencies historical data message with random fluctuations, which also increases the difficulty in using general prediction tools in this paper.
Markov (Markovchain) is generally used in the forecasting model on a random time series, and it is a dynamic system which is based on the transition between the states and all the random factors to reflect the impact [7]. Although it is very suitable for handling random variation problems, it still needs enough raw data to determine the probability of state transitions. Hence, this study uses Markov theory to improve GM(1,1) and reduce the variation when historical data increase and financial crisis happens. By evaluating accuracy, the prediction effect of Markov model can be improved [3].
2. Grey System Theory
Yin [10] has investigated during the time span of 1996 to 2011, the earliest year when the database was available to the best possible full text availability to date. The paper points out that there were 482 publications that met the selection criteria of “grey relational analysis” or “gray relational analysis” topic search. It is seen that the grey system theory is one of the most popular forecasting methods. Ken and Li [6] applied grey system as a module to help predicting sales revenue of convenience chain stores in Taiwan. They do get more accuracy on their prediction. Also, Huang et al. [11] predict the silkworm cocoon yield in china base on the greyMarkov forecasting model. The findings show that greyMarkov has more accuracy than other forecasting model. Meanwhile, Hsu and Chen [12] compared various statistical methods based on the suitability of shortterm forecasting. They got similar findings.
3. Methodology
3.1. GM(1,1) and DGM(1,1) Grey Forecasting Model
GM(1,1) is a type of prediction model, the analysis steps are listed below.
Step 1. The original sequence is defined as
Step 2. Use 1AGO to decrease the randomness of sequence and define the 1AGO as where .
Step 3. Establish a firstorder differential equation
Step 4. From (3), we can get where
Step 5. Use 1AGO to get the real prediction data
Step 6. From Step 5, a trend curve equation can be formed where .
DGM(1,1) discrete grey forecasting model is the dynamic cycle rolling model; it is from original sequence by GM(1,1). Use this result of DGM(1,1) model to become the first stage of forecasting value that is needed. Comparative forecasting value and actual value are examined in order to get error value and accuracy.
3.2. Partition of States by MarkovChain Forecasting Model
The values of are distributed in the region of the trend curve that may be divided into a convenient number of contiguous intervals. When falls in interval , one of such intervals, it may be regarded as corresponding to a state , and can be signified as follows: where and is the amount of states where is time function, is in the state, and are tine varying, and and are fixed value, depending on the original data.
3.3. Calculation of State Transfer Probability
In Markovchain series, the changing probability form to can be presented as follows: where is the transition probability of state transferred from state for steps (in this paper, 1 step stands for 1 quarter), is the number of transition steps each time, is the number of original data of state transferred from state for m steps, and is the number of original data points in state .
The state transition probability matrix is as follows:
The state transition probability reflects the transition rules of a system. The state transition probability matrix describes the probability of transition from state to . Generally, it is necessary to observe the onestep transition matrix . Supposing the object to be forecasted is in state , row in matrix should be considered. If max , then what will most probably happen in the system at the next moment is the transition from state to state . It is difficult to determine the future transition of the state; if two or more transition probabilities in the row of matrix are the same, the transition probability matrix of twostep transition matrix or multistep transition matrix , where , should be considered [4].
3.4. Calculate the Forecasting Data
After the determination of the future state transition of a system, that is, the determination of grey elements , the changing interval of the forecast value is between and . The most probable forecast value, , is considered to be the middle value of the determined state interval; that is,
3.5. Accuracy Inspection Analysis of Forecasting Ability
Numerous methods exist for judging forecasting model accuracy, and no single recognized inspection method exists for forecasting ability. MAPE is often used to measure forecasting accuracy. Smaller MAPE value indicates better forecasting ability
Evaluation of MAPE forecasting ability is divided as forecasting ability is evaluated as follows: (i)10 excellent forecasting ability, (ii)10~20 good forecasting ability, (iii)20~50 reasonable forecasting ability, (iv)50 poor forecasting.
4. Results and Analyses
There are many factors which could influence the Taiwan economy, such as the country’s culture, the government’s policy, and the industry framework. Some factors are clear, and others are not clear. But semiconductor still is the major supporting in Taiwan.
IC Industry Situation in Taiwan. Taiwan Semiconductor Manufacturing Company (TSMC) and United Microelectronics Corporation (UMC) are the first and second in worldwide ranking, with a combined market share exceeding 50%. IC assembly industry total product value is the first of the worldwide [13].
Table 1 shows the historical data series of the IC assembly industry Revenue of Taiwan from 2007 Q1 to 2010 Q3, but fluctuating randomly. So this paper forecasts and analyzes IC assembly industry revenue by GreyMarkov forecasting model.
 
Unit: million NTDs. 
4.1. Build the GM(1,1) Grey Forecasting Model
Using the data of the IC assembly industry revenue of Taiwan from 2007 Q1 to 2010 Q3, in Table 1 and (1)–(7) for model building, the GM(1,1) model of IC assembly industry revenue yield is where is the series number of the quarter and means 2007’Q1.
4.2. Partition of States by MarkovChain Forecasting Model
According to the actual data, four states, that is, four contiguous intervals, are established about the curve of . According to (8)(9), the four states intervals can be got as follows: where denotes the average value of the historical IC assembly industry revenue of Taiwan from 2007 Q1 to 2009 Q4. Figure 1 shows the historical data series and the states intervals.
4.3. Calculation of State Transition
The number of the historical data in every interval can be gotten as follows: where denotes the number of the historical data in the interval and . The numbers of onestep transiting to , , , and from are as follows: Next , where and , can be calculated in the same way. Then calculate the onestep transition probability to every states interval and present them in the transition matrix as follows:
4.4. Calculate the Forecast Value
According to (12), the date of IC assembly industry revenue in 2010 Q1 is calculated; that is,
As Table 1 shows, the actual IC assembly industry revenue in Taiwan in 2010 Q1 is 48204412 million NTD. So the precision of the forecasting is 88.37 percent.
By the same method, the IC assembly industry revenue in 2010 Q2 and Q3 can be calculated, the forecast value of 2010 Q2 is 58,636,460 million NTD, and that of 2010 Q3 is 59105003 million NTD. So the forecasting precisions are 90.27% and 91.13%.
4.5. Comparison of Forecast Value between the GreyMarkov, GM(1,1), and DGM(1,1) Forecasting Model
First build a GM(1,1) Grey forecasting model to forecast the IC assembly industry revenue from 2010 Q1 to Q3. And utilize GM(1,1) to build and construct a 5 DGM(1,1) rolling modeling, forecasting in the same way. Then the forecasting results are presented and the results that are forecasted by GreyMarkov model are presented in Table 2.

The forecast value of GreyMarkov forecasting models is more precise than GM(1,1) grey forecasting model in Figure 2. It is showed that the GreyMarkov forecasting model is better for forecasting the IC assembly industry revenue in Taiwan.Figure 2 is showed that even if it causes the random array to fluctuate too big to face the influence that Financial turmoil, the forecasting curve of GM(1,1), DGM(1,1), and Grey Markov, the GreyMarkov forecasting models are still fit for the predication of the IC assembly industry revenue in Taiwan.
5. Conclusions and Discussions
Grey Markov prediction model is a hybrid model which mixes traditional grey prediction and Markovchain, with fewer limitations which not need to complicate matrix operations. By the empirical analysis, we can prove the ability with the advantages of shortterm forecasts, and by using GM(1,1), most future trends can be reflected. Indeed, the Markov probability matrix can extract the information provided by historical data to response to fluctuation system. Even when facing inconsistent fluctuating events, such as the financial crisis, it still gain a high accuracy [12].
In addition to revenue, the impact of IC packaging and testing industry in Taiwan also includes the country economy, industrial structure, policies of government, and even the climate, culture, and other relevant issues. These indicators are followed by specific historical data trend of random fluctuations which are suitable for grey Markov model. It is best to be applied to big fluctuation timeseries problems, but outliers often appear in the scope of expectations. According to past experience, it is often impossible for anyone to believe the possibility of their emergence. With the tremendous impact, we easily fall into the Black Swan myth because of human nature, and let it becomes hard explained and unpredictable.
When enterprise decisionmakers make decisions for a major future investment, there must be a risk management mechanism; however, the reference data toward decision making always depend on historical data. For future industry prediction, there are many methods which can be used, but most of them can only reach accurate prediction in normal behavior. In the past market which included drastic noncontrollable factor’s changes, instant dropping data could be recognized as special cases or outliers which would not be calculated. But the results may be distorted due to limited or insufficient data. Therefore, this study applied Markov prediction, which includes fewer complicated calculation and higher accuracy, to help the important references and make better investment portfolio decision.
Acknowledgments
The authors gratefully acknowledge Dr. Chung, Yi Shun from National Cheng Kung University (NCKU) and Mr. Chuang, Yun Chung from Metal Industries Research & Development Centre (MIRDC), for their comments on earlier versions of the paper.
References
 N. Nicholas, “The black swan: the impact of the highly improbable,” Journal of the Management Training Institut, vol. 36, no. 3, p. 56, 2008. View at: Google Scholar
 D. JuLong, “Control problems of grey systems,” Systems and Control Letters, vol. 1, no. 5, pp. 288–294, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. Y. B. Yidan, “GreyMarkov forecasting model and its application,” Systems EngineeringTheory & Practice, vol. 4, p. 12, 1992. View at: Google Scholar
 Z. Jun, S. Jianming, and F. Yingjun, “The grey model MGM (1, n) and its application,” Systems EngineeringTheory & Practice, vol. 5, pp. 109–113, 1997. View at: Google Scholar
 Y. He and M. Huang, “A greyMarkov forecasting model for the electric power requirement in China,” in Proceedings of the Advances in Artificial Intelligence: 4th Mexican International Conference on Artificial Intelligence (MICAI '05), pp. 574–582, Springer, 2005. View at: Google Scholar
 M.L. Ken and Y.B. Li, “A study of applying grey predictive module into sales revenue prediction of convenience chain stores an investigation of 7Eleven and familymart,” Journal of Grey System, vol. 12, no. 3, pp. 139–146, 2009. View at: Google Scholar
 K.C. Hung, K.P. Lin, F.Y. Hsu, C.K. Wang, and J.C. Lin, “Grey model with rolling mechanism for radiowave pathloss forecasting in suburban environment,” Journal of Grey System, vol. 13, no. 2, pp. 49–53, 2010. View at: Google Scholar
 G. L. Siqueira, E. A. Vasquez, R. A. Gomes, C. B. Sampaio, and D. A. Borghi, “Comparison between propagation measurements and coverage prediction models for small urban cells,” in Proceedings of the Global Telecommunications Conference—Communications: The Key to Global Prosperity (GLOBECOM '96), vol. 2, pp. 1182–1186, IEEE, London, UK, November 1996. View at: Publisher Site  Google Scholar
 C.T. Lin and S.Y. Yang, “Forecast of the output value of Taiwan's optoelectronics industry using the Grey forecasting model,” Technological Forecasting and Social Change, vol. 70, no. 2, pp. 177–186, 2003. View at: Publisher Site  Google Scholar
 M.S. Yin, “Fifteen years of grey system theory research: a historical review and bibliometric analysis,” Expert Systems with Applications, vol. 40, no. 7, pp. 2767–2775, 2013. View at: Publisher Site  Google Scholar
 L. Huang, P. Jin, Y. He et al., “Prediction of silkworm cocoon yield in china based on greyMarkov forecasting model,” in Proceedings of the Advances in Artificial Intelligence: 5th Mexican International Conference on Artificial Intelligence (MICAI '06), pp. 505–512, Springer, 2006. View at: Google Scholar
 C.C. Hsu and C.Y. Chen, “Applications of improved grey prediction model for power demand forecasting,” Energy Conversion and Management, vol. 44, no. 14, pp. 2241–2249, 2003. View at: Publisher Site  Google Scholar
 C.C. Lou, T.P. Lee, S.C. Gong, and S.L. Lin, “Effects of technical innovation on market value of the U.S. semiconductor industry,” Technological Forecasting and Social Change, vol. 77, no. 8, pp. 1322–1338, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 LeiChuan Lin and ShanYau Wu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.