The Application of Mixed Smoothing Spline and Fourier Series Model in Nonparametric Regression

In daily life, mixed data patterns are often found, namely, those that change at a certain sub-interval or that follow a repeating pattern in a certain trend. To handle this kind of data, a mixed estimator of a Smoothing Spline and a Fourier Series has been developed. This paper describes a simulation study of the estimator in nonparametric regression and its implementation in the case of poor households. The minimum Generalized Cross Validation (GCV) was used in order to select the best model. The simulation study used generation data with a Uniform distribution and a random error with a symmetrical Normal distribution. The result of the simulation study shows that the larger the sample size n, the better the mixed estimator as a model of nonparametric regression for all variances. The smaller the variance, the better the model for all combinations of samples n. Very poor households are characterized predominantly in their consumption of carbohydrates compared to that of fat and protein. The results of this study suggest that the distribution of assistance to poor households is not the same, because in certain groups there are poor households that consume higher carbohydrates, and some households may consume higher fats.


Introduction
The Smoothing Spline has been used in nonparametric regression analysis [1]. The Smoothing Spline can provide better results because the analysis can overcome the patterns of data that show a sharp increase and decrease, resulting in a relatively smooth curve [2]. The advantages of using the Smoothing Spline are its unique statistical properties, it enables visual interpretation, it can handle smooth data and functions, and can readily handle data that change at certain sub-intervals [3][4][5][6]. In addition to the Spline estimator, the Fourier Series estimator is an estimation technique that is quite popular in nonparametric regression. The Fourier Series is a model that provides good statistical and visual interpretation among nonparametric regression models. The advantage of estimating the Fourier Series is that it can handle the characteristics of data that follow a repeating pattern at a certain trend interval, and provides a good statistical interpretation. It can handle data characters that follow a repeating pattern at certain trend intervals [7].
In addition to the development of research on nonparametric regression, mixed estimators in nonparametric regression have been recently developed. A mixed estimator of the Fourier Series and the Truncated Spline was developed by Sudiarsa et al. [8]. The local Polynomial Mixture and Truncated Spline were developed by Suparti et al. [9]. In daily life, mixed data patterns often appear; in particular, data patterns can fluctuate at certain sub-intervals or follow repeating patterns in a certain trend. Thus, a mixed estimator of the Smoothing Spline and the Fourier Series is applied to handle this kind of data. In the current study, to test the performance of the mixed of the Smoothing Spline and the Fourier 2 of 22 Series estimator as a model of nonparametric regression, a simulation study was carried out. In this paper, we discuss the suitability of the sample size and variance in this model of nonparametric regression. After the simulation study was completed, it was applied to the real data of poor households in Bali Province.
The province of Bali consists of nine regencies/cities. Among the nine districts/cities in Bali Province, Karangasem Regency has the highest poverty rate [10]. This is seen in the high percentage of poor people, of 6.28%. One of the most frequently used indicators to determine the level of poverty in a community is the level of household spending. Therefore, a comprehensive approach is needed to alleviate poverty through the identification of the factors that influence the expenditure of poor households in Karangasem Regency. The poor household data used in this study were household expenditure, wage/salary, and carbohydrate, fat, and protein consumption per capita per month. According to the World Bank, poverty can be reflected in household expenditure and income [11]. Household expenditure is influenced by income and consumption factors [12][13][14]. In addition to the income factor, the consumption factor is also very influential on the expenditure of poor households [15]. Based on the factors that influence the consumption of fat per capita per month, it is expected that it can be analyzed using the Fourier Series. This is because the higher the wage/salary, the higher the expenditure; however, at a certain interval, as the wage/salary decreases, the expenditure also decreases. In contrast, the other factors are examined using the Smoothing Spline. Thus, the problem of poor households in Karangasem Regency was examined with a mixed Smoothing Spline and Fourier Series estimator in nonparametric regression.

