EFFECT OF SELECTED INSECTICIDE ON WHITEFLY (Bemisia tabaci) INFESTING BRINJAL PLANTS 298 THE CAUSAL IMPACT OF GOVERNMENT SUPPORT ON THE SMALL STRATEGIC CROP INDUSTRY: MALAYSIA’S EXPERIENCE Kamarudin Othman, Hapiza Omar, Hanani Ahmad Fuad , Jamilah Laidin, Irlisuhayu Mohd Ramli Faculty of Business, Universiti Teknologi MARA, Kedah, Malaysia  hanan5560@uitm.edu.my (Corresponding author) Corresponding author ARTICLE HISTORY: Received: 24-Jan-2020 Accepted: 27-Apr-2020 Online Available: 26-May- 2020 Keywords: ARDL, Fertiliser subsidy, Minimum price guarantee, Paddy production, Rice ABSTRACT Rice is a staple food for countries particularly in Asia, including Malaysia, hence arise concern on its supply. Thus, this study aims to investigate the causal relationship between the incentives (fertiliser subsidy and minimum price guarantee) given and paddy production in Kedah, the rice bowl of Malaysia. A time-series data ranging from the 1996 main season to the main season in 2018 (1996H1 to 2018H1) employed. Using ARDL and Granger approaches, the result shows that any disequilibrium in paddy production is regulated to accomplish the convergence faster. In the short-run, fertiliser subsidy and paddy prices are not significant in the MADA regions. The results of the analysis show that there are various causes of the long-run in the production of rice with subsidized steel and rice prices related to rice yield. The outcome of this study demonstrates the presence of long-run causality between paddy production with fertiliser subsidy and paddy price with respect to paddy yield. Therefore, to boost the paddy yield, the government should increase the number of subsidies, and control the paddy price to prevent it from falling below the market price. Therefore, future research should investigate the causal relationship between fertiliser and paddy price subsidy to the different types of paddy seeds planted. Contribution/ Originality This study significantly contributes to the growing evidence of the relationship between government incentive schemes to paddy production in Malaysia. The research explored how fertiliser subsidy and minimum price guarantee contributed to paddy yield. The evidence provided in this study demonstrates that the presence of long-run causality between paddy production of fertiliser subsidy minimum price guarantee. DOI: 10.18488/journal.1005/2020.10.1/1005.1.298.310 ISSN (P): 2304-1455/ISSN (E):2224-4433 How to cite: Kamarudin Othman, Hapiza Omar, Hanani Ahmad Fuad, Jamilah Laidin, and Irlisuhayu Mohd Ramli (2020). The causal impact of government support on the small strategic crop industry: malaysia’s experience. Asian Journal of Agriculture and Rural Development, 10(1), 298-310. © 2020 Asian Economic and Social Society. All rights reserved. Asian Journal of Agriculture and Rural Development Volume 10, Issue 1 (2020): 298-310 http://www.aessweb.com/journals/5005 mailto:hanan5560@uitm.edu.my http://www.aessweb.com/journals/5005 http://crossmark.crossref.org/dialog/?doi=10.18488/journal.1005/2020.10.1/1005.1.298.310 Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 299 1. INTRODUCTION Paddy that is undoubtedly the principal crop in Malaysia represents the predominant source of the single Malaysians’ staple food that is rice. Paddy planting areas cover almost 6% to 10% of the total agricultural areas in Malaysia (Economic Planning Unit, 2001, and 2006, Fatimah et al., 2020). Although rice is a staple food for Malaysians, by 2015, more than 70% of the agricultural land was occupied by oil palm (Fatimah et al., 2020). Paddy cultivation is initiated on a small scale and the average field size is 1 to 2 hectares. There are two types of paddy crops, namely wetland paddy, and dryland paddy. The average production for the main and off seasons is around 3000 kilograms per hectare. Based on these facts, it seems that the Malaysian paddy yield is below the national paddy production target at a rate of 5% in 2010 (Food and Agriculture Organization, 2002; Othman, 2008; Department of Agriculture, Malaysia, 2011). The actual yield for rice production is between 44% to 68% and under the potential yield (Othman, 2008, Che et al., 2019). This situation causes a shortage of paddy production around 125,552.15 metric tons per annum. The shortage of local paddy production also causes lower local rice supply. On average, rice consumption in Malaysia is roughly 120 kilograms per person every year. This is approximately 26% to 44% of the daily calorie needs (Suthida et al., 2008: USDA, 2012, Che et al., 2019). To solve these problems, the government has to import almost 30% of the rice from neighbouring countries such as Vietnam and Thailand. Therefore, to increase the rice self-sufficiency level (SSL), Malaysia has to increase its local paddy production. The tenth Malaysia plan aims to achieve a minimum self-sufficiency level of 70 percent food security (Saiful and Mohd, 2016). Thus, as to increase paddy production, farmers have to increase their level of productivity. To increase Malaysian paddy production, the government has adopted various series of subsidy policies. The purpose of subsidy policies is basically to reduce the costs of production and to boost farmers’ income. There are also subsidies given to reduce consumers’ burden. Fundamentally, there are three instruments implemented under the subsidy policy, which are the fertiliser scheme, paddy price subsidy scheme, and rice price subsidy scheme. The fertiliser subsidy was introduced in the early 1950s. In the early stage, the government stopped the fertiliser subsidy when the price of fertilisers decreased. Nonetheless, in 1979, the government had to resume this policy when the price of fertilisers increased and until to date, fertilisers have become 100% subsidy. It was only given to the landowners and not to the operators. The average allocation of fertilisers for the years (2003- 2007) was approximately USD 44.4 million. Meanwhile, the government introduced the paddy price subsidy scheme in 1980. The purpose of introducing this scheme was to increase the farmers’ productivity and income (Tan, 1987; Tawang et al., 1999, Fatimah, 2018). The scheme was stemmed from the realisation that the paddy farmers have recorded the highest poverty rates in Malaysia. On average, the paddy price subsidy has increased over the years except in 2006. However, the allocation gap between 2005 and 2006 had a minor difference of about 0.6%. This proves that the paddy price subsidy is the major instrument in boosting the farmers' income. Therefore, the purpose of this study is to analyse the fertiliser and paddy price subsidy schemes on paddy cultivation based on the estimated paddy production function. 2. LITERATURE REVIEW Among the most important inputs in agricultural production is fertiliser. The correct measurement of the fertiliser application may help to increase production output. Generally, fertilisers are classified into two, namely organic fertilisers and chemical fertilisers. Fertilisers are additional nutrients supplied to the plant. They are very important, especially in increasing the level of productivity among paddy farmers. The right application of the fertilisers and the right fertiliser specification help increase the nutrients in the soil, such as nitrogen, phosphorus, potassium, and others. The effective use of fertilisers may increase the output by more than 20 %. Therefore, fertilisers are crucial in Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 300 increasing the fertility of the soil especially in unfertile paddy growing areas (Nurul et al., 2012; Meena et al., 2013). Many of the researchers concluded that fertiliser inputs are complementary to the production of other inputs (Lori et al., 2013; Lenis et al., 2013; Adnan et al., 2017). In other words, to increase the paddy production, farmers have to make sure that they use effective technology and the correct schedule of fertilisers. In a preliminary study of paddy, Jayaraman (1983) found that fertiliser is an important input in the paddy production in Indonesia. Meanwhile, in 2010, Krishna et al. (2010) discovered that fertiliser is an extremely important factor for the family farms in Nepal. Chaudhry et al. (2009) and Oluwatosin (2013) also discovered that fertiliser is significant in agricultural production. In the efficiency study conducted by Sherlund et al. (2002), it was discovered that labour can be substituted by the use of chemical fertilisers. Furthermore, Suresh and Ready (2006) found that the ratio of MPV / MFC land in the district of Kerala, India is high, which is approximately 3:04. This means that land is significant as a determinant factor in paddy production. However, the size of the cultivated areas is small. This causes a low average output. Therefore, to increase paddy production, Ready (2007) suggested that the farmers in Kerala employ fertilisers. Furthermore, the results of the study indicated that one rupee increase in the fertiliser expenditure increases the yield by about 2.83 rupees in the Kerala district of India. The literature above shows strong evidence in supporting that fertilisers are very important in increasing the paddy yield. Similar to other sectors, the agricultural output price also affects the allocation of resources to aggregate production (Liyan and Richard, 1993). Mundlak (1988), in his article, stressed the importance of price as an input to determine the choice of technique used. Furthermore, the price also important in influencing the level of productivity (Mundlak, 1988). When the price of agricultural output is high, farmers are given the incentive to improve their earning and, indirectly, it can help to raise the supply of the agricultural product. Numerous research has shown that price and productivity are positively related. Several studies have shown that price and productivity are positively related to the research in South Asia (Bhavani et al., 2019) According to Timmer et al. (1983), two major roles played by price are, firstly, it shows the level of consumption especially among the poor. Secondly, the price also affects the supply through increased production as an incentive to producers. In reality, there are two different desires among consumers and agricultural producers. Consumers prefer lower prices of goods while farmers want a higher price of output. In the case of Malaysia, this difference is offset by the government through price subsidies for both consumers and agricultural producers. 3. MATERIALS AND METHODS 3.1. Co-integration test: the ARDL approach The underlying regression for this section is based on the log-linear production model as follows: 𝑌 = (𝑁, 𝐾, 𝐿40, 𝐿41, 𝐹, 𝑃) …………………... (1) where Y is the paddy yield. Meanwhile, F and P represent fertiliser and paddy price subsidy schemes respectively. In the meantime, N, K, LL40, and LL41 are selected control variables for land, capital, young farmers, and old farmers. A natural log of both sides of equation 1 is taken and equation 2 will be obtained as follows: 𝑙𝑛𝑌𝑡 = 𝛼0 + 𝛼1𝑙𝑛(𝑁𝑡) + 𝛼2𝑙𝑛(𝐾𝑡) + 𝛼3𝑙𝑛(𝐿40) + 𝛼4𝑙𝑛(𝐿41) + 𝛼5𝑙𝑛(𝐹𝑡) + 𝛼6𝑙𝑛(𝑃𝑡) + 𝜇𝑡 ... (2) Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 301 where α represents the elasticity of paddy production with respect to each independent variable in VAR model. Subscript 𝑡 represents time and µ represents the residual. Equation 2 is the underlying VAR that is used in the empirical analysis of all paddy production regions in MADA. The notation of 𝛼0 is the intercept, 𝜇𝑡 is white-noise random error, and 𝛼1 to 𝛼6 are the coefficients for the long-run model. The F-statistic of the bounds testing method (Pesaran and Pesaran, 1997; Pesaran et al., 2001) is examined for the joint significance of the coefficients of the lagged levels of the variables. Every year has two main seasons. So all together 32 observations. This study employed the Narayan (2004) approach because critical values produced by Narayan is more suitable for small sample sizes. This study employs time-series data from the main season in 1996 (1996H1) to the main season in 2011 (2011H1). This study only focuses on the Muda Agricultural Development Authority (MADA) granary area. The long-run equations which form the production functions model to be analysed in the present study can be portrayed as a general vector autoregressive model of order p in zt as follows: 1 0 1 , 1,2,3,..., p t i t i t i z a t z t T           …………………... (3) where 0a represents a (k+1) vector containing intercepts and β is a (k+1) vector of trend coefficients. Following this, a vector equilibrium error-correction model (VECM) is derived as shown below: 1 0 1 1 p t t i t i t i z a t z z              t = 1,2,3,…, T …………………... (4) where the (k+1) x (k+1) matrices 1 1 p k i i I       and 1 p i j j i      with i = 1,2, …, p-1 have both the long-run multipliers and short-run dynamic coefficients of the VECM. tz contains the vector of variables ty and tx . yt is the dependent variable and xt reflects the ‘forcing’ I (0) and I (1) variables. Assuming there is a unique long-run association between the parameters of interest, the VECM above can be written as: 1 1 0 1 1 1 , 1 0 1,2,..., p p t y yy t xx t i t i i t yt i i y a t y x y x t T                          ….. (5) In equation (5) above, i refers to the paddy sub-sector,  is the long-run multiplier, a is constant, and  is the white-noise error. Adapting the above specification to the variables in the present study, yield equations (6) and (7) below are for the production function model of paddy sub-sector: ∆𝑙𝑛𝑌𝑡 = 𝑐0 + 𝛿1𝑙𝑛𝑌𝑡−1 + 𝛿2𝑙𝑛𝑁𝑡−1 + 𝛿3𝑙𝑛𝐾𝑡−1 + 𝛿4𝑙𝑛𝐿40𝑡−1 + 𝛿5𝑙𝑛𝐿41𝑡−1 + 𝛿6𝑙𝑛𝐹𝑡−1 + 𝛿7𝑙𝑛𝑃𝑡−1 + ∑ 𝜔𝑗 𝑞1 𝑗=0 ∆𝑙𝑛𝑁𝑡−𝑗 + ∑ 𝛾𝑘 𝑞2 𝑘=0 ∆𝑙𝑛𝐾𝑡−𝑘 + ∑ ∅𝑙 𝑞3 𝑙=0 ∆𝑙𝑛𝐿40𝑡−𝑙 + ∑ 𝜗𝑚 𝑞4 𝑚=0 ∆𝑙𝑛𝐿41𝑡−𝑚 + ∑ 𝜇𝑛 𝑞5 𝑛=0 ∆𝑙𝑛𝐹𝑡−𝑛 + ∑ 𝜏𝑜 𝑞6 𝑜=0 ∆𝑙𝑛𝑃𝑡−𝑜 + 휀𝑡 …………………... (6) Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 302 ∆𝑙𝑛𝑁𝑡 = 𝑐0 + 𝛿1𝑙𝑛𝑁𝑡−1 + 𝛿2𝑙𝑛𝑌𝑡−1 + 𝛿3𝑙𝑛𝐾𝑡−1 + 𝛿4𝑙𝑛𝐿40𝑡−1 + 𝛿5𝑙𝑛𝐿41𝑡−1 + 𝛿6𝑙𝑛𝐹𝑡−1 + 𝛿7𝑙𝑛𝑃𝑡−1 + ∑ 𝜔𝑗 𝑞1 𝑗=0 ∆𝑙𝑛𝑁𝑡−𝑗 + ∑ 𝛾𝑘 𝑞2 𝑘=0 ∆𝑙𝑛𝐾𝑡−𝑘 + ∑ ∅𝑙 𝑞3 𝑙=0 ∆𝑙𝑛𝐿40𝑡−𝑙 + ∑ 𝜗𝑚 𝑞4 𝑚=0 ∆𝑙𝑛𝐿41𝑡−𝑚 + ∑ 𝜇𝑛 𝑞5 𝑛=0 ∆𝑙𝑛𝐹𝑡−𝑛 + ∑ 𝜏𝑜 𝑞6 𝑜=0 ∆𝑙𝑛𝑃𝑡−𝑜 + 휀𝑡 …………………... (7) The OLS technique is used to conduct the F-test for the above equation which tests the null hypothesis 0 1 2 3 4 5: 0H          against the alternate hypothesis 1 1 2 3 4 5: 0H          . The rejection of the null hypothesis shows that there exists a long-run relationship between the variables. Therefore, the variables in the model are co-integrated. The F-statistic obtained is a statistic in a generalised Dickey-Fuller type of regression that employed to test the significance of lagged level variables in a conditionally unrestricted equilibrium-correction model (Pesaran et al., 2001). F-statistic is then compared with the two sets of critical values that create a band. The first value presumes all the variables are I (1) while another value presumes all the variables are I (0). Critical values for I (1) series are known as upper bound while the critical values for the I (0) series are known as lower bound. If the computed F-statistic is beyond the upper-bound level, the null hypothesis can be rejected. This shows that there is a co-integration between the parameters or variables. However, if the F-statistic is below the lower-bound value, the null hypothesis cannot be rejected and this demonstrates no co- integration between variables. If the statistic falls between the upper and lower bounds, a conclusive inference cannot be made without knowing the order of integration of the independent variables (Narayan, 2004). If the variables are integrated I (1), then the upper bound is taken as the critical value. If the variables are integrated I (0), the lower bound becomes the critical value. In the meantime, there is also an alternative way to capture the existence of co-integration for the inclusive case. As mentioned by Kremers et al. (1992), the alternative way to capture the co-integration is by using the error-correction term (ECT). Additionally, the ECT value must be negative and significant. This condition indicates that there is a presence of adjustment towards the long-run equilibrium, therefore, the series are co-integrated. The bounds table by Pesaran and Pesaran (1997) and Pesaran et al. (2001) has produced a sample size of 500 and 1,000 observations and 20,000 and 40,000 replications consecutively. However, the present study just employs 31 observations. Therefore, the F-statistic obtained compared by the critical values generated by Narayan (2004), which is more suitable for small sample sizes1. After establishing the presence of co-integration, the long-run relationship between the variables is estimated via the following conditional ARDL (p, q1, q2, q3, q4, q5), a long- run model for paddy sub- sector production function in equations (8) and (9) below: ∆𝑙𝑛𝑌𝑡 = 𝑐0 + ∑ 𝛿1 𝑝 𝑖=1 ∆𝑙𝑛𝑁𝑡−𝑖 + ∑ 𝛿2 𝑞1 𝑖=0 ∆𝑙𝑛𝐾𝑡−𝑖 + ∑ 𝛿3 𝑞2 𝑖=0 ∆𝑙𝑛𝐿40𝑡−𝑖 + ∑ 𝛿4 𝑞3 𝑖=0 ∆𝑙𝑛𝐿41𝑡−𝑖 + ∑ 𝛿5 𝑞4 𝑖=0 ∆𝑙𝑛𝐹𝑡−𝑖 + ∑ 𝛿6 𝑞5 𝑖=0 ∆𝑙𝑛𝑃𝑡−𝑖 + 휀𝑡 …………………... (8) ∆𝑙𝑛𝑁𝑖,𝑡 = 𝑐𝑖,0 + ∑ 𝛿1 𝑝 𝑗=1 ∆𝑙𝑛𝑁𝑖,𝑡−𝑗 + ∑ 𝛿2 𝑞1 𝑙=0 ∆𝑙𝑛𝐾𝑖,𝑡−𝑙 + ∑ 𝛿3 𝑞2 𝑚=0 ∆𝑙𝑛𝐿40𝑖,𝑡−𝑚 + ∑ 𝛿4 𝑞3 𝑛=0 ∆𝑙𝑛𝐿41𝑖,𝑡−𝑛 + ∑ 𝛿5 𝑞4 𝑜=0 ∆𝑙𝑛𝐹𝑖,𝑡−𝑜 + ∑ 𝛿6 𝑞5 𝑝=0 ∆𝑙𝑛𝑃𝑖,𝑡−𝑝 + 휀𝑡 …………………... (9) The optimal lag for each variable (p, q1, q2, q3, q4, q5) is selected based on Schwarz Bayesian Criterion (SBC) or Akaike Information Criterion (AIC). After obtaining the long-run coefficient estimates, the short-run dynamic coefficients are estimated by OLS through the error-correction models associated with the long-run estimates. The error- 1 Narayan (2004) generated a different set of critical values using the same GAUSS code as Pesaran et al. (2001) but tailored his study with 31 observations. These new critical values were calculated using stochastic simulations for T=31 and 40000 replications for the F-statistic. Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 303 correction models corresponding to paddy sub-sector production functions are presented in equations (10) and (11) below: ∆𝑙𝑛𝑌𝑖,𝑡 = 𝛽𝑖 + ∑ ф1 𝑝 𝑗=1 ∆𝑙𝑛𝑁𝑖,𝑡−𝑗 + ∑ 𝜛2 𝑞1 𝑙=0 ∆𝑙𝑛𝐾𝑖,𝑡−𝑙 + ∑ 𝜑3 𝑞2 𝑚=0 ∆𝑙𝑛𝐿40𝑖,𝑡−𝑚 + ∑ 𝛾4 𝑞3 𝑛=0 ∆𝑙𝑛𝐿41𝑖,𝑡−𝑛 + ∑ 𝜂5 𝑞4 𝑜=0 ∆𝑙𝑛𝐹𝑖,𝑡−𝑜 + ∑ 𝜂6 𝑞5 𝑝=0 ∆𝑙𝑛𝑃𝑖,𝑡−𝑝 + 𝜗𝑒𝑐𝑡𝑖,𝑡−1 + 휀𝑖,𝑡 ..…… (10) ∆𝑙𝑛𝑁𝑖,𝑡 = 𝛽𝑖 + ∑ ф1 𝑝 𝑗=1 ∆𝑙𝑛𝑁𝑖,𝑡−𝑗 + ∑ 𝜛2 𝑞1 𝑙=0 ∆𝑙𝑛𝐾𝑖,𝑡−𝑙 + ∑ 𝜑3 𝑞2 𝑚=0 ∆𝑙𝑛𝐿40𝑖,𝑡−𝑚 + ∑ 𝛾4 𝑞3 𝑛=0 ∆𝑙𝑛𝐿41𝑖,𝑡−𝑛 + ∑ 𝜂5 𝑞4 𝑜=0 ∆𝑙𝑛𝐹𝑖,𝑡−𝑜 + ∑ 𝜂6 𝑞5 𝑝=0 ∆𝑙𝑛𝑃𝑖,𝑡−𝑝 + 𝜗𝑒𝑐𝑡𝑖,𝑡−1 + 휀𝑖,𝑡 ……… (11) In equations (10) and (11), the short-run coefficient estimates for the paddy sub-sector production function models are given by , , , , ,      each model converges to its long-run equilibrium. reflects how fast the equilibrium is achieved or how the speed of adjustment is. 3.2. Granger causality If the co-integration test reveals that the variables are co-integrated, we employ the Vector Error- Correction (VEC) Model estimation as in equations 12 to 18y. The Error-Correction Model (ECM) is used to test for Granger-type causality. A benefit of using an ECM specification to test for causality is that it permits the testing of the short-run causality over the lagged different regressors. It also used to test the long-run causality over the lagged Error-Correction Term (ECT). A significant ECT confirms the presence of the long-run causality from the independent variables to the measured variable. Nevertheless, if the variables are not co-integrated, then Vector Autoregressive (VAR) Model is used in the first difference in the estimation given that both variables fertiliser and paddy price is I (1). Thus, ECM is used to examine the Granger causality between variables. ∆𝐿𝑌𝑡 = 𝛼0𝑌 + ∑ 𝛿1𝑌 𝑘 𝑖−1 ∆𝑌𝑡−1 + ∑ 𝛿2𝑌 𝑘 𝑖−1 ∆𝑁𝑡−1 + ∑ 𝛿3𝑌 𝑘 𝑖−1 ∆𝐾𝑡−1 + ∑ 𝛿4𝑌 𝑘 𝑖−1 ∆𝐿𝐿40𝑡−1 + ∑ 𝛿5𝑌 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿6𝑌 𝑘 𝑖−1 ∆𝐿𝐹𝑡−1 + ∑ 𝛿7𝑌 𝑘 𝑖−1 ∆𝐿𝑃𝑡−1 + 𝛾1𝐸𝐶𝑇𝑡−1 + 휀𝑡 ………….. (12) ∆𝐿𝑁𝑡 = 𝛼0𝑁 + ∑ 𝛿1𝑁 𝑘 𝑖−1 ∆𝑁𝑡−1 + ∑ 𝛿2𝑁 𝑘 𝑖−1 ∆𝑌𝑡−1 + ∑ 𝛿3𝑁 𝑘 𝑖−1 ∆𝐾𝑡−1 + ∑ 𝛿4𝑁 𝑘 𝑖−1 ∆𝐿𝐿40𝑡−1 + ∑ 𝛿5𝑁 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿6𝑁 𝑘 𝑖−1 ∆𝐿𝐹𝑡−1 + ∑ 𝛿7𝑁 𝑘 𝑖−1 ∆𝐿𝑃𝑡−1 + 𝛾2𝐸𝐶𝑇𝑡−1 + 휀𝑡 …………. (13) ∆𝐿𝐾𝑡 = 𝛼0𝐾 + ∑ 𝛿1𝐾 𝑘 𝑖−1 ∆𝐾𝑡−1 + ∑ 𝛿2𝐾 𝑘 𝑖−1 ∆𝐾𝑡−1 + ∑ 𝛿3𝐾 𝑘 𝑖−1 𝑌𝑁𝑡−1 + ∑ 𝛿4𝐾 𝑘 𝑖−1 ∆𝐿𝐿40𝑡−1 + ∑ 𝛿5𝐾 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿6𝐾 𝑘 𝑖−1 ∆𝐿𝐹𝑡−1 + ∑ 𝛿7𝐾 𝑘 𝑖−1 ∆𝐿𝑃𝑡−1 + 𝛾3𝐸𝐶𝑇𝑡−1 + 휀𝑡 …………. (14) ∆𝐿𝐿40𝑡 = 𝛼0𝐿𝐿40 + ∑ 𝛿1𝐿𝐿40 𝑘 𝑖−1 ∆𝐿𝐿40𝑡−1 + ∑ 𝛿2𝐿𝐿40 𝑘 𝑖−1 ∆𝑌𝑡−1 + ∑ 𝛿3𝐿𝐿40 𝑘 𝑖−1 ∆𝑁𝑡−1 + ∑ 𝛿4𝐿𝐿40 𝑘 𝑖−1 ∆𝐿𝐾𝑡−1 + ∑ 𝛿5𝐿𝐿40 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿6𝐿𝐿40 𝑘 𝑖−1 ∆𝐿𝐹𝑡−1 + ∑ 𝛿7𝐿𝐿40 𝑘 𝑖−1 ∆𝐿𝑃𝑡−1 + 𝛾4𝐸𝐶𝑇𝑡−1 + 휀𝑡 ………….…………. (15) ∆𝐿𝐿41𝑡 = 𝛼0𝐿𝐿41 + ∑ 𝛿1𝐿𝐿41 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿2𝐿𝐿41 𝑘 𝑖−1 ∆𝑌𝑡−1 + ∑ 𝛿3𝐿𝐿41 𝑘 𝑖−1 ∆𝑁𝑡−1 + ∑ 𝛿4𝐿𝐿41 𝑘 𝑖−1 ∆𝐿𝐾𝑡−1 + ∑ 𝛿5𝐿𝐿41 𝑘 𝑖−1 ∆𝐿𝐿40𝑡−1 + ∑ 𝛿6𝐿𝐿41 𝑘 𝑖−1 ∆𝐿𝐹𝑡−1 + ∑ 𝛿7𝐿𝐿41 𝑘 𝑖−1 ∆𝐿𝑃𝑡−1 + 𝛾5𝐸𝐶𝑇𝑡−1 + 휀𝑡 …………………….. (16) ∆𝐿𝐹𝑡 = 𝛼0𝐿𝐹 + ∑ 𝛿1𝐿𝐹 𝑘 𝑖−1 ∆𝐿𝐹𝑡−1 + ∑ 𝛿2𝐿𝐹 𝑘 𝑖−1 ∆𝑌𝑡−1 + ∑ 𝛿3𝐿𝐹 𝑘 𝑖−1 𝑁𝑌𝑡−1 + ∑ 𝛿4𝐿𝐹 𝑘 𝑖−1 ∆𝐿𝐾𝑡−1 + ∑ 𝛿5𝐿𝐹 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿6𝐿𝐹 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿7𝐿𝐹 𝑘 𝑖−1 ∆𝐿𝑃𝑡−1 + 𝛾6𝐸𝐶𝑇𝑡−1 + 휀𝑡 …….….. (17) ∆𝐿𝑃𝑡 = 𝛼0𝐿𝑃 + ∑ 𝛿1𝐿𝑃 𝑘 𝑖−1 ∆𝐿𝑃𝑡−1 + ∑ 𝛿2𝐿𝑃 𝑘 𝑖−1 ∆𝑌𝑡−1 + ∑ 𝛿3𝐿𝑃 𝑘 𝑖−1 ∆𝑁𝑡−1 + ∑ 𝛿4𝐿𝑃 𝑘 𝑖−1 ∆𝐿𝐾𝑡−1 + ∑ 𝛿5𝐿𝑃 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿6𝐿𝑃 𝑘 𝑖−1 ∆𝐿𝐿41𝑡−1 + ∑ 𝛿7𝐿𝑃 𝑘 𝑖−1 ∆𝐿𝐹𝑡−1 + 𝛾7𝐸𝐶𝑇𝑡−1 + 휀𝑡 ……… (18) Where 𝐸𝐶𝑇𝑡−1is the lagged error-correction term. The Granger causality testing technique includes testing the significance of 𝛿𝑖𝑗 s limited on the optimum lags. The ECM equations 12 to 18 are the alternative tests of causality. The 𝛾𝑖𝑗 is zero, Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 304 which means that the variation in the dependent variable doesn't respond to the deviation in the long- run equilibrium in a period of 𝑡 − 1. For example, equation 1, if 𝛾𝑖𝑗 is zero and all 𝛿𝑖𝑗 is zero, that can be implied that all the variables do not Granger-cause paddy production. The insignificance of t and F-statistic in the Wald test suggests that the variable is weakly exogenous. If there is a presence of the co-integration between variables, we presume that at least one or all of the ECTs should be significantly non-zero. Granger causality between LY and other variables can be revealed by testing the following null hypothesis: a) Short-run Granger causality: (𝐻𝑂: 𝛿2𝑌 = 0, 𝐻𝑂: 𝛿2𝑁 = 0, 𝐻𝑂: 𝛿2𝐾 = 0, 𝐻𝑂: 𝛿2𝐿𝐿40 = 0, 𝐻𝑂: 𝛿2𝐿𝐿41 = 0, 𝐻𝑂: 𝛿2𝐿𝐹 = 0 𝑎𝑛𝑑 𝐻𝑂: 𝛿2𝐿𝑃 = 0, ) b) Long-run Granger causality: (𝐻𝑂: 𝛾1 = 0, 𝐻𝑂: 𝛾2 = 0, 𝐻𝑂: 𝛾3 = 0, 𝐻𝑂: 𝛾4 = 0, 𝐻𝑂: 𝛾5 = 0, 𝐻𝑂: 𝛾6 = 0 𝑎𝑛𝑑 𝐻𝑂: 𝛾7 = 0) c) Joint Granger causality: (𝐻𝑂 = 𝛿2𝑌 = 𝛾1 = 0, 𝐻𝑂=𝛿2𝑁 = 𝛾2 = 0, 𝐻𝑂=𝛿2𝐾 = 𝛾3 = 0, 𝐻𝑂=𝛿2𝐿𝐿40 = 𝛾4 = 0, 𝐻𝑂=𝛿2𝐿𝐿41 = 𝛾5 = 0, 𝐻𝑂=𝛿2𝐿𝐹 = 𝛾6 = 0 𝑎𝑛𝑑 𝐻𝑂=𝛿2𝐿𝑃 = 𝛾7 = 0) 4. EMPIRICAL RESULTS Based on the results of the unit-root test, ARDL co-integration test, vector autoregressive (VAR) Granger-cause test in the framework of vector error-correction (VECM), response functions (IRFs), variance decomposition (VDS), and diagnostic tests shown in table 1 and 2. Table 1 displays the results of Augmented Dickey-Fuller (ADF) and Phillip Perron test (PP) to the paddy yield, land, capital, young farmers, old farmers, fertiliser, and paddy price to MADA paddy production areas. Based on these tests, at level, the absolute t-statistic obtained for all the time-series data fails to reject the null hypothesis of which the data are not stationary time series. This situation describes that, all factors not stationary at 1% and 5% significance level. This research proceeded with the first- difference stationarity test and found out that the t-calculate was larger than the t-critical value at 1% and 5 %. This means that the integration of the variables is at the first difference or I(1). Table 1: Unit-root analysis for MADA Regions Variables ADF Test PP Test Conclusion Constant Without Trend Constant Without Trend LY ΔLY -1.696 -14.489* -2.041 -22.41** I(1) LN ΔLN -1.037 -7.489* -1.455 -7.537* I(1) LK ΔLK 0.254 -4.465* 0.080 -4.465** I(1) LL40 ΔLL40 -2.273 -6.125* -2.407 -6.108* I(1) LL41 ΔLL41 -2.083 -6.402* -2.198 -6.356* I(1) LF ΔLF -1.564 -6.779* -1.565 -6.611* I(1) LP ΔLP -0.115 -5.599** -0.051 -5.599** I(1) Notes: *, and** donated 1% and 5% significant levels respectively To approve the existence of co-integration among the parameters, this study employed the bound test and the Autoregressive Distributed Lag (ARDL) method. The optimal lag length needs to be determined first before the ARDL co-integration test. Based on the AIC suggestion, the optimal lag Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 305 length is 2. Table 2 shows that F-statistic for the MADA area is 3.4337. This F-statistic value is between the lower (2.303) and upper (3.483) bound value at 5% and 10% significant levels. This indicates that the presence of co-integration between variables in the model could not be ascertained. Therefore, the test on the error correction term (ECT) by Kremers et al. (1992) is used to confirm the co-integration between variables. To test the ECT, the long-run equation for the MADA area needs to be generated first. Table 2: Results of bounds test and critical values for case II Case II (intercept and no trend) Level of Significance Lower Upper 1% Significance 5% Significance 10% Significance 3.901 2.780 2.303 5.611 4.084 3.483 Computed F-statistic 3.4337 Notes: The F-statistics are compared with the upper bound i.e. I(1) and lower bound i.e. I(0) critical values for zero restriction on the coefficient of the lagged level variables N = 30 Based on Table 3, the coefficient estimation shows that land, old farmers, fertilisers, and paddy prices are significant in the long-run equation for MADA. The coefficient estimation of land, fertilisers, and paddy price show that paddy yield is elastic towards change in land, fertilisers, and paddy price. 1% increase in land size leads to an increase in output by 28%. Meanwhile, a 1% increase in fertilisers used triggers the increase in output by approximately 172%. This reflects that the use of fertilizers intensively and effectively increases paddy yield. On the other hand, old farmers and paddy prices are inelastic to the change in paddy yield. If there is a 1% increase in the number of old farmers, the output decreases by about 0.4%. Nonetheless, when the paddy price increases by 1%, the output decreases by 5%. Other variables are not significant to determine paddy yield in the long-run. Table 3: Long-run coefficient estimations for the determinants of the paddy yield in MADA Area Variables LN LK L40 L41 LF LP C Coefficient 28.29** -0.06 0.12 -0.04** 172.35** -5.07** -1438.3** t-statistic 2.86 -0.78 1.50 -2.17 2.83 -2.84 -3.31 Notes: *, ** indicates significance at 5% and 10% levels From the information in Table 3, we then computed the error-correction term (ECT). The negative and significant value of the error-correction term (ECT) shows the presence of the long-run relationship among variables. The Output Error-Correction Terms (ECT) for MADA are shown as below: - ΔLY = 0.015 + 46.966ΔLN** + 0.449ΔLK** + 256.334ΔLL40* - 345.454ΔLL41** [1.445] [4.083] [ 2.986)] [6.522)] [-4.233] + 93.228ΔLF - 2.851ΔLP - 1.038ect* [1.297] [-1.288] [-8.355] ……………………. (19) R2: 0.843; Adjusted R-squared: 0.793; AR(2): 0.970(0.396); Heteroscedasticity (1): 0.891(0.559); RESET(2): 1.645(0.217); Normality(2): 1.269(0.530); CUSUM: Stable; CUSUM-SQ: Stable Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 306 Notes: Brackets denote the t-statistics while p-values are in parentheses. Breusch-Godfrey, White, Ramsey RESET, and Jarque Bera tests are used to test for the presence of serial correlation, heteroscedasticity, model miss-specification, and residual non-normality in the model. * and ** denote significant at 1% and 5% level respectively. SR is short-run From the above ECT estimation, we found that the value of ECT is negative and significant. This shows that there exists a long-run relationship among variables. The lagged coefficient of the error- correction term (ECT) depicts the speed of adjustment of paddy yield to shock in the system. The speed of adjustment of paddy production is higher. From the short-run output, the speed of adjustment for MADA is 104%. The result suggested that any disequilibrium in paddy production is adjusted to achieve faster convergence. Additionally, it could be inferred from these results that paddy yield, land, capital, young and old farmers, fertilisers, and local paddy price are the long-run ‘forcing variables’ of paddy yield in MADA paddy production areas. From this result, we could make an initial assumption that at least there is a presence of one-way causality between variables in MADA paddy production areas. The present study also found that in the short-run, the land is significant in the MADA areas. The degree of elasticity of paddy yield with respect to land is that it is elastic and has a positive sign. This means that 1% of the increase in paddy planting areas leads to approximately 46.97 % of the increase in paddy yield. This indicates that in the short-run, the increase in harvested areas may increase the paddy yield (Fatimah et al., 2020). Besides, the present research also found that young and old farmers are also significant in influencing paddy production. The degree of elasticity of paddy yield concerning young and old farmers is elastic. However, both variables have a different sign. Young farmers have a positive effect on paddy yield; a 1% increase in young farmers, paddy yield increases roughly by 611%. The increase in old farmers reduces paddy production by about 611%. This indicates that young farmers are among the most vulnerable target to increases agricultural output (Balezentis et al., 2020). If paddy cultivation is too dependent on old farmers, the productivity may reduce. Nonetheless, this research also found that capital is significant in the short-run. From the above result, the degree of elasticity is inelastic and has a positive sign. Based on the empirical finding, a 1% increase in capital results in the increase of paddy yield by about 0.45 %. This finding is smaller as compared to the study done by Oniah et al. (2008). The diagnostic test was then employed to test the goodness of fit of the model. The goodness of fit of the model is adequate as roughly 84% of the differences in paddy yield clarified by dissimilarities in the regressors. Additionally, there is no autocorrelation, no residual non-normality, and no miss- specification problems associated with the MADA areas model. The present study also found that there is no heteroscedasticity in MADA. The CUSUM and CUSUM-Sq test in Figure 1 shows that the model parameters and residual variance are stable throughout the study. Figure 1: CUSUM and CUSUM-SQ Graphs for MADA paddy production 4.1. Causal relationship The result from Table 4 shows that land and capital are significant at 5% and 1% levels respectively in the paddy yield equation but there are no other explanatory variables that are significant in the Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 307 short-run. Meanwhile, in the land and capital equation, paddy yield is significant at 5% and 1% levels. These clearly show that paddy yield and land, and paddy yield and capital both have a bidirectional Granger causality in the short-run. Furthermore, we also discovered that young farmers and capital, old farmers and capital, old farmers and young farmers, and paddy price and fertilisers also have the presence of bidirectional Granger causality. Table 4: Short-run and long-run granger causality VEC Model Variable s MADA Short-run Causality Long- run Causalit y ∆ 𝑳𝒀𝒕 ∆ 𝑳𝑵𝒕 ∆ 𝑳𝑲𝒕 ∆ 𝑳𝑳𝟒𝟎𝒕 ∆ 𝑳𝑳𝟒𝟏𝒕 ∆ 𝑳𝑭𝒕 ∆ 𝑳𝑷𝒕 ECT Wald F-statistics ∆ 𝐿𝑌𝑡 - 5.96** [0.02] 11.53* [0.00] 0.89 [0.36] 0.20 [0.66] 1.38 [0.25] 1.39 [0.25] 1.03* (-7.85) ∆ 𝐿𝑁𝑡 5.81** [0.03] - 0.50 [0.49] 0.93 [0.35] 1.37 [0.26] 0.39 [0.54] 0.35 [0.56] 0.009** (2.17) ∆ 𝐿𝐾𝑡 24.99* [0.00] 0.57 [0.46] - 11.30* [0.00] 6.00** [0.02] 1.04 [0.32] 1.26 [0.27] 0.86* (4.10) ∆ 𝐿𝐿40𝑡 2.39 [0.14] 0.89 [0.36] 10.71* [0.00] - 101.66* [0.00] 0.09 [0.77] 0.12 [0.73] -0.0006 (-1.58) ∆ 𝐿𝐿41𝑡 -.71 [0.41] 1.29 [0.27] 6.05** [0.02] 100.87* [0.00] - 2.68 [0.12] 2.48 [0.13] 0.0004 (0.82) ∆ 𝐿𝐹𝑡 2.13 [0.16] 0.38 [0.55] 0.82 [0.37] 0.09 [0.76] 2.63 [0.12] - 5213* [0.00] 0.001 (1.53) ∆ 𝐿𝑃𝑡 2.22 [0.15] 0.84 [0.37] 1.23 [0.28] 0.11 [0.74] 2.43 [0.13] 50.80* [0.00] - -0.044 (-1.43) Notes: ( ) is t-statistic, [ ] is p value. *, **, and *** shows 1%, 5%, and 10% level of significance respectively Based on t-statistic, it could be seen that the coefficient of ECT is significant in the paddy yield, land, and capital equations. Meanwhile, ECT is found to be insignificant in fertilisers, young and old farmers, and paddy price equations. This implies that paddy yield, land, and capital carry out the short-run adjustment to the long-run equilibrium. Based on the paddy yield equation, the deviation of paddy yield from the long-run equilibrium, and all the seven variables interact ineffective way to restore the long-run equilibrium. Joint Granger causality in Table 5 suggests that the presence of bidirectional causality in capital and land, young farmers and land, young farmers and capital, old farmers and capital, old farmers and young farmers, fertilisers and old farmers, fertilisers and young farmers, and paddy price and fertilisers have the presence of two-way short-run Granger causality. There are also a few variables having a one-way Granger causality such as young farmers Granger-cause paddy yield, old farmers Granger-cause paddy yield, fertilisers Granger-cause paddy yield, and paddy price Granger-cause paddy yield. Table 5: Joint causality Granger causality VEC model Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 308 Variables MADA Joint Causality ECT ∆ 𝑳𝒀𝒕 ECT ∆ 𝑳𝑵𝒕 ECT ∆ 𝑳𝑲𝒕 ECT ∆ 𝑳𝑳𝟒𝟎𝒕 ECT ∆ 𝑳𝑳𝟒𝟏𝒕 ECT ∆ 𝑳𝑭𝒕 ECT ∆ 𝑳𝑷𝒕 Wald F-statistics ∆ 𝐿𝑌𝑡 - 32.22* [0.00] 30.83* [0.00] 31.95* [0.00] 31.74* [0.00] 30.82* [0.00] 30.83* [0.00] ∆ 𝐿𝑁𝑡 2.91*** [0.08] - 3.14*** [0.07] 4.98** [0.02] 4.30** [0.03] 2.36 [0.12] 2.36 [0.12] ∆ 𝐿𝐾𝑡 12.49* [0.00] 9.47* [0.00] - 16.42* [0.00] 14.09* [0.00] 8.44* [0.00] 8.49* [0.00] ∆ 𝐿𝐿40𝑡 1.29 [0.29] 3.02*** [0.07] 5.49** [0.01] - 63.23* [0.00] 1.27 [0.30] 1.26 [0.31] ∆ 𝐿𝐿41𝑡 0.36 [0.69] 1.64 [0.22] 3.58** [0.05] 57.34* [0.00] - 2.75*** [0.08] 2.61** [0.09] ∆ 𝐿𝐹𝑡 1.18 [0.33] 1.17 [0.33] 1.17 [0.33] 1.18 [0.33] 3.72** [0.04] - 2627* [0.00] ∆ 𝐿𝑃𝑡 1.12 [0.35] 1.08 [0.36] 1.04 [0.37] 1.03 [0.38] 3.33*** [0.06] 2561* [0.00] - Notes: ( ) is t-statistic, [ ] is p value. *, **, and *** shows 1%, 5%, and 10% level of significance respectively 5. CONCLUSION In the short-run, fertiliser subsidy and paddy prices are not significant in the MADA regions. This indicates that fertiliser subsidy and paddy prices do not give any significant effects to the local paddy production in the short run. Nonetheless, in the long-run, fertiliser subsidy and paddy prices have a positive relationship and are significant with respect to paddy yield. Although fertiliser subsidy and paddy prices are significant with paddy yield, the relationship between paddy price and paddy yield is negative. This means that when paddy price increases, paddy yield decreases. Poor soil fertility and the awful weather condition might be the cause of this negative relationship. The present research also discovered the presence of one-way short-run causality between fertilisers and paddy yield in the MADA regions. This clearly shows that any shock in fertilisers will increase paddy yield. Fertiliser is considered as one of the most important non-traditional inputs in paddy production. The evidence from the Granger causality analysis shows that fertilisers play an important role in causing the rise in paddy output in the MADA regions. Therefore, the government should provide information and knowledge related to the use of fertilisers to paddy farmers. The government should also reduce the price of fertilisers by increasing the number of subsidies. Also, paddy farmers themselves should have a better understanding of how to use the fertilisers efficiently. They should also follow the fertilising schedule provided by the government to raise the paddy yield. Therefore, researchers should expand the study further by investigating the causal relationship between fertiliser subsidy and paddy price subsidy to different types of paddy seeds planted. Funding: This study did not receive any specific financial support. Competing Interests: The authors declared that they have no conflict of interests. Contributors/Acknowledgement: The authors would like to express their gratitude to everyone who involved directly or indirectly in completing this study. The causal impact of government support on the small strategic crop industry: Malaysia’s experience. Views and opinions expressed in this study are the views and opinions of the authors, Asian Journal of Agriculture and Rural Development shall not be responsible or answerable for any loss, damage or liability, etc. caused in relation to/arising out of the use of the content. References Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 309 Adnan, N., Nordin, S. M., Rahman, I., & Noor, A. (2017). The impacts and visions of the green fertilizer technologies (GFT). World Journal of Science, Technology and Sustainable Development, 14(4), 336-354. doi.org/10.1108/wjstsd-08-2016-0053. Balezentis, T., Ribasauskiene, E., Morkunas, M., Volkov, A., Streimikiene, D., & Toma, P. (2020). Land use policy young farmers’ support under the common agricultural policy and sustainability of rural regions: evidence from Lithuania. Land Use Policy, 94(February), 104542. https://doi.org/10.1016/j.landusepol.2020.104542. Bhavani, S., Nigel, P., & Frances, A. B. (2019). Agricultural inputs and nutrition in South Asia. Food Policy, Elsevier, 82, 28-38. doi.org/10.1016/j.foodpol.2018.10.011. Chaudhry, I. M., Khan, M. B., & Anwar, M. (2009). Factors affecting cotton production in Pakistan: empirical evidence from Multan district. MPRA, 5(2), 91-100. Che, O. S., Shaharudin, A., & Tumin, S. A. (2019). The status of the paddy and rice industry in Malaysia. in Khazanah research institute. Malaysia http://www.krinstitute.org/assets/contentMS/img/template/editor/20190409_RiceReport_Ful lReport_Final.pdf. Department of Agriculture, Malaysia (2011). Pembangunan industri padi. Retrived on January 2, 2010 from http://www.doa.gov.my/c/document.../get_file? Economic Planning Unit (2001). Eight Malaysia plan (2001-2005). Kuala Lumpur, Prime Minister’s Department. Economic Planning Unit (2006). Ninth Malaysia plan (2006-2010). Kuala Lumpur, Prime Minister’s Department. Fatimah, K. (2018). Evaluation of agricultural subsidies and the welfare of farmers: Malaysia agricultural subsidies report 2018. in IDEAS Policy Research Berhad, Kuala Lumpur. Fatimah, M. A., Kusairi, M. N., Bach, N. L., Illisriyani, I., Abdulla. I., Gregory, H. W. S., Siti, A. T., & Ahmad, A. A. S. (2020). Implications of the dominant shift to industrial crops in Malaysian agriculture. phase i : system dynamics model of the paddy and rice sector, In Khazanah Research Institute. http://www.krinstitute.org/assets/contentMS/img/template/editor/0.Full%20report%2020041 0.pdf. Food and Agriculture Organization (2002). Selected indicators of food and agriculture development in the asia-pacific region 1992-2002. Food and Agriculture Organization of The United Nations, Regional Office for Asia and the Pacific. Bangkok. Retrieved on January 2, 2010 from ftp://ftp.fao.org/docrep/fao/004/AD452e/ad452e00.pdf. Jayaraman, T. K. (1983). Determinants of paddy production in Indonesia: 1972-1981, A Simultaneous Equation Model Approach. Asian Development Bank, Economic Office, Manila Kremers, J., Ericsson, N., & Dolado, J. (1992). The power of cointegration tests. Oxford Bulletin of Economics and Statistics, 54, 349-367. Krishna, M., Yasuhiro, A., & Hideyuki (2010). Estimation of production function and resource use condition of organic coffee cultivation in different farm size and altitude categories in the hill region of Nepal. European Journal of Scientific Research, 45(3), 438-449. Lenis, S., Liverpool, T., & Sheu, S. (2013). Spillover effects of targeted subsidies an assessment of fertilizer and improved seed use in Nigeria. International Food Policy Research Institute (IFPRI) Discussion Paper 01260, Retrieved on Mac 4, 2013 from http://www.ifpri.org/sites/default/files/publications/ifpridp01260.pdf. Liyan, E. F., & Richard, K. P. (1993). Prices and productivity in agriculture. The Review of Economics and Statistics, 75(3), 471-482. Lori, B., Dean, K., Bram, T., & Christopher, U. (2013). Probability of fertilizer: experimental evidence from female rice farmers in mali. Economics Department Working Paper No. 111, Department of Economics, Yale University, Retrieved on Mac 4, from http://www.econ.yale.edu/growth_pdf/cdp1020.pdf. http://www.krinstitute.org/assets/contentMS/img/template/editor/20190409_RiceReport_FullReport_Final.pdf http://www.krinstitute.org/assets/contentMS/img/template/editor/20190409_RiceReport_FullReport_Final.pdf http://www.doa.gov.my/c/document.../get_file http://www.krinstitute.org/assets/contentMS/img/template/editor/0.Full%20report%20200410.pdf http://www.krinstitute.org/assets/contentMS/img/template/editor/0.Full%20report%20200410.pdf ftp://ftp.fao.org/docrep/fao/004/AD452e/ad452e00.pdf http://www.ifpri.org/sites/default/files/publications/ifpridp01260.pdf http://www.econ.yale.edu/growth_pdf/cdp1020.pdf Asian Journal of Agriculture and Rural Development, 10(1)2020: 298-310 310 Meena, D. S., Meena, V. R., & Meena, A. K. (2013). Fertilizer management studies on growth and productivity of hybrid Indian mustard Brassica juncea (L). Journal of Oilseed Brassica, 4(1), 39-42. Mundlak, Y. (1988). Endogenous technology and measurement of productivity. In S.M. Capalbo and J.M. Antle (Ed.), Agricultural productivity: measurement and explanation, Washington D.C. Narayan, P. K. (2004). Reformulating critical values for the bounds F-statistics approach to cointegration: An application to the tourism demand model for Fiji. Department of Economics Discussion Papers No. 02/04, Monash University, Australia. Nurul, N. R., Mad, Nasir, S., Zainalabidin, M., & Alias, R. (2012). The impact of fertilizer subsidy on Malaysia paddy/rice industry using a system dynamics approach. International Journal of Social Science and Humanity, 2(3), 213-219. Oluwatosin, O. F. (2013). Determinants of perceived effectiveness of organic fertilizer used by farmers in Oyo State, Nigeria. Agricultura Tropica Et Subtropica, 46/1, 23-28. Oniah, M. O., Kuye, O. O., & Idiong, I. C. (2008). Efficiency of resource use in small scale seamp rice production in Obubra local government area of cross river state, Nigeria. Middle-East Journal of Scientific Research, 3(3), 145-148. Othman, O. (2008). Rice production and potential for hybrid rice in Malaysia. In International Plantation Industry Conference and Expo - IPiCEX. Shah Alam. Malaysia. Pesaran, M. H., & Pesaran, B. (1997). Working with Microfit 4.0: interactive econometric analysis. Oxford University Press, Oxford. Pesaran, M. H., Shin, Y., & Smith, R. J. (2001). Bounds testing approaches to the analysis of level relationships. Journal of Applied Econometrics, 16, 289-326. Ready, G. P. (2007). Impact of water management on production of rice in Balipatna command area or Orissa, India. Sabaragamnwa University Journal, 1(2), 15-20. Saiful, N. S. B., & Mohd, R. M. (2016), Distributional effects of non-farm incomes in a Malaysian rice bowl. International Journal of Social Economics, 43(2), 205-220. Sherlund, S. M., Christopher, B. B., & Akinwumi A. A. (2002). Smallholder technical controlling for environmental production conditions. Journal of Development Economics, 69, 85-101. Suresh, A., & Ready, T. R. K. (2006). Resource-use efficiency of paddy cultivation in Peechi command area of Thrissur district of Kerala: an economic analysis. Agricultural Economics Research Review, 19, 159-171. Suthida, W., Gil, R., & Jammie, P. (2008). Changing food consumption and import in Malaysia, opportunities for Australian agricultural export, Australia Bureau of Agricultural and Resources Economics, Research Report 08.6 August 2008. Tan, S. H. (1987). Malaysia’s rice policy: a critical analysis. Institute of Strategic and International Studies (ISIS), Malaysia 1987. Tawang, A., Ariffin, T. M., & Ahmad, T. (1999). Effects of trade liberalization on agriculture in Malaysia: commodity aspects. Working Paper No. 46, CGPRT Centre, Bogor, September 1999. Timmer, C., Peter, W., Falcon, P., & Scott, R. P. (1983). Food policy analysis. Baltimore, Md., The Johns Hopkins University Press. USDA Foreign Agricultural Service (2012). Global agricultural information network report. Retrieved from http://static.globaltrade.net/files/pdf/20110406115103794.pdf. http://static.globaltrade.net/files/pdf/20110406115103794.pdf