Tarea 1: Usar un algoritmo de clasificación que impacte directamente tu vida diaria.

library(openxlsx)
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.1     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.0
## ✔ ggplot2   3.4.1     ✔ tibble    3.2.1
## ✔ lubridate 1.9.2     ✔ tidyr     1.3.0
## ✔ purrr     1.0.1     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(factoextra)
## Welcome! Want to learn more? See two factoextra-related books at https://goo.gl/ve3WBa
library(cowplot)
## 
## Attaching package: 'cowplot'
## 
## The following object is masked from 'package:lubridate':
## 
##     stamp
library(ggpubr)
## 
## Attaching package: 'ggpubr'
## 
## The following object is masked from 'package:cowplot':
## 
##     get_legend
library(cluster)
library(purrr)
library(dplyr)
Comentario
Libros <- read.xlsx("/Users/eduardoleyva/Desktop/Evidencia 2/Top_100_Trending_Books.xlsx")
summary(Libros)
##       Rank         book.title          book.price         rating    
##  Min.   :  1.00   Length:100         Min.   : 2.780   Min.   :4.10  
##  1st Qu.: 25.75   Class :character   1st Qu.: 6.303   1st Qu.:4.60  
##  Median : 50.50   Mode  :character   Median :11.480   Median :4.70  
##  Mean   : 50.50                      Mean   :12.709   Mean   :4.69  
##  3rd Qu.: 75.25                      3rd Qu.:16.990   3rd Qu.:4.80  
##  Max.   :100.00                      Max.   :48.770   Max.   :5.00  
##                                                       NA's   :3     
##     author          year.of.publication    genre               url           
##  Length:100         Min.   :1947        Length:100         Length:100        
##  Class :character   1st Qu.:2014        Class :character   Class :character  
##  Mode  :character   Median :2019        Mode  :character   Mode  :character  
##                     Mean   :2014                                             
##                     3rd Qu.:2023                                             
##                     Max.   :2024                                             
## 
Comentario
numeric_columns <- Libros[sapply(Libros, is.numeric)]
Escalar <- scale(numeric_columns, center = TRUE, scale = TRUE)
Escalar <- as.data.frame(Escalar)

Escalar <- Escalar[complete.cases(Escalar), ]
str(Escalar)
## 'data.frame':    97 obs. of  4 variables:
##  $ Rank               : num  -1.71 -1.67 -1.64 -1.6 -1.57 ...
##  $ book.price         : num  0.722 1.039 2.374 1.425 -0.892 ...
##  $ rating             : num  -3.256 -1.047 -1.047 -1.599 0.609 ...
##  $ year.of.publication: num  0.611 0.611 0.611 0.611 0.28 ...
Comentario
total_within = function(n_clusters, data, iter.max=1000, nstart=50){
  
  cluster_means = kmeans(data,centers = n_clusters,
                         iter.max = iter.max,
                         nstart = nstart)
  return(cluster_means$tot.withinss)
}

total_withinss <- map_dbl(.x = 1:15,
                          .f = total_within,
                          data = Escalar)
total_withinss
##  [1] 388.01141 286.64929 224.51635 170.23330 135.40300 114.27131  97.98697
##  [8]  87.25978  77.25084  68.97324  63.33957  58.50453  53.56437  48.97433
## [15]  45.60682
data.frame(n_clusters = 1:15, suma_cuadrados_internos = total_withinss) %>%
  ggplot(aes(x = n_clusters, y = suma_cuadrados_internos)) +
  geom_line() +
  geom_point() +
  scale_x_continuous(breaks = 1:15) +
  labs(title = "Suma total de cuadrados intra-cluster") +
  theme_bw()

Comentario
kmcluster = kmeans(Escalar,centers=4,nstart = 50)
kmcluster
## K-means clustering with 4 clusters of sizes 40, 32, 10, 15
## 
## Cluster means:
##         Rank book.price     rating year.of.publication
## 1 -0.4903229 -0.4272087  0.4847771           0.2356170
## 2  1.0534618  0.2557650 -0.1846194           0.2492484
## 3  0.1240888 -0.9026010  0.7746188          -2.6145882
## 4 -1.0352220  1.0399000 -1.4152967           0.4917224
## 
## Clustering vector:
##   1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19  20 
##   4   4   4   4   1   1   4   1   4   4   1   3   1   1   1   4   4   4   1   1 
##  21  22  23  24  26  27  28  29  30  31  32  33  34  35  36  37  38  39  40  41 
##   1   1   1   4   1   1   1   1   3   3   1   1   4   1   3   1   1   1   1   1 
##  42  43  44  45  47  48  49  50  51  52  53  54  55  56  57  58  59  60  61  62 
##   1   1   3   1   4   1   1   1   4   1   1   1   1   3   1   2   2   1   2   1 
##  63  64  65  66  67  68  69  70  71  72  73  74  75  76  77  78  79  80  81  82 
##   1   4   2   2   1   3   1   2   2   2   2   2   2   3   2   2   2   2   2   2 
##  83  84  85  87  88  89  90  91  92  93  94  95  96  97  98  99 100 
##   2   2   2   2   2   2   2   3   2   2   2   2   2   3   2   2   2 
## 
## Within cluster sum of squares by cluster:
## [1] 47.44112 60.03961 20.42951 42.32306
##  (between_SS / total_SS =  56.1 %)
## 
## Available components:
## 
## [1] "cluster"      "centers"      "totss"        "withinss"     "tot.withinss"
## [6] "betweenss"    "size"         "iter"         "ifault"
Comentario
fviz_cluster(kmcluster, Escalar)+
  theme_minimal()

Comentario
Escalar <- na.omit(Escalar)
fviz_cluster(kmcluster, Escalar, show.clust.cent = T,
             ellipse.type = "euclid", star.plot = T, repel = T) +
  labs(title = "Resultados clustering K-means") +
  theme_bw()

Comentario
Libros2 <- dist(Libros, method = 'euclidean') #Sacamos la distancia euclidiana de los puntos
## Warning in dist(Libros, method = "euclidean"): NAs introduced by coercion
Libros3 <- hclust(Libros2, method = 'average') #Hacemos el análisis de clusters

plot(Libros3, cex=0.5, col="red", hang = -1, 
     main="Dendograma, Distancia Euclídea, Método completo")
rect.hclust(Libros3, k = 4, border = c("blue", "green", "purple", "orange"))

Análisis de Resultados

Tarea 2: Buscar datos del agua en algún estado del país, y entrenar una red neuronal.

library(neuralnet)
## 
## Attaching package: 'neuralnet'
## The following object is masked from 'package:dplyr':
## 
##     compute
library(stats)
library(psych)
## 
## Attaching package: 'psych'
## The following objects are masked from 'package:ggplot2':
## 
##     %+%, alpha
library(MASS)
## 
## Attaching package: 'MASS'
## The following object is masked from 'package:dplyr':
## 
##     select
library(ISLR)
library(fRegression)
library(vcd)
## Loading required package: grid
## 
## Attaching package: 'vcd'
## The following object is masked from 'package:ISLR':
## 
##     Hitters
library(dplyr)
library(mlbench)
library(magrittr)
## 
## Attaching package: 'magrittr'
## The following object is masked from 'package:purrr':
## 
##     set_names
## The following object is masked from 'package:tidyr':
## 
##     extract
library(neuralnet)
library(keras)
library(caret)
## Loading required package: lattice
## 
## Attaching package: 'caret'
## The following object is masked from 'package:purrr':
## 
##     lift

