Introducción a R

Antes de empezar con cualquier tipo de análisis estadístico, se necesita saber la manipulación básica del lenguaje de r.

x <- 5  #Esto es una variable numérica

y <- 10  #Esto también es una variable numérica

nombre <- "Ana"    #Esto es una variable tipo "string "

# Operaciones básicas
suma <- x + y
producto <- x * y #Con las variables numéricas puedes hacer cualquier tipo de operación matemática

saludo <- paste("Hola", nombre) 

# Mostrar resultados
print(suma)
## [1] 15
print(producto)
## [1] 50
print(saludo)
## [1] "Hola Ana"
#print() sirve para mostrar lo que contenga una variable

Ahora, verémos las estructuras que pueden tener los datos que estemos manejando:

# Vectores

edades <- c(20, 25, 30, 22, 20, 32, 28, 25, 24)
nombres <- c("Juan", "María", "Carlos", "Luisa")

# Data frames
datos <- data.frame(
  nombre = c("Ana", "Luis", "Carlos"),
  edad = c(25, 30, 22),
  ciudad = c("Bogotá", "Medellín", "Cali")
)

# Listas
mi_lista <- list(
  vector = edades,
  dataframe = datos,
  texto = "Ejemplo de lista"
)

#Como nota, tanto en las listas como en los dataframes, los vectores que se les asignen todos deben tener la misma canntidad de datos.

Ahora, vamos a ver cómo podemos manipular todos estos datos dependiendo de lo que se requiera.

# Crear un data frame de ejemplo
set.seed(123)
datos_clase <- data.frame(
  id = 1:20,
  edad = sample(18:30, 20, replace = TRUE),
  promedio = round(runif(20, 3.0, 5.0), 2),
  semestre = sample(1:10, 20, replace = TRUE)
)

# Filtrar datos
estudiantes_sobresalientes <- datos_clase %>% 
  filter(promedio > 4.0)

# Ordenar datos
datos_ordenados <- datos_clase %>% 
  arrange(desc(promedio), edad)

# Agregar nueva columna
datos_clase <- datos_clase %>% 
  mutate(categoria = ifelse(promedio >= 4.0, "Alto", "Bajo"))

Por último, vamos a mostrar algunas herramientas útiles que pueden llegar a ser de ayuda dependiendo de la situación:

vector_aleatorio<- runif(5, min = 0, max = 100) ##crea un vector aleaorio con numeros decimales
print(vector_aleatorio)
## [1] 37.98165 61.27710 35.17979 11.11354 24.36195
vector_aleatorio2<-sample(1:100, 6, replace=TRUE) ##crea un vector aleaorio con numeros enteros
print(vector_aleatorio2)
## [1]  6 72 86 86 39 31
vector_aleatorio3<-sample(1:100,size = 6) ##otra forma
vector_aleatorio4<-sample.int(18,45,size=5)   ## crea solo valores enteros
generos<-sample(c("M","F","O"),10,replace = TRUE)##crea un vecto co elementos aleatorios definidos
##"replace = TRUE": indica si el vector puede tener elementos repetidos
##sample: se le puede pedir cialquier tipo de datos
##sample.int: solo se le pueden pedir entros


matriz<- matrix(1:4, nrow = 2, ncol = 2) ##crea una matriz 
print(matriz)
##      [,1] [,2]
## [1,]    1    3
## [2,]    2    4
set.seed(123) ## se escribe antes de crear daos aleatorios para fijarlos
edad<-round(runif(6, min=18, max=65)) ##crea un vector aleatorio de numeros decimales los cuales se redondean con round
print(edad)
## [1] 32 55 37 60 62 20
data_aleatoria<-data.frame()
df<-data.frame(id=sample(1:100,10), edad= sample(18:65, 10, replace=TRUE), ciudad=sample(c("Bogota", "Medellin", "Cali", "b"),10,replace=TRUE)) ##data frame aleatoio

Estadística descriptiva

La estadística descriptiva es el área de la estadística que trata de la organización, resumen, análisis e interpretación de datos mediante un análisis descriptivo.

Las variables pueden ser cuantitativas o cualitativas.

# Datos de ejemplo
notas <- c(3.5, 4.2, 3.8, 4.5, 3.2, 4.0, 3.9, 4.1, 3.7, 4.3)

# Cálculo de medidas
media <- mean(notas)
mediana <- median(notas)
desviacion <- sd(notas)
varianza <- var(notas)
cuartiles <- quantile(notas)

# Mostrar resultados
cat("Media:", media, "\n")
## Media: 3.92
cat("Mediana:", mediana, "\n")
## Mediana: 3.95
cat("Desviación estándar:", desviacion, "\n")
## Desviación estándar: 0.388158
cat("Varianza:", varianza, "\n")
## Varianza: 0.1506667
print("Cuartiles:")
## [1] "Cuartiles:"
print(cuartiles)
##    0%   25%   50%   75%  100% 
## 3.200 3.725 3.950 4.175 4.500

Para hacer cualquier cálculo de las medidas de tendencias centra o medidas de posición, se debe tener en cuenta cuándo la data teine valores perdidos y manejar acorde a la situación.

# Crear dataset con valores NA
set.seed(123)
datos_na <- data.frame(
  id = 1:10,
  edad = c(20, NA, 25, 23, NA, 30, 22, 19, NA, 24),
  ingresos = c(1500, 1800, NA, 2200, 1900, NA, NA, 2100, 2300, 2400),
  educacion = c("Primaria", "Bachillerato", NA, "Universidad", "Técnico", 
                "Universidad", "Bachillerato", NA, "Maestría", "Doctorado")
)

# Identificar NAs por columna
sapply(datos_na, function(x) sum(is.na(x)))
##        id      edad  ingresos educacion 
##         0         3         3         2
# Opción 1: Eliminar filas con NAs
datos_sin_na <- na.omit(datos_na)

# Opción 2: Imputar con media (solo variables numéricas)
datos_imputados <- datos_na %>%
  mutate(
    edad = ifelse(is.na(edad), mean(edad, na.rm = TRUE), edad),
    ingresos = ifelse(is.na(ingresos), median(ingresos, na.rm = TRUE), ingresos)
  )

# Opción 3: Imputar con moda (variables categóricas)
get_mode <- function(v) {
  uniqv <- unique(v)
  uniqv[which.max(tabulate(match(v, uniqv)))]
}

datos_imputados <- datos_imputados %>%
  mutate(
    educacion = ifelse(is.na(educacion), get_mode(educacion), educacion)
  )

Visualización de datos. Tipos de gráficos:

#plot: diagrama de dispersion
altura <- c(160, 165, 170, 175, 180)
peso <- c(55, 60, 65, 72, 80)
plot(altura, peso,
     main = "Altura vs Peso",
     xlab = "Altura (cm)",
     ylab = "Peso (kg)",
     col = "blue",
     pch = 20)# el peso depende de la altura altura: imdependinte, peso: dependiente

#grafico de lineas, añadiendo type="o"
tiempo<-1:5
plot(tiempo, peso, type="o",
     main = "Altura vs Peso",
     xlab = "tiempo",
     ylab = "Peso (kg)",
     col = "blue",
     pch = 20)

#grafico de barras, muestra la frecuencia de los datos
frutas <- c("Manzana", "Banana", "Pera", "Uva")
ventas <- c(30, 45, 25, 50)
barplot(ventas, names.arg = frutas,
        main = "Ventas por tipo de fruta",
        col = "orange",
        ylab = "Cantidad")

#histograma : hist, mustra la distribucion de los datos de forma agrupada
edades <- c(18, 22, 21, 25, 20, 22, 24, 26, 28, 22, 23)
hist(edades,
     main = "Distribución de edades",
     xlab = "Edad",
     col = "skyblue",
     border = "black")

#Diagrama de bigotes, describ el comportamiento de los datos
edades <- c(18, 22, 21, 25, 20, 22, 24, 26, 28, 22, 23, 40)
boxplot(edades,
        main = "Boxplot de edades",
        ylab = "Edad",
        col = "pink")

#grafico de pastel / grafico de sectores
categorias <- c("Aprobados", "Reprobados", "Retirados")
valores <- c(50, 30, 20)
pie(valores, labels = categorias,
    main = "Resultados de los estudiantes",
    col = c("green", "red", "yellow"))

#barras apiladas
ventas <- matrix(c(10, 15, 20, 8, 12, 18), nrow = 2, byrow = TRUE)
colnames(ventas) <- c("Ene", "Feb", "Mar")
rownames(ventas) <- c("Producto A", "Producto B")
barplot(ventas,
        beside = FALSE,
        main = "Ventas mensuales por producto",
        col = c("blue", "red"),
        legend = rownames(ventas))

#barras lado a lado
ventas <- matrix(c(10, 15, 20, 8, 12, 18), nrow = 2, byrow = TRUE)
colnames(ventas) <- c("Ene", "Feb", "Mar")
rownames(ventas) <- c("Producto A", "Producto B")
barplot(ventas,
        beside = TRUE,
        main = "Ventas mensuales por producto",
        col = c("blue", "red"),
        legend = rownames(ventas))

#ggplot2
#dispersion
datos <- data.frame(
  altura = c(160, 165, 170, 175, 180),
  peso = c(55, 60, 65, 72, 80)
)

ggplot(datos, aes(x = altura, y = peso))+
  geom_point(color = "blue") +
  ggtitle("Altura vs Peso") +
  xlab("Altura (cm)") +
  ylab("Peso (kg)")

# Cargar datos de ejemplo (mtcars)
data(mtcars)

# Resumen estadístico
summary(mtcars)
##       mpg             cyl             disp             hp       
##  Min.   :10.40   Min.   :4.000   Min.   : 71.1   Min.   : 52.0  
##  1st Qu.:15.43   1st Qu.:4.000   1st Qu.:120.8   1st Qu.: 96.5  
##  Median :19.20   Median :6.000   Median :196.3   Median :123.0  
##  Mean   :20.09   Mean   :6.188   Mean   :230.7   Mean   :146.7  
##  3rd Qu.:22.80   3rd Qu.:8.000   3rd Qu.:326.0   3rd Qu.:180.0  
##  Max.   :33.90   Max.   :8.000   Max.   :472.0   Max.   :335.0  
##       drat             wt             qsec             vs        
##  Min.   :2.760   Min.   :1.513   Min.   :14.50   Min.   :0.0000  
##  1st Qu.:3.080   1st Qu.:2.581   1st Qu.:16.89   1st Qu.:0.0000  
##  Median :3.695   Median :3.325   Median :17.71   Median :0.0000  
##  Mean   :3.597   Mean   :3.217   Mean   :17.85   Mean   :0.4375  
##  3rd Qu.:3.920   3rd Qu.:3.610   3rd Qu.:18.90   3rd Qu.:1.0000  
##  Max.   :4.930   Max.   :5.424   Max.   :22.90   Max.   :1.0000  
##        am              gear            carb      
##  Min.   :0.0000   Min.   :3.000   Min.   :1.000  
##  1st Qu.:0.0000   1st Qu.:3.000   1st Qu.:2.000  
##  Median :0.0000   Median :4.000   Median :2.000  
##  Mean   :0.4062   Mean   :3.688   Mean   :2.812  
##  3rd Qu.:1.0000   3rd Qu.:4.000   3rd Qu.:4.000  
##  Max.   :1.0000   Max.   :5.000   Max.   :8.000
# Análisis por grupos
mtcars %>%
  group_by(cyl) %>%
  summarise(
    media_mpg = mean(mpg),
    mediana_mpg = median(mpg),
    sd_mpg = sd(mpg),
    n = n()
  )
