Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- bddendro<-read.table("bddendro.csv", h=T, sep = ";", dec = ",")
- bddendro$dbh_new<-(bddendro$cbh_init + (bddendro$dendro_length * 0.2))/pi
- bddendro<- bddendro %>%
- filter(med != 0) %>%
- group_by(parc, tree) %>%
- mutate(dbh_increase = ifelse(dendro_length < lag(dendro_length), 0 ,dbh_new - lag(dbh_new))) %>%
- mutate(dbh_cumsum = cumsum(dbh_increase))
- mutate(dbh_increase = ifelse(is.na(dbh_increase), 0, dbh_increase))
- bddendro2 <- bddendro %>%
- filter(med != 0) %>%
- group_by(parc, tree) %>%
- arrange(dendro_length) %>%
- mutate(dbh_increase = ifelse(is.na(lag(dbh_new)), 0, dbh_new - lag(dbh_new))) %>%
- mutate(dbh_cumsum = cumsum(dbh_increase)) %>%
- ungroup()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement