Clase del día jueves 28
Profesor: Alejandro Acevedo A., Ph.D. Clase del día jueves 28 Establecer directorio de trabajo:
setwd(“~/clase_28_mayo_2020”)
library(tidyverse)
library(GGally)
library(ggsci)
library(ggplot2)
library(gridExtra)
library(autoplotly)
library(ggmuller)
library(gapminder)
library(gganimate)
library(gifski)
library(av)
library(webshot)data('psychademic')
data('diamonds')
#data('iris')
sets.de.datos <- list()
sets.de.datos[['sicoacademico']] <- psychademic
sets.de.datos[['diamantes']] <- diamonds
#sets.de.datos[['flores']] <- iris
sets.de.datos$sicoacademico %>% head()## locus_of_control self_concept motivation read write math science sex
## 1 -0.84 -0.24 4 54.8 64.5 44.5 52.6 female
## 2 -0.38 -0.47 3 62.7 43.7 44.7 52.6 female
## 3 0.89 0.59 3 60.6 56.7 70.5 58.0 male
## 4 0.71 0.28 3 62.7 56.7 54.7 58.0 male
## 5 -0.64 0.03 4 41.6 46.3 38.4 36.3 female
## 6 1.11 0.90 2 62.7 64.5 61.4 58.0 femalesets.de.datos$diamantes %>% head()## # A tibble: 6 x 10
## carat cut color clarity depth table price x y z
## <dbl> <ord> <ord> <ord> <dbl> <dbl> <int> <dbl> <dbl> <dbl>
## 1 0.23 Ideal E SI2 61.5 55 326 3.95 3.98 2.43
## 2 0.21 Premium E SI1 59.8 61 326 3.89 3.84 2.31
## 3 0.23 Good E VS1 56.9 65 327 4.05 4.07 2.31
## 4 0.290 Premium I VS2 62.4 58 334 4.2 4.23 2.63
## 5 0.31 Good J SI2 63.3 58 335 4.34 4.35 2.75
## 6 0.24 Very Good J VVS2 62.8 57 336 3.94 3.96 2.48#sets.de.datos$flores %>% head()mi.set = sets.de.datos$sicoacademico
mi.set %>% class()## [1] "data.frame"mi.set %>% dim()## [1] 600 8mi.set %>% length()## [1] 8mi.set %>% head()## locus_of_control self_concept motivation read write math science sex
## 1 -0.84 -0.24 4 54.8 64.5 44.5 52.6 female
## 2 -0.38 -0.47 3 62.7 43.7 44.7 52.6 female
## 3 0.89 0.59 3 60.6 56.7 70.5 58.0 male
## 4 0.71 0.28 3 62.7 56.7 54.7 58.0 male
## 5 -0.64 0.03 4 41.6 46.3 38.4 36.3 female
## 6 1.11 0.90 2 62.7 64.5 61.4 58.0 femalemi.set %>% as.matrix()## locus_of_control self_concept motivation read write math science
## [1,] "-0.84" "-0.24" "4" "54.8" "64.5" "44.5" "52.6"
## [2,] "-0.38" "-0.47" "3" "62.7" "43.7" "44.7" "52.6"
## [3,] " 0.89" " 0.59" "3" "60.6" "56.7" "70.5" "58.0"
## [4,] " 0.71" " 0.28" "3" "62.7" "56.7" "54.7" "58.0"
## [5,] "-0.64" " 0.03" "4" "41.6" "46.3" "38.4" "36.3"
## [6,] " 1.11" " 0.90" "2" "62.7" "64.5" "61.4" "58.0"
## [7,] " 0.06" " 0.03" "3" "41.6" "39.1" "56.3" "45.0"
## [8,] "-0.91" "-0.59" "3" "44.2" "39.1" "46.3" "36.3"
## [9,] " 0.45" " 0.03" "4" "62.7" "51.5" "54.4" "49.8"
## [10,] " 0.00" " 0.03" "3" "62.7" "64.5" "38.3" "55.8"
## [11,] " 0.24" "-0.43" "2" "70.7" "43.7" "58.8" "66.1"
## [12,] "-1.09" "-0.26" "2" "44.2" "41.1" "45.1" "47.1"
## [13,] " 0.46" " 0.03" "3" "57.4" "59.3" "53.9" "49.8"
## [14,] " 0.68" " 0.06" "3" "49.5" "51.5" "41.2" "41.7"
## [15,] "-0.14" "-1.05" "4" "70.7" "65.1" "66.4" "63.4"
## [16,] " 0.10" "-0.16" "2" "49.5" "59.3" "51.0" "47.1"
## [17,] " 0.45" " 0.65" "4" "57.4" "56.7" "46.9" "52.6"
## [18,] "-1.56" " 0.03" "2" "62.7" "54.1" "53.0" "55.3"
## [19,] "-1.02" "-1.67" "3" "45.3" "43.0" "52.3" "60.7"
## [20,] " 0.94" " 0.03" "4" "68.0" "56.7" "59.6" "58.0"
## [21,] "-0.40" " 0.03" "4" "44.2" "54.1" "59.3" "58.0"
## [22,] "-0.19" "-1.73" "4" "44.7" "47.6" "50.3" "39.0"
## [23,] " 0.26" " 0.03" "2" "48.4" "48.9" "52.2" "43.9"
## [24,] " 0.21" " 0.94" "1" "38.9" "48.9" "41.8" "58.5"
## [25,] " 0.70" " 0.87" "4" "46.9" "54.1" "46.4" "49.8"
## [26,] " 1.36" " 0.94" "4" "57.4" "64.5" "49.9" "55.3"
## [27,] "-0.84" "-0.57" "2" "33.6" "33.3" "41.0" "36.3"
## [28,] "-0.68" " 0.03" "2" "38.9" "38.5" "41.2" "39.0"
## [29,] "-0.17" " 0.31" "4" "65.4" "67.1" "75.5" "71.5"
## [30,] " 0.06" " 0.28" "4" "60.1" "59.3" "56.2" "60.7"
## [31,] "-0.14" " 0.56" "4" "54.8" "61.9" "54.6" "47.1"
## [32,] " 0.96" " 0.03" "3" "44.2" "51.5" "43.2" "44.4"
## [33,] " 0.23" "-0.26" "1" "52.1" "61.9" "56.3" "58.0"
## [34,] " 0.66" " 0.00" "4" "52.1" "48.9" "48.5" "44.4"
## [35,] " 0.08" " 0.59" "4" "57.4" "56.7" "53.9" "60.7"
## [36,] "-0.90" " 0.03" "3" "36.3" "44.3" "36.1" "33.6"
## [37,] " 0.06" " 0.03" "3" "62.7" "61.9" "58.8" "60.7"
## [38,] " 0.04" "-0.47" "3" "33.6" "59.3" "44.7" "47.1"
## [39,] " 0.66" " 0.34" "4" "60.1" "51.5" "53.0" "63.4"
## [40,] " 0.46" " 0.03" "4" "57.4" "67.1" "47.7" "55.8"
## [41,] "-0.14" "-0.24" "3" "54.8" "38.5" "57.0" "52.6"
## [42,] " 0.26" "-0.57" "4" "41.6" "58.6" "54.6" "55.3"
## [43,] " 0.02" "-0.14" "1" "60.1" "64.5" "67.0" "58.5"
## [44,] " 1.11" " 0.34" "4" "65.4" "64.5" "64.1" "66.1"
## [45,] " 1.16" "-0.81" "3" "65.4" "64.5" "48.0" "63.4"
## [46,] "-0.40" " 0.65" "3" "44.2" "43.7" "38.6" "33.6"
## [47,] "-0.38" " 0.34" "3" "38.9" "28.1" "35.3" "39.0"
## [48,] " 1.13" " 0.03" "4" "38.9" "54.1" "47.5" "49.8"
## [49,] " 0.08" " 0.94" "2" "49.5" "52.8" "50.6" "48.8"
## [50,] " 0.06" " 0.94" "4" "52.1" "51.5" "57.9" "60.7"
## [51,] "-0.64" "-0.51" "3" "49.5" "56.7" "47.7" "44.4"
## [52,] " 0.46" " 0.03" "4" "70.7" "56.7" "51.3" "71.5"
## [53,] "-1.13" "-0.55" "2" "44.2" "46.9" "45.5" "39.0"
## [54,] "-1.05" "-1.65" "1" "46.9" "46.3" "38.6" "47.1"
## [55,] "-0.84" " 0.03" "1" "36.3" "59.3" "47.7" "39.0"
## [56,] " 0.06" " 0.65" "2" "47.4" "64.5" "53.9" "55.8"
## [57,] "-1.30" " 0.13" "2" "44.2" "41.1" "51.8" "47.1"
## [58,] " 0.91" "-0.47" "3" "44.2" "41.1" "50.5" "58.0"
## [59,] " 0.28" " 0.03" "4" "54.8" "59.3" "49.5" "63.4"
## [60,] "-0.40" " 0.03" "1" "54.8" "49.5" "55.3" "47.1"
## [61,] "-0.86" " 0.28" "4" "36.3" "48.9" "54.4" "60.7"
## [62,] " 0.46" " 0.94" "4" "65.4" "61.9" "67.9" "66.1"
## [63,] " 0.50" "-0.22" "4" "54.8" "61.9" "61.3" "60.7"
## [64,] " 0.27" "-1.05" "2" "60.1" "54.1" "56.3" "55.3"
## [65,] "-0.86" " 1.19" "2" "33.6" "28.1" "31.8" "39.6"
## [66,] "-0.89" " 0.56" "2" "44.2" "41.1" "40.3" "49.8"
## [67,] "-1.28" " 1.19" "4" "41.6" "41.1" "45.9" "44.4"
## [68,] " 1.36" " 0.87" "4" "36.3" "43.7" "37.2" "41.7"
## [69,] " 0.46" " 0.65" "4" "49.5" "48.9" "60.5" "55.3"
## [70,] " 0.25" "-0.26" "2" "49.5" "30.7" "40.3" "33.6"
## [71,] " 0.23" " 0.03" "4" "62.7" "61.9" "55.5" "55.3"
## [72,] "-1.33" "-0.60" "2" "34.7" "35.2" "40.2" "50.9"
## [73,] " 0.32" "-1.17" "2" "52.1" "54.1" "58.1" "47.1"
## [74,] "-0.44" "-0.31" "3" "49.5" "59.3" "42.1" "52.6"
## [75,] " 0.25" " 0.65" "1" "49.5" "48.9" "50.4" "63.4"
## [76,] "-0.37" "-1.90" "3" "54.8" "36.5" "37.7" "49.8"
## [77,] " 0.06" " 1.19" "4" "65.4" "64.5" "67.9" "63.4"
## [78,] "-0.60" " 0.32" "4" "38.9" "56.7" "41.2" "33.6"
## [79,] " 0.71" " 0.03" "2" "46.9" "59.3" "63.0" "52.6"
## [80,] "-1.78" " 0.56" "4" "28.3" "46.3" "42.8" "44.4"
## [81,] " 0.04" " 0.03" "3" "65.4" "51.5" "61.2" "68.8"
## [82,] " 0.66" "-0.60" "1" "42.6" "56.7" "40.3" "49.8"
## [83,] "-1.09" "-0.90" "3" "57.4" "59.9" "50.5" "52.6"
## [84,] " 0.06" " 0.03" "4" "62.7" "64.5" "71.3" "55.3"
## [85,] "-0.66" "-0.08" "4" "44.7" "33.3" "33.7" "39.0"
## [86,] " 0.71" " 1.19" "3" "54.8" "48.9" "52.4" "58.0"
## [87,] " 0.75" " 1.19" "4" "68.0" "61.9" "52.2" "60.7"
## [88,] " 0.06" "-0.60" "3" "62.7" "56.0" "57.3" "60.7"
## [89,] " 0.43" "-1.09" "3" "63.3" "64.5" "52.1" "58.0"
## [90,] " 0.30" " 0.56" "4" "44.2" "33.3" "37.7" "30.9"
## [91,] " 0.26" " 0.94" "1" "54.8" "38.5" "46.8" "36.3"
## [92,] " 0.46" "-0.30" "4" "53.2" "60.6" "61.2" "56.9"
## [93,] "-0.90" "-2.54" "1" "44.2" "43.7" "40.2" "39.6"
## [94,] "-0.19" "-0.73" "3" "73.3" "60.6" "64.7" "66.1"
## [95,] " 0.45" " 0.56" "4" "52.1" "43.7" "41.9" "47.1"
## [96,] " 0.30" "-0.60" "4" "49.5" "59.3" "53.1" "60.7"
## [97,] "-1.78" " 0.56" "2" "37.3" "43.7" "45.4" "39.0"
## [98,] " 0.91" " 0.59" "4" "65.4" "67.1" "67.1" "66.1"
## [99,] "-0.60" "-0.47" "3" "46.9" "54.1" "49.0" "52.6"
## [100,] " 0.68" " 0.03" "1" "44.2" "35.9" "43.6" "47.1"
## [101,] " 0.06" "-1.07" "1" "38.9" "48.9" "45.2" "44.4"
## [102,] " 0.00" " 0.34" "2" "46.9" "59.3" "53.7" "58.0"
## [103,] " 0.46" " 0.03" "4" "60.1" "64.5" "56.3" "63.4"
## [104,] "-0.63" " 0.44" "2" "41.6" "43.7" "46.8" "36.3"
## [105,] "-0.40" "-0.26" "2" "52.1" "54.1" "55.3" "52.6"
## [106,] " 0.44" " 0.03" "4" "55.3" "51.5" "48.0" "58.0"
## [107,] " 0.91" " 0.65" "2" "54.8" "48.9" "41.9" "52.6"
## [108,] "-1.15" "-0.47" "3" "44.2" "56.7" "52.2" "44.4"
## [109,] "-0.03" " 0.63" "3" "52.1" "54.1" "48.2" "55.3"
## [110,] "-0.40" "-0.47" "4" "54.8" "59.3" "51.8" "41.7"
## [111,] " 0.02" " 0.13" "2" "73.3" "59.9" "70.5" "60.7"
## [112,] " 1.36" " 0.44" "3" "70.7" "58.0" "65.4" "63.4"
## [113,] "-0.14" " 0.13" "3" "41.6" "56.7" "51.3" "47.1"
## [114,] "-0.17" "-0.84" "3" "62.7" "56.7" "54.3" "58.0"
## [115,] " 0.93" "-1.63" "3" "38.9" "41.1" "40.3" "34.1"
## [116,] "-1.99" " 0.03" "1" "39.4" "54.1" "38.7" "47.1"
## [117,] " 0.31" " 1.19" "4" "40.5" "46.9" "40.4" "39.0"
## [118,] " 0.71" " 0.03" "4" "54.8" "61.2" "53.7" "48.8"
## [119,] " 0.52" " 0.34" "4" "65.4" "54.1" "61.4" "58.0"
## [120,] " 0.42" " 0.56" "3" "70.7" "56.0" "62.0" "67.7"
## [121,] "-1.74" " 0.61" "2" "45.8" "34.6" "37.9" "39.0"
## [122,] " 0.06" " 0.03" "4" "41.6" "46.3" "46.2" "39.0"
## [123,] " 1.13" " 0.87" "4" "57.4" "54.7" "51.7" "49.8"
## [124,] " 0.75" " 1.19" "4" "60.1" "61.9" "67.1" "49.8"
## [125,] " 1.16" " 1.19" "4" "70.7" "64.5" "72.2" "66.1"
## [126,] " 0.46" " 0.63" "3" "46.9" "52.8" "49.3" "53.1"
## [127,] "-0.45" "-0.60" "3" "52.1" "59.9" "49.5" "48.8"
## [128,] "-1.50" " 0.03" "3" "33.6" "48.9" "38.6" "42.3"
## [129,] " 0.71" " 0.34" "1" "62.7" "64.5" "57.4" "60.7"
## [130,] " 0.46" "-0.55" "2" "41.6" "54.1" "50.3" "49.8"
## [131,] " 0.71" " 1.19" "2" "54.8" "59.3" "62.5" "68.8"
## [132,] " 1.36" " 1.19" "4" "65.4" "48.9" "66.3" "58.0"
## [133,] " 0.30" " 0.03" "3" "36.8" "59.3" "40.7" "49.8"
## [134,] " 0.93" " 0.65" "4" "62.7" "59.3" "71.3" "68.8"
## [135,] " 0.55" " 0.90" "4" "62.7" "61.9" "59.6" "60.7"
## [136,] " 0.27" " 0.03" "3" "57.4" "41.1" "57.0" "55.3"
## [137,] " 0.44" "-0.16" "3" "62.7" "48.9" "48.8" "66.1"
## [138,] "-1.34" " 1.19" "4" "44.2" "41.1" "44.2" "44.4"
## [139,] " 0.06" " 0.94" "3" "45.3" "54.7" "44.3" "33.6"
## [140,] " 0.26" " 0.03" "4" "60.1" "64.5" "55.7" "63.4"
## [141,] "-1.33" " 0.65" "1" "40.5" "38.5" "47.1" "44.4"
## [142,] "-1.10" "-0.28" "3" "38.9" "41.1" "45.9" "47.1"
## [143,] " 0.96" " 0.63" "4" "65.4" "64.5" "70.3" "66.1"
## [144,] " 0.06" " 0.90" "3" "38.9" "52.8" "54.4" "49.8"
## [145,] " 0.46" " 0.03" "1" "73.3" "61.9" "73.1" "68.8"
## [146,] " 0.68" " 0.03" "3" "71.2" "63.2" "60.2" "65.5"
## [147,] " 0.68" " 0.32" "4" "36.3" "56.7" "41.9" "49.8"
## [148,] "-0.44" "-1.13" "4" "54.8" "61.9" "69.6" "60.7"
## [149,] " 0.55" "-0.60" "4" "54.8" "54.1" "51.3" "41.7"
## [150,] "-0.33" " 0.38" "3" "62.7" "59.3" "56.5" "55.3"
## [151,] " 0.50" " 0.03" "3" "68.0" "59.3" "58.8" "66.1"
## [152,] " 1.11" "-0.09" "1" "44.2" "48.9" "48.0" "49.8"
## [153,] "-0.43" " 0.03" "4" "52.1" "44.3" "53.1" "58.0"
## [154,] " 0.00" "-0.42" "3" "49.5" "51.5" "52.9" "39.0"
## [155,] " 0.71" " 0.34" "3" "68.0" "61.9" "64.5" "68.8"
## [156,] " 0.00" "-1.09" "3" "62.7" "51.5" "45.9" "47.1"
## [157,] " 0.91" "-0.28" "4" "60.1" "67.1" "56.2" "37.4"
## [158,] " 0.06" " 0.03" "4" "41.6" "54.1" "41.2" "41.7"
## [159,] "-0.59" " 0.03" "2" "33.6" "54.1" "41.0" "41.7"
## [160,] " 0.69" " 0.34" "4" "52.1" "56.7" "53.4" "60.7"
## [161,] " 0.68" " 0.03" "4" "65.4" "64.5" "51.3" "66.1"
## [162,] "-0.82" " 0.63" "4" "36.8" "36.5" "41.5" "33.1"
## [163,] "-0.66" "-0.47" "4" "52.7" "41.1" "40.2" "58.0"
## [164,] " 0.03" " 0.32" "4" "44.7" "56.7" "50.3" "30.9"
## [165,] " 0.71" "-0.29" "3" "57.4" "61.9" "55.5" "58.0"
## [166,] " 0.47" "-0.09" "2" "33.6" "33.9" "38.8" "39.6"
## [167,] "-1.10" " 0.03" "4" "38.9" "38.5" "42.8" "41.7"
## [168,] "-0.40" "-0.76" "3" "52.1" "35.9" "50.3" "47.1"
## [169,] " 0.23" " 0.44" "4" "49.5" "51.5" "45.5" "49.8"
## [170,] " 0.68" "-0.47" "2" "46.9" "61.9" "60.5" "68.8"
## [171,] " 1.36" " 0.59" "4" "62.7" "67.1" "69.6" "68.8"
## [172,] "-0.21" "-1.38" "1" "34.2" "46.3" "44.5" "39.0"
## [173,] " 0.06" " 0.56" "2" "46.9" "41.1" "45.3" "47.1"
## [174,] " 0.91" "-1.67" "2" "57.4" "54.1" "54.8" "60.7"
## [175,] "-0.40" " 0.88" "2" "54.8" "45.0" "45.9" "58.0"
## [176,] " 0.29" " 0.39" "2" "41.6" "33.3" "37.6" "28.2"
## [177,] "-0.39" " 1.19" "4" "40.5" "59.3" "41.9" "33.6"
## [178,] " 0.30" "-0.59" "4" "54.8" "46.3" "45.5" "58.0"
## [179,] " 0.31" " 0.63" "4" "38.9" "56.7" "46.3" "58.0"
## [180,] " 0.69" " 0.03" "3" "49.5" "39.8" "38.6" "49.3"
## [181,] " 0.52" "-0.28" "3" "38.9" "41.7" "33.7" "30.9"
## [182,] "-0.38" " 0.03" "1" "52.1" "43.7" "49.4" "55.3"
## [183,] " 0.10" " 1.19" "4" "52.1" "56.7" "62.3" "66.1"
## [184,] " 0.10" " 0.03" "2" "49.5" "56.7" "48.0" "47.1"
## [185,] "-0.42" "-0.60" "4" "54.8" "56.7" "64.7" "58.0"
## [186,] " 0.33" "-0.26" "4" "52.1" "67.1" "57.0" "63.4"
## [187,] " 0.06" "-0.59" "1" "44.2" "41.1" "59.8" "44.4"
## [188,] " 1.14" " 1.19" "4" "65.9" "67.1" "67.1" "60.7"
## [189,] " 0.46" " 0.03" "3" "52.1" "56.7" "53.0" "47.1"
## [190,] " 0.46" " 0.34" "3" "62.2" "56.7" "51.9" "54.7"
## [191,] " 0.68" " 0.94" "4" "62.7" "52.1" "66.1" "68.8"
## [192,] " 0.10" " 0.56" "3" "54.8" "61.9" "60.4" "55.3"
## [193,] " 0.05" " 0.15" "4" "44.2" "49.5" "40.5" "39.0"
## [194,] " 0.00" " 0.65" "4" "68.0" "64.5" "58.3" "58.5"
## [195,] "-0.60" " 0.34" "1" "49.5" "41.7" "50.3" "36.3"
## [196,] "-1.57" " 0.65" "2" "52.1" "56.7" "45.2" "55.3"
## [197,] " 1.36" " 0.94" "4" "52.1" "48.9" "51.3" "41.7"
## [198,] " 0.00" " 0.65" "3" "46.9" "33.3" "50.4" "47.1"
## [199,] "-0.93" "-0.80" "2" "52.1" "51.5" "49.5" "55.3"
## [200,] " 0.32" " 0.90" "3" "52.1" "59.3" "58.1" "47.1"
## [201,] " 0.26" " 0.03" "4" "57.4" "51.5" "55.3" "60.7"
## [202,] " 0.20" "-0.47" "3" "56.4" "49.5" "51.9" "47.1"
## [203,] " 0.23" " 0.15" "2" "49.5" "51.5" "48.0" "52.6"
## [204,] "-0.16" " 0.65" "4" "68.0" "61.9" "69.7" "71.5"
## [205,] "-0.61" " 0.36" "1" "44.2" "48.9" "46.1" "47.1"
## [206,] " 0.66" " 0.34" "3" "54.8" "56.7" "61.9" "63.4"
## [207,] " 0.46" " 0.03" "3" "54.8" "54.7" "56.9" "58.0"
## [208,] " 0.29" " 0.03" "4" "49.5" "48.9" "56.2" "47.1"
## [209,] "-0.66" "-1.07" "3" "49.5" "61.9" "60.4" "47.1"
## [210,] " 0.26" " 1.19" "4" "52.1" "52.8" "37.6" "52.6"
## [211,] " 0.44" "-0.47" "2" "62.7" "64.5" "48.0" "63.4"
## [212,] "-0.19" "-0.85" "1" "54.8" "60.6" "55.5" "58.0"
## [213,] " 0.52" " 0.65" "2" "54.3" "62.5" "56.6" "54.7"
## [214,] " 0.91" " 0.34" "4" "52.1" "54.1" "58.1" "55.8"
## [215,] " 0.06" " 0.32" "3" "49.5" "43.7" "55.5" "68.8"
## [216,] "-1.55" " 0.03" "1" "36.3" "41.1" "43.5" "33.6"
## [217,] " 0.42" "-0.47" "2" "73.3" "61.2" "57.4" "54.7"
## [218,] " 0.67" " 0.03" "3" "52.1" "56.7" "51.1" "55.3"
## [219,] " 0.06" " 0.59" "4" "52.1" "51.5" "45.4" "39.0"
## [220,] "-0.16" " 0.03" "2" "44.2" "59.3" "49.5" "55.3"
## [221,] " 0.49" " 1.19" "4" "46.9" "51.5" "50.6" "49.8"
## [222,] "-0.80" " 0.15" "2" "41.6" "41.1" "39.5" "47.1"
## [223,] "-0.47" " 0.28" "4" "65.4" "61.9" "47.1" "60.7"
## [224,] " 0.46" " 0.03" "4" "41.6" "64.5" "47.1" "53.1"
## [225,] "-0.34" " 0.59" "4" "38.9" "33.9" "35.1" "44.4"
## [226,] " 0.10" " 1.19" "4" "38.9" "35.9" "46.8" "39.0"
## [227,] " 0.52" "-0.59" "3" "41.6" "59.3" "38.6" "36.3"
## [228,] " 0.23" " 0.03" "4" "52.1" "59.3" "52.9" "60.7"
## [229,] " 0.49" " 0.94" "4" "49.5" "39.8" "44.5" "55.3"
## [230,] " 0.08" " 0.03" "4" "54.8" "67.1" "47.4" "49.8"
## [231,] " 0.48" " 0.69" "3" "62.7" "43.0" "58.9" "65.0"
## [232,] " 0.57" "-2.62" "1" "38.9" "28.1" "38.4" "44.4"
## [233,] "-0.38" " 0.37" "3" "44.2" "43.7" "46.1" "47.1"
## [234,] " 1.11" "-1.05" "2" "70.7" "67.1" "63.0" "63.4"
## [235,] "-0.04" "-0.29" "3" "41.6" "38.5" "40.2" "33.6"
## [236,] " 0.46" " 0.59" "4" "52.1" "62.5" "53.6" "56.3"
## [237,] "-0.83" " 1.19" "3" "38.9" "38.5" "44.4" "26.0"
## [238,] " 0.00" " 0.34" "4" "62.7" "59.9" "65.4" "53.6"
## [239,] " 0.28" " 0.32" "4" "52.1" "54.1" "56.5" "55.3"
## [240,] "-0.22" "-1.34" "3" "46.9" "31.3" "43.6" "36.3"
## [241,] " 0.00" "-0.18" "1" "62.7" "59.3" "55.5" "45.5"
## [242,] "-0.44" "-0.55" "2" "46.9" "38.5" "47.1" "41.7"
## [243,] " 0.26" " 0.88" "4" "65.4" "59.3" "62.3" "60.7"
## [244,] " 0.46" "-0.47" "1" "49.5" "61.9" "41.4" "55.3"
## [245,] " 1.16" " 1.19" "2" "55.3" "46.9" "49.1" "54.2"
## [246,] " 0.01" " 0.32" "1" "52.1" "41.1" "50.6" "52.6"
## [247,] " 0.93" " 0.03" "4" "65.9" "67.1" "70.5" "52.6"
## [248,] " 1.16" "-0.72" "3" "56.9" "64.5" "72.2" "54.2"
## [249,] " 0.46" " 0.65" "4" "68.0" "59.3" "57.9" "74.2"
## [250,] " 0.89" " 0.32" "3" "46.9" "44.3" "48.7" "53.1"
## [251,] " 0.91" " 0.03" "4" "57.4" "59.3" "55.5" "55.3"
## [252,] " 0.73" " 0.03" "4" "62.7" "62.5" "68.9" "60.7"
## [253,] " 0.46" " 0.03" "1" "49.5" "46.3" "46.2" "41.7"
## [254,] "-0.22" " 0.32" "2" "35.2" "38.5" "39.9" "34.7"
## [255,] " 0.00" " 0.15" "2" "73.3" "64.5" "75.5" "60.7"
## [256,] " 0.73" "-1.42" "2" "38.9" "41.1" "41.0" "36.3"
## [257,] " 0.25" " 0.34" "4" "49.5" "61.9" "42.9" "41.7"
## [258,] " 0.31" " 0.03" "4" "54.8" "48.9" "61.2" "60.7"
## [259,] " 0.93" " 0.34" "4" "46.9" "54.1" "54.6" "55.3"
## [260,] " 0.51" "-1.29" "2" "39.4" "41.7" "46.9" "35.2"
## [261,] "-0.49" " 0.03" "1" "46.9" "37.2" "42.8" "41.7"
## [262,] "-0.59" " 1.19" "3" "36.3" "36.5" "35.1" "33.6"
## [263,] "-0.59" " 0.28" "4" "46.9" "51.5" "48.5" "49.8"
## [264,] " 0.29" " 0.88" "3" "33.6" "43.7" "40.2" "39.0"
## [265,] "-0.23" " 0.69" "3" "46.9" "61.9" "48.0" "39.0"
## [266,] " 0.45" " 0.03" "3" "60.1" "61.9" "51.9" "53.1"
## [267,] " 0.28" "-0.06" "2" "44.2" "30.7" "35.3" "47.1"
## [268,] "-0.49" " 0.03" "2" "65.4" "59.3" "56.8" "66.1"
## [269,] " 0.44" " 0.65" "4" "52.1" "54.1" "54.6" "41.7"
## [270,] " 0.53" " 0.81" "3" "54.8" "59.3" "61.4" "47.1"
## [271,] " 0.46" " 0.03" "2" "44.7" "38.5" "45.9" "44.4"
## [272,] " 0.48" "-0.47" "2" "52.1" "61.9" "55.5" "60.7"
## [273,] "-0.39" "-0.28" "3" "44.2" "51.5" "51.1" "63.4"
## [274,] "-0.43" " 0.13" "2" "38.9" "44.3" "51.6" "43.9"
## [275,] " 0.75" "-0.72" "3" "42.6" "46.3" "41.2" "28.2"
## [276,] " 0.91" " 0.34" "4" "68.0" "59.3" "53.0" "63.4"
## [277,] "-1.76" "-0.93" "3" "36.3" "38.5" "36.4" "36.3"
## [278,] " 0.46" "-0.86" "4" "68.0" "59.3" "60.5" "55.3"
## [279,] " 0.66" " 0.03" "4" "52.1" "56.7" "41.9" "52.6"
## [280,] "-1.58" "-0.26" "2" "62.7" "48.9" "35.3" "66.1"
## [281,] "-0.19" " 0.03" "2" "54.8" "51.5" "42.8" "60.7"
## [282,] " 0.46" " 1.19" "3" "65.9" "64.5" "61.8" "60.7"
## [283,] "-0.41" " 0.28" "2" "46.9" "35.9" "42.6" "58.0"
## [284,] " 0.00" " 0.65" "4" "52.1" "61.9" "62.1" "58.0"
## [285,] " 0.22" " 0.03" "2" "52.1" "54.7" "49.5" "53.6"
## [286,] " 0.06" "-0.29" "3" "57.4" "64.5" "51.1" "63.4"
## [287,] " 0.46" " 0.34" "4" "57.4" "54.1" "59.6" "60.7"
## [288,] " 0.00" " 0.03" "2" "52.1" "59.3" "50.4" "60.7"
## [289,] "-0.65" " 0.03" "2" "41.6" "39.1" "42.3" "41.7"
## [290,] "-0.60" "-1.18" "3" "54.8" "59.3" "68.0" "49.3"
## [291,] "-0.40" "-1.80" "2" "50.0" "45.0" "42.9" "41.7"
## [292,] " 0.71" " 0.34" "3" "52.1" "41.1" "51.3" "60.7"
## [293,] "-0.45" "-0.26" "1" "38.9" "35.9" "44.2" "55.3"
## [294,] " 0.52" " 0.03" "4" "46.9" "51.5" "53.7" "41.7"
## [295,] "-1.15" "-1.38" "4" "49.5" "56.7" "46.1" "53.6"
## [296,] " 1.36" " 0.63" "4" "68.0" "64.5" "58.9" "63.4"
## [297,] " 0.68" " 1.19" "2" "41.6" "41.1" "57.2" "71.5"
## [298,] " 0.49" " 0.44" "3" "53.7" "43.7" "51.1" "52.0"
## [299,] "-0.11" " 0.25" "4" "44.2" "44.3" "45.6" "39.0"
## [300,] " 0.71" " 0.03" "4" "46.9" "61.9" "46.2" "60.7"
## [301,] "-1.58" " 0.56" "3" "49.5" "33.3" "48.7" "44.4"
## [302,] "-0.21" "-1.13" "1" "38.9" "41.1" "43.6" "55.3"
## [303,] " 0.48" "-0.82" "2" "62.7" "64.5" "58.2" "61.8"
## [304,] " 0.26" " 0.03" "4" "62.7" "59.3" "53.2" "63.4"
## [305,] " 0.75" "-0.49" "2" "41.6" "41.1" "41.2" "44.4"
## [306,] " 0.06" " 0.03" "4" "54.8" "41.1" "42.8" "47.1"
## [307,] "-0.84" "-0.18" "1" "33.6" "46.3" "38.4" "36.3"
## [308,] " 1.36" " 0.65" "4" "46.9" "51.5" "57.0" "49.8"
## [309,] " 0.23" " 0.94" "4" "57.4" "67.1" "57.9" "60.7"
## [310,] "-0.44" "-0.78" "4" "58.0" "48.9" "43.8" "44.4"
## [311,] "-0.38" " 0.03" "4" "60.1" "54.1" "58.0" "58.0"
## [312,] "-0.18" " 0.03" "3" "44.2" "48.9" "56.3" "41.7"
## [313,] " 1.13" " 0.63" "3" "60.1" "64.5" "57.9" "60.7"
## [314,] "-0.40" " 0.03" "1" "44.2" "33.3" "47.8" "44.4"
## [315,] " 1.16" " 1.19" "4" "60.1" "59.3" "55.3" "58.0"
## [316,] " 0.06" " 0.28" "4" "65.9" "63.2" "50.3" "60.1"
## [317,] " 0.46" " 0.03" "3" "54.8" "64.5" "58.6" "55.3"
## [318,] " 0.28" " 1.19" "4" "56.4" "59.3" "38.5" "55.3"
## [319,] " 1.13" " 0.03" "2" "43.2" "54.1" "54.8" "55.3"
## [320,] " 0.30" "-1.09" "1" "51.1" "56.7" "55.3" "47.7"
## [321,] " 0.27" "-2.52" "2" "38.9" "32.0" "35.3" "46.6"
## [322,] "-1.05" "-0.34" "3" "44.2" "38.5" "41.4" "42.3"
## [323,] " 1.13" " 0.03" "3" "62.7" "64.5" "64.5" "52.6"
## [324,] " 0.06" " 0.03" "1" "46.9" "51.5" "57.2" "52.6"
## [325,] "-0.40" "-1.34" "4" "57.4" "61.9" "55.5" "49.8"
## [326,] " 0.91" " 0.03" "4" "60.6" "59.3" "49.4" "44.4"
## [327,] "-0.17" " 0.03" "3" "44.2" "59.3" "45.7" "55.3"
## [328,] " 0.68" " 0.25" "4" "44.2" "51.5" "36.9" "49.8"
## [329,] " 0.33" " 0.28" "3" "52.1" "46.3" "48.1" "55.3"
## [330,] "-0.16" "-0.53" "3" "62.7" "61.9" "67.0" "66.1"
## [331,] " 0.68" " 0.87" "4" "54.8" "56.7" "56.9" "58.0"
## [332,] "-1.09" " 0.03" "2" "46.9" "43.7" "41.9" "41.7"
## [333,] " 0.28" " 0.03" "4" "73.3" "64.5" "68.0" "63.4"
## [334,] " 0.47" " 0.34" "3" "33.6" "43.0" "41.0" "49.8"
## [335,] " 0.46" " 0.34" "4" "52.1" "46.3" "50.4" "52.6"
## [336,] "-0.40" "-0.89" "2" "41.6" "56.7" "45.0" "50.4"
## [337,] "-1.28" " 0.34" "2" "31.0" "35.9" "46.1" "39.0"
## [338,] " 0.46" "-0.47" "3" "68.0" "64.5" "70.6" "68.8"
## [339,] "-0.38" " 0.56" "4" "33.6" "36.5" "45.7" "39.0"
## [340,] "-0.03" " 0.87" "4" "53.2" "46.3" "43.0" "41.7"
## [341,] " 0.49" " 0.94" "4" "54.8" "56.7" "63.8" "60.7"
## [342,] " 0.11" "-1.42" "4" "60.1" "59.3" "62.1" "60.7"
## [343,] "-2.23" " 1.19" "2" "44.2" "38.5" "41.2" "49.8"
## [344,] " 1.13" " 0.03" "2" "45.8" "41.7" "43.1" "53.6"
## [345,] " 0.47" " 0.34" "3" "44.2" "48.9" "48.0" "39.0"
## [346,] " 0.22" "-0.28" "1" "36.3" "48.9" "39.5" "41.7"
## [347,] " 0.68" "-0.26" "4" "62.7" "61.9" "56.2" "47.1"
## [348,] "-0.71" "-2.29" "2" "65.4" "64.5" "63.7" "58.0"
## [349,] " 0.20" "-0.47" "2" "38.9" "38.5" "53.0" "52.6"
## [350,] " 0.48" " 0.32" "4" "49.5" "54.1" "38.7" "49.8"
## [351,] " 0.91" "-0.47" "3" "57.4" "59.3" "57.2" "59.1"
## [352,] " 0.06" " 0.03" "3" "54.8" "54.1" "47.1" "49.8"
## [353,] " 0.02" " 0.03" "2" "41.6" "35.9" "37.7" "41.7"
## [354,] " 0.46" " 0.94" "2" "36.3" "30.7" "44.8" "39.6"
## [355,] "-0.89" "-1.23" "4" "44.2" "48.9" "43.8" "35.2"
## [356,] " 0.46" " 0.32" "3" "49.5" "54.1" "47.8" "60.7"
## [357,] " 0.26" " 0.03" "4" "49.5" "35.9" "47.8" "52.6"
## [358,] "-0.27" " 0.88" "4" "52.1" "64.5" "60.6" "60.7"
## [359,] "-0.93" " 0.65" "4" "47.4" "45.6" "49.4" "42.3"
## [360,] " 0.46" "-1.11" "2" "65.4" "67.1" "67.0" "66.1"
## [361,] " 0.23" "-0.76" "4" "65.4" "64.5" "59.6" "66.1"
## [362,] " 0.96" " 0.87" "4" "57.4" "64.5" "48.0" "47.1"
## [363,] "-0.36" "-1.67" "4" "44.2" "43.7" "56.4" "58.0"
## [364,] " 0.06" " 0.32" "1" "60.1" "54.1" "49.9" "49.8"
## [365,] " 0.22" "-0.76" "4" "31.0" "41.1" "49.2" "33.6"
## [366,] "-1.05" " 0.28" "3" "62.7" "64.5" "66.3" "58.0"
## [367,] "-0.19" "-1.67" "1" "42.1" "39.8" "42.2" "38.5"
## [368,] "-1.33" " 0.03" "3" "41.6" "30.7" "56.9" "47.1"
## [369,] " 0.21" "-0.28" "3" "46.9" "41.1" "53.6" "41.7"
## [370,] "-0.44" " 1.19" "4" "60.1" "54.1" "59.6" "55.3"
## [371,] "-0.23" " 0.44" "2" "33.6" "35.2" "40.9" "29.3"
## [372,] " 0.32" "-0.88" "2" "57.4" "54.1" "42.2" "44.4"
## [373,] "-0.93" "-0.22" "3" "65.4" "64.5" "65.5" "60.7"
## [374,] " 0.48" " 0.03" "4" "60.1" "46.3" "51.4" "52.6"
## [375,] "-1.81" "-0.64" "4" "41.6" "54.1" "46.9" "47.1"
## [376,] " 0.71" "-0.05" "3" "52.7" "47.6" "52.2" "45.0"
## [377,] "-0.54" "-0.59" "4" "38.9" "35.9" "38.4" "33.6"
## [378,] " 0.10" " 0.03" "1" "62.7" "48.9" "42.1" "66.1"
## [379,] " 0.26" " 0.94" "4" "60.1" "59.3" "51.8" "53.1"
## [380,] " 1.11" " 0.34" "4" "73.3" "67.1" "62.3" "58.0"
## [381,] "-0.84" "-2.60" "2" "44.2" "64.5" "53.0" "55.3"
## [382,] " 0.02" "-2.03" "3" "38.9" "59.3" "39.5" "39.0"
## [383,] " 0.66" "-0.47" "3" "54.8" "59.3" "50.6" "60.7"
## [384,] " 0.27" " 0.03" "3" "46.9" "56.7" "63.8" "55.3"
## [385,] " 0.45" " 0.03" "3" "41.6" "25.5" "44.3" "36.3"
## [386,] " 0.55" " 0.34" "4" "62.7" "61.9" "72.9" "63.4"
## [387,] " 0.46" " 0.65" "4" "68.0" "54.1" "63.0" "63.4"
## [388,] "-0.37" "-0.47" "3" "41.6" "48.9" "42.6" "49.8"
## [389,] "-0.40" "-1.09" "3" "62.7" "43.7" "47.1" "52.6"
## [390,] " 0.55" " 0.63" "4" "54.8" "59.3" "62.3" "58.0"
## [391,] " 0.87" " 0.87" "3" "57.4" "61.9" "62.8" "55.3"
## [392,] " 0.68" " 0.34" "4" "57.4" "56.7" "62.8" "58.0"
## [393,] " 0.71" " 0.03" "3" "41.6" "59.3" "41.9" "44.4"
## [394,] "-1.33" " 0.03" "4" "41.6" "48.9" "33.7" "44.4"
## [395,] " 0.24" "-0.35" "4" "30.5" "35.9" "36.9" "33.6"
## [396,] " 0.43" "-0.47" "4" "42.6" "50.8" "44.6" "39.0"
## [397,] " 0.65" "-0.30" "4" "60.1" "61.9" "74.6" "58.0"
## [398,] " 0.75" "-0.26" "3" "41.6" "35.9" "42.0" "31.4"
## [399,] " 0.00" "-0.82" "3" "54.8" "67.1" "63.8" "49.8"
## [400,] " 0.06" " 0.28" "4" "44.2" "54.1" "47.1" "58.0"
## [401,] " 0.02" "-1.05" "2" "60.1" "64.5" "61.2" "60.7"
## [402,] " 0.32" " 0.28" "3" "62.7" "63.2" "68.9" "74.2"
## [403,] "-0.88" "-0.29" "2" "44.2" "54.1" "41.0" "49.8"
## [404,] " 0.43" "-0.05" "3" "54.8" "45.6" "51.8" "60.7"
## [405,] "-1.28" "-1.05" "2" "57.4" "43.7" "57.8" "60.7"
## [406,] "-0.35" "-0.89" "1" "49.5" "67.1" "66.4" "66.1"
## [407,] "-0.16" "-2.58" "4" "52.1" "59.3" "58.1" "52.6"
## [408,] "-0.16" "-0.47" "3" "36.3" "61.9" "52.2" "39.0"
## [409,] " 0.70" "-0.16" "4" "49.5" "61.9" "50.9" "46.1"
## [410,] " 0.55" "-0.24" "3" "46.9" "56.7" "54.6" "47.1"
## [411,] " 1.12" "-0.74" "3" "31.0" "41.1" "36.0" "36.9"
## [412,] " 0.00" " 0.65" "4" "54.8" "64.5" "51.1" "49.8"
## [413,] "-0.65" "-1.17" "1" "36.3" "30.7" "45.9" "36.3"
## [414,] " 0.66" " 0.03" "2" "65.4" "59.3" "67.1" "66.1"
## [415,] "-0.47" "-0.47" "1" "52.1" "57.3" "63.7" "58.0"
## [416,] " 0.52" " 0.34" "4" "52.1" "54.1" "51.4" "49.8"
## [417,] "-0.38" " 0.03" "4" "45.3" "38.5" "43.5" "40.1"
## [418,] "-0.83" " 0.13" "4" "41.6" "35.9" "36.7" "36.3"
## [419,] " 0.67" "-0.59" "3" "54.8" "64.5" "55.3" "49.8"
## [420,] " 0.06" " 0.03" "4" "68.0" "61.9" "55.5" "49.8"
## [421,] " 0.27" " 0.34" "4" "49.5" "48.9" "48.0" "41.7"
## [422,] "-0.16" "-0.45" "3" "49.5" "54.1" "61.9" "52.6"
## [423,] "-0.21" " 0.03" "4" "41.6" "43.7" "51.1" "49.8"
## [424,] "-0.68" " 0.32" "3" "36.3" "44.3" "34.4" "47.1"
## [425,] " 1.16" " 0.59" "4" "36.3" "45.0" "44.3" "39.0"
## [426,] " 1.13" " 0.87" "4" "76.0" "62.5" "60.0" "67.2"
## [427,] "-0.89" "-0.18" "3" "44.2" "43.7" "44.5" "41.7"
## [428,] " 0.00" " 0.03" "3" "61.1" "61.2" "62.2" "57.4"
## [429,] " 0.22" "-0.76" "3" "45.8" "52.1" "54.9" "44.4"
## [430,] " 0.46" "-0.22" "4" "70.7" "64.5" "68.7" "58.0"
## [431,] " 0.48" " 0.63" "4" "54.8" "48.2" "58.7" "51.5"
## [432,] " 0.06" " 0.03" "4" "57.4" "54.1" "58.9" "47.1"
## [433,] "-0.89" "-1.29" "2" "60.1" "64.5" "41.0" "52.6"
## [434,] "-0.69" "-1.13" "3" "46.9" "54.1" "45.4" "44.4"
## [435,] "-0.14" "-0.86" "1" "38.9" "46.3" "50.1" "36.3"
## [436,] " 0.07" " 0.34" "4" "65.4" "59.3" "51.2" "60.7"
## [437,] " 0.68" " 0.03" "4" "44.2" "48.9" "46.2" "47.7"
## [438,] " 0.20" "-2.54" "1" "54.8" "61.9" "63.6" "63.4"
## [439,] "-0.04" " 0.65" "4" "44.2" "44.3" "61.1" "49.8"
## [440,] "-0.44" "-0.84" "2" "46.9" "46.3" "45.9" "37.9"
## [441,] " 0.73" " 0.03" "4" "65.4" "59.3" "69.6" "63.4"
## [442,] " 0.06" "-0.28" "4" "65.9" "64.5" "54.6" "60.7"
## [443,] " 0.20" " 0.56" "3" "60.1" "52.8" "51.5" "58.0"
## [444,] " 0.68" " 0.03" "4" "65.4" "67.1" "64.5" "60.7"
## [445,] " 0.68" "-0.60" "2" "52.1" "64.5" "59.7" "55.8"
## [446,] " 0.02" " 0.28" "3" "44.2" "51.5" "52.9" "52.6"
## [447,] " 0.68" "-0.53" "4" "48.4" "48.2" "60.5" "60.1"
## [448,] "-0.18" "-1.09" "3" "62.7" "56.7" "58.8" "71.5"
## [449,] " 1.36" " 0.65" "3" "59.0" "56.7" "52.9" "48.8"
## [450,] "-1.33" "-0.26" "3" "38.4" "37.2" "33.4" "30.9"
## [451,] " 0.06" " 0.44" "3" "52.1" "61.9" "66.3" "47.1"
## [452,] " 0.91" " 0.15" "1" "54.8" "56.7" "55.5" "58.0"
## [453,] " 0.06" "-0.29" "3" "49.5" "56.7" "68.1" "49.8"
## [454,] "-0.42" " 0.03" "2" "46.9" "35.9" "41.9" "36.3"
## [455,] " 0.02" "-0.32" "4" "46.9" "46.3" "52.0" "47.7"
## [456,] " 0.32" " 1.19" "4" "57.4" "61.9" "64.0" "55.8"
## [457,] " 0.22" " 0.03" "3" "60.1" "41.7" "52.0" "52.6"
## [458,] "-0.67" "-0.47" "1" "44.2" "48.9" "43.1" "55.3"
## [459,] " 0.48" " 0.03" "4" "65.4" "61.9" "60.4" "63.4"
## [460,] " 0.89" " 0.34" "2" "57.4" "61.9" "72.2" "60.7"
## [461,] "-1.33" " 0.03" "3" "57.4" "56.7" "59.8" "58.0"
## [462,] "-1.06" " 0.03" "3" "46.9" "56.7" "48.0" "44.4"
## [463,] " 0.89" " 0.15" "3" "68.0" "59.3" "56.6" "50.4"
## [464,] "-0.19" "-2.58" "1" "41.6" "54.1" "41.0" "41.7"
## [465,] "-1.33" " 0.34" "1" "44.2" "41.1" "36.9" "41.7"
## [466,] " 0.51" " 0.03" "2" "60.1" "51.5" "53.9" "63.4"
## [467,] " 1.36" " 0.63" "4" "62.7" "61.9" "68.6" "68.8"
## [468,] "-0.27" " 0.94" "1" "44.7" "54.1" "50.4" "45.0"
## [469,] "-0.22" "-2.27" "1" "57.4" "51.5" "40.3" "60.7"
## [470,] " 0.69" " 0.05" "3" "60.1" "59.3" "57.9" "40.6"
## [471,] " 1.13" " 0.56" "4" "60.1" "48.9" "50.4" "55.3"
## [472,] "-0.44" "-0.47" "2" "38.9" "46.3" "41.9" "41.7"
## [473,] " 1.36" "-1.18" "4" "62.7" "61.9" "63.1" "60.7"
## [474,] " 0.44" " 0.65" "4" "46.9" "46.3" "48.5" "32.5"
## [475,] " 0.06" " 0.65" "1" "65.4" "59.3" "56.8" "46.1"
## [476,] "-0.24" "-0.16" "2" "62.7" "61.9" "52.4" "47.1"
## [477,] " 0.27" " 0.03" "4" "60.1" "51.5" "48.1" "52.6"
## [478,] " 0.46" " 0.03" "3" "62.7" "54.1" "54.8" "55.3"
## [479,] "-1.14" "-0.06" "4" "36.3" "46.3" "36.9" "41.7"
## [480,] " 0.46" "-0.28" "4" "57.4" "56.7" "56.4" "63.4"
## [481,] " 0.22" " 1.19" "2" "38.9" "33.3" "37.7" "47.1"
## [482,] " 1.16" " 0.63" "3" "58.5" "54.1" "46.3" "48.2"
## [483,] " 1.16" " 0.03" "4" "70.7" "67.1" "61.2" "60.7"
## [484,] " 1.36" "-0.80" "3" "52.1" "59.3" "59.6" "59.1"
## [485,] " 0.08" "-0.24" "2" "52.1" "59.3" "48.0" "55.3"
## [486,] "-0.40" " 0.65" "1" "52.1" "46.3" "51.3" "58.0"
## [487,] " 0.31" " 1.19" "4" "46.9" "44.3" "41.9" "39.0"
## [488,] "-0.61" " 0.34" "3" "68.0" "54.1" "74.6" "66.1"
## [489,] "-1.33" " 0.03" "2" "38.9" "30.7" "40.4" "39.0"
## [490,] " 0.71" " 0.28" "4" "60.1" "64.5" "63.0" "68.8"
## [491,] " 0.28" " 0.94" "4" "54.8" "56.0" "43.5" "49.8"
## [492,] " 0.91" "-0.47" "4" "46.9" "64.5" "55.8" "61.8"
## [493,] "-1.10" " 0.03" "1" "46.9" "40.4" "42.8" "45.0"
## [494,] "-0.40" "-0.22" "3" "52.1" "38.5" "36.1" "44.4"
## [495,] " 0.06" " 0.03" "4" "46.9" "41.1" "46.2" "39.6"
## [496,] " 1.16" " 0.32" "4" "60.1" "54.1" "48.6" "58.0"
## [497,] " 0.05" "-0.88" "3" "62.7" "63.2" "74.6" "71.5"
## [498,] " 1.36" " 0.03" "3" "46.9" "59.3" "56.3" "66.1"
## [499,] " 0.26" " 0.87" "4" "49.5" "61.9" "61.4" "63.4"
## [500,] " 0.06" " 0.03" "1" "44.2" "33.3" "53.9" "58.0"
## [501,] " 0.52" " 0.59" "4" "68.0" "59.9" "64.0" "68.8"
## [502,] "-0.43" "-1.22" "1" "41.6" "46.9" "55.5" "44.4"
## [503,] " 0.89" " 0.88" "4" "52.1" "59.3" "64.0" "61.2"
## [504,] " 0.06" " 0.03" "4" "44.2" "54.1" "51.4" "39.0"
## [505,] " 0.02" " 0.03" "3" "65.4" "56.7" "62.3" "58.5"
## [506,] "-0.60" " 0.28" "4" "33.6" "41.1" "42.9" "41.7"
## [507,] " 0.03" " 0.28" "2" "52.1" "54.1" "52.0" "55.3"
## [508,] " 1.14" " 0.63" "4" "55.3" "53.4" "44.5" "53.6"
## [509,] " 0.00" " 0.34" "2" "46.9" "56.7" "57.2" "58.0"
## [510,] "-0.44" " 0.03" "2" "41.6" "38.5" "38.7" "56.3"
## [511,] " 0.71" "-0.55" "4" "60.1" "59.3" "64.3" "60.7"
## [512,] " 0.25" "-0.47" "2" "76.0" "67.1" "72.2" "66.1"
## [513,] " 0.06" " 0.32" "4" "62.7" "64.5" "70.5" "68.8"
## [514,] " 0.06" " 0.56" "4" "49.5" "56.7" "46.2" "55.3"
## [515,] "-0.18" " 0.15" "1" "38.9" "28.1" "36.0" "36.3"
## [516,] "-0.83" " 0.65" "2" "44.2" "59.3" "55.5" "36.3"
## [517,] "-0.84" " 0.34" "1" "41.6" "46.3" "59.6" "63.4"
## [518,] " 0.06" " 0.03" "2" "68.0" "59.3" "71.3" "66.1"
## [519,] " 0.70" "-0.16" "2" "68.0" "59.3" "55.7" "63.4"
## [520,] " 0.22" " 1.19" "4" "57.4" "54.1" "48.8" "55.3"
## [521,] " 0.32" " 0.90" "3" "33.6" "51.5" "41.9" "49.8"
## [522,] "-0.60" " 0.56" "2" "36.3" "43.7" "44.3" "47.1"
## [523,] " 1.16" " 0.09" "4" "46.9" "52.1" "42.5" "47.7"
## [524,] " 0.45" "-0.91" "2" "41.6" "41.1" "42.8" "49.8"
## [525,] "-0.20" " 1.19" "4" "41.6" "46.9" "41.7" "47.1"
## [526,] " 0.22" "-0.57" "4" "61.1" "67.1" "64.5" "50.4"
## [527,] " 0.03" " 0.03" "4" "41.6" "43.0" "44.5" "39.0"
## [528,] " 0.52" " 0.34" "4" "46.9" "64.5" "59.6" "49.8"
## [529,] " 0.32" " 0.65" "4" "54.8" "61.9" "58.1" "58.0"
## [530,] " 0.68" " 0.03" "2" "73.3" "64.5" "64.0" "60.7"
## [531,] " 0.96" " 0.94" "3" "46.9" "50.2" "46.5" "47.1"
## [532,] "-0.60" " 0.90" "3" "46.9" "61.9" "45.4" "33.6"
## [533,] " 0.71" " 0.03" "4" "68.0" "61.9" "72.2" "66.1"
## [534,] " 0.75" " 0.34" "4" "47.4" "41.7" "51.6" "38.5"
## [535,] " 0.46" "-0.18" "4" "52.1" "46.3" "48.3" "58.0"
## [536,] " 0.46" " 0.34" "3" "62.7" "61.9" "52.9" "44.4"
## [537,] "-0.45" "-0.64" "2" "65.4" "56.7" "55.2" "52.6"
## [538,] "-0.19" " 1.19" "4" "60.1" "61.9" "54.6" "60.7"
## [539,] "-1.33" " 0.03" "4" "36.3" "28.1" "53.7" "33.6"
## [540,] "-0.38" " 0.03" "4" "38.9" "41.7" "37.7" "44.4"
## [541,] " 0.45" " 0.03" "4" "49.5" "46.3" "52.9" "52.6"
## [542,] " 0.06" " 0.03" "1" "38.9" "41.1" "32.7" "41.7"
## [543,] "-0.18" " 0.28" "2" "44.2" "51.5" "46.1" "47.1"
## [544,] " 0.66" " 0.34" "2" "33.6" "38.5" "42.9" "47.1"
## [545,] " 0.71" " 0.03" "2" "46.9" "46.9" "41.0" "41.7"
## [546,] " 0.06" " 0.03" "3" "46.9" "60.6" "51.3" "63.4"
## [547,] " 0.23" "-0.60" "4" "57.4" "67.1" "68.0" "59.1"
## [548,] "-0.24" " 0.03" "3" "45.8" "54.1" "49.5" "47.1"
## [549,] " 1.11" " 0.90" "2" "55.3" "50.2" "41.7" "58.5"
## [550,] " 0.68" " 0.88" "4" "46.9" "61.9" "53.0" "52.6"
## [551,] " 0.66" " 0.34" "4" "56.9" "49.5" "49.5" "50.9"
## [552,] " 1.16" " 0.65" "4" "65.4" "67.1" "63.1" "55.3"
## [553,] "-0.44" "-0.47" "4" "54.8" "65.1" "66.1" "49.8"
## [554,] "-0.14" "-0.60" "2" "44.2" "48.9" "54.3" "49.8"
## [555,] " 0.46" " 0.56" "3" "44.2" "54.1" "52.2" "44.4"
## [556,] "-0.68" "-0.26" "1" "60.1" "51.5" "56.8" "60.7"
## [557,] " 0.46" " 0.34" "4" "76.0" "52.1" "64.1" "63.9"
## [558,] " 0.46" "-0.47" "3" "50.0" "41.1" "45.1" "44.4"
## [559,] "-0.41" " 0.03" "3" "49.5" "51.5" "45.1" "58.0"
## [560,] "-0.60" "-0.26" "3" "49.5" "41.1" "38.4" "52.6"
## [561,] "-0.51" "-1.98" "1" "60.1" "57.3" "58.8" "58.0"
## [562,] " 1.13" "-0.60" "3" "73.3" "67.1" "71.3" "63.4"
## [563,] " 1.13" " 0.88" "2" "49.5" "51.5" "52.7" "55.3"
## [564,] " 0.89" " 0.65" "3" "68.0" "48.9" "56.3" "63.4"
## [565,] " 0.06" " 0.03" "4" "43.2" "47.6" "45.6" "47.1"
## [566,] "-0.60" "-1.43" "4" "38.9" "51.5" "38.6" "39.0"
## [567,] " 0.46" "-1.17" "1" "68.0" "62.5" "65.4" "60.7"
## [568,] " 0.06" " 0.03" "3" "46.9" "56.7" "56.1" "60.7"
## [569,] " 0.26" "-0.26" "1" "49.5" "56.7" "43.8" "47.1"
## [570,] " 0.05" "-0.60" "1" "57.4" "51.5" "40.6" "46.6"
## [571,] " 0.47" " 0.03" "3" "62.7" "59.3" "48.8" "55.3"
## [572,] "-0.03" "-0.16" "1" "33.6" "43.7" "42.6" "47.1"
## [573,] " 0.20" "-0.55" "3" "60.1" "61.9" "48.7" "49.8"
## [574,] "-0.14" " 0.34" "2" "49.5" "46.3" "44.8" "58.0"
## [575,] "-0.41" " 0.34" "3" "52.1" "38.5" "57.9" "52.6"
## [576,] "-0.44" "-0.22" "1" "44.2" "33.3" "40.9" "47.1"
## [577,] " 0.28" " 0.03" "3" "62.7" "67.1" "65.6" "61.8"
## [578,] "-1.30" " 0.03" "4" "54.8" "41.1" "40.2" "44.4"
## [579,] "-0.16" " 0.03" "4" "57.4" "57.3" "64.0" "58.0"
## [580,] " 0.46" "-1.17" "2" "46.9" "30.7" "40.3" "58.0"
## [581,] " 0.71" "-0.06" "3" "44.2" "37.8" "49.2" "39.0"
## [582,] " 0.20" "-0.47" "2" "38.9" "54.1" "53.7" "52.6"
## [583,] "-0.84" "-1.09" "2" "41.6" "41.7" "45.1" "36.3"
## [584,] " 0.22" " 0.32" "4" "49.5" "56.7" "62.8" "49.8"
## [585,] " 0.47" " 0.01" "4" "68.0" "64.5" "61.5" "55.3"
## [586,] "-0.82" "-0.76" "1" "57.4" "43.7" "59.6" "52.6"
## [587,] "-2.23" " 1.19" "4" "36.3" "38.5" "39.3" "39.0"
## [588,] " 0.09" "-1.42" "2" "52.1" "61.9" "59.6" "58.0"
## [589,] " 0.10" " 0.32" "3" "40.0" "40.4" "42.4" "48.8"
## [590,] " 0.08" " 0.94" "3" "73.3" "61.9" "60.5" "71.5"
## [591,] " 0.06" " 0.03" "4" "54.8" "59.3" "46.2" "47.1"
## [592,] " 0.03" " 0.56" "3" "65.4" "67.1" "63.0" "66.1"
## [593,] "-0.61" " 0.03" "2" "44.2" "54.1" "40.3" "52.6"
## [594,] " 0.02" " 0.03" "2" "62.7" "54.1" "64.7" "58.0"
## [595,] " 0.46" " 0.03" "4" "52.1" "56.7" "62.8" "47.1"
## [596,] " 0.94" "-0.30" "4" "60.1" "67.1" "52.4" "55.3"
## [597,] " 0.23" " 0.03" "4" "65.4" "56.7" "65.4" "58.0"
## [598,] " 0.46" " 0.03" "4" "65.4" "51.5" "61.4" "60.7"
## [599,] " 0.51" " 0.03" "4" "54.8" "54.1" "66.4" "41.7"
## [600,] " 0.25" " 0.03" "3" "49.5" "51.5" "55.5" "44.4"
## sex
## [1,] "female"
## [2,] "female"
## [3,] "male"
## [4,] "male"
## [5,] "female"
## [6,] "female"
## [7,] "male"
## [8,] "male"
## [9,] "female"
## [10,] "female"
## [11,] "male"
## [12,] "male"
## [13,] "female"
## [14,] "female"
## [15,] "female"
## [16,] "male"
## [17,] "female"
## [18,] "female"
## [19,] "female"
## [20,] "female"
## [21,] "female"
## [22,] "female"
## [23,] "male"
## [24,] "male"
## [25,] "female"
## [26,] "female"
## [27,] "male"
## [28,] "female"
## [29,] "male"
## [30,] "female"
## [31,] "female"
## [32,] "female"
## [33,] "female"
## [34,] "female"
## [35,] "female"
## [36,] "female"
## [37,] "female"
## [38,] "female"
## [39,] "male"
## [40,] "female"
## [41,] "male"
## [42,] "female"
## [43,] "female"
## [44,] "male"
## [45,] "male"
## [46,] "male"
## [47,] "male"
## [48,] "male"
## [49,] "male"
## [50,] "male"
## [51,] "female"
## [52,] "male"
## [53,] "male"
## [54,] "female"
## [55,] "female"
## [56,] "female"
## [57,] "male"
## [58,] "male"
## [59,] "male"
## [60,] "female"
## [61,] "male"
## [62,] "male"
## [63,] "female"
## [64,] "male"
## [65,] "male"
## [66,] "male"
## [67,] "male"
## [68,] "female"
## [69,] "male"
## [70,] "male"
## [71,] "female"
## [72,] "female"
## [73,] "female"
## [74,] "male"
## [75,] "male"
## [76,] "male"
## [77,] "male"
## [78,] "female"
## [79,] "male"
## [80,] "female"
## [81,] "male"
## [82,] "female"
## [83,] "female"
## [84,] "male"
## [85,] "male"
## [86,] "male"
## [87,] "male"
## [88,] "female"
## [89,] "female"
## [90,] "male"
## [91,] "female"
## [92,] "male"
## [93,] "female"
## [94,] "male"
## [95,] "female"
## [96,] "female"
## [97,] "male"
## [98,] "male"
## [99,] "female"
## [100,] "female"
## [101,] "female"
## [102,] "female"
## [103,] "female"
## [104,] "female"
## [105,] "female"
## [106,] "female"
## [107,] "female"
## [108,] "female"
## [109,] "male"
## [110,] "female"
## [111,] "male"
## [112,] "male"
## [113,] "female"
## [114,] "male"
## [115,] "male"
## [116,] "female"
## [117,] "female"
## [118,] "male"
## [119,] "female"
## [120,] "male"
## [121,] "male"
## [122,] "female"
## [123,] "male"
## [124,] "female"
## [125,] "female"
## [126,] "male"
## [127,] "male"
## [128,] "male"
## [129,] "female"
## [130,] "female"
## [131,] "male"
## [132,] "male"
## [133,] "female"
## [134,] "female"
## [135,] "female"
## [136,] "female"
## [137,] "male"
## [138,] "female"
## [139,] "female"
## [140,] "female"
## [141,] "male"
## [142,] "female"
## [143,] "male"
## [144,] "female"
## [145,] "male"
## [146,] "male"
## [147,] "female"
## [148,] "female"
## [149,] "female"
## [150,] "male"
## [151,] "female"
## [152,] "male"
## [153,] "male"
## [154,] "female"
## [155,] "female"
## [156,] "female"
## [157,] "female"
## [158,] "female"
## [159,] "female"
## [160,] "female"
## [161,] "male"
## [162,] "male"
## [163,] "male"
## [164,] "male"
## [165,] "female"
## [166,] "female"
## [167,] "male"
## [168,] "male"
## [169,] "female"
## [170,] "male"
## [171,] "male"
## [172,] "male"
## [173,] "male"
## [174,] "female"
## [175,] "male"
## [176,] "male"
## [177,] "female"
## [178,] "male"
## [179,] "female"
## [180,] "male"
## [181,] "female"
## [182,] "male"
## [183,] "female"
## [184,] "female"
## [185,] "male"
## [186,] "female"
## [187,] "male"
## [188,] "female"
## [189,] "female"
## [190,] "male"
## [191,] "male"
## [192,] "male"
## [193,] "female"
## [194,] "female"
## [195,] "male"
## [196,] "male"
## [197,] "female"
## [198,] "male"
## [199,] "male"
## [200,] "female"
## [201,] "male"
## [202,] "male"
## [203,] "male"
## [204,] "male"
## [205,] "female"
## [206,] "female"
## [207,] "male"
## [208,] "female"
## [209,] "female"
## [210,] "male"
## [211,] "male"
## [212,] "female"
## [213,] "male"
## [214,] "male"
## [215,] "male"
## [216,] "male"
## [217,] "female"
## [218,] "female"
## [219,] "female"
## [220,] "male"
## [221,] "male"
## [222,] "female"
## [223,] "female"
## [224,] "female"
## [225,] "female"
## [226,] "female"
## [227,] "female"
## [228,] "female"
## [229,] "male"
## [230,] "female"
## [231,] "male"
## [232,] "female"
## [233,] "female"
## [234,] "female"
## [235,] "female"
## [236,] "female"
## [237,] "male"
## [238,] "female"
## [239,] "male"
## [240,] "male"
## [241,] "female"
## [242,] "male"
## [243,] "male"
## [244,] "female"
## [245,] "male"
## [246,] "male"
## [247,] "female"
## [248,] "female"
## [249,] "male"
## [250,] "male"
## [251,] "female"
## [252,] "female"
## [253,] "male"
## [254,] "female"
## [255,] "male"
## [256,] "female"
## [257,] "female"
## [258,] "male"
## [259,] "male"
## [260,] "female"
## [261,] "female"
## [262,] "male"
## [263,] "female"
## [264,] "female"
## [265,] "female"
## [266,] "female"
## [267,] "male"
## [268,] "male"
## [269,] "male"
## [270,] "female"
## [271,] "female"
## [272,] "female"
## [273,] "male"
## [274,] "female"
## [275,] "female"
## [276,] "female"
## [277,] "female"
## [278,] "female"
## [279,] "female"
## [280,] "male"
## [281,] "male"
## [282,] "female"
## [283,] "male"
## [284,] "male"
## [285,] "female"
## [286,] "male"
## [287,] "female"
## [288,] "male"
## [289,] "female"
## [290,] "female"
## [291,] "female"
## [292,] "male"
## [293,] "male"
## [294,] "female"
## [295,] "female"
## [296,] "female"
## [297,] "male"
## [298,] "male"
## [299,] "male"
## [300,] "male"
## [301,] "male"
## [302,] "male"
## [303,] "female"
## [304,] "male"
## [305,] "male"
## [306,] "male"
## [307,] "female"
## [308,] "female"
## [309,] "female"
## [310,] "female"
## [311,] "male"
## [312,] "male"
## [313,] "male"
## [314,] "male"
## [315,] "female"
## [316,] "female"
## [317,] "male"
## [318,] "female"
## [319,] "female"
## [320,] "female"
## [321,] "female"
## [322,] "male"
## [323,] "female"
## [324,] "male"
## [325,] "female"
## [326,] "female"
## [327,] "female"
## [328,] "female"
## [329,] "male"
## [330,] "male"
## [331,] "male"
## [332,] "female"
## [333,] "male"
## [334,] "female"
## [335,] "male"
## [336,] "female"
## [337,] "male"
## [338,] "male"
## [339,] "male"
## [340,] "male"
## [341,] "male"
## [342,] "female"
## [343,] "male"
## [344,] "male"
## [345,] "female"
## [346,] "female"
## [347,] "male"
## [348,] "female"
## [349,] "male"
## [350,] "female"
## [351,] "female"
## [352,] "female"
## [353,] "female"
## [354,] "male"
## [355,] "female"
## [356,] "female"
## [357,] "male"
## [358,] "male"
## [359,] "female"
## [360,] "female"
## [361,] "female"
## [362,] "female"
## [363,] "male"
## [364,] "female"
## [365,] "female"
## [366,] "male"
## [367,] "female"
## [368,] "male"
## [369,] "male"
## [370,] "male"
## [371,] "female"
## [372,] "male"
## [373,] "male"
## [374,] "male"
## [375,] "female"
## [376,] "female"
## [377,] "female"
## [378,] "male"
## [379,] "female"
## [380,] "female"
## [381,] "female"
## [382,] "female"
## [383,] "female"
## [384,] "female"
## [385,] "male"
## [386,] "male"
## [387,] "male"
## [388,] "male"
## [389,] "male"
## [390,] "male"
## [391,] "female"
## [392,] "male"
## [393,] "female"
## [394,] "male"
## [395,] "male"
## [396,] "female"
## [397,] "male"
## [398,] "female"
## [399,] "female"
## [400,] "female"
## [401,] "female"
## [402,] "male"
## [403,] "male"
## [404,] "male"
## [405,] "male"
## [406,] "female"
## [407,] "female"
## [408,] "female"
## [409,] "male"
## [410,] "female"
## [411,] "female"
## [412,] "male"
## [413,] "female"
## [414,] "male"
## [415,] "female"
## [416,] "male"
## [417,] "female"
## [418,] "male"
## [419,] "female"
## [420,] "female"
## [421,] "female"
## [422,] "female"
## [423,] "male"
## [424,] "male"
## [425,] "female"
## [426,] "female"
## [427,] "female"
## [428,] "female"
## [429,] "female"
## [430,] "female"
## [431,] "female"
## [432,] "male"
## [433,] "female"
## [434,] "female"
## [435,] "female"
## [436,] "female"
## [437,] "female"
## [438,] "male"
## [439,] "male"
## [440,] "male"
## [441,] "male"
## [442,] "female"
## [443,] "male"
## [444,] "female"
## [445,] "female"
## [446,] "male"
## [447,] "male"
## [448,] "female"
## [449,] "female"
## [450,] "female"
## [451,] "male"
## [452,] "female"
## [453,] "female"
## [454,] "male"
## [455,] "female"
## [456,] "male"
## [457,] "male"
## [458,] "male"
## [459,] "female"
## [460,] "female"
## [461,] "male"
## [462,] "female"
## [463,] "male"
## [464,] "female"
## [465,] "female"
## [466,] "female"
## [467,] "male"
## [468,] "male"
## [469,] "female"
## [470,] "male"
## [471,] "male"
## [472,] "female"
## [473,] "female"
## [474,] "female"
## [475,] "female"
## [476,] "female"
## [477,] "male"
## [478,] "male"
## [479,] "male"
## [480,] "female"
## [481,] "male"
## [482,] "female"
## [483,] "female"
## [484,] "male"
## [485,] "female"
## [486,] "male"
## [487,] "female"
## [488,] "male"
## [489,] "male"
## [490,] "female"
## [491,] "female"
## [492,] "female"
## [493,] "male"
## [494,] "female"
## [495,] "female"
## [496,] "male"
## [497,] "male"
## [498,] "female"
## [499,] "female"
## [500,] "male"
## [501,] "female"
## [502,] "male"
## [503,] "female"
## [504,] "male"
## [505,] "female"
## [506,] "female"
## [507,] "male"
## [508,] "female"
## [509,] "male"
## [510,] "male"
## [511,] "female"
## [512,] "female"
## [513,] "female"
## [514,] "male"
## [515,] "male"
## [516,] "female"
## [517,] "male"
## [518,] "male"
## [519,] "male"
## [520,] "male"
## [521,] "female"
## [522,] "male"
## [523,] "female"
## [524,] "male"
## [525,] "female"
## [526,] "female"
## [527,] "female"
## [528,] "female"
## [529,] "female"
## [530,] "male"
## [531,] "female"
## [532,] "female"
## [533,] "female"
## [534,] "female"
## [535,] "male"
## [536,] "female"
## [537,] "male"
## [538,] "female"
## [539,] "male"
## [540,] "male"
## [541,] "female"
## [542,] "female"
## [543,] "female"
## [544,] "male"
## [545,] "male"
## [546,] "female"
## [547,] "female"
## [548,] "female"
## [549,] "female"
## [550,] "female"
## [551,] "female"
## [552,] "female"
## [553,] "female"
## [554,] "female"
## [555,] "female"
## [556,] "male"
## [557,] "male"
## [558,] "female"
## [559,] "female"
## [560,] "female"
## [561,] "male"
## [562,] "male"
## [563,] "female"
## [564,] "female"
## [565,] "female"
## [566,] "male"
## [567,] "female"
## [568,] "male"
## [569,] "female"
## [570,] "male"
## [571,] "male"
## [572,] "male"
## [573,] "female"
## [574,] "male"
## [575,] "male"
## [576,] "male"
## [577,] "female"
## [578,] "male"
## [579,] "female"
## [580,] "male"
## [581,] "male"
## [582,] "female"
## [583,] "male"
## [584,] "male"
## [585,] "male"
## [586,] "male"
## [587,] "male"
## [588,] "male"
## [589,] "female"
## [590,] "male"
## [591,] "female"
## [592,] "male"
## [593,] "female"
## [594,] "male"
## [595,] "female"
## [596,] "female"
## [597,] "female"
## [598,] "female"
## [599,] "female"
## [600,] "female"mi.set %>% names()## [1] "locus_of_control" "self_concept" "motivation" "read"
## [5] "write" "math" "science" "sex"mi.set %>% rownames()## [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12"
## [13] "13" "14" "15" "16" "17" "18" "19" "20" "21" "22" "23" "24"
## [25] "25" "26" "27" "28" "29" "30" "31" "32" "33" "34" "35" "36"
## [37] "37" "38" "39" "40" "41" "42" "43" "44" "45" "46" "47" "48"
## [49] "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60"
## [61] "61" "62" "63" "64" "65" "66" "67" "68" "69" "70" "71" "72"
## [73] "73" "74" "75" "76" "77" "78" "79" "80" "81" "82" "83" "84"
## [85] "85" "86" "87" "88" "89" "90" "91" "92" "93" "94" "95" "96"
## [97] "97" "98" "99" "100" "101" "102" "103" "104" "105" "106" "107" "108"
## [109] "109" "110" "111" "112" "113" "114" "115" "116" "117" "118" "119" "120"
## [121] "121" "122" "123" "124" "125" "126" "127" "128" "129" "130" "131" "132"
## [133] "133" "134" "135" "136" "137" "138" "139" "140" "141" "142" "143" "144"
## [145] "145" "146" "147" "148" "149" "150" "151" "152" "153" "154" "155" "156"
## [157] "157" "158" "159" "160" "161" "162" "163" "164" "165" "166" "167" "168"
## [169] "169" "170" "171" "172" "173" "174" "175" "176" "177" "178" "179" "180"
## [181] "181" "182" "183" "184" "185" "186" "187" "188" "189" "190" "191" "192"
## [193] "193" "194" "195" "196" "197" "198" "199" "200" "201" "202" "203" "204"
## [205] "205" "206" "207" "208" "209" "210" "211" "212" "213" "214" "215" "216"
## [217] "217" "218" "219" "220" "221" "222" "223" "224" "225" "226" "227" "228"
## [229] "229" "230" "231" "232" "233" "234" "235" "236" "237" "238" "239" "240"
## [241] "241" "242" "243" "244" "245" "246" "247" "248" "249" "250" "251" "252"
## [253] "253" "254" "255" "256" "257" "258" "259" "260" "261" "262" "263" "264"
## [265] "265" "266" "267" "268" "269" "270" "271" "272" "273" "274" "275" "276"
## [277] "277" "278" "279" "280" "281" "282" "283" "284" "285" "286" "287" "288"
## [289] "289" "290" "291" "292" "293" "294" "295" "296" "297" "298" "299" "300"
## [301] "301" "302" "303" "304" "305" "306" "307" "308" "309" "310" "311" "312"
## [313] "313" "314" "315" "316" "317" "318" "319" "320" "321" "322" "323" "324"
## [325] "325" "326" "327" "328" "329" "330" "331" "332" "333" "334" "335" "336"
## [337] "337" "338" "339" "340" "341" "342" "343" "344" "345" "346" "347" "348"
## [349] "349" "350" "351" "352" "353" "354" "355" "356" "357" "358" "359" "360"
## [361] "361" "362" "363" "364" "365" "366" "367" "368" "369" "370" "371" "372"
## [373] "373" "374" "375" "376" "377" "378" "379" "380" "381" "382" "383" "384"
## [385] "385" "386" "387" "388" "389" "390" "391" "392" "393" "394" "395" "396"
## [397] "397" "398" "399" "400" "401" "402" "403" "404" "405" "406" "407" "408"
## [409] "409" "410" "411" "412" "413" "414" "415" "416" "417" "418" "419" "420"
## [421] "421" "422" "423" "424" "425" "426" "427" "428" "429" "430" "431" "432"
## [433] "433" "434" "435" "436" "437" "438" "439" "440" "441" "442" "443" "444"
## [445] "445" "446" "447" "448" "449" "450" "451" "452" "453" "454" "455" "456"
## [457] "457" "458" "459" "460" "461" "462" "463" "464" "465" "466" "467" "468"
## [469] "469" "470" "471" "472" "473" "474" "475" "476" "477" "478" "479" "480"
## [481] "481" "482" "483" "484" "485" "486" "487" "488" "489" "490" "491" "492"
## [493] "493" "494" "495" "496" "497" "498" "499" "500" "501" "502" "503" "504"
## [505] "505" "506" "507" "508" "509" "510" "511" "512" "513" "514" "515" "516"
## [517] "517" "518" "519" "520" "521" "522" "523" "524" "525" "526" "527" "528"
## [529] "529" "530" "531" "532" "533" "534" "535" "536" "537" "538" "539" "540"
## [541] "541" "542" "543" "544" "545" "546" "547" "548" "549" "550" "551" "552"
## [553] "553" "554" "555" "556" "557" "558" "559" "560" "561" "562" "563" "564"
## [565] "565" "566" "567" "568" "569" "570" "571" "572" "573" "574" "575" "576"
## [577] "577" "578" "579" "580" "581" "582" "583" "584" "585" "586" "587" "588"
## [589] "589" "590" "591" "592" "593" "594" "595" "596" "597" "598" "599" "600"mi.set %>% colnames()## [1] "locus_of_control" "self_concept" "motivation" "read"
## [5] "write" "math" "science" "sex"mi.set %>% attributes()## $names
## [1] "locus_of_control" "self_concept" "motivation" "read"
## [5] "write" "math" "science" "sex"
##
## $row.names
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
## [19] 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
## [37] 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
## [55] 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72
## [73] 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90
## [91] 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108
## [109] 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126
## [127] 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144
## [145] 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162
## [163] 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180
## [181] 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198
## [199] 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216
## [217] 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234
## [235] 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252
## [253] 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270
## [271] 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288
## [289] 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306
## [307] 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324
## [325] 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342
## [343] 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360
## [361] 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378
## [379] 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396
## [397] 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414
## [415] 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432
## [433] 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450
## [451] 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468
## [469] 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486
## [487] 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504
## [505] 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522
## [523] 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540
## [541] 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558
## [559] 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576
## [577] 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594
## [595] 595 596 597 598 599 600
##
## $academic
## [1] "read" "write" "math" "science" "sex"
##
## $psychology
## [1] "locus_of_control" "self_concept" "motivation"
##
## $class
## [1] "data.frame"mi.set %>% str()## 'data.frame': 600 obs. of 8 variables:
## $ locus_of_control: num -0.84 -0.38 0.89 0.71 -0.64 1.11 0.06 -0.91 0.45 0 ...
## $ self_concept : num -0.24 -0.47 0.59 0.28 0.03 0.9 0.03 -0.59 0.03 0.03 ...
## $ motivation : chr "4" "3" "3" "3" ...
## $ read : num 54.8 62.7 60.6 62.7 41.6 62.7 41.6 44.2 62.7 62.7 ...
## $ write : num 64.5 43.7 56.7 56.7 46.3 64.5 39.1 39.1 51.5 64.5 ...
## $ math : num 44.5 44.7 70.5 54.7 38.4 61.4 56.3 46.3 54.4 38.3 ...
## $ science : num 52.6 52.6 58 58 36.3 58 45 36.3 49.8 55.8 ...
## $ sex : chr "female" "female" "male" "male" ...
## - attr(*, "academic")= chr "read" "write" "math" "science" ...
## - attr(*, "psychology")= chr "locus_of_control" "self_concept" "motivation"mi.set %>% glimpse()## Rows: 600
## Columns: 8
## $ locus_of_control <dbl> -0.84, -0.38, 0.89, 0.71, -0.64, 1.11, 0.06, -0.91, …
## $ self_concept <dbl> -0.24, -0.47, 0.59, 0.28, 0.03, 0.90, 0.03, -0.59, 0…
## $ motivation <chr> "4", "3", "3", "3", "4", "2", "3", "3", "4", "3", "2…
## $ read <dbl> 54.8, 62.7, 60.6, 62.7, 41.6, 62.7, 41.6, 44.2, 62.7…
## $ write <dbl> 64.5, 43.7, 56.7, 56.7, 46.3, 64.5, 39.1, 39.1, 51.5…
## $ math <dbl> 44.5, 44.7, 70.5, 54.7, 38.4, 61.4, 56.3, 46.3, 54.4…
## $ science <dbl> 52.6, 52.6, 58.0, 58.0, 36.3, 58.0, 45.0, 36.3, 49.8…
## $ sex <chr> "female", "female", "male", "male", "female", "femal…mi.set %>% colnames() %>% class()## [1] "character"nombres.filas <- mi.set %>% rownames()
mi.set %>% colnames() -> nombres.columnas
nombres.columnas %>% str_replace_all('_',' ')## [1] "locus of control" "self concept" "motivation" "read"
## [5] "write" "math" "science" "sex"nombres.columnas %>%
str_replace_all(c("locus_of_control" = "control", "self_concept" = "consciencia", 'motivation'='motivación',
'read'='lectura','write'='escritura','math'='matemáticas', 'science'='ciencias', 'sex'='sexo')) -> nuevas.colscolnames(mi.set) <- nuevas.cols
attr(mi.set, "psychology") <- nuevas.cols[c(1,2,3)]
attr(mi.set, "academic") <- nuevas.cols[c(4,5,6,7,8)]
attributes(mi.set)## $names
## [1] "control" "consciencia" "motivación" "lectura" "escritura"
## [6] "matemáticas" "ciencias" "sexo"
##
## $row.names
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
## [19] 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
## [37] 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
## [55] 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72
## [73] 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90
## [91] 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108
## [109] 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126
## [127] 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144
## [145] 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162
## [163] 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180
## [181] 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198
## [199] 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216
## [217] 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234
## [235] 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252
## [253] 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270
## [271] 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288
## [289] 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306
## [307] 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324
## [325] 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342
## [343] 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360
## [361] 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378
## [379] 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396
## [397] 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414
## [415] 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432
## [433] 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450
## [451] 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468
## [469] 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486
## [487] 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504
## [505] 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522
## [523] 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540
## [541] 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558
## [559] 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576
## [577] 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594
## [595] 595 596 597 598 599 600
##
## $academic
## [1] "lectura" "escritura" "matemáticas" "ciencias" "sexo"
##
## $psychology
## [1] "control" "consciencia" "motivación"
##
## $class
## [1] "data.frame"mi.set$motivación <- as.factor(mi.set$motivación)
dp <- ggplot(mi.set, aes(x=motivación, y=matemáticas, fill=motivación)) +
geom_violin(trim=FALSE)+
geom_boxplot(width=0.1, fill="white")+
labs(title="Plot of length by dose",x="motivación", y = "ciencias")
dp + scale_fill_brewer(palette="Blues") + theme_classic()
psych_variables <- attr(mi.set, "psychology")
academic_variables <- attr(mi.set, "academic")
ggpairs(mi.set, academic_variables, title = "Within Academic Variables")
mi.set2 = sets.de.datos$diamantes
p1 <- ggplot(
subset(mi.set2, carat >= 2.2),
aes(x = table, y = price, colour = cut)
) +
geom_point(alpha = 0.7) +
geom_smooth(method = "loess", alpha = 0.05, size = 1, span = 1) +
theme_bw()
p2 <- ggplot(
subset(mi.set2, carat > 2.2 & depth > 55 & depth < 70),
aes(x = depth, fill = cut)
) +
geom_histogram(colour = "black", binwidth = 1, position = "dodge") +
theme_bw()p1_lancet <- p1 + scale_color_lancet()
p2_lancet <- p2 + scale_fill_lancet()
grid.arrange(p1_lancet, p2_lancet, ncol = 2)
# Automatically generate interactive plot for results produced by `stats::prcomp`
p <- autoplotly(prcomp(iris[c(1, 2, 3, 4)]), data = iris,
colour = 'Species', label = TRUE, label.size = 3, frame = TRUE)
# You can apply additional ggplot2 elements to the generated interactive plot
p +
ggplot2::ggtitle("Principal Components Analysis") +
ggplot2::labs(y = "Second Principal Components", x = "First Principal Components")# Or apply additional plotly elements to the generated interactive plot
p %>% plotly::layout(annotations = list(
text = "Example Text",
font = list(
family = "Courier New, monospace",
size = 18,
color = "black"),
x = 0,
y = 0,
showarrow = TRUE))edges3 <- data.frame(Parent = paste0("clone_",
LETTERS[c(rep(1:3, each = 2), 2, 5)]),
Identity = paste0("clone_", LETTERS[2:9]))
# a function for generating exponential growth curves:
pop_seq <- function(gens, lambda, start_gen) c(rep(0, start_gen),
exp(lambda * gens[0:(length(gens) - start_gen)]))
lambda <- 0.1 # baseline fitness
gens <- 0:150 # time points
fitnesses <- c(1, 2, 2.2, 2.5, 3, 3.2, 3.5, 3.5, 3.8) # relative fitnesses of genotypes
pop3 <- data.frame(Generation = rep(gens, 9),
Identity = paste0("clone_", LETTERS[rep(1:9, each = length(gens))]),
Population = c(1E2 * pop_seq(gens, fitnesses[1]*lambda, 0),
pop_seq(gens, fitnesses[2]*lambda, 0),
pop_seq(gens, fitnesses[3]*lambda, 10),
pop_seq(gens, fitnesses[4]*lambda, 20),
pop_seq(gens, fitnesses[5]*lambda, 30),
pop_seq(gens, fitnesses[6]*lambda, 40),
pop_seq(gens, fitnesses[7]*lambda, 50),
pop_seq(gens, fitnesses[8]*lambda, 50),
pop_seq(gens, fitnesses[9]*lambda, 60)),
Fitness = rep(fitnesses, each = length(gens)))
Muller_df3 <- get_Muller_df(edges3, pop3)Muller_plot(Muller_df3, add_legend = TRUE, xlab = "Time", ylab = "Proportion")
Muller_plot(Muller_df3, colour_by = "Fitness", add_legend = TRUE)
p <- ggplot(Muller_df3, aes(Generation,Population, size = Frequency , colour = Identity)) +
geom_point(alpha = 0.7, show.legend = T) +
scale_size(range = c(2, 12)) +
scale_y_log10() +
labs(title = 'Generation: {frame_time}', x = 'Generation', y = 'Population') +
transition_time(Generation) +
ease_aes('linear')
p
# Render and show the video
#q <- 2
#df <- animate(p, renderer = av_renderer('animation.mp4'), width = 720*q, height = 480*q, res = 72*q, fps = 25)
#utils::browseURL('animation.mp4')library(readxl)
download.file(url='https://www.pnas.org/highwire/filestream/794560/field_highwire_adjunct_files/2/pnas.1800165115.sd03.xlsx', destfile='pnas.1800165115.sd03.xlsx', method='wget')
clark_2018 <- read_excel("pnas.1800165115.sd03.xlsx")
clark_2018 %>% head()## # A tibble: 6 x 11
## GeneID `Entrez Gene Na… P7MeanFPKM P32MeanFPKM `10wkMeanFPKM` `9.5mMeanFPKM`
## <chr> <chr> <dbl> <dbl> <dbl> <dbl>
## 1 Oasl2 2'-5' oligoaden… 1.11 2.03 1.36 3.39
## 2 Ctss cathepsin S 0.73 0.613 5.92 7.88
## 3 C4b complement comp… 1.19 2.46 2.62 3.95
## 4 Lyz2 lysozyme 2.03 1.65 3.01 4.06
## 5 Serpi… serpin peptidas… 1.29 2 1.71 4.06
## 6 Heyl hes-related fam… 15.1 5.89 3.06 9.37
## # … with 5 more variables: `2yrMeanFPKM` <dbl>, `log2 FC` <dbl>, logCPM <dbl>,
## # PValue <dbl>, FDR <dbl>