## # A tibble: 3 × 5
##     cyl media_mpg mediana_mpg sd_mpg     n
##   <dbl>     <dbl>       <dbl>  <dbl> <int>
## 1     4      26.7        26     4.51    11
## 2     6      19.7        19.7   1.45     7
## 3     8      15.1        15.2   2.56    14
# Visualización
ggplot(mtcars, aes(x = factor(cyl), y = mpg)) +
  geom_boxplot(fill = "lightcoral") +
  labs(title = "Consumo de combustible por cilindrada",
       x = "Número de cilindros",
       y = "Millas por galón") +
  theme_minimal()

Por último, podemos tener en cuenta los siguientes códigos si queremos estudiar una data de algún lugar en específico:

#archivos de excel library(readxl)
DATA <- read_excel(“ejemp.xlsx”) DATA

#archivos csv library(readr) dataset <- read_csv(NULL) View(dataset)

#para enlaces .txt url <- ‘https://raw.githubusercontent.com/fhernanb/datos/master/babies.txt’ dt_w <- read.table(url, header=TRUE, sep=’).

Estadística inferencial

La correlación mide si dos variables numéricas están relacionadas. Puede ser positiva (ambas aumentan juntas) o negativa (una sube cuando la otra baja). Pero ojo: que estén correlacionadas no significa que una cause la otra. Usamos coeficientes como Pearson o Spearman para medir esta relación, donde valores cercanos a cero indican poca conexión. En la vida real casi nunca hay correlaciones perfectas. Lo clave es recordar que correlación ≠ causalidad.

  1. Se hizo un estudio de automoviles de diferentes marcas y modelos que están dados por la siguiente tabla:
library(readxl)

Automoviles <- read_excel("C:\\Users\\Ryuma\\OneDrive\\Escritorio\\Base de datos automoviles.xlsx")

print(Automoviles)
## # A tibble: 2,151 × 14
##      No. Marca Modelo Consumo Cilindraje `Tipo de alimentación` Potencia
##    <dbl> <dbl>  <dbl>   <dbl>      <dbl>                  <dbl>    <dbl>
##  1     1     6     19    57         1368                      3     99  
##  2     2     2      7    50         1496                      1    109. 
##  3     3     6     21    50         1999                      3    155  
##  4     4     6     19    57         1368                      3     99  
##  5     5     1      1    59         1000                      2    115  
##  6     6     1      3    63         1200                      3     80.5
##  7     7     6     21    48         1999                      3    155  
##  8     8     5     17    63.5       1598                      3    106  
##  9     9     3     10    55         1598                      3    111  
## 10    10     4     16    59.1       1591                      3    122  
## # ℹ 2,141 more rows
## # ℹ 7 more variables: `Tipo de caja` <dbl>, Peso <dbl>, Torque <dbl>,
## #   Válvulas <dbl>, `Tipo de combustible` <dbl>, `Tipos de dirección` <dbl>,
## #   Satisfacción <dbl>

Un comerciante se pregunta si el cilindraje tiene algo que ver con el peso del automovil y así saber qué componentes manipular para ahorrar costos de pruducción y maximizar la eficiencia.

Para esto, primero vamos a manipular las variables que queremos analizar:

Automoviles_imputados <- Automoviles %>%
  mutate(
    Cilindraje = ifelse(is.na(Cilindraje), mean(Cilindraje, na.rm = TRUE), Cilindraje),
    Peso = ifelse(is.na(Peso), median(Peso, na.rm = TRUE), Peso)
  )

Cilindrajes <- Automoviles_imputados$Cilindraje

Pesos <- Automoviles_imputados$Peso


C_vs_P <- data.frame(Cilindrajes, Pesos)