Mixed Smoothing Spline and Fourier Series Estimator in Nonparametric Multivariable
The mixed Smoothing Spline and Fourier Series estimator in multivariable nonparametric regression is an estimator that combines two methods that have different characteristics, i.e., data that change at sub-intervals are characteristic of the Smoothing Spline and data that have repeating patterns are characteristic of the Fourier Series. The following describes the mixed Smoothing Spline and Fourier Series estimator. Paired data is provided as x 1i , x 2i , . . . , x pi , t 1i , t 2i , . . . , t qi , y i with y i , i = 1, 2, . . . , n, which indicates many observations are response variables. The relationship between the response variables y i and predictor variables x 1i , x 2i , . . . , x pi , t 1i , t 2i , . . . , t qi follows a multivariable nonparametric regression model: y i = µ x 1i , x 2i , . . . , x pi , t 1i , t 2i , . . . , t qi + ε i , i = 1, 2, . . . , n.
where µ is the regression curve and ε i is random error, which is assumed to be independent and identical with a symmetrical Normal distribution having zero mean and variance σ 2 . Assume that the shape of the regression curve µ is unknown and additive in nature, so the regression model is obtained as follows: h k (t ki ) is a component of the Fourier Series. The component regression curve is g j x ji , j = 1, 2, . . . p and assumed to be smooth in the sense that it is contained in the Sobolev space, namely: Function g is a function of unknown shape and is assumed to be smooth in the sense that it fits within Sobolev space W. The Sobolev space W can be decomposed into a direct sum of two spaces W 0 and W 1 that are mutually perpendicular, namely: α uk cos ut ki , k = 1, 2, . . . , q.
The mixed estimator Smoothing Spline and Fourier Series in nonparametric regression can be obtained by optimization Penalized Least Squares (PLS) as follows [16]: Equation (1) becomes: The explanation of the transformation from Equation (1) to Equation (2) is described in Appendix A.
The estimator of the mixture is highly dependent on the smoothing parameter and the oscillation parameter. These parameters are selected using the Generalized Cross Validation method. The selection of the optimal smoothing and oscillation parameters in the mixed Smoothing Spline and Fourier Series estimator requires the following lemma:   Proof. The GCV function can be written as: Based on Lemma 2, we obtain: The smoothing parameter (λ) and the optimal oscillation parameter (K) are obtained from the minimization function GCV(λ , K ), i.e., GCV λ optimal , K optimal = Min

