Я нашел этот пост и решение о ежегодных темпах роста с использованием dplyr действительно полезным:
R годовая скорость изменения (темп роста) с лагом dplyr
Вместо Рассчитывая годовой темп изменений, я хочу рассчитать ежемесячные темпы роста для раз в два месяца данных за часть года, с мая по сентябрь. Я думаю, что я понял это, используя пост, упомянутый выше в качестве руководства. Вот воспроизводимый пример:
#build toy dataset. In contains the Plant ID, leaf length, Sampling Month, and Sampling Time within the month (either T1 or T2)
Plant_ID <- c("365","365","365","365","365","365","365","365","365","365")
Leaf_length <- c(4, 10, 15, 17, 20, 25, 30, 50, 45, 47)
Month <- c(5,5,6,6,7,7,8,8,9,9)
Period <- c("T1","T2","T1","T2","T1","T2","T1","T2","T1","T2")
toy_growthrate <- data.frame(Plant_ID, Leaf_length, Month, Period)
#look at dataset
toy_growthrate
#try to calculate bimonthly percentage change
toy_growthrate <- toy_growthrate %>% mutate(change=(Leaf_length-lag(Leaf_length,2))/lag(Leaf_length,2)*100)
#the new column "change" is filled with month-to-month growth
toy_growthrate
Тем не менее, я все еще застрял на том, как сделать этот ежемесячный расчет роста, если в определенные месяцы у меня есть данные за два месяца, но в другие месяцы у меня есть 4 точки данных (т.е. еженедельные данные)? Нужно ли мне брать средние значения еженедельных данных, чтобы преобразовать их в данные за два месяца, чтобы все месяцы имели одинаковое количество точек данных?
Вот еще один пример кода с этим новым поворотом:
#build toy dataset. In contains the Plant ID, leaf length, Sampling Month, and Sampling Time within the month (T1, T2, T3, T4 in May, T1 and T2 in remaining months)
Plant_ID <- c("365","365","365","365","365","365","365","365","365","365","365","365")
Leaf_length <- c(1,2,4, 10, 15, 17, 20, 25, 30, 50, 45, 47)
Month <- c(5,5,5,5,6,6,7,7,8,8,9,9)
Period <- c("T1","T2","T3","T4","T1","T2","T1","T2","T1","T2","T1","T2")
toy_growthrate_with_twist <- data.frame(Plant_ID, Leaf_length, Month, Period)
#look at dataset
toy_growthrate_with_twist
В этом новом наборе данных имеется 4 измерения длины листа в мае, но только 2 в оставшиеся месяцы. Как я могу сделать подсчет роста в этом месяце?
Спасибо!