library(openintro)
library("tidyverse")
library(palmerpenguins)
sum(is.na(penguins))
## [1] 19
penguins %>%
drop_na() -> OKpenguins
sum(is.na(OKpenguins))
## [1] 0
OKpenguins %>%
group_by(species) %>%
summarize(mean.bl=mean(bill_length_mm), sd.bl=sd(bill_length_mm),
mean.bd=mean(bill_depth_mm), sd.bd=sd(bill_depth_mm),
n=n()) %>%
print.data.frame(., digits=3)
## species mean.bl sd.bl mean.bd sd.bd n
## 1 Adelie 38.8 2.66 18.3 1.219 146
## 2 Chinstrap 48.8 3.34 18.4 1.135 68
## 3 Gentoo 47.6 3.11 15.0 0.986 119
ggplot(OKpenguins, aes(x=species, y=bill_depth_mm, fill=species)) +
geom_boxplot(outlier.colour = "blue", outlier.size=2) +
geom_jitter(width=0.2) +
labs(title="Species Bill Depth") +
theme_light()
## Exercies 4
OKpenguins %>%
group_by(species) %>%
summarize(r.bill =cor(bill_length_mm, bill_depth_mm), r.flipper.bodymass=cor(flipper_length_mm, body_mass_g)) %>%
print.data.frame(., digits=3)
## species r.bill r.flipper.bodymass
## 1 Adelie 0.386 0.465
## 2 Chinstrap 0.654 0.642
## 3 Gentoo 0.654 0.711
ggplot(OKpenguins, aes(x=bill_length_mm, y=bill_depth_mm)) +
geom_point(aes(colour=species), size= 2.5)+
labs(title = "Bill Dimensions") +
theme_light()
## Excerise 6