R – adding group value into result of auto.arima
If you want to use arima with data that are divided into groups, you can do something like this.
data <- read.csv("some_input.txt", header=TRUE, sep="\t") split_data <- split( data, data$group ) for (name in names(split_data)) { pred <- c(spt1[[name]][['value']]) indata <- stats::ts( pred ) result<-forecast(auto.arima(indata,seasonal=FALSE),seasonal=FALSE,h=12) print(result) }
Result will be grouped in sections, like this
[1] "1" Point Forecast Lo 80 Hi 80 Lo 95 Hi 95 [1] "2" Point Forecast Lo 80 Hi 80 Lo 95 Hi 95 [1] "3" Point Forecast Lo 80 Hi 80 Lo 95 Hi 95 [1] "4" Point Forecast Lo 80 Hi 80 Lo 95 Hi 95 [1] "5" Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
If you want to add group No. into result, you can easily achieve that with awk
cat result | awk 'BEGIN{regex="^\\[";} \ { if (match($0, regex)) { split($0,arr,"\""); number = arr[2]} \ else \ { if ($0 ~ /Point/) {} \ else \ { print $0,number } }}'
June 29th, 2015 in
main entries