Simulation Study for Mixed Smoothing Spline and Fourier Series
A simulation study for the mixed Smoothing Spline and Fourier Series estimator in multivariable nonparametric regression was conducted. The simulation was carried out to describe the ability of the Smoothing Spline estimator and the Fourier Series in estimating the multivariable nonparametric regression curve. The steps taken are as follows.
Using the mixed Smoothing Spline and Fourier Series estimator in Equation (1), this simulation uses two conditions, such that the variance of the observation size is n = 25, 50, 100, and 200, and the error variance is σ 2 = 0.25, 0.5, and 1. The sample size n was used to show if increasing the sample size improved the model estimate. The sample sizes n = 25, 50, 100, and 200 are considered sufficient to represent small and large sample sizes. Similarly, the variance measure was used to determine if a smaller variance resulted in a better model estimate. The measures of the variance error, σ 2 = 0.25, 0.5, and 1, were deemed sufficient to represent small and large variance sizes. The combination for the simulation model can be seen in Table 1 below. All simulation experiments in Table 1 will then look for the GCV value. The regression equation designed for this simulation study is as follows: Furthermore, the function error ε i is generated from the distribution N(0, 1), x i ∼ U(0, 1), t i ∼ U(0, 1). The scatterplot simulation data are shown in Figure 1. 14) Compare the results in each scenario using the GCV, MSE, and R 2 criteria. 15) Form a conclusion about the ability of the mixed Smoothing Spline and Fourier Serie estimator to estimate the multivariable nonparametric regression curve.
Point 10 states the selection of the smoothing parameters ( λ ) and oscillation param eters (K) are optimal using the Generalized Cross Validation (GCV) method.
Using the mixed Smoothing Spline and Fourier Series estimator in Equation (1), thi simulation uses two conditions, such that the variance of the observation size is n = 25, 50 100, and 200, and the error variance is 2 σ = 0.25, 0.5, and 1. The sample size n was used to show if increasing the sample size improved the model estimate. The sample sizes n = 25 50, 100, and 200 are considered sufficient to represent small and large sample sizes. Simi larly, the variance measure was used to determine if a smaller variance resulted in a bette model estimate. The measures of the variance error, 2 σ = 0.25, 0.5, and 1, were deemed suf ficient to represent small and large variance sizes. The combination for the simulation mode can be seen in Table 1 below. All simulation experiments in Table 1 will then look for the GCV value. The regres sion equation designed for this simulation study is as follows: Furthermore, the function error is generated from the distribution . The scatterplot simulation data are shown in Figure 1.   Figure 1 shows the scatterplot with a sample size of n = 100 and an error variance σ 2 = 0.25. Figure 1 indicates that the pattern of the relationship between y and x tends to change in certain sub-intervals. The change in the pattern of the data can be seen at several points, namely, points 0.6, 0.7, and 0.8. At intervals of 0 to 0.6, the data pattern tends to decrease, between 0.6 and 0.7 the data pattern tends to increase, between 0.7 and 0.8 the data pattern increases, and at intervals of 0.8 to 1 the data pattern once again decreased. This data pattern is in accordance with the characteristics of the Smoothing Spline. The pattern of the relationship between y and t shows that the plot indicates an unclear pattern or leads to a change in the shape of the repeating data behavior pattern following an upward trend line. Based on the results of this description, this pattern is characteristic of the Fourier Series.
Based on Figure 1, it can be said that, in the use of the two predictors, there are two different data patterns. The pattern of the relationship between y and x tends to change at certain sub-intervals, whereas the relationship between y and t follows a repeating pattern on a trend. Thus, there is a different pattern between the predictors, so that the data was modeled using a mixed model of the Smoothing Spline and Fourier Series. The best model was selected based on the optimal value of the smoothing parameter and the oscillation parameter. The optimal values of the smoothing parameter and the oscillation parameter were selected based on the minimum GCV value. Table 2 presents the selection of the smoothing parameters and optimal oscillation parameters based on the minimum GCV value.  Table 2 shows that, for the various combinations of smoothing parameters and oscillation parameters used for the Smoothing Spline and Fourier Series mixed models, the optimal smoothing parameter value obtained was 0.042 and the oscillation parameter was 1. These results provide the minimum GCV value of 0.0000184. Based on this, the value of R 2 = 86.439% was obtained and the MSE value was 0.257. The obtained values of R 2 and MSE demonstrate the suitable performance of the Smoothing Spline and Fourier Series mixed model. Figure 2 presents a scatter plot between the value of y and the valueŷ. In Figure 2, the estimated data plot is very close to the original data, so this model can be used to make accurate predictions.
Furthermore, with the same procedure, one hundred repetitions were carried out for scenario 4. The minimum GCV, R 2 , and MSE values for one hundred repetitions are sum-  In Figure 2, the estimated data plot is very close to the original data, so this model can be used to make accurate predictions.
Furthermore, with the same procedure, one hundred repetitions were carried out for scenario 4. The minimum GCV, R 2 , and MSE values for one hundred repetitions are summarized in Table 3. In order to make it easier to compare the effect of the error variance, Table 4 presents a comparison of the R 2 value for one hundred replications based on n = 100. Determining whether there is a significant difference based on the different error variance values was tested based on the average difference between the groups using the hypothesis:

Hypothesis 1 (H1):
There is no significant difference between the average,

