Advertisement
Guest User

Untitled

a guest
Jun 19th, 2019
73
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.67 KB | None | 0 0
  1. data = read.table(text = "
  2. start end av y1 y2 y3
  3. 719000 721000 720000 1 1 2
  4. 1199000 1201000 1200000 1 1 2
  5. 3039000 3041000 3040000 1 1 2
  6. 3679000 3681000 3680000 1 1 2
  7. 4119000 4121000 4120000 1 1 2
  8. 4999000 5001000 5000000 1 1 2", sep = " ", stringsAsFactors = F, header = T)
  9.  
  10. ids <- factor(paste(1:nrow(data), 1:nrow(data), sep = "."))
  11.  
  12. positions <- data.frame(
  13. id = rep(ids, times = nrow(data)/2),
  14. x = c(data$start, data$end, data$av),
  15. y = c(data$y1, data$y2, data$y3)
  16. )
  17.  
  18. values <- data.frame(
  19. id = ids,
  20. value = 1:length(ids)
  21. )
  22.  
  23. datapoly <- merge(values, positions, by = c("id"))
  24.  
  25. ggplot(datapoly, aes(x = x, y = y)) + geom_polygon(aes(fill = value, group = id))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement