Descriptives

source("0_helpers.R")

Load data

msleep = readRDS("msleep_processed.rds")

Descriptives

msleep %>% 
    select(sleep_total, sleep_rem, brainwt, bodywt) %>% 
    summarise_all_vars()

variable	mean	sd	min	max	n_miss	n_nonmiss	n_distinct
sleep_total	10.43	4.45	1.9	19.9	0	83	65
sleep_rem	1.875	1.298	0.1	6.6	22	61	33
brainwt	0.2816	0.9764	0.00014	5.712	27	56	54
bodywt	166.1	786.8	0.005	6654	0	83	82

Distributions

qplot(eats, data = msleep)

qplot(sleep_total, data = msleep, binwidth = 0.5)

qplot(sleep_rem, data = msleep, binwidth = 0.5)

## Warning: Removed 22 rows containing non-finite values (stat_bin).

qplot(sleep_total, sleep_rem, data = msleep)

## Warning: Removed 22 rows containing missing values (geom_point).

Look at the data

library(DT)
datatable(msleep, 
    filter = 'top', options = list(
  pageLength = 50, autoWidth = TRUE
), rownames = F)

IyBEZXNjcmlwdGl2ZXMgey50YWJzZXQgLnRhYnNldC1zdGlja3l9CgpgYGB7ciBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQpzb3VyY2UoIjBfaGVscGVycy5SIikKYGBgCgojIyBMb2FkIGRhdGEKCmBgYHtyfQptc2xlZXAgPSByZWFkUkRTKCJtc2xlZXBfcHJvY2Vzc2VkLnJkcyIpCmBgYAoKIyMgRGVzY3JpcHRpdmVzIHsuYWN0aXZlfQpgYGB7ciB3YXJuaW5nPUZBTFNFfQptc2xlZXAgJT4lIAoJc2VsZWN0KHNsZWVwX3RvdGFsLCBzbGVlcF9yZW0sIGJyYWlud3QsIGJvZHl3dCkgJT4lIAoJc3VtbWFyaXNlX2FsbF92YXJzKCkKYGBgCgojIyBEaXN0cmlidXRpb25zCgpgYGB7cn0KcXBsb3QoZWF0cywgZGF0YSA9IG1zbGVlcCkKcXBsb3Qoc2xlZXBfdG90YWwsIGRhdGEgPSBtc2xlZXAsIGJpbndpZHRoID0gMC41KQpxcGxvdChzbGVlcF9yZW0sIGRhdGEgPSBtc2xlZXAsIGJpbndpZHRoID0gMC41KQpxcGxvdChzbGVlcF90b3RhbCwgc2xlZXBfcmVtLCBkYXRhID0gbXNsZWVwKQpgYGAKCiMjIExvb2sgYXQgdGhlIGRhdGEKCmBgYHtyfQpsaWJyYXJ5KERUKQpkYXRhdGFibGUobXNsZWVwLCAKCWZpbHRlciA9ICd0b3AnLCBvcHRpb25zID0gbGlzdCgKICBwYWdlTGVuZ3RoID0gNTAsIGF1dG9XaWR0aCA9IFRSVUUKKSwgcm93bmFtZXMgPSBGKQpgYGAKCg==