print(C_vs_P)
##      Cilindrajes  Pesos
## 1       1368.000 1092.0
## 2       1496.000 1087.0
## 3       1999.000 1460.0
## 4       1368.000 1062.0
## 5       1000.000 1122.0
## 6       1200.000  940.0
## 7       1999.000 1460.0
## 8       1598.000  959.7
## 9       1598.000 1070.0
## 10      1591.000 1035.0
## 11      1496.000 1087.0
## 12      1368.000 1092.0
## 13      1598.000 1072.0
## 14      1368.000 1062.0
## 15      1200.000  940.0
## 16      1200.000  940.0
## 17      1248.000  912.0
## 18      1000.000 1122.0
## 19      1998.000 1370.0
## 20      1368.000 1092.0
## 21      1598.000 1100.0
## 22      1248.000  912.0
## 23       998.000  885.0
## 24      1248.000  894.0
## 25      1998.000 1386.0
## 26      1000.000 1122.0
## 27      1998.000 1330.0
## 28      1368.000 1092.0
## 29      1598.000  974.9
## 30       999.000  785.0
## 31      1197.000 1003.0
## 32      1198.000  975.0
## 33      1496.000 1060.0
## 34      1248.000  951.0
## 35      1000.000 1122.0
## 36      1399.000 1336.0
## 37      1399.000 1336.0
## 38      1248.000  932.0
## 39      1999.000 1460.0
## 40      1496.000 1082.0
## 41      1998.000 1398.0
## 42      1998.000 1386.0
## 43      1598.000  945.0
## 44      1796.000 1345.0
## 45      1998.000 1432.0
## 46      1368.000 1062.0
## 47      1000.000 1122.0
## 48      1248.000  951.0
## 49       998.000  885.0
## 50      1998.000 1320.0
## 51      1600.000 1109.0
## 52      1998.000 1386.0
## 53      1999.000 1439.0
## 54      1000.000 1122.0
## 55      1200.000  940.0
## 56      1591.000 1035.0
## 57      1368.000 1092.0
## 58      1200.000  940.0
## 59      1368.000 1062.0
## 60      1496.000 1082.0
## 61      1398.000 1014.0
## 62      1998.000 1432.0
## 63      1368.000 1015.0
## 64      1598.000 1070.0
## 65      1998.000 1320.0
## 66      1368.000 1092.0
## 67      1600.000 1109.0
## 68      1197.000 1003.0
## 69      1600.000 1109.0
## 70      1598.000 1072.0
## 71      1999.000 1439.0
## 72      1248.000  912.0
## 73      1197.000 1003.0
## 74      1368.000 1092.0
## 75      1198.000  975.0
## 76      1999.000 1460.0
## 77      1591.000 1035.0
## 78       998.000  885.0
## 79      1200.000  940.0
## 80      1000.000 1122.0
## 81      1368.000 1092.0
## 82      1000.000 1122.0
## 83      1200.000  940.0
## 84      1197.000 1003.0
## 85       998.000  885.0
## 86      1598.000  974.9
## 87      1598.000 1100.0
## 88       998.000  885.0
## 89      1598.000  974.9
## 90      1398.000 1012.0
## 91      1600.000 1109.0
## 92      1198.000  975.0
## 93      1000.000 1122.0
## 94      1248.000  913.0
## 95      1999.000 1439.0
## 96      1496.000 1060.0
## 97      1496.000 1082.0
## 98      1248.000  894.0
## 99      1368.000 1062.0
## 100     1598.000 1082.0
## 101     1598.000 1082.0
## 102     1598.000  959.7
## 103     1496.000 1087.0
## 104     1398.000 1014.0
## 105      999.000  776.0
## 106      998.000  885.0
## 107     1200.000  940.0
## 108     1598.000  945.0
## 109     1598.000 1100.0
## 110     1368.000 1092.0
## 111     1368.000 1062.0
## 112     1600.000 1109.0
## 113     1591.000 1035.0
## 114     1999.000 1439.0
## 115     1200.000  940.0
## 116     1000.000 1122.0
## 117     1000.000 1122.0
## 118     1248.000  894.0
## 119     1198.000  975.0
## 120     1000.000 1122.0
## 121     1999.000 1460.0
## 122     1598.000 1072.0
## 123     1496.000 1082.0
## 124      998.000  885.0
## 125     1598.000 1071.0
## 126     1368.000 1062.0
## 127     1368.000 1092.0
## 128     1598.000 1082.0
## 129     1598.000 1082.0
## 130     1368.000 1015.0
## 131     1600.000 1090.0
## 132     1000.000 1122.0
## 133     1368.000 1062.0
## 134     1600.000 1109.0
## 135     1248.000  894.0
## 136     1598.000 1100.0
## 137      998.000  885.0
## 138     1600.000 1109.0
## 139     1998.000 1320.0
## 140     1248.000  912.0
## 141     1999.000 1460.0
## 142     1248.000  912.0
## 143     1597.000 1169.0
## 144     1200.000  940.0
## 145     1398.000 1012.0
## 146     1598.000  974.9
## 147      998.000  885.0
## 148     1598.000  959.7
## 149     1600.000 1109.0
## 150     1597.000 1169.0
## 151     1999.000 1448.0
## 152     1597.000 1169.0
## 153     1598.000  974.9
## 154     1600.000 1109.0
## 155     1368.000 1092.0
## 156     1200.000  940.0
## 157     1248.000  912.0
## 158     1597.000 1169.0
## 159     1598.000 1078.0
## 160     1999.000 1460.0
## 161     1399.000 1336.0
## 162     1598.000 1072.0
## 163     1591.000 1035.0
## 164     1248.000  932.0
## 165      999.000  776.0
## 166      999.000  760.0
## 167      998.000  885.0
## 168     1598.000  959.7
## 169     1998.000 1386.0
## 170     1368.000 1062.0
## 171     1200.000  940.0
## 172     1197.000 1003.0
## 173     1796.000 1345.0
## 174     1796.000 1345.0
## 175     1198.000  975.0
## 176     1368.000 1035.0
## 177     1598.000 1070.0
## 178     1368.000 1035.0
## 179     1368.000 1092.0
## 180     1597.000 1169.0
## 181     1591.000 1035.0
## 182     1999.000 1465.0
## 183     1598.000  959.7
## 184     1200.000  940.0
## 185     1368.000 1062.0
## 186      998.000  885.0
## 187     1000.000 1122.0
## 188     1597.000 1169.0
## 189     1200.000  940.0
## 190     1000.000 1122.0
## 191     1000.000 1122.0
## 192     1368.000 1092.0
## 193     1496.000 1087.0
## 194     1000.000 1122.0
## 195     1000.000 1122.0
## 196      999.000  760.0
## 197     1248.000  912.0
## 198     1248.000  912.0
## 199     1999.000 1460.0
## 200     1398.000 1014.0
## 201     1368.000 1035.0
## 202     1198.000  975.0
## 203     1200.000  940.0
## 204     1598.000 1072.0
## 205      998.000  885.0
## 206     1368.000 1092.0
## 207     1598.000 1082.0
## 208     1368.000 1092.0
## 209     1796.000 1345.0
## 210     1248.000  932.0
## 211     1198.000  975.0
## 212     1368.000 1092.0
## 213     1368.000 1062.0
## 214     1496.000 1087.0
## 215     1200.000  940.0
## 216     1368.000 1015.0
## 217     1999.000 1439.0
## 218     1368.000 1062.0
## 219     1496.000 1087.0
## 220     1598.000 1071.0
## 221     1591.000 1035.0
## 222     1600.000 1109.0
## 223     1399.000 1290.0
## 224     1999.000 1439.0
## 225      998.000  885.0
## 226     1999.000 1460.0
## 227     1598.000 1082.0
## 228     1598.000 1071.0
## 229     1591.000 1035.0
## 230     1198.000  975.0
## 231     1999.000 1439.0
## 232     1999.000 1439.0
## 233     1998.000 1386.0
## 234     1998.000 1386.0
## 235      999.000  760.0
## 236     1200.000  940.0
## 237      998.000  885.0
## 238     1998.000 1398.0
## 239     1597.000 1169.0
## 240     1598.000 1094.0
## 241     1600.000 1109.0
## 242     1600.000 1090.0
## 243     1597.000 1169.0
## 244     1598.000 1078.0
## 245     1398.000 1014.0
## 246     1248.000  932.0
## 247      999.000  785.0
## 248     1399.000 1336.0
## 249      999.000  760.0
## 250     1600.000 1109.0
## 251     1368.000 1062.0
## 252     1368.000 1035.0
## 253     1000.000 1122.0
## 254     1998.000 1320.0
## 255     1398.000 1012.0
## 256      999.000  785.0
## 257     1368.000 1062.0
## 258     1368.000 1062.0
## 259     1368.000 1092.0
## 260     1248.000  932.0
## 261     1598.000  959.7
## 262     1248.000  894.0
## 263     1368.000 1062.0
## 264      998.000  885.0
## 265     1598.000  959.7
## 266     1998.000 1386.0
## 267      998.000  885.0
## 268     1398.000 1012.0
## 269     1398.000 1013.0
## 270     1368.000 1092.0
## 271      999.000  776.0
## 272     1600.000 1090.0
## 273     1000.000 1122.0
## 274     1197.000 1003.0
## 275     1496.000 1090.0
## 276     1598.000  959.7
## 277     1200.000  940.0
## 278     1368.000 1015.0
## 279     1999.000 1439.0
## 280     1198.000  975.0
## 281     1368.000 1062.0
## 282     1368.000 1062.0
## 283     1198.000  975.0
## 284     1597.000 1169.0
## 285     1248.000  912.0
## 286     1198.000  975.0
## 287     1999.000 1460.0
## 288     1200.000  940.0
## 289     1598.000 1100.0
## 290     1368.000 1062.0
## 291     1197.000 1003.0
## 292     1999.000 1439.0
## 293      998.000  885.0
## 294     1598.000 1072.0
## 295     1248.000  951.0
## 296     1368.000 1015.0
## 297     1248.000  913.0
## 298     1000.000 1122.0
## 299     1998.000 1432.0
## 300     1598.000 1072.0
## 301     1999.000 1448.0
## 302     1998.000 1432.0
## 303     1496.000 1082.0
## 304     1368.000 1035.0
## 305     1200.000  940.0
## 306     1200.000  940.0
## 307     1368.000 1062.0
## 308      999.000  776.0
## 309     1000.000 1122.0
## 310     1000.000 1122.0
## 311     1496.000 1082.0
## 312     1598.000  959.7
## 313     1591.000 1035.0
## 314     1200.000  940.0
## 315     1368.000 1062.0
## 316     1598.000 1100.0
## 317     1248.000  912.0
## 318     1368.000 1092.0
## 319      999.000  760.0
## 320     1398.000 1012.0
## 321     1600.000 1109.0
## 322     1398.000 1012.0
## 323     1598.000 1094.0
## 324     1200.000  940.0
## 325     1597.000 1169.0
## 326     1597.000 1169.0
## 327     1796.000 1345.0
## 328     1999.000 1460.0
## 329     1998.000 1432.0
## 330      999.000  785.0
## 331     1597.000 1169.0
## 332     1796.000 1345.0
## 333      998.000  885.0
## 334     1368.000 1092.0
## 335     1598.000 1078.0
## 336      999.000  760.0
## 337      998.000  885.0
## 338     1999.000 1448.0
## 339     1600.000 1090.0
## 340      999.000  760.0
## 341     1368.000 1092.0
## 342     1198.000  975.0
## 343     1368.000 1035.0
## 344     1368.000 1062.0
## 345     1600.000 1109.0
## 346     1796.000 1345.0
## 347     1368.000 1092.0
## 348     1200.000  940.0
## 349     1000.000 1122.0
## 350     1248.000  912.0
## 351     1248.000  932.0
## 352     1600.000 1109.0
## 353     1496.000 1087.0
## 354      999.000  760.0
## 355     1598.000  959.7
## 356      998.000  885.0
## 357      998.000  885.0
## 358     1999.000 1439.0
## 359     1598.000 1072.0
## 360      999.000  785.0
## 361     1248.000  912.0
## 362     1368.000 1062.0
## 363     1600.000 1090.0
## 364     1248.000  932.0
## 365     1999.000 1460.0
## 366     1598.000 1078.0
## 367     1368.000 1092.0
## 368     1198.000  975.0
## 369     1999.000 1439.0
## 370     1496.000 1087.0
## 371     1398.000 1013.0
## 372     1000.000 1122.0
## 373     1999.000 1460.0
## 374     1598.000 1100.0
## 375      998.000  885.0
## 376     1998.000 1370.0
## 377     1598.000  945.0
## 378     1600.000 1090.0
## 379     1368.000 1062.0
## 380     1597.000 1169.0
## 381     1598.000 1071.0
## 382     1600.000 1090.0
## 383     1248.000  912.0
## 384     1248.000  932.0
## 385     1598.000 1071.0
## 386     1200.000  940.0
## 387     1598.000 1078.0
## 388     1998.000 1320.0
## 389     1000.000 1122.0
## 390     1000.000 1122.0
## 391     1591.000 1035.0
## 392     1598.000  974.9
## 393     1200.000  940.0
## 394     1399.000 1290.0
## 395     1600.000 1109.0
## 396     1598.000  959.7
## 397      998.000  885.0
## 398     1200.000  940.0
## 399     1598.000 1071.0
## 400     1496.000 1087.0
## 401     1368.000 1092.0
## 402     1600.000 1109.0
## 403     1600.000 1109.0
## 404     1998.000 1398.0
## 405     1368.000 1092.0
## 406     1600.000 1109.0
## 407     1398.000 1014.0
## 408     1399.000 1290.0
## 409     1998.000 1330.0
## 410     1368.000 1092.0
## 411     1200.000  940.0
## 412     1368.000 1035.0
## 413     1597.000 1169.0
## 414     1368.000 1062.0
## 415     1600.000 1090.0
## 416     1398.000 1014.0
## 417     1999.000 1448.0
## 418     1999.000 1460.0
## 419     1999.000 1465.0
## 420     1368.000 1015.0
## 421      999.000  760.0
## 422     1999.000 1439.0
## 423     1368.000 1092.0
## 424     1496.000 1060.0
## 425      999.000  785.0
## 426     1600.000 1109.0
## 427     1998.000 1320.0
## 428     1368.000 1062.0
## 429     1600.000 1090.0
## 430     1600.000 1109.0
## 431     1248.000  912.0
## 432     1197.000 1003.0
## 433     1598.000  974.9
## 434     1398.000 1013.0
## 435     1248.000  912.0
## 436     1198.000  975.0
## 437     1796.000 1345.0
## 438     1398.000 1013.0
## 439      998.000  885.0
## 440     1998.000 1330.0
## 441     1368.000 1035.0
## 442     1998.000 1432.0
## 443     1000.000 1122.0
## 444     1496.000 1087.0
## 445     1368.000 1062.0
## 446     1597.000 1169.0
## 447     1368.000 1035.0
## 448     1598.000  974.9
## 449     1198.000  975.0
## 450     1368.000 1092.0
## 451     1000.000 1122.0
## 452     1368.000 1092.0
## 453     1999.000 1460.0
## 454     1598.000 1100.0
## 455     1591.000 1035.0
## 456     1598.000 1078.0
## 457      999.000  760.0
## 458     1368.000 1092.0
## 459     1368.000 1062.0
## 460     1598.000 1070.0
## 461     1000.000 1122.0
## 462     1368.000 1092.0
## 463     1000.000 1122.0
## 464     1200.000  940.0
## 465     1597.000 1169.0
## 466     1598.000  959.7
## 467     1600.000 1109.0
## 468     1598.000  974.9
## 469     1200.000  940.0
## 470     1999.000 1448.0
## 471      999.000  760.0
## 472     1999.000 1439.0
## 473     1999.000 1460.0
## 474     1000.000 1122.0
## 475     1598.000 1078.0
## 476     1597.000 1169.0
## 477     1999.000 1448.0
## 478     1200.000  940.0
## 479     1597.000 1169.0
## 480     1398.000 1013.0
## 481      999.000  785.0
## 482     1368.000 1062.0
## 483     1999.000 1460.0
## 484     1598.000  959.7
## 485     1598.000 1070.0
## 486     1600.000 1109.0
## 487      998.000  885.0
## 488      998.000  885.0
## 489     1598.000 1078.0
## 490     1600.000 1090.0
## 491     1999.000 1465.0
## 492     1998.000 1320.0
## 493     1398.000 1014.0
## 494     1999.000 1465.0
## 495     1200.000  940.0
## 496     1398.000 1012.0
## 497     1248.000  951.0
## 498     1600.000 1109.0
## 499     1399.000 1290.0
## 500     1598.000  974.9
## 501     1796.000 1345.0
## 502     1496.000 1060.0
## 503     1398.000 1013.0
## 504     1496.000 1082.0
## 505     1398.000 1012.0
## 506      999.000  776.0
## 507     1591.000 1035.0
## 508     1248.000  912.0
## 509     1248.000  951.0
## 510     1248.000  932.0
## 511     1998.000 1320.0
## 512     1591.000 1035.0
## 513     1998.000 1432.0
## 514     1999.000 1439.0
## 515     1368.000 1092.0
## 516     1496.000 1090.0
## 517     1368.000 1035.0
## 518     1496.000 1090.0
## 519     1000.000 1122.0
## 520     1999.000 1439.0
## 521     1399.000 1336.0
## 522     1000.000 1122.0
## 523     1200.000  940.0
## 524     1598.000 1100.0
## 525     1999.000 1448.0
## 526     1197.000 1003.0
## 527     1999.000 1460.0
## 528     1999.000 1465.0
## 529     1198.000  975.0
## 530      999.000  785.0
## 531      999.000  760.0
## 532     1598.000 1072.0
## 533     1999.000 1448.0
## 534     1368.000 1092.0
## 535     1600.000 1090.0
## 536     1600.000 1109.0
## 537     1598.000  974.9
## 538     1600.000 1109.0
## 539     1368.000 1092.0
## 540     1198.000  975.0
## 541     1998.000 1320.0
## 542     1368.000 1062.0
## 543     1000.000 1122.0
## 544     1368.000 1035.0
## 545     1248.000  932.0
## 546      999.000  776.0
## 547     1598.000 1070.0
## 548     1000.000 1122.0
## 549     1598.000 1072.0
## 550     1200.000  940.0
## 551     1248.000  932.0
## 552      999.000  776.0
## 553     1368.000 1062.0
## 554     1598.000 1071.0
## 555     1368.000 1062.0
## 556     1999.000 1439.0
## 557     1198.000  975.0
## 558     1248.000  932.0
## 559     1398.000 1013.0
## 560     1598.000  945.0
## 561     1368.000 1015.0
## 562     1999.000 1465.0
## 563     1600.000 1109.0
## 564     1999.000 1460.0
## 565     1998.000 1432.0
## 566     1598.000  974.9
## 567     1368.000 1092.0
## 568     1600.000 1109.0
## 569     1248.000  913.0
## 570     1368.000 1062.0
## 571     1398.000 1014.0
## 572     1998.000 1330.0
## 573     1368.000 1035.0
## 574     1598.000 1071.0
## 575     1999.000 1439.0
## 576     1598.000  959.7
## 577      998.000  885.0
## 578     1368.000 1062.0
## 579     1248.000  894.0
## 580     1496.000 1090.0
## 581     1600.000 1109.0
## 582     1496.000 1060.0
## 583     1368.000 1062.0
## 584     1368.000 1035.0
## 585     1999.000 1460.0
## 586     1598.000 1100.0
## 587      999.000  785.0
## 588      999.000  760.0
## 589     1000.000 1122.0
## 590     1368.000 1092.0
## 591      999.000  776.0
## 592     1368.000 1092.0
## 593      999.000  776.0
## 594     1999.000 1460.0
## 595     1248.000  894.0
## 596     1000.000 1122.0
## 597     1398.000 1014.0
## 598     1600.000 1109.0
## 599      999.000  785.0
## 600     1496.000 1090.0
## 601     1796.000 1345.0
## 602     1398.000 1014.0
## 603     1198.000  975.0
## 604     1368.000 1062.0
## 605     1248.000  912.0
## 606     1598.000  974.9
## 607     1368.000 1092.0
## 608      999.000  785.0
## 609     1999.000 1460.0
## 610     1198.000  975.0
## 611     1600.000 1109.0
## 612      999.000  760.0
## 613      998.000  885.0
## 614     1200.000  940.0
## 615     1368.000 1035.0
## 616     1368.000 1062.0
## 617     1398.000 1014.0
## 618     1598.000 1070.0
## 619     1999.000 1460.0
## 620     1999.000 1439.0
## 621     1398.000 1012.0
## 622      998.000  885.0
## 623     1368.000 1035.0
## 624     1368.000 1062.0
## 625     1248.000  951.0
## 626     1368.000 1035.0
## 627     1598.000  959.7
## 628      998.000  885.0
## 629     1248.000  912.0
## 630     1598.000 1100.0
## 631     1591.000 1035.0
## 632     1000.000 1122.0
## 633     1999.000 1439.0
## 634     1368.000 1092.0
## 635     1999.000 1460.0
## 636     1368.000 1092.0
## 637     1598.000 1100.0
## 638     1999.000 1439.0
## 639     1600.000 1109.0
## 640     1198.000  975.0
## 641     1998.000 1330.0
## 642     1398.000 1013.0
## 643     1598.000  959.7
## 644     1198.000  975.0
## 645     1591.000 1035.0
## 646     1598.000  974.9
## 647      999.000  776.0
## 648     1368.000 1092.0
## 649     1598.000 1072.0
## 650      998.000  885.0
## 651     1000.000 1122.0
## 652      999.000  785.0
## 653     1598.000  974.9
## 654     1000.000 1122.0
## 655     1398.000 1014.0
## 656     1368.000 1062.0
## 657     1000.000 1122.0
## 658     1598.000 1072.0
## 659     1200.000  940.0
## 660     1598.000  945.0
## 661     1600.000 1090.0
## 662     1368.000 1092.0
## 663     1998.000 1386.0
## 664     1598.000 1070.0
## 665     1000.000 1122.0
## 666     1496.000 1090.0
## 667     1248.000  912.0
## 668     1600.000 1090.0
## 669     1197.000 1003.0
## 670     1198.000  975.0
## 671     1368.000 1092.0
## 672     1200.000  940.0
## 673     1597.000 1169.0
## 674     1998.000 1320.0
## 675     1398.000 1012.0
## 676     1000.000 1122.0
## 677      998.000  885.0
## 678     1248.000  912.0
## 679     1591.000 1035.0
## 680     1398.000 1012.0
## 681     1248.000  912.0
## 682     1600.000 1109.0
## 683     1597.000 1169.0
## 684     1796.000 1345.0
## 685     1999.000 1439.0
## 686     1598.000 1071.0
## 687     1198.000  975.0
## 688     1591.000 1035.0
## 689      999.000  776.0
## 690     1000.000 1122.0
## 691     1999.000 1460.0
## 692     1000.000 1122.0
## 693     1598.000 1082.0
## 694     1999.000 1439.0
## 695     1597.000 1169.0
## 696     1399.000 1290.0
## 697     1000.000 1122.0
## 698     1368.000 1035.0
## 699     1598.000 1071.0
## 700      999.000  785.0
## 701     1999.000 1460.0
## 702     1200.000  940.0
## 703     1998.000 1398.0
## 704     1998.000 1432.0
## 705     1598.000  959.7
## 706     1248.000  912.0
## 707     1368.000 1062.0
## 708     1368.000 1015.0
## 709     1368.000 1092.0
## 710     1999.000 1439.0
## 711     1200.000  940.0
## 712     1368.000 1062.0
## 713     1598.000 1094.0
## 714     1999.000 1460.0
## 715     1398.000 1014.0
## 716     1198.000  975.0
## 717     1598.000  974.9
## 718     1598.000 1072.0
## 719     1598.000  974.9
## 720     1368.000 1092.0
## 721     1248.000  913.0
## 722     1368.000 1062.0
## 723     1496.000 1087.0
## 724     1399.000 1336.0
## 725     1600.000 1090.0
## 726     1998.000 1370.0
## 727      998.000  885.0
## 728     1998.000 1432.0
## 729     1000.000 1122.0
## 730     1598.000 1100.0
## 731     1999.000 1460.0
## 732     1197.000 1003.0
## 733     1198.000  975.0
## 734     1598.000 1071.0
## 735     1598.000 1078.0
## 736     1398.000 1012.0
## 737     1197.000 1003.0
## 738     1197.000 1003.0
## 739      998.000  885.0
## 740     1999.000 1439.0
## 741     1200.000  940.0
## 742     1368.000 1035.0
## 743     1598.000  974.9
## 744     1200.000  940.0
## 745     1198.000  975.0
## 746     1496.000 1082.0
## 747     1368.000 1092.0
## 748      998.000  885.0
## 749     1598.000 1078.0
## 750     1398.000 1014.0
## 751     1598.000  959.7
## 752      999.000  776.0
## 753      998.000  885.0
## 754     1600.000 1109.0
## 755     1496.000 1087.0
## 756     1399.000 1336.0
## 757      999.000  760.0
## 758     1398.000 1014.0
## 759     1598.000  974.9
## 760     1368.000 1062.0
## 761     1368.000 1062.0
## 762     1368.000 1015.0
## 763     1000.000 1122.0
## 764     1496.000 1060.0
## 765     1197.000 1003.0
## 766     1796.000 1345.0
## 767      999.000  785.0
## 768     1496.000 1090.0
## 769     1368.000 1062.0
## 770     1598.000  974.9
## 771     1200.000  940.0
## 772     1600.000 1090.0
## 773     1000.000 1122.0
## 774     1198.000  975.0
## 775     1198.000  975.0
## 776     1368.000 1062.0
## 777     1598.000 1094.0
## 778     1999.000 1439.0
## 779     1999.000 1460.0
## 780      999.000  785.0
## 781     1398.000 1012.0
## 782     1368.000 1015.0
## 783     1200.000  940.0
## 784     1998.000 1432.0
## 785     1598.000 1070.0
## 786     1598.000  945.0
## 787     1399.000 1290.0
## 788     1368.000 1092.0
## 789      998.000  885.0
## 790     1198.000  975.0
## 791     1200.000  940.0
## 792     1368.000 1015.0
## 793     1998.000 1320.0
## 794      999.000  785.0
## 795     1200.000  940.0
## 796      998.000  885.0
## 797     1998.000 1398.0
## 798     1197.000 1003.0
## 799     1398.000 1013.0
## 800     1248.000  912.0
## 801      999.000  785.0
## 802     1248.000  913.0
## 803     1591.000 1035.0
## 804     1000.000 1122.0
## 805     1998.000 1320.0
## 806     1598.000 1070.0
## 807      999.000  785.0
## 808     1600.000 1109.0
## 809     1198.000  975.0
## 810     1368.000 1092.0
## 811     1598.000  959.7
## 812     1598.000 1100.0
## 813     1198.000  975.0
## 814     1496.000 1090.0
## 815     1248.000  951.0
## 816     1368.000 1092.0
## 817     1200.000  940.0
## 818     1998.000 1386.0
## 819     1248.000  894.0
## 820     1398.000 1014.0
## 821     1368.000 1062.0
## 822     1248.000  912.0
## 823     1000.000 1122.0
## 824     1368.000 1062.0
## 825     1368.000 1092.0
## 826     1598.000 1082.0
## 827      998.000  885.0
## 828      998.000  885.0
## 829     1591.000 1035.0
## 830     1248.000  951.0
## 831     1598.000 1072.0
## 832     1597.000 1169.0
## 833     1368.000 1092.0
## 834     1496.000 1060.0
## 835     1999.000 1460.0
## 836     1591.000 1035.0
## 837     1200.000  940.0
## 838     1248.000  951.0
## 839     1999.000 1460.0
## 840     1368.000 1062.0
## 841     1999.000 1465.0
## 842     1496.000 1087.0
## 843      998.000  885.0
## 844     1368.000 1015.0
## 845     1600.000 1109.0
## 846     1200.000  940.0
## 847     1998.000 1398.0
## 848     1999.000 1460.0
## 849     1368.000 1062.0
## 850      999.000  760.0
## 851     1198.000  975.0
## 852     1368.000 1062.0
## 853     1398.000 1014.0
## 854     1398.000 1013.0
## 855     1368.000 1092.0
## 856     1600.000 1090.0
## 857     1598.000 1100.0
## 858     1200.000  940.0
## 859     1598.000  959.7
## 860     1368.000 1062.0
## 861     1999.000 1439.0
## 862     1368.000 1092.0
## 863     1496.000 1060.0
## 864     1200.000  940.0
## 865     1398.000 1012.0
## 866     1496.000 1060.0
## 867     1248.000  951.0
## 868     1600.000 1090.0
## 869     1496.000 1090.0
## 870     1999.000 1439.0
## 871     1248.000  913.0
## 872     1796.000 1345.0
## 873     1200.000  940.0
## 874     1000.000 1122.0
## 875     1198.000  975.0
## 876     1197.000 1003.0
## 877     1496.000 1082.0
## 878     1399.000 1290.0
## 879     1248.000  951.0
## 880     1368.000 1092.0
## 881     1600.000 1109.0
## 882     1368.000 1092.0
## 883      998.000  885.0
## 884     1398.000 1013.0
## 885     1398.000 1012.0
## 886     1200.000  940.0
## 887     1496.000 1082.0
## 888     1368.000 1035.0
## 889     1248.000  912.0
## 890     1600.000 1109.0
## 891     1000.000 1122.0
## 892     1368.000 1092.0
## 893     1998.000 1320.0
## 894     1248.000  912.0
## 895     1597.000 1169.0
## 896     1598.000 1082.0
## 897     1000.000 1122.0
## 898     1198.000  975.0
## 899     1198.000  975.0
## 900      999.000  760.0
## 901     1399.000 1336.0
## 902     1368.000 1092.0
## 903      999.000  776.0
## 904     1198.000  975.0
## 905      999.000  785.0
## 906     1368.000 1092.0
## 907     1496.000 1087.0
## 908     1248.000  913.0
## 909     1999.000 1460.0
## 910     1000.000 1122.0
## 911     1496.000 1087.0
## 912     1999.000 1460.0
## 913     1600.000 1090.0
## 914     1598.000 1072.0
## 915     1598.000 1070.0
## 916     1398.000 1012.0
## 917     1999.000 1460.0
## 918     1597.000 1169.0
## 919     1600.000 1109.0
## 920     1000.000 1122.0
## 921     1197.000 1003.0
## 922     1200.000  940.0
## 923     1248.000  894.0
## 924     1248.000  912.0
## 925      998.000  885.0
## 926     1368.000 1092.0
## 927     1597.000 1169.0
## 928     1600.000 1090.0
## 929     1200.000  940.0
## 930     1398.000 1013.0
## 931     1200.000  940.0
## 932     1496.000 1060.0
## 933     1999.000 1439.0
## 934     1248.000  912.0
## 935      998.000  885.0
## 936     1198.000  975.0
## 937      998.000  885.0
## 938     1999.000 1460.0
## 939     1999.000 1448.0
## 940     1399.000 1290.0
## 941     1496.000 1082.0
## 942     1200.000  940.0
## 943     1598.000 1094.0
## 944     1998.000 1320.0
## 945     1368.000 1015.0
## 946     1598.000 1082.0
## 947     1197.000 1003.0
## 948     1598.000  974.9
## 949     1368.000 1062.0
## 950     1496.000 1090.0
## 951     1600.000 1090.0
## 952     1999.000 1465.0
## 953     1368.000 1062.0
## 954     1998.000 1398.0
## 955     1998.000 1398.0
## 956     1998.000 1370.0
## 957     1796.000 1345.0
## 958     1248.000  951.0
## 959     1000.000 1122.0
## 960     1598.000  959.7
## 961     1600.000 1109.0
## 962     1368.000 1062.0
## 963     1000.000 1122.0
## 964     1368.000 1092.0
## 965     1999.000 1460.0
## 966      999.000  785.0
## 967     1598.000 1078.0
## 968     1999.000 1448.0
## 969     1398.000 1014.0
## 970     1248.000  912.0
## 971     1368.000 1062.0
## 972     1200.000  940.0
## 973     1999.000 1460.0
## 974     1197.000 1003.0
## 975     1600.000 1090.0
## 976     1998.000 1398.0
## 977     1598.000  974.9
## 978     1197.000 1003.0
## 979     1000.000 1122.0
## 980     1999.000 1448.0
## 981     1368.000 1035.0
## 982     1598.000 1070.0
## 983     1600.000 1109.0
## 984     1999.000 1448.0
## 985     1999.000 1439.0
## 986     1600.000 1090.0
## 987     1368.000 1092.0
## 988     1598.000 1094.0
## 989     1796.000 1345.0
## 990      999.000  785.0
## 991     1368.000 1092.0
## 992     1598.000 1071.0
## 993     1200.000  940.0
## 994     1368.000 1092.0
## 995     1368.000 1092.0
## 996      999.000  760.0
## 997     1368.000 1062.0
## 998     1198.000  975.0
## 999     1197.000 1003.0
## 1000    1000.000 1122.0
## 1001    1398.000 1012.0
## 1002    1000.000 1122.0
## 1003    1000.000 1122.0
## 1004    1598.000 1070.0
## 1005    1598.000  959.7
## 1006    1368.000 1062.0
## 1007     998.000  885.0
## 1008    1198.000  975.0
## 1009    1796.000 1345.0
## 1010    1398.000 1014.0
## 1011    1999.000 1460.0
## 1012    1368.000 1015.0
## 1013    1999.000 1460.0
## 1014    1000.000 1122.0
## 1015    1999.000 1460.0
## 1016    1200.000  940.0
## 1017     999.000  776.0
## 1018    1368.000 1062.0
## 1019    1368.000 1092.0
## 1020    1200.000  940.0
## 1021     999.000  760.0
## 1022    1591.000 1035.0
## 1023    1368.000 1015.0
## 1024    1399.000 1336.0
## 1025    1496.000 1090.0
## 1026    1368.000 1092.0
## 1027    1998.000 1370.0
## 1028    1197.000 1003.0
## 1029     999.000  760.0
## 1030    1200.000  940.0
## 1031    1248.000  912.0
## 1032     999.000  760.0
## 1033    1248.000  951.0
## 1034    1598.000  959.7
## 1035    1000.000 1122.0
## 1036    1999.000 1448.0
## 1037    1999.000 1465.0
## 1038    1248.000  913.0
## 1039    1000.000 1122.0
## 1040     999.000  776.0
## 1041    1197.000 1003.0
## 1042    1998.000 1398.0
## 1043    1399.000 1336.0
## 1044    1398.000 1013.0
## 1045    1368.000 1035.0
## 1046    1398.000 1012.0
## 1047    1398.000 1012.0
## 1048    1197.000 1003.0
## 1049    1248.000  932.0
## 1050    1999.000 1448.0
## 1051    1000.000 1122.0
## 1052    1600.000 1109.0
## 1053     998.000  885.0
## 1054    1398.000 1014.0
## 1055    1200.000  940.0
## 1056    1999.000 1460.0
## 1057    1598.000  945.0
## 1058    1368.000 1062.0
## 1059    1200.000  940.0
## 1060     998.000  885.0
## 1061    1598.000  974.9
## 1062    1200.000  940.0
## 1063    1591.000 1035.0
## 1064    1999.000 1460.0
## 1065    1998.000 1398.0
## 1066    1248.000  913.0
## 1067    1600.000 1109.0
## 1068     999.000  785.0
## 1069    1200.000  940.0
## 1070    1598.000  945.0
## 1071    1368.000 1062.0
## 1072    1598.000 1072.0
## 1073    1398.000 1013.0
## 1074    1598.000 1082.0
## 1075    1248.000  951.0
## 1076    1200.000  940.0
## 1077    1368.000 1062.0
## 1078    1600.000 1109.0
## 1079    1597.000 1169.0
## 1080    1200.000  940.0
## 1081    1598.000  959.7
## 1082    1600.000 1109.0
## 1083    1398.000 1014.0
## 1084    1796.000 1345.0
## 1085    1998.000 1398.0
## 1086    1998.000 1330.0
## 1087    1200.000  940.0
## 1088    1248.000  913.0
## 1089    1999.000 1439.0
## 1090    1368.000 1092.0
## 1091    1998.000 1370.0
## 1092    1998.000 1330.0
## 1093    1398.000 1013.0
## 1094    1000.000 1122.0
## 1095    1598.000 1072.0
## 1096    1197.000 1003.0
## 1097    1591.000 1035.0
## 1098    1198.000  975.0
## 1099    1598.000 1070.0
## 1100    1198.000  975.0
## 1101    1999.000 1448.0
## 1102    1198.000  975.0
## 1103    1368.000 1062.0
## 1104    1591.000 1035.0
## 1105    1198.000  975.0
## 1106    1399.000 1290.0
## 1107    1368.000 1062.0
## 1108    1368.000 1062.0
## 1109    1598.000 1071.0
## 1110    1368.000 1062.0
## 1111    1368.000 1062.0
## 1112    1998.000 1432.0
## 1113    1600.000 1109.0
## 1114     999.000  760.0
## 1115    1998.000 1370.0
## 1116    1399.000 1290.0
## 1117    1796.000 1345.0
## 1118    1496.000 1090.0
## 1119    1999.000 1460.0
## 1120    1999.000 1460.0
## 1121    1796.000 1345.0
## 1122     998.000  885.0
## 1123    1999.000 1460.0
## 1124    1000.000 1122.0
## 1125     999.000  776.0
## 1126    1368.000 1062.0
## 1127    1598.000 1082.0
## 1128    1368.000 1062.0
## 1129    1398.000 1013.0
## 1130    1598.000  974.9
## 1131    1398.000 1012.0
## 1132    1597.000 1169.0
## 1133    1248.000  912.0
## 1134    1496.000 1082.0
## 1135    1598.000  974.9
## 1136    1591.000 1035.0
## 1137    1398.000 1013.0
## 1138    1496.000 1082.0
## 1139    1999.000 1460.0
## 1140    1368.000 1092.0
## 1141    1496.000 1082.0
## 1142    1368.000 1092.0
## 1143    1597.000 1169.0
## 1144    1999.000 1439.0
## 1145     999.000  785.0
## 1146    1796.000 1345.0
## 1147    1398.000 1012.0
## 1148    1796.000 1345.0
## 1149    1998.000 1386.0
## 1150    1600.000 1109.0
## 1151    1398.000 1014.0
## 1152    1200.000  940.0
## 1153    1200.000  940.0
## 1154    1368.000 1092.0
## 1155    1000.000 1122.0
## 1156    1248.000  951.0
## 1157    1999.000 1439.0
## 1158    1368.000 1092.0
## 1159    1368.000 1092.0
## 1160    1398.000 1013.0
## 1161    1597.000 1169.0
## 1162    1600.000 1109.0
## 1163    1999.000 1460.0
## 1164    1598.000  959.7
## 1165    1999.000 1460.0
## 1166    1368.000 1092.0
## 1167    1398.000 1013.0
## 1168    1368.000 1062.0
## 1169    1399.000 1290.0
## 1170    1597.000 1169.0
## 1171    1999.000 1460.0
## 1172    1000.000 1122.0
## 1173    1598.000  959.7
## 1174    1999.000 1448.0
## 1175    1598.000 1072.0
## 1176    1368.000 1062.0
## 1177     998.000  885.0
## 1178    1796.000 1345.0
## 1179    1998.000 1432.0
## 1180    1368.000 1035.0
## 1181    1496.000 1090.0
## 1182    1368.000 1015.0
## 1183    1591.000 1035.0
## 1184    1999.000 1460.0
## 1185    1398.000 1012.0
## 1186    1200.000  940.0
## 1187    1368.000 1015.0
## 1188    1598.000  945.0
## 1189    1368.000 1092.0
## 1190    1591.000 1035.0
## 1191    1398.000 1012.0
## 1192    1368.000 1062.0
## 1193    1198.000  975.0
## 1194    1598.000  974.9
## 1195    1600.000 1109.0
## 1196     998.000  885.0
## 1197     998.000  885.0
## 1198    1000.000 1122.0
## 1199    1598.000  974.9
## 1200    1998.000 1330.0
## 1201    1600.000 1109.0
## 1202    1198.000  975.0
## 1203    1597.000 1169.0
## 1204    1998.000 1386.0
## 1205    1998.000 1370.0
## 1206    1600.000 1090.0
## 1207     999.000  776.0
## 1208    1598.000  945.0
## 1209    1600.000 1090.0
## 1210    1598.000  959.7
## 1211    1368.000 1092.0
## 1212    1999.000 1439.0
## 1213    1000.000 1122.0
## 1214    1198.000  975.0
## 1215    1998.000 1330.0
## 1216    1598.000 1071.0
## 1217    1368.000 1092.0
## 1218     999.000  760.0
## 1219    1398.000 1013.0
## 1220    1398.000 1012.0
## 1221    1591.000 1035.0
## 1222    1368.000 1062.0
## 1223    1368.000 1015.0
## 1224    1600.000 1109.0
## 1225    1598.000 1100.0
## 1226    1998.000 1320.0
## 1227    1368.000 1062.0
## 1228    1999.000 1439.0
## 1229    1598.000 1072.0
## 1230    1368.000 1062.0
## 1231    1496.000 1090.0
## 1232    1398.000 1012.0
## 1233    1000.000 1122.0
## 1234    1598.000 1082.0
## 1235    1398.000 1012.0
## 1236    1000.000 1122.0
## 1237     998.000  885.0
## 1238    1200.000  940.0
## 1239    1368.000 1062.0
## 1240    1597.000 1169.0
## 1241    1368.000 1062.0
## 1242    1598.000 1082.0
## 1243    1600.000 1109.0
## 1244    1600.000 1109.0
## 1245    1200.000  940.0
## 1246    1368.000 1035.0
## 1247    1998.000 1386.0
## 1248     999.000  776.0
## 1249    1597.000 1169.0
## 1250    1998.000 1320.0
## 1251     998.000  885.0
## 1252    1248.000  894.0
## 1253    1368.000 1062.0
## 1254    1796.000 1345.0
## 1255    1998.000 1432.0
## 1256    1368.000 1092.0
## 1257    1600.000 1109.0
## 1258    1368.000 1062.0
## 1259    1999.000 1460.0
## 1260    1200.000  940.0
## 1261    1597.000 1169.0
## 1262     999.000  776.0
## 1263    1000.000 1122.0
## 1264    1999.000 1439.0
## 1265    1368.000 1062.0
## 1266    1198.000  975.0
## 1267    1368.000 1062.0
## 1268    1591.000 1035.0
## 1269    1597.000 1169.0
## 1270     998.000  885.0
## 1271    1368.000 1062.0
## 1272    1597.000 1169.0
## 1273    1998.000 1386.0
## 1274    1496.000 1082.0
## 1275     998.000  885.0
## 1276    1368.000 1062.0
## 1277    1999.000 1460.0
## 1278    1198.000  975.0
## 1279    1998.000 1370.0
## 1280    1368.000 1062.0
## 1281    1598.000 1082.0
## 1282    1200.000  940.0
## 1283    1598.000  974.9
## 1284    1496.000 1060.0
## 1285     999.000  785.0
## 1286    1000.000 1122.0
## 1287    1999.000 1460.0
## 1288    1998.000 1330.0
## 1289     999.000  776.0
## 1290    1368.000 1062.0
## 1291    1598.000 1072.0
## 1292    1598.000 1078.0
## 1293    1200.000  940.0
## 1294    1600.000 1109.0
## 1295    1368.000 1092.0
## 1296    1598.000 1078.0
## 1297     999.000  760.0
## 1298    1598.000  974.9
## 1299    1398.000 1012.0
## 1300     999.000  776.0
## 1301    1598.000  959.7
## 1302    1796.000 1345.0
## 1303    1598.000 1072.0
## 1304    1999.000 1460.0
## 1305     998.000  885.0
## 1306    1368.000 1062.0
## 1307    1198.000  975.0
## 1308    1998.000 1398.0
## 1309    1591.000 1035.0
## 1310    1598.000 1082.0
## 1311    1198.000  975.0
## 1312    1998.000 1398.0
## 1313    1598.000 1078.0
## 1314    1398.000 1013.0
## 1315    1598.000  974.9
## 1316    1598.000  974.9
## 1317    1368.000 1015.0
## 1318    1999.000 1460.0
## 1319    1368.000 1062.0
## 1320     999.000  785.0
## 1321    1597.000 1169.0
## 1322    1368.000 1035.0
## 1323    1998.000 1370.0
## 1324    1598.000  974.9
## 1325    1399.000 1336.0
## 1326    1368.000 1062.0
## 1327    1591.000 1035.0
## 1328    1398.000 1012.0
## 1329    1998.000 1370.0
## 1330    1200.000  940.0
## 1331    1496.000 1082.0
## 1332    1598.000  945.0
## 1333    1999.000 1439.0
## 1334    1000.000 1122.0
## 1335    1248.000  913.0
## 1336    1198.000  975.0
## 1337    1999.000 1460.0
## 1338    1399.000 1290.0
## 1339    1200.000  940.0
## 1340    1200.000  940.0
## 1341    1368.000 1092.0
## 1342     999.000  776.0
## 1343    1600.000 1109.0
## 1344    1598.000  974.9
## 1345    1598.000 1072.0
## 1346     999.000  760.0
## 1347    1368.000 1062.0
## 1348    1368.000 1062.0
## 1349     998.000  885.0
## 1350    1598.000 1071.0
## 1351     999.000  776.0
## 1352    1248.000  951.0
## 1353    1597.000 1169.0
## 1354    1000.000 1122.0
## 1355     998.000  885.0
## 1356    1598.000 1071.0
## 1357    1598.000  959.7
## 1358    1198.000  975.0
## 1359    1998.000 1432.0
## 1360    1398.000 1013.0
## 1361    1000.000 1122.0
## 1362    1398.000 1012.0
## 1363    1248.000  894.0
## 1364    1368.000 1062.0
## 1365    1999.000 1460.0
## 1366    1600.000 1109.0
## 1367    1999.000 1439.0
## 1368    1368.000 1092.0
## 1369    1198.000  975.0
## 1370    1999.000 1460.0
## 1371    1999.000 1460.0
## 1372    1999.000 1460.0
## 1373    1600.000 1109.0
## 1374    1796.000 1345.0
## 1375    1999.000 1460.0
## 1376    1368.000 1015.0
## 1377    1248.000  913.0
## 1378    1398.000 1014.0
## 1379    1598.000 1070.0
## 1380    1398.000 1013.0
## 1381    1597.000 1169.0
## 1382    1368.000 1062.0
## 1383    1496.000 1087.0
## 1384    1368.000 1092.0
## 1385    1796.000 1345.0
## 1386    1999.000 1448.0
## 1387    1198.000  975.0
## 1388    1198.000  975.0
## 1389    1999.000 1460.0
## 1390    1600.000 1109.0
## 1391    1496.000 1087.0
## 1392     999.000  760.0
## 1393    1248.000  913.0
## 1394    1197.000 1003.0
## 1395    1600.000 1090.0
## 1396    1368.000 1062.0
## 1397    1597.000 1169.0
## 1398    1198.000  975.0
## 1399    1999.000 1448.0
## 1400    1000.000 1122.0
## 1401    1598.000 1078.0
## 1402    1398.000 1014.0
## 1403    1368.000 1092.0
## 1404    1597.000 1169.0
## 1405    1000.000 1122.0
## 1406    1597.000 1169.0
## 1407    1598.000 1100.0
## 1408    1999.000 1460.0
## 1409    1368.000 1062.0
## 1410    1496.000 1060.0
## 1411    1598.000 1078.0
## 1412     998.000  885.0
## 1413     998.000  885.0
## 1414    1368.000 1092.0
## 1415    1248.000  951.0
## 1416    1200.000  940.0
## 1417    1248.000  951.0
## 1418    1248.000  912.0
## 1419    1998.000 1330.0
## 1420    1200.000  940.0
## 1421     998.000  885.0
## 1422     999.000  776.0
## 1423    1999.000 1439.0
## 1424    1200.000  940.0
## 1425    1000.000 1122.0
## 1426    1198.000  975.0
## 1427    1598.000  974.9
## 1428     999.000  776.0
## 1429    1600.000 1090.0
## 1430    1598.000  959.7
## 1431     999.000  785.0
## 1432    1598.000  945.0
## 1433    1198.000  975.0
## 1434    1368.000 1062.0
## 1435    1591.000 1035.0
## 1436    1198.000  975.0
## 1437    1598.000  945.0
## 1438    1999.000 1439.0
## 1439    1368.000 1092.0
## 1440    1368.000 1062.0
## 1441     998.000  885.0
## 1442     999.000  760.0
## 1443    1198.000  975.0
## 1444    1368.000 1092.0
## 1445     999.000  785.0
## 1446    1198.000  975.0
## 1447    1000.000 1122.0
## 1448    1598.000 1100.0
## 1449    1200.000  940.0
## 1450    1597.000 1169.0
## 1451     999.000  785.0
## 1452    1398.000 1013.0
## 1453    1200.000  940.0
## 1454    1598.000 1078.0
## 1455    1197.000 1003.0
## 1456    1398.000 1012.0
## 1457    1200.000  940.0
## 1458    1200.000  940.0
## 1459    1598.000 1070.0
## 1460    1598.000  974.9
## 1461    1998.000 1432.0
## 1462    1368.000 1062.0
## 1463    1999.000 1460.0
## 1464    1197.000 1003.0
## 1465    1998.000 1320.0
## 1466    1248.000  894.0
## 1467    1368.000 1062.0
## 1468    1598.000 1094.0
## 1469    1000.000 1122.0
## 1470    1368.000 1015.0
## 1471    1597.000 1169.0
## 1472    1198.000  975.0
## 1473    1999.000 1460.0
## 1474    1999.000 1465.0
## 1475    1600.000 1109.0
## 1476    1398.000 1013.0
## 1477     999.000  760.0
## 1478    1368.000 1092.0
## 1479    1368.000 1015.0
## 1480    1368.000 1062.0
## 1481    1200.000  940.0
## 1482    1368.000 1092.0
## 1483    1598.000 1072.0
## 1484    1796.000 1345.0
## 1485    1368.000 1062.0
## 1486    1796.000 1345.0
## 1487    1598.000 1100.0
## 1488    1000.000 1122.0
## 1489    1598.000 1094.0
## 1490    1597.000 1169.0
## 1491    1999.000 1460.0
## 1492     998.000  885.0
## 1493     998.000  885.0
## 1494    1496.000 1060.0
## 1495    1368.000 1035.0
## 1496    1598.000  974.9
## 1497    1000.000 1122.0
## 1498    1248.000  932.0
## 1499    1200.000  940.0
## 1500    1796.000 1345.0
## 1501    1598.000 1070.0
## 1502    1796.000 1345.0
## 1503    1496.000 1087.0
## 1504    1600.000 1109.0
## 1505     999.000  776.0
## 1506    1999.000 1439.0
## 1507    1598.000  959.7
## 1508    1200.000  940.0
## 1509    1591.000 1035.0
## 1510    1248.000  912.0
## 1511     999.000  760.0
## 1512    1600.000 1109.0
## 1513    1200.000  940.0
## 1514    1248.000  912.0
## 1515    1248.000  912.0
## 1516    1999.000 1465.0
## 1517    1598.000 1100.0
## 1518    1398.000 1012.0
## 1519    1368.000 1062.0
## 1520    1796.000 1345.0
## 1521     999.000  785.0
## 1522    1598.000 1100.0
## 1523    1999.000 1460.0
## 1524    1998.000 1330.0
## 1525    1999.000 1460.0
## 1526    1999.000 1448.0
## 1527     999.000  760.0
## 1528    1248.000  912.0
## 1529    1796.000 1345.0
## 1530     998.000  885.0
## 1531    1999.000 1460.0
## 1532    1368.000 1035.0
## 1533    1598.000  974.9
## 1534     999.000  776.0
## 1535    1598.000  945.0
## 1536    1999.000 1465.0
## 1537    1998.000 1398.0
## 1538     999.000  776.0
## 1539    1198.000  975.0
## 1540    1999.000 1460.0
## 1541    1496.000 1090.0
## 1542    1496.000 1087.0
## 1543    1496.000 1082.0
## 1544    1598.000  974.9
## 1545    1368.000 1062.0
## 1546    1496.000 1087.0
## 1547    1999.000 1460.0
## 1548    1600.000 1109.0
## 1549    1600.000 1109.0
## 1550    1248.000  912.0
## 1551    1368.000 1062.0
## 1552    1398.000 1013.0
## 1553    1999.000 1448.0
## 1554    1197.000 1003.0
## 1555     998.000  885.0
## 1556    1368.000 1062.0
## 1557    1796.000 1345.0
## 1558    1200.000  940.0
## 1559    1496.000 1060.0
## 1560    1398.000 1014.0
## 1561    1999.000 1439.0
## 1562    1999.000 1448.0
## 1563    1398.000 1012.0
## 1564    1200.000  940.0
## 1565    1248.000  912.0
## 1566    1999.000 1439.0
## 1567    1368.000 1092.0
## 1568    1597.000 1169.0
## 1569    1597.000 1169.0
## 1570    1398.000 1014.0
## 1571    1248.000  932.0
## 1572    1000.000 1122.0
## 1573    1496.000 1060.0
## 1574    1598.000  959.7
## 1575    1368.000 1062.0
## 1576    1600.000 1090.0
## 1577     999.000  760.0
## 1578    1368.000 1092.0
## 1579    1998.000 1370.0
## 1580    1398.000 1013.0
## 1581    1597.000 1169.0
## 1582    1398.000 1012.0
## 1583    1999.000 1439.0
## 1584    1598.000  959.7
## 1585    1398.000 1013.0
## 1586    1496.000 1060.0
## 1587    1398.000 1012.0
## 1588    1600.000 1109.0
## 1589    1398.000 1012.0
## 1590    1398.000 1013.0
## 1591    1198.000  975.0
## 1592     999.000  776.0
## 1593     999.000  760.0
## 1594    1496.000 1082.0
## 1595    1198.000  975.0
## 1596    1200.000  940.0
## 1597    1368.000 1092.0
## 1598    1796.000 1345.0
## 1599    1597.000 1169.0
## 1600    1200.000  940.0
## 1601    1598.000 1071.0
## 1602    1248.000  913.0
## 1603    1998.000 1386.0
## 1604    1999.000 1465.0
## 1605    1368.000 1092.0
## 1606    1368.000 1092.0
## 1607     999.000  776.0
## 1608     999.000  760.0
## 1609    1598.000 1094.0
## 1610    1197.000 1003.0
## 1611    1398.000 1013.0
## 1612    1999.000 1448.0
## 1613    1000.000 1122.0
## 1614    1598.000 1078.0
## 1615     999.000  760.0
## 1616    1200.000  940.0
## 1617     999.000  785.0
## 1618    1496.000 1090.0
## 1619    1496.000 1090.0
## 1620    1368.000 1092.0
## 1621    1999.000 1460.0
## 1622    1398.000 1013.0
## 1623    1200.000  940.0
## 1624    1999.000 1439.0
## 1625    1398.000 1014.0
## 1626    1000.000 1122.0
## 1627    1591.000 1035.0
## 1628    1000.000 1122.0
## 1629    1198.000  975.0
## 1630    1597.000 1169.0
## 1631    1597.000 1169.0
## 1632    1591.000 1035.0
## 1633    1000.000 1122.0
## 1634    1398.000 1012.0
## 1635    1200.000  940.0
## 1636    1598.000  959.7
## 1637    1368.000 1062.0
## 1638    1198.000  975.0
## 1639     998.000  885.0
## 1640    1598.000 1082.0
## 1641    1598.000 1078.0
## 1642    1597.000 1169.0
## 1643    1368.000 1092.0
## 1644    1000.000 1122.0
## 1645    1496.000 1082.0
## 1646     998.000  885.0
## 1647    1998.000 1432.0
## 1648    1197.000 1003.0
## 1649    1999.000 1448.0
## 1650    1399.000 1336.0
## 1651    1368.000 1092.0
## 1652    1248.000  894.0
## 1653    1597.000 1169.0
## 1654    1496.000 1090.0
## 1655    1200.000  940.0
## 1656    1198.000  975.0
## 1657    1248.000  932.0
## 1658    1597.000 1169.0
## 1659    1598.000 1078.0
## 1660    1248.000  894.0
## 1661    1998.000 1370.0
## 1662    1496.000 1082.0
## 1663    1198.000  975.0
## 1664    1368.000 1062.0
## 1665     999.000  785.0
## 1666    1000.000 1122.0
## 1667    1999.000 1460.0
## 1668    1198.000  975.0
## 1669    1197.000 1003.0
## 1670     998.000  885.0
## 1671    1197.000 1003.0
## 1672    1598.000 1070.0
## 1673    1368.000 1062.0
## 1674    1600.000 1109.0
## 1675     998.000  885.0
## 1676    1998.000 1330.0
## 1677    1496.000 1090.0
## 1678    1198.000  975.0
## 1679    1796.000 1345.0
## 1680    1598.000 1072.0
## 1681    1000.000 1122.0
## 1682    1998.000 1386.0
## 1683     999.000  776.0
## 1684    1198.000  975.0
## 1685    1248.000  951.0
## 1686    1200.000  940.0
## 1687    1999.000 1460.0
## 1688    1368.000 1092.0
## 1689    1399.000 1336.0
## 1690    1200.000  940.0
## 1691    1999.000 1448.0
## 1692    1496.000 1087.0
## 1693    1368.000 1092.0
## 1694    1496.000 1082.0
## 1695     999.000  785.0
## 1696    1600.000 1109.0
## 1697    1597.000 1169.0
## 1698    1200.000  940.0
## 1699    1248.000  894.0
## 1700    1398.000 1014.0
## 1701     999.000  776.0
## 1702    1398.000 1013.0
## 1703    1000.000 1122.0
## 1704    1598.000 1071.0
## 1705    1998.000 1432.0
## 1706    1000.000 1122.0
## 1707    1600.000 1090.0
## 1708    1999.000 1460.0
## 1709     998.000  885.0
## 1710    1598.000 1072.0
## 1711    1598.000  959.7
## 1712    1998.000 1330.0
## 1713    1998.000 1398.0
## 1714    1368.000 1062.0
## 1715    1998.000 1370.0
## 1716    1598.000 1071.0
## 1717    1248.000  912.0
## 1718    1598.000 1094.0
## 1719    1248.000  894.0
## 1720    1398.000 1013.0
## 1721    1000.000 1122.0
## 1722    1000.000 1122.0
## 1723    1999.000 1460.0
## 1724    1000.000 1122.0
## 1725    1398.000 1012.0
## 1726    1496.000 1060.0
## 1727    1998.000 1398.0
## 1728    1796.000 1345.0
## 1729    1368.000 1035.0
## 1730    1000.000 1122.0
## 1731    1368.000 1092.0
## 1732    1598.000 1071.0
## 1733    1999.000 1439.0
## 1734    1368.000 1092.0
## 1735     998.000  885.0
## 1736    1248.000  894.0
## 1737    1248.000  894.0
## 1738    1600.000 1109.0
## 1739     999.000  785.0
## 1740    1591.000 1035.0
## 1741    1598.000 1070.0
## 1742    1598.000  974.9
## 1743     999.000  760.0
## 1744    1368.000 1092.0
## 1745    1999.000 1460.0
## 1746    1000.000 1122.0
## 1747    1368.000 1092.0
## 1748    1999.000 1439.0
## 1749    1999.000 1460.0
## 1750    1198.000  975.0
## 1751    1398.000 1014.0
## 1752     999.000  760.0
## 1753    1000.000 1122.0
## 1754    1248.000  932.0
## 1755    1248.000  912.0
## 1756    1368.000 1062.0
## 1757    1399.000 1336.0
## 1758    1200.000  940.0
## 1759    1200.000  940.0
## 1760    1600.000 1109.0
## 1761    1368.000 1092.0
## 1762    1000.000 1122.0
## 1763    1398.000 1012.0
## 1764    1998.000 1370.0
## 1765    1398.000 1012.0
## 1766    1368.000 1015.0
## 1767    1368.000 1092.0
## 1768     999.000  776.0
## 1769    1598.000 1078.0
## 1770     999.000  776.0
## 1771     998.000  885.0
## 1772    1598.000 1072.0
## 1773    1200.000  940.0
## 1774    1597.000 1169.0
## 1775    1598.000 1071.0
## 1776    1598.000  974.9
## 1777    1198.000  975.0
## 1778    1598.000 1100.0
## 1779    1598.000 1094.0
## 1780    1197.000 1003.0
## 1781    1598.000 1072.0
## 1782    1198.000  975.0
## 1783    1999.000 1460.0
## 1784    1600.000 1109.0
## 1785    1398.000 1014.0
## 1786    1496.000 1087.0
## 1787    1248.000  913.0
## 1788    1496.000 1082.0
## 1789    1999.000 1460.0
## 1790    1600.000 1109.0
## 1791    1248.000  912.0
## 1792    1368.000 1092.0
## 1793    1496.000 1090.0
## 1794    1597.000 1169.0
## 1795    1200.000  940.0
## 1796    1368.000 1062.0
## 1797    1368.000 1092.0
## 1798    1248.000  932.0
## 1799    1998.000 1370.0
## 1800    1591.000 1035.0
## 1801     998.000  885.0
## 1802    1368.000 1062.0
## 1803     999.000  760.0
## 1804    1998.000 1432.0
## 1805    1598.000 1082.0
## 1806     998.000  885.0
## 1807    1000.000 1122.0
## 1808    1598.000  945.0
## 1809     998.000  885.0
## 1810    1598.000  959.7
## 1811    1200.000  940.0
## 1812    1398.000 1014.0
## 1813     999.000  760.0
## 1814    1598.000 1070.0
## 1815    1398.000 1014.0
## 1816    1600.000 1090.0
## 1817    1198.000  975.0
## 1818     998.000  885.0
## 1819    1200.000  940.0
## 1820    1200.000  940.0
## 1821    1248.000  932.0
## 1822    1600.000 1090.0
## 1823    1600.000 1109.0
## 1824    1000.000 1122.0
## 1825    1598.000  974.9
## 1826    1398.000 1013.0
## 1827    1398.000 1014.0
## 1828    1598.000 1082.0
## 1829    1496.000 1090.0
## 1830    1999.000 1465.0
## 1831    1998.000 1386.0
## 1832    1368.000 1062.0
## 1833    1399.000 1290.0
## 1834    1999.000 1460.0
## 1835     998.000  885.0
## 1836    1198.000  975.0
## 1837    1598.000 1082.0
## 1838    1398.000 1012.0
## 1839     999.000  785.0
## 1840     999.000  785.0
## 1841    1368.000 1015.0
## 1842    1598.000  974.9
## 1843    1999.000 1439.0
## 1844    1998.000 1330.0
## 1845    1000.000 1122.0
## 1846    1598.000 1070.0
## 1847    1496.000 1082.0
## 1848    1598.000 1100.0
## 1849    1000.000 1122.0
## 1850    1000.000 1122.0
## 1851    1600.000 1109.0
## 1852    1598.000  974.9
## 1853     999.000  776.0
## 1854    1598.000 1070.0
## 1855    1368.000 1035.0
## 1856     998.000  885.0
## 1857    1368.000 1092.0
## 1858    1198.000  975.0
## 1859     999.000  760.0
## 1860    1398.000 1012.0
## 1861    1999.000 1460.0
## 1862     998.000  885.0
## 1863     999.000  776.0
## 1864    1398.000 1014.0
## 1865    1398.000 1014.0
## 1866    1998.000 1398.0
## 1867    1197.000 1003.0
## 1868     999.000  760.0
## 1869    1598.000 1082.0
## 1870    1200.000  940.0
## 1871     998.000  885.0
## 1872    1998.000 1398.0
## 1873    1999.000 1460.0
## 1874    1999.000 1460.0
## 1875    1591.000 1035.0
## 1876    1998.000 1320.0
## 1877    1598.000 1094.0
## 1878    1197.000 1003.0
## 1879    1368.000 1092.0
## 1880    1200.000  940.0
## 1881    1000.000 1122.0
## 1882    1598.000 1100.0
## 1883     999.000  785.0
## 1884    1399.000 1336.0
## 1885    1999.000 1439.0
## 1886    1999.000 1460.0
## 1887    1398.000 1012.0
## 1888    1600.000 1109.0
## 1889    1999.000 1448.0
## 1890    1368.000 1062.0
## 1891    1200.000  940.0
## 1892    1496.000 1087.0
## 1893    1598.000 1072.0
## 1894    1000.000 1122.0
## 1895    1200.000  940.0
## 1896    1598.000 1071.0
## 1897    1368.000 1035.0
## 1898    1200.000  940.0
## 1899    1496.000 1060.0
## 1900    1398.000 1014.0
## 1901     998.000  885.0
## 1902    1600.000 1109.0
## 1903    1368.000 1035.0
## 1904     999.000  785.0
## 1905    1248.000  894.0
## 1906    1598.000 1094.0
## 1907    1600.000 1109.0
## 1908     999.000  760.0
## 1909    1248.000  894.0
## 1910    1398.000 1013.0
## 1911    1398.000 1013.0
## 1912    1598.000  959.7
## 1913    1597.000 1169.0
## 1914    1368.000 1062.0
## 1915    1000.000 1122.0
## 1916    1248.000  913.0
## 1917    1998.000 1320.0
## 1918    1496.000 1087.0
## 1919    1600.000 1109.0
## 1920     998.000  885.0
## 1921    1496.000 1060.0
## 1922    1000.000 1122.0
## 1923    1368.000 1015.0
## 1924    1248.000  951.0
## 1925     998.000  885.0
## 1926    1398.000 1014.0
## 1927    1200.000  940.0
## 1928    1998.000 1432.0
## 1929    1796.000 1345.0
## 1930    1999.000 1460.0
## 1931     998.000  885.0
## 1932    1399.000 1290.0
## 1933    1600.000 1109.0
## 1934    1600.000 1090.0
## 1935    1398.000 1014.0
## 1936    1600.000 1109.0
## 1937    1999.000 1460.0
## 1938    1598.000 1071.0
## 1939    1368.000 1035.0
## 1940    1597.000 1169.0
## 1941    1398.000 1014.0
## 1942    1197.000 1003.0
## 1943    1368.000 1062.0
## 1944    1000.000 1122.0
## 1945    1248.000  913.0
## 1946     998.000  885.0
## 1947    1598.000 1072.0
## 1948    1598.000  945.0
## 1949    1368.000 1062.0
## 1950    1598.000 1094.0
## 1951    1368.000 1062.0
## 1952    1197.000 1003.0
## 1953    1399.000 1336.0
## 1954    1368.000 1035.0
## 1955     998.000  885.0
## 1956    1368.000 1035.0
## 1957    1198.000  975.0
## 1958    1591.000 1035.0
## 1959    1248.000  912.0
## 1960    1000.000 1122.0
## 1961    1398.000 1012.0
## 1962    1998.000 1320.0
## 1963    1368.000 1092.0
## 1964    1368.000 1062.0
## 1965    1398.000 1014.0
## 1966    1998.000 1330.0
## 1967    1598.000  959.7
## 1968    1398.000 1012.0
## 1969    1368.000 1092.0
## 1970     999.000  760.0
## 1971    1597.000 1169.0
## 1972    1496.000 1090.0
## 1973     998.000  885.0
## 1974    1000.000 1122.0
## 1975    1368.000 1092.0
## 1976    1197.000 1003.0
## 1977    1399.000 1336.0
## 1978    1597.000 1169.0
## 1979    1998.000 1386.0
## 1980     999.000  785.0
## 1981    1598.000  974.9
## 1982    1248.000  951.0
## 1983    1600.000 1109.0
## 1984    1398.000 1014.0
## 1985    1496.000 1087.0
## 1986    1197.000 1003.0
## 1987    1368.000 1062.0
## 1988    1998.000 1386.0
## 1989    1000.000 1122.0
## 1990    1999.000 1439.0
## 1991     999.000  776.0
## 1992    1999.000 1448.0
## 1993    1496.000 1060.0
## 1994    1998.000 1370.0
## 1995    1998.000 1386.0
## 1996    1998.000 1330.0
## 1997    1368.000 1035.0
## 1998    1368.000 1062.0
## 1999    1796.000 1345.0
## 2000    1368.000 1092.0
## 2001    1368.000 1092.0
## 2002    1999.000 1465.0
## 2003    1368.000 1015.0
## 2004    1597.000 1169.0
## 2005    1999.000 1460.0
## 2006    1200.000  940.0
## 2007    1597.000 1169.0
## 2008    1200.000  940.0
## 2009    1000.000 1122.0
## 2010    1368.000 1062.0
## 2011    1598.000 1078.0
## 2012    1600.000 1109.0
## 2013    1796.000 1345.0
## 2014    1200.000  940.0
## 2015    1368.000 1062.0
## 2016    1200.000  940.0
## 2017    1598.000  974.9
## 2018    1368.000 1092.0
## 2019    1398.000 1013.0
## 2020    1248.000  894.0
## 2021    1999.000 1460.0
## 2022    1998.000 1432.0
## 2023    1598.000  945.0
## 2024    1200.000  940.0
## 2025    1796.000 1345.0
## 2026    1398.000 1012.0
## 2027    1399.000 1290.0
## 2028    1598.000 1071.0
## 2029    1398.000 1012.0
## 2030    1796.000 1345.0
## 2031    1200.000  940.0
## 2032    1368.000 1062.0
## 2033    1368.000 1092.0
## 2034    1000.000 1122.0
## 2035    1496.000 1090.0
## 2036    1398.000 1013.0
## 2037    1248.000  951.0
## 2038    1598.000  945.0
## 2039    1999.000 1460.0
## 2040    1999.000 1460.0
## 2041    1591.000 1035.0
## 2042    1598.000 1100.0
## 2043    1368.000 1015.0
## 2044    1998.000 1370.0
## 2045    1000.000 1122.0
## 2046    1368.000 1092.0
## 2047    1998.000 1432.0
## 2048    1999.000 1465.0
## 2049    1598.000 1082.0
## 2050    1999.000 1439.0
## 2051    1200.000  940.0
## 2052    1598.000  959.7
## 2053    1998.000 1330.0
## 2054    1999.000 1439.0
## 2055    1398.000 1013.0
## 2056    1598.000 1082.0
## 2057    1368.000 1062.0
## 2058    1198.000  975.0
## 2059    1999.000 1460.0
## 2060    1999.000 1460.0
## 2061    1598.000 1082.0
## 2062    1368.000 1062.0
## 2063    1368.000 1092.0
## 2064    1598.000 1070.0
## 2065    1248.000  912.0
## 2066    1998.000 1320.0
## 2067    1591.000 1035.0
## 2068    1000.000 1122.0
## 2069    1368.000 1035.0
## 2070    1368.000 1092.0
## 2071    1368.000 1035.0
## 2072    1368.000 1035.0
## 2073    1198.000  975.0
## 2074    1197.000 1003.0
## 2075     999.000  760.0
## 2076     999.000  776.0
## 2077    1398.000 1013.0
## 2078     998.000  885.0
## 2079    1248.000  912.0
## 2080    1368.000 1092.0
## 2081    1597.000 1169.0
## 2082    1999.000 1460.0
## 2083    1200.000  940.0
## 2084    1000.000 1122.0
## 2085    1998.000 1370.0
## 2086    1200.000  940.0
## 2087    1398.000 1014.0
## 2088    1368.000 1062.0
## 2089    1398.000 1013.0
## 2090     999.000  785.0
## 2091    1368.000 1062.0
## 2092    1368.000 1035.0
## 2093    1998.000 1320.0
## 2094    1200.000  940.0
## 2095    1496.000 1060.0
## 2096    1600.000 1109.0
## 2097     999.000  785.0
## 2098    1496.000 1060.0
## 2099    1398.000 1013.0
## 2100    1496.000 1060.0
## 2101     998.000  885.0
## 2102    1600.000 1109.0
## 2103    1000.000 1122.0
## 2104     999.000  785.0
## 2105    1999.000 1460.0
## 2106    1368.000 1015.0
## 2107    1368.000 1062.0
## 2108     998.000  885.0
## 2109    1368.000 1092.0
## 2110    1496.000 1090.0
## 2111    1200.000  940.0
## 2112    1368.000 1062.0
## 2113    1496.000 1060.0
## 2114     998.000  885.0
## 2115    1248.000  932.0
## 2116    1368.000 1092.0
## 2117     999.000  760.0
## 2118    1200.000  940.0
## 2119    1368.000 1062.0
## 2120    1999.000 1460.0
## 2121    1197.000 1003.0
## 2122    1200.000  940.0
## 2123    1998.000 1398.0
## 2124    1999.000 1465.0
## 2125    1200.000  940.0
## 2126     998.000  885.0
## 2127    1999.000 1439.0
## 2128    1600.000 1109.0
## 2129    1198.000  975.0
## 2130    1248.000  932.0
## 2131     999.000  760.0
## 2132    1200.000  940.0
## 2133    1200.000  940.0
## 2134     998.000  885.0
## 2135    1398.000 1014.0
## 2136    1598.000  974.9
## 2137    1591.000 1035.0
## 2138    1999.000 1460.0
## 2139    1000.000 1122.0
## 2140     999.000  785.0
## 2141     998.000  885.0
## 2142    1200.000  940.0
## 2143     999.000  785.0
## 2144    1398.000 1014.0
## 2145    1368.000 1062.0
## 2146    1368.000 1092.0
## 2147    1000.000 1122.0
## 2148    1796.000 1345.0
## 2149    1598.000 1070.0
## 2150    1368.000 1062.0
## 2151    1441.753 1062.0

