Box plot command ggplot geom package basic ggplot

# 1. START OF SESSION '
# 1a. Preliminaries'
# setwd( ) to set working directory'
# Tips: 1) R wants forward slashes and 2) must enclose in quotes' setwd("/Users/cbigelow/Desktop/")'
# 1b. Clear workspace '
rm(list=ls())'
# 1c. Turn OFF scientific notation'
options(scipen=1000)'

library(ggplot2)

## id matage hyp gestwks sex bweight' ## 636 636 34 0 41.15 male 2972' ## 637 637 28 0 38.58 female 2850' ## 638 638 38 1 38.44 male 3182' ## 639 639 26 0 38.94 female 3048' ## 640 640 31 0 40.43 female 3183' ## 641 641 31 0 38.15 male 2920'

# With aesthetics: I recommend building step by step, then displaying!' p <-ggplot(data=ivf, aes(x=factor(hyp)))'
p1 <-p +geom_bar(color="black", fill="blue")'
p2 <-p1 +xlab("Hypertension") '
p3 <-p2 +ylab("Frequency")'
p4 <-p3 +ggtitle("Bar Graph of Hypertension")'
p5 <-p4 +theme_bw()'
p5'

# 4a.2. Quantiles'
# Command quantile() in package=mosaic' quantile(~bweight, data=ivf)'

## 0% 25% 50% 75% 100% '
## 630 2850 3200 3550 4650'

# 4b. Graphical'
# 4b.1. Box plot'
# Command ggplot() + geom_boxplot() in package=ggplot2' # Basic'
ggplot(data=ivf, aes(x=1, y=matage)) +geom_boxplot() '

# 4b.2 Histogram'
# Command ggplot() + geom_histogram() in package=ggplot2' # Basic'
ggplot(data=ivf, aes(x=matage)) +geom_histogram() '

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.'

# With aesthetics - step by step'
p18 <-ggplot(data=ivf, aes(x=matage)) +geom_histogram(color="black", fill="blue", binwidth=2)' p19 <-p18 +xlab("Maternal Age (years)")'
p20 <-p19 +ylab("Frequency")'
p21 <-p20 +ggtitle("Histogram of Maternal Age")'
p22 <-p21 +theme_bw()'
p22'

# 5. TWO VARIABLES - ONE CONTINUOUS, ONE DISCRETE'
# 5a. Numerical'
# 5a.1. Detailed descriptives, by group'
# Command describeBy() in package=psych'
# Tip - Declare your discrete group variable to be factor' library(psych)'
ivf$sex <-as.factor(ivf$sex)'
describeBy(ivf$bweight, group = ivf$sex,digits=4)'

# 5b.3. Histogram, by group - Overlay '
# Command ggplot() + geom_histogram() in package=ggplot2' # With aesthetics - step by step'

p32