Штриховая диаграмма: изменение оси x и добавление строки

#r #plot #zoo

#r #график #зоопарк

Вопрос:

У меня есть зоопарк с ежедневными данными, который выглядит следующим образом: > head(almorol) 1973-10-02 1973-10-03 1973-10-04 1973-10-05 1973-10-06 1973-10-07 183.9 208.2 153.7 84.8 52.5 35.5

Я хочу построить годовые итоги и скользящую среднюю, что я и сделал:

y<-apply.yearly(almorol, FUN=sum) plot(y, main="Annual totals - Tagus (Almorol)",ylab="Q (m3/s)") lines(rollapply(y, 10, mean, na.rm=TRUE), col="red", lwd=2)

Который работает нормально, но поскольку данные не являются непрерывными, линейный график некорректен. Если я делаю это с точками, за этим трудно уследить, поэтому я хотел сделать штриховую диаграмму. Однако, когда я это делаю barplot(y) , ось x отображается в формате yy-mm-dd вместо просто лет, и я не могу добавить строку со скользящим средним (это не выдает ошибку, это просто не отображается).

Спасибо за помощь!

Редактировать: > y 1973-12-31 1974-12-31 1975-12-31 1976-12-31 19224.20 103766.30 72180.90 55939.80 1977-12-31 1978-12-31 1979-12-31 1980-12-31 215905.00 231014.21 319481.02 58979.84 1981-12-31 1982-12-31 1983-12-31 1984-12-31 32931.17 67989.06 83920.62 99431.75 1985-12-31 1986-12-31 1987-12-31 1988-12-31 161357.10 82910.87 101154.81 147541.80 1989-12-31 1990-12-31 1991-12-31 1992-12-31 137684.21 134974.39 89039.02 20774.72 1993-12-31 1994-12-31 1995-12-31 1996-12-31 58916.95 55187.38 52361.77 271064.34 1997-12-31 1998-12-31 1999-12-31 2000-12-31 220510.88 125116.62 42170.95 103915.99 2001-12-31 2002-12-31 2003-12-31 2004-12-31 244513.59 66811.15 167458.93 67223.66 2005-12-31 2006-12-31 2007-12-31 2008-12-31 8906.31 116874.33 79105.30 40142.61 2009-12-31 2010-12-31 2011-05-31 43835.34 161491.77 73093.06

1. Вместо head(almorol) не могли бы вы, пожалуйста, показать dput(almorol) , чтобы мы могли подыграть?

2. dput (almorol) выдает сотни строк, но выглядит примерно так: structure(c(183.9, 208.2, 153.7, 84.8, 52.5,......), index = structure(c(1370, 1371, 1372, 1373, 1374, 1375....), class = "Date"), class = "zoo")

3. Пожалуйста, сократите это до осмысленного размера для воспроизводимого примера, например, «как задать хороший вопрос».

Ответ №1:

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

Функция updateusr в пакете TeachingDemos может использоваться для изменения масштаба координат в соответствии с данными, которые вы хотите добавить.

Другой альтернативой является использование функции plot с type='h' (и посмотрите на lwd и lend опции в ?par ), чтобы создать свой собственный график, подобный barplot, используя интересующие координаты, затем добавьте свою опорную линию.