R/grandistrib.R

In G2Sd: Grain-Size Statistics and Description of Sediment

grandistrib <-
  function(
    x,
    main = "",
    scale = "fine",
    xlab = "Stations",
    ylab = "Percentage"
  ) {
    if (scale == "fine") {
      Descript <- .sedim.descript(x)
      Descript <- Descript |>
        pivot_longer(cols = -samples, names_to = "class", values_to = "value")
      Descript$class <- factor(
        Descript$class,
        c(
          "boulder",
          "vcgravel",
          "cgravel",
          "mgravel",
          "fgravel",
          "vfgravel",
          "vcsand",
          "csand",
          "msand",
          "fsand",
          "vfsand",
          "vcsilt",
          "silt"
        ),
        ordered = T
      )
    }

    if (scale == "large") {
      Descript <- .texture.sedim(x)
      Descript <- Descript |>
        pivot_longer(
          cols = -c(samples, texture),
          names_to = "class",
          values_to = "value"
        )
      Descript$class <- factor(
        Descript$class,
        c("Boulder", "Gravel", "Sand", "Mud"),
        ordered = T
      )
    }

    p <- ggplot(Descript, aes(x = samples, y = value, fill = class)) +
      geom_bar(stat = "identity", position = "stack") +
      theme_bw() +
      scale_fill_viridis_d() +
      xlab(xlab) +
      ylab(ylab) +
      ggtitle(main) +
      guides(fill = guide_legend(title = "Classes")) +
      theme(
        axis.title = element_text(face = "bold", size = 13),
        axis.text = element_text(size = 11),
        legend.text = element_text(size = 11),
        legend.title = element_text(face = "bold", size = 11)
      )

    print(p)
  }