---
title: "SUMMARY AND GRAPH REGARDING PENGUINS IN PALMERPENGUINS PACKAGE"
output: 
  flexdashboard::flex_dashboard:
    orientation: columns
    vertical_layout: scroll
    social: ["menu"]
    source_code: embed
---

```{r setup, include=FALSE}
library(flexdashboard)
library(tidyverse)
library(palmerpenguins)
library(dplyr)
library(plotly)
library(gt)
library(flextable)
data("penguins")

```
Graph
======================
Column {data-width=650}
-----------------------------------------------------------------------

### scatter-plot of flipper depth vs body mass

```{r}
a<-penguins%>%ggplot(aes(x=bill_depth_mm, y=body_mass_g, fill = body_mass_g, colour = species))+geom_point()+theme_classic()

plotly::ggplotly(a)
```

Column {data-width=350}
-----------------------------------------------------------------------

### Box-plot of flipper length by sex

```{r}
penguins%>%ggplot(aes(x=flipper_length_mm, y=sex, fill = sex))+geom_boxplot()+theme_gray()
```

### Histogram of body mass by island

```{r}
penguins%>%ggplot(aes(x=body_mass_g,fill=island))+geom_histogram()+facet_grid(~island)+theme_light()
```

Data
==========================================================================
### palmer penguins data
```{r}
DT::datatable(penguins)
```
Summary
================
>Mean and standard deviation that was separated by each nominal and ordinal data.

### Data summary by species
```{r}
species1 <- na.omit(penguins) %>%
  group_by(species) %>%
  summarise(
    mean_bill_length = mean(bill_length_mm, na.rm = TRUE),
    sd_bill_length = sd(bill_length_mm, na.rm = TRUE),
    mean_bill_depth = mean(bill_depth_mm, na.rm = TRUE),
    sd_bill_depth = sd(bill_depth_mm, na.rm = TRUE),
    mean_flipper_length = mean(flipper_length_mm, na.rm = TRUE),
    sd_flipper_length = sd(flipper_length_mm, na.rm = TRUE),
    mean_body_mass = mean(body_mass_g, na.rm = TRUE),
    sd_body_mass = sd(body_mass_g, na.rm = TRUE)
  )

species1_flextable <- flextable(species1) %>%
  theme_vanilla() %>%
  set_caption("Summary Statistics by Penguin Species")

species1_flextable

```

### Data summary by sex
```{r}
library(dplyr)
library(gt)

Sex1 <- na.omit(penguins) %>%
  group_by(sex) %>%
  summarise(
    mean_bill_length = mean(bill_length_mm, na.rm = TRUE),
    sd_bill_length = sd(bill_length_mm, na.rm = TRUE),
    mean_bill_depth = mean(bill_depth_mm, na.rm = TRUE),
    sd_bill_depth = sd(bill_depth_mm, na.rm = TRUE),
    mean_flipper_length = mean(flipper_length_mm, na.rm = TRUE),
    sd_flipper_length = sd(flipper_length_mm, na.rm = TRUE),
    mean_body_mass = mean(body_mass_g, na.rm = TRUE),
    sd_body_mass = sd(body_mass_g, na.rm = TRUE)
  )

Sex1_flextable <- flextable(Sex1) %>%
  theme_vanilla() %>%
  set_caption("Summary Statistics by Penguin Sex")

Sex1_flextable
```

### Data summary by islands
```{r}
island1 <- na.omit(penguins) %>%
  group_by(island) %>%
  summarise(
    mean_bill_length = mean(bill_length_mm, na.rm = TRUE),
    sd_bill_length = sd(bill_length_mm, na.rm = TRUE),
    mean_bill_depth = mean(bill_depth_mm, na.rm = TRUE),
    sd_bill_depth = sd(bill_depth_mm, na.rm = TRUE),
    mean_flipper_length = mean(flipper_length_mm, na.rm = TRUE),
    sd_flipper_length = sd(flipper_length_mm, na.rm = TRUE),
    mean_body_mass = mean(body_mass_g, na.rm = TRUE),
    sd_body_mass = sd(body_mass_g, na.rm = TRUE)
  )

island1_flextable <- flextable(island1) %>%
  theme_vanilla() %>%
  set_caption("Summary Statistics by Penguin Species")

island1_flextable
```
### Data summary by year
```{r}
year1 <- na.omit(penguins) %>%
  group_by(year) %>%
  summarise(
    mean_bill_length = mean(bill_length_mm, na.rm = TRUE),
    sd_bill_length = sd(bill_length_mm, na.rm = TRUE),
    mean_bill_depth = mean(bill_depth_mm, na.rm = TRUE),
    sd_bill_depth = sd(bill_depth_mm, na.rm = TRUE),
    mean_flipper_length = mean(flipper_length_mm, na.rm = TRUE),
    sd_flipper_length = sd(flipper_length_mm, na.rm = TRUE),
    mean_body_mass = mean(body_mass_g, na.rm = TRUE),
    sd_body_mass = sd(body_mass_g, na.rm = TRUE)
  )

year1_flextable <- flextable(year1) %>%
  theme_vanilla() %>%
  set_caption("Summary Statistics by Penguin Species")

year1_flextable
```