Hypothesis 2 (H2):
There is a real difference between the average of the average, The output presented in the Appendix C shows that the p-value 0.000 < 0.05. Therefore, it was concluded that H 0 should be rejected, meaning that there is a difference on average between R 2 σ 2 = 0.25 , R 2 σ 2 = 0.5 , and R 2 σ 2 = 1 .
Furthermore, the Tukey test was used to compare the average of all couples. The output of presented in Appendix B shows that all the error variance values are in different groups, which means that, on average, there is a difference between the values of R 2 σ 2 = 0.25 , R 2 σ 2 = 0.5 , and R 2 σ 2 = 1 . Based on the output above, it can also be seen that R 2 σ 2 = 0.25 gives better results than the other options, followed by R 2 σ 2 = 0.5 and R 2 σ 2 = 1 . Thus, for n = 100 it can be seen that, the smaller the value of σ 2 , the better the results that were obtained. The same steps were carried out for 12 trials, and the whole experiment is summarized in Table 5. In order to make it easier to compare the simulation results, a scatter diagram is presented in Figures 3 and 4 for the GCV and R 2 values of the Polynomial-Trigonometric experimental function for each variation of the number n.    Table 5 shows that the obtained results improve as the sample size increases. This is indicated by the resulting GCV value. It can be seen that the larger the sample size, the smaller the result of the GCV value. The GCV value of the Smoothing Spline and Fourier  Table 5 shows that the obtained results improve as the sample size increases. This is indicated by the resulting GCV value. It can be seen that the larger the sample size, the smaller the result of the GCV value. The GCV value of the Smoothing Spline and Fourier Series mixed model for the number of samples n = 200 is smaller than that for the sample size n = 100. Similarly, for the number of samples n = 100, the resulting GCV value was smaller than that for the sample size n = 50, etc. In addition to using GCV criteria, we also used MSE criteria and R 2 to determine the suitability of the model. The simulation results show that a sample size of n = 200 gives better results based on the highest value R 2 and lowest MSE for all given error variance values.   Table 5 shows that the obtained results improve as the sample size increases. This is indicated by the resulting GCV value. It can be seen that the larger the sample size, the smaller the result of the GCV value. The GCV value of the Smoothing Spline and Fourier Series mixed model for the number of samples n = 200 is smaller than that for the sample size n = 100. Similarly, for the number of samples n = 100, the resulting GCV value was smaller than that for the sample size n = 50, etc. In addition to using GCV criteria, we also used MSE criteria and R 2 to determine the suitability of the model. The simulation results show that a sample size of n = 200 gives better results based on the highest value R 2 and lowest MSE for all given error variance values. Table 5 also shows that the smaller the sample error variance, the better the model estimate obtained for all sample size variations. This is indicated by the smaller the size of the error variance, the smaller the GCV value. This condition is as shown in Table 5: at a sample size of n = 25, the error variance 2  is presented in Table 6.   Table 5 also shows that the smaller the sample error variance, the better the model estimate obtained for all sample size variations. This is indicated by the smaller the size of the error variance, the smaller the GCV value. This condition is as shown in Table 5: at a sample size of n = 25, the error variance σ 2 = 1 gives a GCV value of 0.01741598, then decreases at σ 2 = 0.5 to 0.00378041, and decreases again at σ 2 = 0.25 to 0.00185812. This can also be seen in the sample size n = 50, where the error variance is σ 2 = 1 gives a GCV value of 0.00046407, then decreases at σ 2 = 0.5 to 0.00019175, and decreases again at σ 2 = 0.25 to 0.00012196. Similar results were seen for n = 100 and n = 200; for the sample size n = 100, for variance values of σ 2 = 1, σ 2 = 0.5, and σ 2 = 0.25, the GCV results were 0.00006150, 0.00002668, and 0.00001073, respectively, whereas for the sample size n = 200, the GCV results were respectively 0.00000629, 0.00000243, and 0.0000040. Furthermore, to determine the power of the mixed Smoothing Spline and Fourier Series, a comparison was made between the nonparametric regression of the Smoothing Spline, the Fourier Series, and the mixed Smoothing Spline and Fourier Series. In this simulation, the observed size variance was n = 200 and the error variance was σ 2 = 0.25, and 100 repetitions were undertaken. The selection of the best model was determined by the minimum GCV value. A comparison table for n = 200 and the variance of error σ 2 = 0.25 is presented in Table 6. The results show that the mixed Smoothing Spline and Fourier Series in nonparametric regression performed better than the individual Smoothing Spline or Fourier Series. This can be seen based on the minimum GCV value in Table 6. After performing the simulation, real-life applications were carried out regarding the use of the mixed method of the Smoothing Spline and Fourier Series in multivariable nonparametric regression.