Leer base de datos:

Comentario
Agua <- read.xlsx("/Users/eduardoleyva/Desktop/Evidencia 2/El_Mojarral_20082.xlsx")
str(Agua)
## 'data.frame':    4322 obs. of  4 variables:
##  $ Fecha      : num  39448 39448 39448 39448 39448 ...
##  $ Hora       : num  0 0.0417 0.0833 0.125 0.1667 ...
##  $ Nivel      : num  0.393 0.393 0.396 0.399 0.397 ...
##  $ Temperatura: num  28.7 28.5 28.4 28.1 27.6 ...
colSums(is.na(Agua))
##       Fecha        Hora       Nivel Temperatura 
##           0           0           0           0

Ajustar y Visualizar Redes Neuronales:

Comentario
set.seed(13)
train = createDataPartition(y = Agua$Temperatura, p=0.7,list=FALSE, times = 1)
Agua_train = Agua[train,]
Agua_test = Agua[-train,]
RN <- neuralnet(Temperatura ~ Fecha+Hora+Nivel,
                data = Agua_train,
                hidden = c(8,6),
                linear.output = T,
                lifesign = 'full',
                threshold = 0.05,
                rep=1)
## hidden: 8, 6    thresh: 0.05    rep: 1/1    steps:      56   error: 921.67135    time: 0.17 secs
Comentario
#train/test split en matrices y separando variable a predecir
training <- as.matrix(Agua_train[, c(1, 2, 3)])
trainingtarget <- as.matrix(Agua_train[,4])
test <- as.matrix(Agua_test[, c(1, 2, 3)])
testtarget <- as.matrix(Agua_test[,4])

#Estandarización de variables
m <- colMeans(training) #Obtener medias por columna
s <- apply(training, 2, sd) #Calcular StandDev por columna (por ello el apply lleva el 2, si pusieran 1 sería por renglón)
training <- scale(training, center = m, scale = s)
test <- scale(test, center = m, scale = s)
Comentario
Agua_train_S <- as.data.frame(cbind(training,(trainingtarget - mean(trainingtarget))/sd(trainingtarget)))
colnames(Agua_train_S) <- colnames(Agua_train)

RNS <- neuralnet(Temperatura ~ Fecha+Hora+Nivel,
               data = Agua_train_S,
               hidden = c(8,6),
               linear.output = T, 
               #linear.output se debe poner como T en modelos de regresion y como F en modelos de clasificación
               lifesign = 'full',
               rep=1,
               # threshold=0.02,
               stepmax=4000000)
