Приложения
Приложение А — исходный код.
Приложение а
library(forecast) #прогнозирование
library("TTR") #сглаживание
library("ggplot2") #графики
library("sophisthse") #загрузка
с sophist.hse.ru
wag <- sophisthse("WAG_M")
wag
gov <- sophisthse("GOV_M")
gov
ggplot(data = as.data.frame(wag)) +
geom_point(mapping = aes(x = time(wag), y =
WAG_C_M)) +
labs(title = 'Заработная плата',
x = 'Время',
y = 'Заработная плата, руб.')
ggplot(data = as.data.frame(gov)) +
geom_point(mapping = aes(x = time(gov), y =
FBREV_M)) +
labs(title = 'Доходы федерального бюджета',
x = 'Время',
y = 'Доход, млрд. руб.')
ggplot(data = as.data.frame(wag)) +
geom_line(mapping = aes(x = time(wag), y =
WAG_C_M)) +
labs(title = 'Заработная плата',
x = 'Время',
y = 'Заработная плата, руб.')
ggplot(data = as.data.frame(gov)) +
geom_line(mapping = aes(x = time(gov), y =
FBREV_M)) +
labs(title = 'Доходы федерального бюджета',
x = 'Время',
y = 'Доход, млрд.руб.')
ssw <- SMA(wag[,3], n=10)
ssw
ssg <- SMA(gov[,1], n=10)
ssg
ggplot(data = as.data.frame(wag)) +
geom_line(mapping = aes(x = time(wag), y =
WAG_C_M), color = 'blue') +
geom_line(mapping = aes(x = time(wag), y =
ssw), color = 'red') +
labs(title = 'Заработная плата',
x = 'Время',
y = 'Заработная плата, руб.')
ggplot(data = as.data.frame(gov)) +
geom_line(mapping = aes(x = time(gov), y =
FBREV_M), color = 'blue') +
geom_line(mapping = aes(x = time(gov), y =
ssg), color = 'red') +
labs(title = 'Доходы федерального бюджета',
x = 'Время',
y = 'Доход, млрд.руб.')
res <- auto.arima (ssw)
forecasted <- forecast(res, 35)
forecasted
plot(forecasted,
xlab = "Время",
ylab = "Заработная плата",
main = "Прогноз доходов населения до
2025 года")
