Exploratory Data Analysis

UNG MATH 3350 (online)

I. Initialization Block

Initializing RStudio

If you execute the code block below, RStudio will prompt you to install Mosaic if needed.

library(mosaic)

II. Exercises

1. Using the built-in R data frame ChickWeight which examines weight vs age of chicks on different diets, conduct an analysis of the variable weight. Include all relevant plots, and include titles and axis labels for your histogram and density plot.

Built-in R Datasets

R has dozens of built-in data frames like iris.
Type variable name only and execute code block.
Output shows headers, variable names and
observations. Search variable name in help tab
of bottom-right panel for details.

2. Using the built-in R data frame InsectSprays which examines the number of insects in an area depending upon type of insecticide used, conduct an analysis of the variable count. Include all relevant plots, and include titles and axis labels for your histogram and density plot.

3. Use the built-in R data frame cars which compares the numeric variables speed measured in miles per hours (mph) and stopping distance (variable is dist, measured in feet). The data is from the 1920’s, by the way. Create an xyplot for the two variables, and include axis labels on your plot.

4. Using the built-in R data frame mtcars which includes miles per gallon and ten other variables for 32 different models in 1974, create an xyplot for miles per gallon (mpg variable) vs. displacement (disp variable) using a grouping variable of transmission type, automatic vs. manual (am).

5. Using the built-in R data frame Dimes which compares the mass of a dime to the year in which is was minted, conduct an analysis of the variable mass. Include all relevant plots, and include titles and axis labels for your histogram and density plot.

III. Code Blocks

W = c( 28, 11, 18, 35, 36, 6, -38, 14, -19, 43, -14, -30, -16, -25, 0, 40, 16, -79, 3, 11)

favstats(W)

right = mean(W) + 2 * sd(W)
left = mean(W) - 2 * sd(W)
left
right

histogram(W)

histogram(W, width = 20, center = 10 , type = "count")

histogram(W,
     width = 10,
     center = 0,
     fit = "normal",
     main = "Histogram: Mandy's Winnings",
     xlab = "Dollars Won or Lost",
     ylab = "Sessions")

boxplot(W)

boxplot(W, horizontal = TRUE,
     main = "Boxplot: Mandy's Winnings",
     xlab = "Dollars Won or Lost")

densityplot(W)

boxplot(iris$Petal.Length, horizontal = TRUE)
stem(iris$Petal.Length)

faithful

xyplot(eruptions ~ waiting, data = faithful)

cor(eruptions ~ waiting, data = faithful)

favstats(~ waiting, data = faithful)

histogram(~ waiting, data = faithful)

boxplot(faithful, horizontal = TRUE)

densityplot(~ waiting, data = faithful)

favstats(~ eruptions, data = faithful)

histogram(~ eruptions, data = faithful)

densityplot(~ waiting, data = faithful)

Births78

xyplot(births ~ day_of_year, data=Births78)

xyplot(births ~ day_of_year, data = Births78, 
       groups = wday,
       main = "Births vs. Day of Year for 1978",
       xlab = "Day of Year",
       ylab = "Births"
       )