Application on Poor Household Data in Karangasem, Bali
Linearity tests and scatter plots can be used for modeling the pattern of relationships between response variables and predictor variables, and for examining the behavior of data patterns. This is explained in the following. A non-linearity test was undertaken based on the Terasvirta Neural Network Test to identify data patterns using a significance level of 5%. This identification was used to compare the data patterns between using a parametric regression or a nonparametric regression. The non-linearity test results are presented in Table 7. Based on Table 7, it can be concluded that the relationship between the variables of wages/salaries per capita per month, consumption of carbohydrates per capita per month, consumption of fat per capita per month, and consumption of protein per capita per month, with household expenditure variables per capita per month, was non-linear. Because the form of this non-linear function was not known, it was modeled using multivariable nonparametric regression. This was clarified using a scatter plot, as explained below.
The presentation of data in the form of a scatter plot aims to identify the behavioral pattern of the distribution of the data, and whether the data has a pattern that changes in sub-intervals or follows a repeating pattern. Figure 5a shows that the relationship between the response variable of household expenditure per capita per month and the predictor variable of wages/salaries does not follow a certain pattern. The data plot indicates that the relationship does not show a clear pattern or that there is a tendency to exhibit a change in behavior at a certain interval. Wages/salaries increased and household expenditure per capita per month also increased; this can be seen in the interval of 100,000 IDR until 250,000 IDR, in which there was an increase in household expenditure per capita per month. By comparison, for wages/salaries in the interval of 250,000 IDR until 292,500 IDR, monthly household expenditure per capita experienced a slight decrease from the previous interval. Based on the results of this description analysis, the pattern of the relationship between the variable wage/salary and household expenditure per capita a month followed the Smoothing Spline approach. Figure 5b shows that the pattern of the relationship between the response variable of household expenditure per capita per month and the predictor variable of carbohydrate consumption per capita per month changes at certain sub-intervals. Carbohydrate consumption at the interval of 1200-1300 g increased, and carbohydrate consumption at the interval of 1300-1400 g decreased slightly, and increased in the interval of 1400-1700 g. If the consumption of carbohydrates is high, the expenditure will also increase. Carbohydrate consumption in the interval of 1700-1800 g decreased. Based on this illustration, the relationship pattern between the variable consumption of carbohydrates per capita per month and household expenditure per capita per month was approached using the Smoothing Spline. The third variable that is thought to influence household expenditure per capita per month is the consumption of fat per capita per month.
interval of 1300-1400 g decreased slightly, and increased in the interval of 1400-1700 g. If the consumption of carbohydrates is high, the expenditure will also increase. Carbohydrate consumption in the interval of 1700-1800 g decreased. Based on this illustration, the relationship pattern between the variable consumption of carbohydrates per capita per month and household expenditure per capita per month was approached using the Smoothing Spline. The third variable that is thought to influence household expenditure per capita per month is the consumption of fat per capita per month.  Figure 5c shows that the pattern of the relationship between the response variable of household expenditure per capita per month and the predictor variable of fat consumption per capita per month does not follow a certain pattern or lead to changes in the form of repetitive data behavior patterns that follow an upward trend line. Based on the results of this description, the pattern of the relationship between the consumption of fat per capita per month with household expenditure per capita per month was approached using the Fourier Series. Figure 5d shows that the pattern of the relationship between the response variable of household expenditure per capita per month with the predictor variable of protein consumption per capita per month does not follow a certain pattern. The data plot indicates that it does not show a clear pattern or a tendency to change behavior at a certain interval. Consumption of protein per capita a month in the interval of 80-110 g increases slowly. Then, the consumption of protein per capita per month in the interval  Figure 5c shows that the pattern of the relationship between the response variable of household expenditure per capita per month and the predictor variable of fat consumption per capita per month does not follow a certain pattern or lead to changes in the form of repetitive data behavior patterns that follow an upward trend line. Based on the results of this description, the pattern of the relationship between the consumption of fat per capita per month with household expenditure per capita per month was approached using the Fourier Series. Figure 5d shows that the pattern of the relationship between the response variable of household expenditure per capita per month with the predictor variable of protein consumption per capita per month does not follow a certain pattern. The data plot indicates that it does not show a clear pattern or a tendency to change behavior at a certain interval. Consumption of protein per capita a month in the interval of 80-110 g increases slowly. Then, the consumption of protein per capita per month in the interval of 110-125 g experienced an increase with that in household expenditure per capita per month. Based on the results of this descriptive analysis, the pattern of the relationship between the response variable of household expenditure per capita per month with the predictor variable of protein consumption per capita per month used the Smoothing Spline approach.
Based on the data and information obtained from the data exploration, the nonlinearity test and scatter plot showed that the Smoothing Spline Mixed Estimator and Fourier Series are the appropriate methods in nonparametric regression. This is because the pattern of the relationship between the response variable and the predictor variable tends to change its behavior in sub-intervals, and there is also a pattern of relationship between the response variable and the predictor variable that repeats itself following an upward trend. The general model of the mixed regression Smoothing Spline and Fourier Series for a case study of monthly household expenditure per capita is as follows: In matrix form this can be written as ỹ = μ + ε.
3.705 × 10 9 · · · 2.84 × 10 9   , d = 4.49 × 10 −7 ; 4.06 × 10 −7 ; · · · ; −6.61 × 10 −7 Oscillation parameters for K = 1 and the smoothing parameter were λ 1 = 0.03; λ 2 = 0.03; λ 3 = 0.03. From the model above, we obtained a plot of the data estimation results using the Smoothing Spline and Fourier Series mixed estimators. The plot of the goodness of the mixed Smoothing Spline and Fourier Series model can be seen in Figure 6. Figure 6 shows that the data pattern approach using a mixed estimate of the Smoothing Spline and Fourier Series in nonparametric regression is good. This is because the estimation curve of mixed Smoothing Spline and Fourier Series estimator is in visual accordance with the data, or close to the actual data. Furthermore, based on the mixed model of the Smoothing Spline and the Fourier Series shown above, we obtained the scatter plot between the y value and the valueŷ shown in Table 7. Figure 7 shows that the predicted value of monthly household expenditure per capita (ŷ) was relatively close to the actual monthly expenses (y) for each poor household in Karangasem Regency, Bali Province. In addition, the mixed model of the Smoothing Spline and Fourier Series had a relatively large coefficient of determination, namely 98.99%. This shows that the mixed model of the Smoothing Spline and Fourier Series is suitable to be used in modeling the monthly per capita household expenditure of poor households in Karangasem Regency, Bali Province.
Oscillation parameters for K = 1 and the smoothing parameter were     Based on Figure 6 and the best model, the following can be concluded: 1. At wages/salaries above Rp. 250,000, there is a decrease in the monthly expenditure of poor households, whereas the characteristics of these households are dominant in consumption of carbohydrates and fat compared to protein consumption. 2. Consumption of fat below 90 g has decreased the expenditure of poor households; in this position the household consumes more carbohydrates than fat and protein, and the household has a wage/salary income below Rp. 120,000 per month. 3. Consumption of carbohydrates above 1700 g has decreased the expenditure of poor households per month; these households have the characteristic of consuming protein. These households have an income above Rp. 270,000. 4. To ensure household consumption is nutritious, safe, and balanced, and that the consumption of carbohydrates, fats, and proteins is in accordance with the recommendations of the PUGS of the Ministry of Health, poor households must increase their carbohydrate consumption to 5280.75 g, fat consumption to 1429.89 g, and protein Based on Figure 6 and the best model, the following can be concluded: 1.
At wages/salaries above Rp. 250,000, there is a decrease in the monthly expenditure of poor households, whereas the characteristics of these households are dominant in consumption of carbohydrates and fat compared to protein consumption.

