Переназначить начальную юлианскую дату (1 июля = юлианская дата 1, южное полушарие)

Question

Мой набор данных включает в себя различные наблюдения на разных этапах в течение года.

• year когда образцы были собраны.
• site место измерения
• Class физическая стадия во время измерения
• date дата измерения
• Julian Юлианская дата

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

year site Class date Julian 1 2009 10C Early 2008-09-15 259 2 2009 10C L2 2008-09-29 273 3 2009 10C L3 2008-12-15 350 4 2010 10C Early 2009-08-31 243 5 2010 10C L2 2009-09-14 257 6 2010 10C L3 2009-12-11 345 7 2012 10C Early 2011-08-23 235 8 2012 10C L2 2011-09-22 265 9 2012 10C L3 2011-12-03 337 10 2012 10C LSample 2012-03-26 86 11 2013 10C Early 2012-09-07 251 12 2013 10C L2 2012-09-30 274 13 2013 10C L3 2012-12-17 352 14 2014 10C Early 2013-09-02 245 15 2014 10C L2 2013-09-16 259 16 2014 10C L3 2013-12-16 350 17 2014 10C LMid 2014-01-07 7 18 2015 10C Early 2014-09-08 251 19 2015 10C L2 2014-09-30 273 20 2015 10C L3 2014-12-01 335

Мне трудно конвертировать / переназначать юлианскую дату начала на 1 июля вместо 1 января. Точечный график ниже иллюстрирует окончательную выборку, которая происходит в начале года (февраль-март).chron Пакет имеет возможность изменить порядок origin но я не могу заставить его работать должным образом с моими данными.

library(chron) library(dplyr) data.date <- data %>% mutate(July.Julian = chron(date,format = c(dates = "ymd"), options(chron.origin = c(month=7, day=1, year=2008)))) Error in chron(c("2008-09-15", "2008-09-29", "2008-12-15", "2009-08-31", : misspecified chron format(s) length

или же

July.Julian = chron(data$date, format = c(dates = "ymd"), options(chron.origin = c(month=7, day=1, year=2008))) Error in chron(c("2008-09-15", "2008-09-29", "2008-12-15", "2009-08-31", : misspecified chron format(s) length

Я пытаюсь начать юлианскую дату с 1 вместо 182.

Мысли или предложения приветствуются.

Answer 1

При условии, что July.Julian должен быть Julian дни после 1 июля:

transform(date.data, July.Julian = as.chron(sprintf("%d-07-01", year)) + Julian)

или же

date.data %>% mutate(July.Julian = as.chron(sprintf("%d-07-01", year)) + Julian)

Обратите внимание, что здесь на самом деле не нужен хрон. Просто замени as.chron с as.Date и любой из этих работ.