Advertisement
Guest User

Untitled

a guest
Aug 20th, 2017
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.67 KB | None | 0 0
  1. bddendro<-read.table("bddendro.csv", h=T, sep = ";", dec = ",")
  2. bddendro$dbh_new<-(bddendro$cbh_init + (bddendro$dendro_length * 0.2))/pi
  3.  
  4. bddendro<- bddendro %>%
  5. filter(med != 0) %>%
  6. group_by(parc, tree) %>%
  7. mutate(dbh_increase = ifelse(dendro_length < lag(dendro_length), 0 ,dbh_new - lag(dbh_new))) %>%
  8. mutate(dbh_cumsum = cumsum(dbh_increase))
  9.  
  10. mutate(dbh_increase = ifelse(is.na(dbh_increase), 0, dbh_increase))
  11.  
  12. bddendro2 <- bddendro %>%
  13. filter(med != 0) %>%
  14. group_by(parc, tree) %>%
  15. arrange(dendro_length) %>%
  16. mutate(dbh_increase = ifelse(is.na(lag(dbh_new)), 0, dbh_new - lag(dbh_new))) %>%
  17. mutate(dbh_cumsum = cumsum(dbh_increase)) %>%
  18. ungroup()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement