Выписка Дата, прогноз, пакет R

Question

Я использую прогнозный пакет и сохранил результаты

t <- data.frame(forecast(prod.arima, h=4))

к фрейму данных. Названия строк представляют мои даты, которые при печати выглядят как

> t
         Point.Forecast    Lo.80    Hi.80    Lo.95     Hi.95
2010.856       812.9849 630.5707 995.3992 534.0064 1091.9634
2010.876       670.3363 485.1885 855.4840 387.1772  953.4953
2010.895       769.4848 584.2552 954.7145 486.2005 1052.7691
2010.914       692.8951 507.6630 878.1272 409.6070  976.1832

Я не близок (пока) к тому, чтобы завершить это, но я надеюсь автоматизировать отчет LaTex, где я буду печатать прогнозы на странице отчета и графике прогноза. Я полагаю, что могу передать кадр данных в таблицу LaTex, но мне нужно, чтобы даты были отформатированы таким образом, чтобы конечный пользователь понял, что они представляют выходные дни (установлено по воскресеньям).

Я преобразовал свой исходный временной ряд с «читаемыми» датами в объект временного ряда, используя

ts(prod, start = 2009 +(31+28+31+5)/365, f=52 )

, где начало моих временных рядов было 5 апреля 2009 года и представляло еженедельные точки данных.

Буду очень признателен за любую помощь, которую вы сможете оказать. Излишне говорить, что я довольно новичок в R.

Answer 1

Я могу придумать два пути: один из них - использовать функцию as.Date.cal.yr из пакета Epi.Другой вариант - умножить эти даты на количество секунд в году, а затем передать их как .POSIXct (но я этого не проверял.) Из этих двух первый кажется наименее болезненным:

require(Epi)
as.Date.cal.yr(as.numeric(row.names(t)))
# [1] "2010-11-10" "2010-11-17" "2010-11-24" "2010-12-01"
as.POSIXlt(as.Date.cal.yr(as.numeric(row.names(t))))$wday
#  [1] 3 3 3 3
paste(weekdays(as.Date.cal.yr(as.numeric(row.names(t)))) , 
         as.Date.cal.yr(as.numeric(row.names(t))), sep=",")
 #   [1] "Wednesday,2010-11-10" "Wednesday,2010-11-17" "Wednesday,2010-11-24"
 #   [4] "Wednesday,2010-12-01"

Так что я не получаю воскресенья, так как это будет wday == 0 в системе POSIXlt.См. Справку для DateTimeClasses.

Что касается части вопроса о LaTeX, установите пакет Hmisc и используйте функцию latex ().