## hidden: 8, 6    thresh: 0.01    rep: 1/1    steps:    1000   min thresh: 1.45716721363483
##                                                       2000   min thresh: 0.987372860026007
##                                                       3000   min thresh: 0.4278997409004
##                                                       4000   min thresh: 0.397328373352213
##                                                       5000   min thresh: 0.317897572474276
##                                                       6000   min thresh: 0.244005529679809
##                                                       7000   min thresh: 0.223498894845103
##                                                       8000   min thresh: 0.21607577665216
##                                                       9000   min thresh: 0.21607577665216
##                                                      10000   min thresh: 0.21607577665216
##                                                      11000   min thresh: 0.21607577665216
##                                                      12000   min thresh: 0.194131654692289
##                                                      13000   min thresh: 0.194131654692289
##                                                      14000   min thresh: 0.194131654692289
##                                                      15000   min thresh: 0.193754805303248
##                                                      16000   min thresh: 0.190566650517719
##                                                      17000   min thresh: 0.172340473288804
##                                                      18000   min thresh: 0.154847574608854
##                                                      19000   min thresh: 0.138599661518144
##                                                      20000   min thresh: 0.138599661518144
##                                                      21000   min thresh: 0.138599661518144
##                                                      22000   min thresh: 0.138599661518144
##                                                      23000   min thresh: 0.138599661518144
##                                                      24000   min thresh: 0.138599661518144
##                                                      25000   min thresh: 0.0973981775147962
##                                                      26000   min thresh: 0.0950049912111781
##                                                      27000   min thresh: 0.0950049912111781
##                                                      28000   min thresh: 0.0950049912111781
##                                                      29000   min thresh: 0.0950049912111781
##                                                      30000   min thresh: 0.0950049912111781
##                                                      31000   min thresh: 0.0950049912111781
##                                                      32000   min thresh: 0.0950049912111781
##                                                      33000   min thresh: 0.0950049912111781
##                                                      34000   min thresh: 0.0950049912111781
##                                                      35000   min thresh: 0.0950049912111781
##                                                      36000   min thresh: 0.0950049912111781
##                                                      37000   min thresh: 0.0950049912111781
##                                                      38000   min thresh: 0.0950049912111781
##                                                      39000   min thresh: 0.0950049912111781
##                                                      40000   min thresh: 0.0950049912111781
##                                                      41000   min thresh: 0.0950049912111781
##                                                      42000   min thresh: 0.0950049912111781
##                                                      43000   min thresh: 0.0950049912111781
##                                                      44000   min thresh: 0.0950049912111781
##                                                      45000   min thresh: 0.0950049912111781
##                                                      46000   min thresh: 0.0950049912111781
##                                                      47000   min thresh: 0.0871072649873084
##                                                      48000   min thresh: 0.0871072649873084
##                                                      49000   min thresh: 0.0871072649873084
##                                                      50000   min thresh: 0.0835197739399279
##                                                      51000   min thresh: 0.0835197739399279
##                                                      52000   min thresh: 0.0835197739399279
##                                                      53000   min thresh: 0.0835197739399279
##                                                      54000   min thresh: 0.0835197739399279
##                                                      55000   min thresh: 0.0835197739399279
##                                                      56000   min thresh: 0.0835197739399279
##                                                      57000   min thresh: 0.0835197739399279
##                                                      58000   min thresh: 0.0835197739399279
##                                                      59000   min thresh: 0.0830505697769267
##                                                      60000   min thresh: 0.0830505697769267
##                                                      61000   min thresh: 0.0830505697769267
##                                                      62000   min thresh: 0.0806155955728976
##                                                      63000   min thresh: 0.0774870426314016
##                                                      64000   min thresh: 0.0699481233159423
##                                                      65000   min thresh: 0.0699481233159423
##                                                      66000   min thresh: 0.0699481233159423
##                                                      67000   min thresh: 0.0699481233159423
##                                                      68000   min thresh: 0.0699481233159423
##                                                      69000   min thresh: 0.0699481233159423
##                                                      70000   min thresh: 0.0699481233159423
##                                                      71000   min thresh: 0.0699481233159423
##                                                      72000   min thresh: 0.0699481233159423
##                                                      73000   min thresh: 0.0699481233159423
##                                                      74000   min thresh: 0.0699481233159423
##                                                      75000   min thresh: 0.0699481233159423
##                                                      76000   min thresh: 0.0699481233159423
##                                                      77000   min thresh: 0.0699481233159423
##                                                      78000   min thresh: 0.0699481233159423
##                                                      79000   min thresh: 0.0699481233159423
##                                                      80000   min thresh: 0.0699481233159423
##                                                      81000   min thresh: 0.0699481233159423
##                                                      82000   min thresh: 0.0689794595911323
##                                                      83000   min thresh: 0.0618888029649585
##                                                      84000   min thresh: 0.0618888029649585
##                                                      85000   min thresh: 0.0618888029649585
##                                                      86000   min thresh: 0.0618888029649585
##                                                      87000   min thresh: 0.0618888029649585
##                                                      88000   min thresh: 0.0618888029649585
##                                                      89000   min thresh: 0.0618888029649585
##                                                      90000   min thresh: 0.0618888029649585
##                                                      91000   min thresh: 0.0618888029649585
##                                                      92000   min thresh: 0.0618888029649585
##                                                      93000   min thresh: 0.0615997512135286
##                                                      94000   min thresh: 0.0615997512135286
##                                                      95000   min thresh: 0.0615997512135286
##                                                      96000   min thresh: 0.0562842989084955
##                                                      97000   min thresh: 0.0562842989084955
##                                                      98000   min thresh: 0.0562842989084955
##                                                      99000   min thresh: 0.0504368327217668
##                                                      1e+05   min thresh: 0.0504368327217668
##                                                     101000   min thresh: 0.0504368327217668
##                                                     102000   min thresh: 0.0485489121376166
##                                                     103000   min thresh: 0.0485489121376166
##                                                     104000   min thresh: 0.0485489121376166
##                                                     105000   min thresh: 0.0485489121376166
##                                                     106000   min thresh: 0.0485489121376166
##                                                     107000   min thresh: 0.0485489121376166
##                                                     108000   min thresh: 0.0485489121376166
##                                                     109000   min thresh: 0.0485489121376166
##                                                     110000   min thresh: 0.0485489121376166
##                                                     111000   min thresh: 0.0482135670447367
##                                                     112000   min thresh: 0.0482135670447367
##                                                     113000   min thresh: 0.0482135670447367
##                                                     114000   min thresh: 0.0476906770286839
##                                                     115000   min thresh: 0.0476906770286839
##                                                     116000   min thresh: 0.0476906770286839
##                                                     117000   min thresh: 0.0476906770286839
##                                                     118000   min thresh: 0.0458577599686919
##                                                     119000   min thresh: 0.0439749183828209
##                                                     120000   min thresh: 0.0439749183828209
##                                                     121000   min thresh: 0.0435758136664992
##                                                     122000   min thresh: 0.0435758136664992
##                                                     123000   min thresh: 0.0435758136664992
##                                                     124000   min thresh: 0.0435758136664992
##                                                     125000   min thresh: 0.0435758136664992
##                                                     126000   min thresh: 0.0435758136664992
##                                                     127000   min thresh: 0.0435758136664992
##                                                     128000   min thresh: 0.0435758136664992
##                                                     129000   min thresh: 0.0401619757439649
##                                                     130000   min thresh: 0.0401619757439649
##                                                     131000   min thresh: 0.0401619757439649
##                                                     132000   min thresh: 0.0401619757439649
##                                                     133000   min thresh: 0.0401619757439649
##                                                     134000   min thresh: 0.0401619757439649
##                                                     135000   min thresh: 0.0401619757439649
##                                                     136000   min thresh: 0.0401619757439649
##                                                     137000   min thresh: 0.0401619757439649
##                                                     138000   min thresh: 0.0401619757439649
##                                                     139000   min thresh: 0.0401619757439649
##                                                     140000   min thresh: 0.0401619757439649
##                                                     141000   min thresh: 0.0401619757439649
##                                                     142000   min thresh: 0.0401619757439649
##                                                     143000   min thresh: 0.0401619757439649
##                                                     144000   min thresh: 0.0401619757439649
##                                                     145000   min thresh: 0.0401619757439649
##                                                     146000   min thresh: 0.0401619757439649
##                                                     147000   min thresh: 0.0401619757439649
##                                                     148000   min thresh: 0.0401619757439649
##                                                     149000   min thresh: 0.0401619757439649
##                                                     150000   min thresh: 0.0401619757439649
##                                                     151000   min thresh: 0.0401619757439649
##                                                     152000   min thresh: 0.0401619757439649
##                                                     153000   min thresh: 0.0401619757439649
##                                                     154000   min thresh: 0.0401619757439649
##                                                     155000   min thresh: 0.0401619757439649
##                                                     156000   min thresh: 0.0401619757439649
##                                                     157000   min thresh: 0.0401619757439649
##                                                     158000   min thresh: 0.0384054705450031
##                                                     159000   min thresh: 0.0384054705450031
##                                                     160000   min thresh: 0.0384054705450031
##                                                     161000   min thresh: 0.0384054705450031
##                                                     162000   min thresh: 0.0384054705450031
##                                                     163000   min thresh: 0.0384054705450031
##                                                     164000   min thresh: 0.0384054705450031
##                                                     165000   min thresh: 0.0384054705450031
##                                                     166000   min thresh: 0.0384054705450031
##                                                     167000   min thresh: 0.0375953111595492
##                                                     168000   min thresh: 0.0375953111595492
##                                                     169000   min thresh: 0.0375953111595492
##                                                     170000   min thresh: 0.0375953111595492
##                                                     171000   min thresh: 0.0375953111595492
##                                                     172000   min thresh: 0.0354241355438344
##                                                     173000   min thresh: 0.0354241355438344
##                                                     174000   min thresh: 0.0335440990138856
##                                                     175000   min thresh: 0.0335440990138856
##                                                     176000   min thresh: 0.0335440990138856
##                                                     177000   min thresh: 0.0335440990138856
##                                                     178000   min thresh: 0.0335440990138856
##                                                     179000   min thresh: 0.0335440990138856
##                                                     180000   min thresh: 0.0335440990138856
##                                                     181000   min thresh: 0.0335440990138856
##                                                     182000   min thresh: 0.0335440990138856
##                                                     183000   min thresh: 0.0335440990138856
##                                                     184000   min thresh: 0.0335440990138856
##                                                     185000   min thresh: 0.0335440990138856
##                                                     186000   min thresh: 0.0335440990138856
##                                                     187000   min thresh: 0.0327874666383128
##                                                     188000   min thresh: 0.0327874666383128
##                                                     189000   min thresh: 0.0327874666383128
##                                                     190000   min thresh: 0.0327874666383128
##                                                     191000   min thresh: 0.0327874666383128
##                                                     192000   min thresh: 0.0327874666383128
##                                                     193000   min thresh: 0.0327874666383128
##                                                     194000   min thresh: 0.0327874666383128
##                                                     195000   min thresh: 0.0327874666383128
##                                                     196000   min thresh: 0.0313985103728857
##                                                     197000   min thresh: 0.0313985103728857
##                                                     198000   min thresh: 0.0313985103728857
##                                                     199000   min thresh: 0.0313985103728857
##                                                      2e+05   min thresh: 0.0312864439799706
##                                                     201000   min thresh: 0.0312864439799706
##                                                     202000   min thresh: 0.0312864439799706
##                                                     203000   min thresh: 0.0312864439799706
##                                                     204000   min thresh: 0.0312864439799706
##                                                     205000   min thresh: 0.0312864439799706
##                                                     206000   min thresh: 0.0312864439799706
##                                                     207000   min thresh: 0.030495036643212
##                                                     208000   min thresh: 0.029404528002537
##                                                     209000   min thresh: 0.029404528002537
##                                                     210000   min thresh: 0.029404528002537
##                                                     211000   min thresh: 0.029404528002537
##                                                     212000   min thresh: 0.029404528002537
##                                                     213000   min thresh: 0.029404528002537
##                                                     214000   min thresh: 0.029404528002537
##                                                     215000   min thresh: 0.029404528002537
##                                                     216000   min thresh: 0.029404528002537
##                                                     217000   min thresh: 0.029404528002537
##                                                     218000   min thresh: 0.029404528002537
##                                                     219000   min thresh: 0.029404528002537
##                                                     220000   min thresh: 0.029404528002537
##                                                     221000   min thresh: 0.029404528002537
##                                                     222000   min thresh: 0.0293090059218773
##                                                     223000   min thresh: 0.028718355065831
##                                                     224000   min thresh: 0.028718355065831
##                                                     225000   min thresh: 0.028718355065831
##                                                     226000   min thresh: 0.028718355065831
##                                                     227000   min thresh: 0.028718355065831
##                                                     228000   min thresh: 0.028718355065831
##                                                     229000   min thresh: 0.028718355065831
##                                                     230000   min thresh: 0.028718355065831
##                                                     231000   min thresh: 0.028718355065831
##                                                     232000   min thresh: 0.028718355065831
##                                                     233000   min thresh: 0.028718355065831
##                                                     234000   min thresh: 0.028718355065831
##                                                     235000   min thresh: 0.028718355065831
##                                                     236000   min thresh: 0.028718355065831
##                                                     237000   min thresh: 0.028718355065831
##                                                     238000   min thresh: 0.028718355065831
##                                                     239000   min thresh: 0.028718355065831
##                                                     240000   min thresh: 0.028718355065831
##                                                     241000   min thresh: 0.028718355065831
##                                                     242000   min thresh: 0.028718355065831
##                                                     243000   min thresh: 0.028718355065831
##                                                     244000   min thresh: 0.028718355065831
##                                                     245000   min thresh: 0.028718355065831
##                                                     246000   min thresh: 0.028718355065831
##                                                     247000   min thresh: 0.028718355065831
##                                                     248000   min thresh: 0.028718355065831
##                                                     249000   min thresh: 0.028718355065831
##                                                     250000   min thresh: 0.028718355065831
##                                                     251000   min thresh: 0.028718355065831
##                                                     252000   min thresh: 0.028718355065831
##                                                     253000   min thresh: 0.028718355065831
##                                                     254000   min thresh: 0.028718355065831
##                                                     255000   min thresh: 0.028718355065831
##                                                     256000   min thresh: 0.028718355065831
##                                                     257000   min thresh: 0.028718355065831
##                                                     258000   min thresh: 0.028718355065831
##                                                     259000   min thresh: 0.028718355065831
##                                                     260000   min thresh: 0.028718355065831
##                                                     261000   min thresh: 0.028718355065831
##                                                     262000   min thresh: 0.028718355065831
##                                                     263000   min thresh: 0.028718355065831
##                                                     264000   min thresh: 0.028718355065831
##                                                     265000   min thresh: 0.028718355065831
##                                                     266000   min thresh: 0.028718355065831
##                                                     267000   min thresh: 0.028718355065831
##                                                     268000   min thresh: 0.028718355065831
##                                                     269000   min thresh: 0.028718355065831
##                                                     270000   min thresh: 0.028718355065831
##                                                     271000   min thresh: 0.028718355065831
##                                                     272000   min thresh: 0.028718355065831
##                                                     273000   min thresh: 0.028718355065831
##                                                     274000   min thresh: 0.028718355065831
##                                                     275000   min thresh: 0.028718355065831
##                                                     276000   min thresh: 0.028718355065831
##                                                     277000   min thresh: 0.028718355065831
##                                                     278000   min thresh: 0.028718355065831
##                                                     279000   min thresh: 0.028718355065831
##                                                     280000   min thresh: 0.028718355065831
##                                                     281000   min thresh: 0.028718355065831
##                                                     282000   min thresh: 0.028718355065831
##                                                     283000   min thresh: 0.028718355065831
##                                                     284000   min thresh: 0.028718355065831
##                                                     285000   min thresh: 0.028718355065831
##                                                     286000   min thresh: 0.028718355065831
##                                                     287000   min thresh: 0.028718355065831
##                                                     288000   min thresh: 0.028718355065831
##                                                     289000   min thresh: 0.028718355065831
##                                                     290000   min thresh: 0.028718355065831
##                                                     291000   min thresh: 0.028718355065831
##                                                     292000   min thresh: 0.028718355065831
##                                                     293000   min thresh: 0.028718355065831
##                                                     294000   min thresh: 0.028718355065831
##                                                     295000   min thresh: 0.028718355065831
##                                                     296000   min thresh: 0.028718355065831
##                                                     297000   min thresh: 0.028718355065831
##                                                     298000   min thresh: 0.028718355065831
##                                                     299000   min thresh: 0.028718355065831
##                                                      3e+05   min thresh: 0.028718355065831
##                                                     301000   min thresh: 0.028718355065831
##                                                     302000   min thresh: 0.028718355065831
##                                                     303000   min thresh: 0.028718355065831
##                                                     304000   min thresh: 0.028718355065831
##                                                     305000   min thresh: 0.028718355065831
##                                                     306000   min thresh: 0.028718355065831
##                                                     307000   min thresh: 0.028718355065831
##                                                     308000   min thresh: 0.028718355065831
##                                                     309000   min thresh: 0.028718355065831
##                                                     310000   min thresh: 0.028718355065831
##                                                     311000   min thresh: 0.028718355065831
##                                                     312000   min thresh: 0.028718355065831
##                                                     313000   min thresh: 0.028718355065831
##                                                     314000   min thresh: 0.028718355065831
##                                                     315000   min thresh: 0.028718355065831
##                                                     316000   min thresh: 0.028718355065831
##                                                     317000   min thresh: 0.028718355065831
##                                                     318000   min thresh: 0.028718355065831
##                                                     319000   min thresh: 0.028718355065831
##                                                     320000   min thresh: 0.028718355065831
##                                                     321000   min thresh: 0.028718355065831
##                                                     322000   min thresh: 0.028718355065831
##                                                     323000   min thresh: 0.028718355065831
##                                                     324000   min thresh: 0.028718355065831
##                                                     325000   min thresh: 0.028718355065831
##                                                     326000   min thresh: 0.028718355065831
##                                                     327000   min thresh: 0.028718355065831
##                                                     328000   min thresh: 0.028718355065831
##                                                     329000   min thresh: 0.028718355065831
##                                                     330000   min thresh: 0.028718355065831
##                                                     331000   min thresh: 0.028718355065831
##                                                     332000   min thresh: 0.028718355065831
##                                                     333000   min thresh: 0.028718355065831
##                                                     334000   min thresh: 0.028718355065831
##                                                     335000   min thresh: 0.028718355065831
##                                                     336000   min thresh: 0.028718355065831
##                                                     337000   min thresh: 0.028718355065831
##                                                     338000   min thresh: 0.028718355065831
##                                                     339000   min thresh: 0.028718355065831
##                                                     340000   min thresh: 0.028718355065831
##                                                     341000   min thresh: 0.028718355065831
##                                                     342000   min thresh: 0.028718355065831
##                                                     343000   min thresh: 0.028718355065831
##                                                     344000   min thresh: 0.028718355065831
##                                                     345000   min thresh: 0.028718355065831
##                                                     346000   min thresh: 0.028718355065831
##                                                     347000   min thresh: 0.028718355065831
##                                                     348000   min thresh: 0.028718355065831
##                                                     349000   min thresh: 0.028718355065831
##                                                     350000   min thresh: 0.028718355065831
##                                                     351000   min thresh: 0.028718355065831
##                                                     352000   min thresh: 0.028718355065831
##                                                     353000   min thresh: 0.028718355065831
##                                                     354000   min thresh: 0.028718355065831
##                                                     355000   min thresh: 0.028718355065831
##                                                     356000   min thresh: 0.028718355065831
##                                                     357000   min thresh: 0.028718355065831
##                                                     358000   min thresh: 0.0286874785545854
##                                                     359000   min thresh: 0.028640661724615
##                                                     360000   min thresh: 0.028640661724615
##                                                     361000   min thresh: 0.028640661724615
##                                                     362000   min thresh: 0.028640661724615
##                                                     363000   min thresh: 0.028640661724615
##                                                     364000   min thresh: 0.028640661724615
##                                                     365000   min thresh: 0.028640661724615
##                                                     366000   min thresh: 0.028640661724615
##                                                     367000   min thresh: 0.028640661724615
##                                                     368000   min thresh: 0.0272449814177336
##                                                     369000   min thresh: 0.0272449814177336
##                                                     370000   min thresh: 0.0272449814177336
##                                                     371000   min thresh: 0.0272449814177336
##                                                     372000   min thresh: 0.0272449814177336
##                                                     373000   min thresh: 0.0272449814177336
##                                                     374000   min thresh: 0.0272285501945671
##                                                     375000   min thresh: 0.0272285501945671
##                                                     376000   min thresh: 0.0256338789811313
##                                                     377000   min thresh: 0.0256338789811313
##                                                     378000   min thresh: 0.0256338789811313
##                                                     379000   min thresh: 0.0256338789811313
##                                                     380000   min thresh: 0.0256338789811313
##                                                     381000   min thresh: 0.0256338789811313
##                                                     382000   min thresh: 0.0256338789811313
##                                                     383000   min thresh: 0.0256338789811313
##                                                     384000   min thresh: 0.0256338789811313
##                                                     385000   min thresh: 0.0254560840199696
##                                                     386000   min thresh: 0.0241917491991012
##                                                     387000   min thresh: 0.0241917491991012
##                                                     388000   min thresh: 0.0241917491991012
##                                                     389000   min thresh: 0.0241917491991012
##                                                     390000   min thresh: 0.0241917491991012
##                                                     391000   min thresh: 0.0241917491991012
##                                                     392000   min thresh: 0.0241917491991012
##                                                     393000   min thresh: 0.0241917491991012
##                                                     394000   min thresh: 0.0241917491991012
##                                                     395000   min thresh: 0.0241917491991012
##                                                     396000   min thresh: 0.0241468104681154
##                                                     397000   min thresh: 0.0235531404948609
##                                                     398000   min thresh: 0.0235531404948609
##                                                     399000   min thresh: 0.0235531404948609
##                                                      4e+05   min thresh: 0.0229995900981518
##                                                     401000   min thresh: 0.0229995900981518
##                                                     402000   min thresh: 0.0229995900981518
##                                                     403000   min thresh: 0.0229995900981518
##                                                     404000   min thresh: 0.0229995900981518
##                                                     405000   min thresh: 0.0229995900981518
##                                                     406000   min thresh: 0.0228457262511518
##                                                     407000   min thresh: 0.0228457262511518
##                                                     408000   min thresh: 0.0213438994876207
##                                                     409000   min thresh: 0.0213438994876207
##                                                     410000   min thresh: 0.0213438994876207
##                                                     411000   min thresh: 0.0213438994876207
##                                                     412000   min thresh: 0.0213438994876207
##                                                     413000   min thresh: 0.0213438994876207
##                                                     414000   min thresh: 0.0213438994876207
##                                                     415000   min thresh: 0.0213438994876207
##                                                     416000   min thresh: 0.0209253026101019
##                                                     417000   min thresh: 0.0199210707896106
##                                                     418000   min thresh: 0.0199210707896106
##                                                     419000   min thresh: 0.0199210707896106
##                                                     420000   min thresh: 0.0199210707896106
##                                                     421000   min thresh: 0.0199210707896106
##                                                     422000   min thresh: 0.0199210707896106
##                                                     423000   min thresh: 0.0197704601100825
##                                                     424000   min thresh: 0.0197704601100825
##                                                     425000   min thresh: 0.0197704601100825
##                                                     426000   min thresh: 0.0197704601100825
##                                                     427000   min thresh: 0.0197017530733182
##                                                     428000   min thresh: 0.0197017530733182
##                                                     429000   min thresh: 0.0190040168388245
##                                                     430000   min thresh: 0.0190040168388245
##                                                     431000   min thresh: 0.0190040168388245
##                                                     432000   min thresh: 0.0188435319571201
##                                                     433000   min thresh: 0.0188435319571201
##                                                     434000   min thresh: 0.0188435319571201
##                                                     435000   min thresh: 0.0183672030075905
##                                                     436000   min thresh: 0.0183672030075905
##                                                     437000   min thresh: 0.0183672030075905
##                                                     438000   min thresh: 0.0183672030075905
##                                                     439000   min thresh: 0.0183672030075905
##                                                     440000   min thresh: 0.0183672030075905
##                                                     441000   min thresh: 0.0177022706600465
##                                                     442000   min thresh: 0.0177022706600465
##                                                     443000   min thresh: 0.0177022706600465
##                                                     444000   min thresh: 0.0177022706600465
##                                                     445000   min thresh: 0.0177022706600465
##                                                     446000   min thresh: 0.0165103106872605
##                                                     447000   min thresh: 0.0165103106872605
##                                                     448000   min thresh: 0.0165103106872605
##                                                     449000   min thresh: 0.0165103106872605
##                                                     450000   min thresh: 0.0165103106872605
##                                                     451000   min thresh: 0.0165103106872605
##                                                     452000   min thresh: 0.0165103106872605
##                                                     453000   min thresh: 0.0165103106872605
##                                                     454000   min thresh: 0.0165103106872605
##                                                     455000   min thresh: 0.0165103106872605
##                                                     456000   min thresh: 0.0165103106872605
##                                                     457000   min thresh: 0.0165103106872605
##                                                     458000   min thresh: 0.0165103106872605
##                                                     459000   min thresh: 0.0165103106872605
##                                                     460000   min thresh: 0.0165103106872605
##                                                     461000   min thresh: 0.0165103106872605
##                                                     462000   min thresh: 0.0165103106872605
##                                                     463000   min thresh: 0.0165103106872605
##                                                     464000   min thresh: 0.0165103106872605
##                                                     465000   min thresh: 0.0165103106872605
##                                                     466000   min thresh: 0.0165103106872605
##                                                     467000   min thresh: 0.0165103106872605
##                                                     468000   min thresh: 0.0165103106872605
##                                                     469000   min thresh: 0.0165103106872605
##                                                     470000   min thresh: 0.0163238931996648
##                                                     471000   min thresh: 0.0163238931996648
##                                                     472000   min thresh: 0.0160968497507938
##                                                     473000   min thresh: 0.0160968497507938
##                                                     474000   min thresh: 0.0160968497507938
##                                                     475000   min thresh: 0.0160968497507938
##                                                     476000   min thresh: 0.0160968497507938
##                                                     477000   min thresh: 0.0160968497507938
##                                                     478000   min thresh: 0.0160968497507938
##                                                     479000   min thresh: 0.0158015925581235
##                                                     480000   min thresh: 0.0158015925581235
##                                                     481000   min thresh: 0.0158015925581235
##                                                     482000   min thresh: 0.0158015925581235
##                                                     483000   min thresh: 0.0158015925581235
##                                                     484000   min thresh: 0.0158015925581235
##                                                     485000   min thresh: 0.0158015925581235
##                                                     486000   min thresh: 0.0158015925581235
##                                                     487000   min thresh: 0.0158015925581235
##                                                     488000   min thresh: 0.0158015925581235
##                                                     489000   min thresh: 0.0158015925581235
##                                                     490000   min thresh: 0.0158015925581235
##                                                     491000   min thresh: 0.0158015925581235
##                                                     492000   min thresh: 0.0158015925581235
##                                                     493000   min thresh: 0.0158015925581235
##                                                     494000   min thresh: 0.0158015925581235
##                                                     495000   min thresh: 0.0158015925581235
##                                                     496000   min thresh: 0.0158015925581235
##                                                     497000   min thresh: 0.0158015925581235
##                                                     498000   min thresh: 0.0158015925581235
##                                                     499000   min thresh: 0.0158015925581235
##                                                      5e+05   min thresh: 0.0158015925581235
##                                                     501000   min thresh: 0.0158015925581235
##                                                     502000   min thresh: 0.0158015925581235
##                                                     503000   min thresh: 0.0158015925581235
##                                                     504000   min thresh: 0.0158015925581235
##                                                     505000   min thresh: 0.0158015925581235
##                                                     506000   min thresh: 0.0158015925581235
##                                                     507000   min thresh: 0.0158015925581235
##                                                     508000   min thresh: 0.0158015925581235
##                                                     509000   min thresh: 0.0158015925581235
##                                                     510000   min thresh: 0.0158015925581235
##                                                     511000   min thresh: 0.0158015925581235
##                                                     512000   min thresh: 0.0158015925581235
##                                                     513000   min thresh: 0.0158015925581235
##                                                     514000   min thresh: 0.0158015925581235
##                                                     515000   min thresh: 0.0158015925581235
##                                                     516000   min thresh: 0.0158015925581235
##                                                     517000   min thresh: 0.0158015925581235
##                                                     518000   min thresh: 0.0158015925581235
##                                                     519000   min thresh: 0.0158015925581235
##                                                     520000   min thresh: 0.0158015925581235
##                                                     521000   min thresh: 0.0158015925581235
##                                                     522000   min thresh: 0.0158015925581235
##                                                     523000   min thresh: 0.0158015925581235
##                                                     524000   min thresh: 0.0158015925581235
##                                                     525000   min thresh: 0.0158015925581235
##                                                     526000   min thresh: 0.0158015925581235
##                                                     527000   min thresh: 0.0158015925581235
##                                                     528000   min thresh: 0.0158015925581235
##                                                     529000   min thresh: 0.0158015925581235
##                                                     530000   min thresh: 0.0158015925581235
##                                                     531000   min thresh: 0.0158015925581235
##                                                     532000   min thresh: 0.0158015925581235
##                                                     533000   min thresh: 0.0158015925581235
##                                                     534000   min thresh: 0.0158015925581235
##                                                     535000   min thresh: 0.0158015925581235
##                                                     536000   min thresh: 0.0158015925581235
##                                                     537000   min thresh: 0.0158015925581235
##                                                     538000   min thresh: 0.0158015925581235
##                                                     539000   min thresh: 0.0158015925581235
##                                                     540000   min thresh: 0.0158015925581235
##                                                     541000   min thresh: 0.0158015925581235
##                                                     542000   min thresh: 0.0158015925581235
##                                                     543000   min thresh: 0.0158015925581235
##                                                     544000   min thresh: 0.0158015925581235
##                                                     545000   min thresh: 0.0158015925581235
##                                                     546000   min thresh: 0.0158015925581235
##                                                     547000   min thresh: 0.0158015925581235
##                                                     548000   min thresh: 0.0158015925581235
##                                                     549000   min thresh: 0.0158015925581235
##                                                     550000   min thresh: 0.0158015925581235
##                                                     551000   min thresh: 0.0158015925581235
##                                                     552000   min thresh: 0.0158015925581235
##                                                     553000   min thresh: 0.0158015925581235
##                                                     554000   min thresh: 0.0158015925581235
##                                                     555000   min thresh: 0.0158015925581235
##                                                     556000   min thresh: 0.0158015925581235
##                                                     557000   min thresh: 0.0158015925581235
##                                                     558000   min thresh: 0.0158015925581235
##                                                     559000   min thresh: 0.0158015925581235
##                                                     560000   min thresh: 0.0158015925581235
##                                                     561000   min thresh: 0.0158015925581235
##                                                     562000   min thresh: 0.0158015925581235
##                                                     563000   min thresh: 0.0158015925581235
##                                                     564000   min thresh: 0.0158015925581235
##                                                     565000   min thresh: 0.0158015925581235
##                                                     566000   min thresh: 0.0158015925581235
##                                                     567000   min thresh: 0.0158015925581235
##                                                     568000   min thresh: 0.0158015925581235
##                                                     569000   min thresh: 0.0158015925581235
##                                                     570000   min thresh: 0.0158015925581235
##                                                     571000   min thresh: 0.0158015925581235
##                                                     572000   min thresh: 0.0158015925581235
##                                                     573000   min thresh: 0.0158015925581235
##                                                     574000   min thresh: 0.0158015925581235
##                                                     575000   min thresh: 0.0158015925581235
##                                                     576000   min thresh: 0.0158015925581235
##                                                     577000   min thresh: 0.0158015925581235
##                                                     578000   min thresh: 0.0158015925581235
##                                                     579000   min thresh: 0.0158015925581235
##                                                     580000   min thresh: 0.0158015925581235
##                                                     581000   min thresh: 0.0158015925581235
##                                                     582000   min thresh: 0.0158015925581235
##                                                     583000   min thresh: 0.0158015925581235
##                                                     584000   min thresh: 0.0158015925581235
##                                                     585000   min thresh: 0.0158015925581235
##                                                     586000   min thresh: 0.0158015925581235
##                                                     587000   min thresh: 0.0158015925581235
##                                                     588000   min thresh: 0.0158015925581235
##                                                     589000   min thresh: 0.0158015925581235
##                                                     590000   min thresh: 0.0158015925581235
##                                                     591000   min thresh: 0.0158015925581235
##                                                     592000   min thresh: 0.0158015925581235
##                                                     593000   min thresh: 0.0158015925581235
##                                                     594000   min thresh: 0.0158015925581235
##                                                     595000   min thresh: 0.0158015925581235
##                                                     596000   min thresh: 0.0158015925581235
##                                                     597000   min thresh: 0.0158015925581235
##                                                     598000   min thresh: 0.0158015925581235
##                                                     599000   min thresh: 0.0158015925581235
##                                                      6e+05   min thresh: 0.0158015925581235
##                                                     601000   min thresh: 0.0158015925581235
##                                                     602000   min thresh: 0.0158015925581235
##                                                     603000   min thresh: 0.0158015925581235
##                                                     604000   min thresh: 0.0158015925581235
##                                                     605000   min thresh: 0.0158015925581235
##                                                     606000   min thresh: 0.0158015925581235
##                                                     607000   min thresh: 0.0158015925581235
##                                                     608000   min thresh: 0.0158015925581235
##                                                     609000   min thresh: 0.0154552140870816
##                                                     610000   min thresh: 0.0154552140870816
##                                                     611000   min thresh: 0.0154552140870816
##                                                     612000   min thresh: 0.0154552140870816
##                                                     613000   min thresh: 0.0154552140870816
##                                                     614000   min thresh: 0.0154552140870816
##                                                     615000   min thresh: 0.0154552140870816
##                                                     616000   min thresh: 0.0154552140870816
##                                                     617000   min thresh: 0.0154552140870816
##                                                     618000   min thresh: 0.0154552140870816
##                                                     619000   min thresh: 0.0154552140870816
##                                                     620000   min thresh: 0.0154552140870816
##                                                     621000   min thresh: 0.0154552140870816
##                                                     622000   min thresh: 0.0151832676823584
##                                                     623000   min thresh: 0.0150551275284734
##                                                     624000   min thresh: 0.0146191325772651
##                                                     625000   min thresh: 0.0133434385107595
##                                                     626000   min thresh: 0.0133434385107595
##                                                     627000   min thresh: 0.0133434385107595
##                                                     628000   min thresh: 0.0133434385107595
##                                                     629000   min thresh: 0.0133434385107595
##                                                     630000   min thresh: 0.0133434385107595
##                                                     631000   min thresh: 0.0133434385107595
##                                                     632000   min thresh: 0.0133434385107595
##                                                     633000   min thresh: 0.0133434385107595
##                                                     634000   min thresh: 0.0133434385107595
##                                                     635000   min thresh: 0.0133434385107595
##                                                     636000   min thresh: 0.0133434385107595
##                                                     637000   min thresh: 0.0133434385107595
##                                                     638000   min thresh: 0.0133434385107595
##                                                     639000   min thresh: 0.0133434385107595
##                                                     640000   min thresh: 0.0133434385107595
##                                                     641000   min thresh: 0.0133434385107595
##                                                     642000   min thresh: 0.0133434385107595
##                                                     643000   min thresh: 0.0133434385107595
##                                                     644000   min thresh: 0.0133434385107595
##                                                     645000   min thresh: 0.0133434385107595
##                                                     646000   min thresh: 0.0133434385107595
##                                                     647000   min thresh: 0.0133434385107595
##                                                     648000   min thresh: 0.0133434385107595
##                                                     649000   min thresh: 0.0131598904887637
##                                                     650000   min thresh: 0.0131598904887637
##                                                     651000   min thresh: 0.0131598904887637
##                                                     652000   min thresh: 0.0131598904887637
##                                                     653000   min thresh: 0.0131598904887637
##                                                     654000   min thresh: 0.0131598904887637
##                                                     655000   min thresh: 0.0131598904887637
##                                                     656000   min thresh: 0.0131598904887637
##                                                     657000   min thresh: 0.0131598904887637
##                                                     658000   min thresh: 0.0131017481909182
##                                                     659000   min thresh: 0.0130772779105495
##                                                     660000   min thresh: 0.0130772779105495
##                                                     661000   min thresh: 0.0130772779105495
##                                                     662000   min thresh: 0.0129969739191129
##                                                     663000   min thresh: 0.0129969739191129
##                                                     664000   min thresh: 0.0125573917270332
##                                                     665000   min thresh: 0.0125573917270332
##                                                     666000   min thresh: 0.0125573917270332
##                                                     667000   min thresh: 0.0125573917270332
##                                                     668000   min thresh: 0.0125573917270332
##                                                     669000   min thresh: 0.0125573917270332
##                                                     670000   min thresh: 0.0125573917270332
##                                                     671000   min thresh: 0.0125573917270332
##                                                     672000   min thresh: 0.0125573917270332
##                                                     673000   min thresh: 0.0125573917270332
##                                                     674000   min thresh: 0.0125573917270332
##                                                     675000   min thresh: 0.0125573917270332
##                                                     676000   min thresh: 0.0125573917270332
##                                                     677000   min thresh: 0.0125573917270332
##                                                     678000   min thresh: 0.0125573917270332
##                                                     679000   min thresh: 0.0125573917270332
##                                                     680000   min thresh: 0.0125573917270332
##                                                     681000   min thresh: 0.0125573917270332
##                                                     682000   min thresh: 0.0125573917270332
##                                                     683000   min thresh: 0.0125573917270332
##                                                     684000   min thresh: 0.0124366742778832
##                                                     685000   min thresh: 0.0124366742778832
##                                                     686000   min thresh: 0.0124366742778832
##                                                     687000   min thresh: 0.0124366742778832
##                                                     688000   min thresh: 0.012078609985964
##                                                     689000   min thresh: 0.012078609985964
##                                                     690000   min thresh: 0.012078609985964
##                                                     691000   min thresh: 0.012078609985964
##                                                     692000   min thresh: 0.012078609985964
##                                                     693000   min thresh: 0.0118578246123828
##                                                     694000   min thresh: 0.0118578246123828
##                                                     695000   min thresh: 0.0118578246123828
##                                                     696000   min thresh: 0.0118578246123828
##                                                     697000   min thresh: 0.0118578246123828
##                                                     698000   min thresh: 0.0118578246123828
##                                                     699000   min thresh: 0.0118578246123828
##                                                      7e+05   min thresh: 0.0118578246123828
##                                                     701000   min thresh: 0.0117662253566247
##                                                     702000   min thresh: 0.0117662253566247
##                                                     703000   min thresh: 0.0117662253566247
##                                                     704000   min thresh: 0.0117662253566247
##                                                     705000   min thresh: 0.0117565974398038
##                                                     706000   min thresh: 0.0116277007661352
##                                                     707000   min thresh: 0.0116277007661352
##                                                     708000   min thresh: 0.0116277007661352
##                                                     709000   min thresh: 0.0116277007661352
##                                                     710000   min thresh: 0.0116277007661352
##                                                     711000   min thresh: 0.0116277007661352
##                                                     712000   min thresh: 0.0112592258345485
##                                                     713000   min thresh: 0.0112592258345485
##                                                     714000   min thresh: 0.0112592258345485
##                                                     715000   min thresh: 0.0112592258345485
##                                                     716000   min thresh: 0.0112592258345485
##                                                     717000   min thresh: 0.0112592258345485
##                                                     718000   min thresh: 0.0112592258345485
##                                                     719000   min thresh: 0.0112592258345485
##                                                     720000   min thresh: 0.0112592258345485
##                                                     721000   min thresh: 0.0112592258345485
##                                                     722000   min thresh: 0.0106879962460397
##                                                     723000   min thresh: 0.0106879962460397
##                                                     724000   min thresh: 0.0106879962460397
##                                                     725000   min thresh: 0.0106879962460397
##                                                     726000   min thresh: 0.0106879962460397
##                                                     727000   min thresh: 0.0106879962460397
##                                                     728000   min thresh: 0.0106879962460397
##                                                     729000   min thresh: 0.0106879962460397
##                                                     730000   min thresh: 0.0106879962460397
##                                                     731000   min thresh: 0.0106879962460397
##                                                     732000   min thresh: 0.0106879962460397
##                                                     733000   min thresh: 0.0106879962460397
##                                                     734000   min thresh: 0.0106879962460397
##                                                     735000   min thresh: 0.0106879962460397
##                                                     736000   min thresh: 0.0106879962460397
##                                                     737000   min thresh: 0.0106879962460397
##                                                     738000   min thresh: 0.0106879962460397
##                                                     739000   min thresh: 0.0106879962460397
##                                                     740000   min thresh: 0.0106879962460397
##                                                     741000   min thresh: 0.0106879962460397
##                                                     742000   min thresh: 0.0106879962460397
##                                                     743000   min thresh: 0.0106879962460397
##                                                     744000   min thresh: 0.0106879962460397
##                                                     745000   min thresh: 0.0106879962460397
##                                                     746000   min thresh: 0.0106438318106519
##                                                     747000   min thresh: 0.0106438318106519
##                                                     748000   min thresh: 0.0106438318106519
##                                                     749000   min thresh: 0.0106438318106519
##                                                     750000   min thresh: 0.0106438318106519
##                                                     751000   min thresh: 0.0106438318106519
##                                                     752000   min thresh: 0.0106438318106519
##                                                     753000   min thresh: 0.0106438318106519
##                                                     754000   min thresh: 0.0106438318106519
##                                                     755000   min thresh: 0.0106438318106519
##                                                     756000   min thresh: 0.0106438318106519
##                                                     757000   min thresh: 0.0106438318106519
##                                                     758000   min thresh: 0.0106438318106519
##                                                     759000   min thresh: 0.0104524272514834
##                                                     760000   min thresh: 0.0104524272514834
##                                                     761000   min thresh: 0.0104524272514834
##                                                     762000   min thresh: 0.0104511861405935
##                                                     763000   min thresh: 0.0104511861405935
##                                                     764000   min thresh: 0.0104511861405935
##                                                     765000   min thresh: 0.0103793680183831
##                                                     766000   min thresh: 0.0103793680183831
##                                                     767000   min thresh: 0.0103793680183831
##                                                     768000   min thresh: 0.0103793680183831
##                                                     769000   min thresh: 0.0102070261739132
##                                                     770000   min thresh: 0.0102070261739132
##                                                     771000   min thresh: 0.0102070261739132
##                                                     772000   min thresh: 0.0100205114943799
##                                                     773000   min thresh: 0.0100205114943799
##                                                     774000   min thresh: 0.0100205114943799
##                                                     775000   min thresh: 0.0100205114943799
##                                                     776000   min thresh: 0.0100205114943799
##                                                     777000   min thresh: 0.0100205114943799
##                                                     778000   min thresh: 0.0100205114943799
##                                                     779000   min thresh: 0.0100205114943799
##                                                     780000   min thresh: 0.0100205114943799
##                                                     781000   min thresh: 0.0100205114943799
##                                                     782000   min thresh: 0.0100205114943799
##                                                     783000   min thresh: 0.0100205114943799
##                                                     784000   min thresh: 0.0100205114943799
##                                                     785000   min thresh: 0.0100205114943799
##                                                     786000   min thresh: 0.0100205114943799
##                                                     787000   min thresh: 0.0100205114943799
##                                                     788000   min thresh: 0.0100205114943799
##                                                     789000   min thresh: 0.0100205114943799
##                                                     790000   min thresh: 0.0100205114943799
##                                                     790492   error: 32.94333 time: 13.24 mins
Comentario
plot(RNS,col.hidden = 'darkgreen',     
col.hidden.synapse = 'darkgreen',
     show.weights = T,
     information = F,
     fill = 'lightblue')