2.
Consumption of fat below 90 g has decreased the expenditure of poor households; in this position the household consumes more carbohydrates than fat and protein, and the household has a wage/salary income below Rp. 120,000 per month.

3.
Consumption of carbohydrates above 1700 g has decreased the expenditure of poor households per month; these households have the characteristic of consuming protein.
These households have an income above Rp. 270,000.

4.
To ensure household consumption is nutritious, safe, and balanced, and that the consumption of carbohydrates, fats, and proteins is in accordance with the recommendations of the PUGS of the Ministry of Health, poor households must increase their carbohydrate consumption to 5280.75 g, fat consumption to 1429.89 g, and protein consumption to 1542.32 g. Income (wages/salaries) of 441,000 IDR per month will support the expenditure of poor households of 436,113 IDR, thus positioning the poor households above the poverty line. In efforts to meet the nutritional intake of poor households, the government must provide assistance in the form of basic food items that can meet the community's nutrition needs, in terms of carbohydrates, proteins, and fats.
Based on the results of this study, the best model was obtained using the smoothing parameters of λ 1 = 0.03; λ 2 = 0.03; λ 3 = 0.03, and the oscillation parameter of K = 1. It was found that, in certain positions, greater consumption of fat does not always increase the expenditure because, in this position, very poor households receive assistance from the government. Very poor households are characterized by consuming more carbohydrates than fat and protein. The results of this study suggest that the level of assistance provided to poor households should not be the same because in certain groups there are poor households that consume higher carbohydrates, higher fat, or higher protein. In addition, in an effort to increase the income of poor households, the government can provide assistance to create jobs by opening small and medium businesses, and provide counseling to the community by employing scholars to construct villages. The Karangasem district government can create community empowerment programs, such as the independent National Community Empowerment Program, the Simantri Program, and the Joint Business Group Program. These schemes also include the Savings and Loans program implemented by women funded by the independent National Community Empowerment Program, which is often referred to as the Women's Savings and Loans Program. Poverty alleviation programs are referred to as the Integrated Village Development Movement. Furthermore, programs classified as productive include People's Business Credit, Unsecured Credit, and the Bali Mandara Credit Guarantee.
A poverty alleviation program can not only reduce the number of poor people, but also reduce the number of unemployed people around places of business. This program is said to be a truly productive program that aims to increase the income of the poor. Several government policies that relate to poverty alleviation programs can be implemented to address poverty in Karangasem Regency, Bali Province. These include reducing the expenditures of the poor, as implemented by other parties such as the government or other communities, and increasing the income of the poor so that they can escape from poverty. Several policies aimed at reducing expenditures for the poor include the Raskin Program, the Direct Cash Assistance Program, basic food assistance, the Bali Mandara Health Insurance Program, the Community Health Insurance, and poor scholarships, namely, scholarships provided to children or students who have dropped out of school.

Conclusions
The The best model was obtained using the smoothing parameters of 0.03, 0.03, and 0.03, and the oscillation parameter of K = 1. It was found that, in a certain position, the greater consumption of fat does not always increase the household's expenditure because, in this position, very poor households receive assistance from the government. Very poor households are characterized by consuming more carbohydrates than fat and protein.
Based on the results of this study, it is suggested that the distribution of assistance to poor households should not be the same, because some poor households consume a higher level of carbohydrates, and others consume a higher level of fats. and for each function v j ∈ W 1 as follows: v j = ς j1 ψ j1 + ς j2 ψ j2 + . . . + ς jn ψ jn , where c j and d j are constants. Thus, for each function g j ∈ W can be described as follows: Based on Riesz's Representation Theorem, by describing L x as a linear function in the space W and g j ∈ W, then Equation (A1) can be written as follows: = L x u j + L x v j , = u j (x ji ) + v j (x ji ), = g j (x ji ). L x is a linear function limited to the space W, so that a single value is obtained; η i ∈ W is a representation of L x and can satisfy the equation: where ., . is the inner product. Based on properties of the inner product, then Equation (A2) can be written as: Based on Equation (A3), for i = 1 it can be stated as follows: = c j1 η j1 , θ j1 + c j2 η j1 , θ j2 + . . . + c jm η j1 , θ jm + d 1 τ j η j1 , ψ j1 + d 2 τ j η j1 , ψ j2 + . . . + d 2 τ j η j1 , ψ jn .