R Assignment

load(“parenthood.Rdata”)

describe(parenthood)

           vars   n  mean    sd median trimmed   mad   min
dan.sleep     1 100  6.97  1.02   7.03    7.00  1.09  4.84
baby.sleep    2 100  8.05  2.07   7.95    8.05  2.33  3.25
dan.grump     3 100 63.71 10.05  62.00   63.16  9.64 41.00
day           4 100 50.50 29.01  50.50   50.50 37.06  1.00
              max range  skew kurtosis   se
dan.sleep    9.00  4.16 -0.29    -0.72 0.10
baby.sleep  12.07  8.82 -0.02    -0.69 0.21
dan.grump   91.00 50.00  0.43    -0.16 1.00
day        100.00 99.00  0.00    -1.24 2.90

The mean hours show some variance, ranging from 6.97 to 63.71. The standard deviation of dan.sleep and baby. sleep are fairly small while dan.grump is fairly big. The max value dan.sleep reaches is 9.00. The max value for baby.sleep is 12.07. The max value for dan.grump is 91.00. The minimum value for dan.sleep is 4.84. The minimum value for baby.sleep is 3.25. The minimum value of dan.grump is 41.00. The range for dan.sleep is 4.16. The range for baby.sleep is 12.07. The range for dan.grump is 91.00. Dan.sleep and baby.sleep is negatively skewed, skewed to the left. Dan.grump is positively skewed, skewed to the right.

par(mfrow=c(1,3))
hist(parenthood$dan.sleep, 
     xlim=c(2,14), 
     breaks=seq(0,14,by=1),
     main = "# hours father sleeping", 
     xlab = "# hours",
     ylab = "Frequency",
     col= "Blue"
)
hist(parenthood$baby.sleep,
     xlim=c(2,14), 
     breaks=seq(0,14,by=1),
     main = "# hours baby sleep", 
     xlab = "# hours",
     ylab = "Frequency",
     col= "Green"
)

hist(parenthood$dan.grump,
    xlim=c(0,100), 
     breaks=seq(0,100,by=1),
     main = "How much father is grumpy", 
     xlab = "How much Dan grumpy",
     ylab = "Frequency",
     col= "Red"
  )

The histogram for dan.sleep looks unimodal and skewed to the left. The histogram for baby.sleep looks bimodal and skewed to the left. The histogram for dan.grump looks to be multimodal with a skew to the right.

par(mfrow=c(1,3))
boxplot(parenthood$dan.sleep, 
        ylim=c(0,14),
        main = "# hours father is sleeping", 
        ylab = "#hours",
        col= "Blue"
)
boxplot(parenthood$baby.sleep,
        ylim=c(0,14),
        main = "# hours baby is sleeping", 
        ylab = "# hours",
        col= "Green"
)

boxplot(parenthood$dan.grump,
    ylim=c(0,100), 
     main = "How much father grumpy", 
     ylab = "Frequency",
     col= "Red")

The boxplot for dan.sleep shows a median around 7 hours. No outliers are present for this boxplot. The boxplot for baby.sleep shows a median around 8 hours. No outliers are present for this boxplot. The boxplot for dan.grump shows a median of around 60 with no outliers present.