Comentario
Agua_test_S <- as.data.frame(test)
colnames(Agua_test_S) <- colnames(Agua_test)[1:3]

RNSPredictions <- predict(RNS,Agua_test_S)
cor(RNSPredictions,(testtarget-mean(trainingtarget))/sd(trainingtarget))
##           [,1]
## [1,] 0.9872836
Comentario
RNSPred <- RNSPredictions*sd(trainingtarget) + mean(trainingtarget)
plot(RNSPred,testtarget)
abline(a=0, b=1)

RSSnn <- (RNSPred - testtarget)^2
sum(RSSnn)/nrow(testtarget)
## [1] 0.01570494
1 - sum(RSSnn)/sum((testtarget - mean(trainingtarget))^2)
## [1] 0.9747231
Comentario
LRM <- lm(Temperatura ~ Fecha+Hora+Nivel, data=Agua_train)
LRMPred <- predict(LRM, Agua_test[, c(1, 2, 3)])
cor(LRMPred,Agua_test[,4])
## [1] 0.7921928
plot(LRMPred,Agua_test[,4])
abline(a=0,b=1)

LRMRSS <- (LRMPred - Agua_test[,4])^2
sum(LRMRSS)/nrow(testtarget)
## [1] 0.2319171
1 - sum(LRMRSS)/sum((Agua_test[,4] - mean(Agua_train[, c(1, 2, 3)]))^2)
## Warning in mean.default(Agua_train[, c(1, 2, 3)]): argument is not numeric or
## logical: returning NA
## [1] NA
Comentario
par(mfrow=c(1,2))
plot(Agua_test$Temperatura,RNSPred,col='red',main='Real vs predicted NN',pch=19,cex=1)
abline(0,1,lwd=2)
legend('bottomright',legend='NN',pch=18,col='red', bty='n')
plot(Agua_test$Temperatura,LRMPred,col='blue',main='Real vs predicted LM',pch=15, cex=1)
abline(0,1,lwd=2)
legend('bottomright',legend='LM',pch=18,col='blue', bty='n', cex=.95)

Comentario
plot(Agua_test$Temperatura,RNSPred,col='red',main='Real vs predicted NN',pch=19,cex=1)
points(Agua_test$Temperatura,LRMPred,col='blue',pch=15,cex=1)
abline(0,1,lwd=2)
legend('bottomright',legend=c('NN','LM'),pch=c(19,15),col=c('red','blue'))

#### Análisis de Resultados