Ahora, nos interesa saber si estos datos tienen un comportamiento normal o no, para esto primero vamos a ver el gráfico de dispersión y analizar si pueden llegar a ser datos con distribución normal y luego harémos la respectiva prueba. Luego, verémos si estos tienen correlación:

\[H_0= \text{El peso no tiene correlación con el cilindraje } p>0.05\\ H_A= \text{El peso sí tiene una correlación on el cilindraje } p<0.05\]

plot(Cilindrajes, Pesos)

ks.test(Cilindrajes, "pnorm") #Se hace esta prueba ya que son muchos datos y, no parece que tenga una distribución normal.
## 
##  Asymptotic one-sample Kolmogorov-Smirnov test
## 
## data:  Cilindrajes
## D = 1, p-value < 2.2e-16
## alternative hypothesis: two-sided
ks.test(Pesos,"pnorm")
## 
##  Asymptotic one-sample Kolmogorov-Smirnov test
## 
## data:  Pesos
## D = 1, p-value < 2.2e-16
## alternative hypothesis: two-sided
cor.test(Cilindrajes, Pesos, method = "spearman")
## 
##  Spearman's rank correlation rho
## 
## data:  Cilindrajes and Pesos
## S = 474934983, p-value < 2.2e-16
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.7136717

Como podemos observar, efectivamente los datos no presentan una distribución normal, pero, por otro lado si presentan una correlación fuerte positiva. Con esto, podemos construir un modelo de regresión Theil-Sen.

library(mblm)

Regresión_TS <- mblm(Pesos~Cilindrajes) #Primero se coloca la variable dependiente y luego la independiente

summary(Regresión_TS)
## 
## Call:
## mblm(formula = Pesos ~ Cilindrajes)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -183.75  -34.76    4.34   60.42  291.32 
## 
## Coefficients:
##             Estimate      MAD V value Pr(>|V|)    
## (Intercept) 453.6000 146.1495 2314476   <2e-16 ***
## Cilindrajes   0.4225   0.1031 2293555   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 102.7 on 2149 degrees of freedom

Donde Estimate es el valor esperado cuando el cilindraje = 0 y, por cada aumento de una unidad de cilindraje, el peso aumenta en 0.4225. En otras palabras, nuestro modelo sería: 0.4225x + 453.6000

Para mayor visualización, veamoslo en un gráfico:

plot(Cilindrajes, Pesos, main = "Regresión de Theil-Sen")
abline(Regresión_TS, col = "red", lwd = 2)

Claramente existen más pruebas estadísticas dependiendo si son variables cualitativas, cuantitativas, con distribución normal o sin distribución normal, se debe ver cuál se usa según la situación.