# Script to generate pictures of simulated density plots of normal distributions

# Wipes any pictures you generated
close.screen(,all.screens=T)

# Makes 4 Panels for plotting
split.screen(c(2,2))

#Assigns to Screen 1
screen(1)

# Plot.density plots a density object. Density evaluates your data (similar to histograms)
# rnorm randomly creates data following a normal distribution with default mean and stand. dev.
# the only argument is the size of simulation 100

plot(density(rnorm(100)),xlab="Trait Value",ylab="Density",main="100 observations")

screen(2)

plot(density(rnorm(1000)),xlab="Trait Value",ylab="Density",main="1,000 observations")

screen(3)

plot(density(rnorm(10000)),xlab="Trait Value",ylab="Density",main="10,000 observations")

screen(4)

plot(density(rnorm(1000000)),xlab="Trait Value",ylab="Density",main="1,000,000 observations")