DATA 609: Final Presentation

Dan Brooks, Keith Folsom
May 18, 2017

1989: The Midge Classification Problem (Keith)

Two species of midges, Af and Apf, have been identified by biologists Grogan and Wirth (1981) on the basis of antenna and wing length (Figure A.8). Each of nine Af midges is denoted by a circle, and each of six Apf midges is denoted by a square. It is important to be able to classify a specimen as Af or Apf, given the antenna and wing length

  • The Apf midge is a carrier of a debilitating disease that causes swelling of the brain when a human is bitten by an infected midge
  • The Af midge is quite harmless and a valuable pollinator

Midge Data

The midge data provided by species appears as shown below:

Species Ant.Length Wing.Length
Af 1.38 1.64
Af 1.40 1.70
Af 1.24 1.72
Af 1.36 1.74
Af 1.38 1.82
Af 1.48 1.82
Af 1.54 1.82
Af 1.38 1.90
Af 1.56 2.08
Apf 1.14 1.78
Apf 1.20 1.86
Apf 1.18 1.96
Apf 1.30 1.96
Apf 1.26 2.00
Apf 1.28 2.00

Apf & Af Midge Scatterplot

plot of chunk unnamed-chunk-3

Given a midge that you know is species Af or Apf, how would you go about classifying it?

  • Leverage the existence of the pointless region between the species to establish boundaries.

plot of chunk unnamed-chunk-4

Problem Statement

Apply your method to three specimens with (antenna, wing) lengths (1.24, 1.80), (1.28, 1.84), (1.40, 2.04)

Problem : The three specimens fall within the pointless region between the two midge borders.

plot of chunk unnamed-chunk-5

Assume that species Af is a valuable pollinator and that species Apf is a carrier of a debilitating disease. Would you modify your classification scheme and if so, how?

  • Yes, the classification scheme needs to be modified to classifiy the species falling in between the existing Apf and Af boundaries.

Classification Method 1

  • Consider the ratio of midge antenna to wing length as a means of classification

plot of chunk unnamed-chunk-6

Classification Method 1 - Continued

Using the mean ratio for each species, we arrive at the following equations:

\[ \text{Af Midge: }Antenna.Length = 0.7846034 * Wing.Length \]

\[ \text{Apf Midge:} Antenna.Length = 0.6368195 * Wing.Length \]

Classification Method 1 - Continued

Superimposing on top of the midge graph:

plot of chunk unnamed-chunk-7

Using the mean of the ratios between the two species, we establish a border between the two groups:

plot of chunk unnamed-chunk-8

Classification Method 1 - Classification Scheme:

\[ \text{if Antenna.Length > 0.7107115 * Wing.Length then Af Midge} \]

\[ \text{if Antenna.Length < 0.7107115 * Wing.Length then Apf Midge} \]

Classification Method 2 - Linear Regression

  • Create linear regression models for the Af and Apf species.

plot of chunk unnamed-chunk-9

Applying predicted values based on the model with an 85% confidence level:

plot of chunk unnamed-chunk-10

The predicted lower and upper bands for each species do not cover the unclassified midge at (1.8, 1.24).

Classification Method 2 - Linear Regression Continued

Given this missed midge, use a similar approach using the mid-way points between the two boundaries established by the linear models:

\[ \text{Af Midge:} Antenna.Length = 0.4789 * Wing.Length + 0.5491 \] \[ \text{Apf Midge:} Antenna.Length = 0.5582 * Wing.Length + 0.1512 \]

plot of chunk unnamed-chunk-11

Classification Method 2 - Classification Scheme:

  • \[ \text{if Antenna.Length > 0.5185763 * Wing.Length + 0.3501444 then Af Midge} \]

  • \[ \text{if Antenna.Length < 0.5185763 * Wing.Length + 0.3501444 then Apf Midge} \]

1986: The Emergency-Facilities Location Problem (Dan)

The township of Rio Rancho has hitherto not had its own emergency facilities. It has secured funds to erect two emergency facilities in 1986, each of which will combine ambulance, fire, and police services. Figure A.4 indicates the demand, or number of emergencies per square block, for 1985. The L region in th north is an obstacle, whereas the rectangle in the south is a park with a shallow pond. It takes an emergency vehicle an average of 15 seconds to go one block in the N-S direction and 20 seconds in the E-W direction. Your task is to locate the two facilities so as to minimize the total response time.

Assumptions

  • Assume that the demand is concentrated at the center of the block and that the facilities will be located on corners.

  • Assume that the demand is uniformly distributed on the streets bordering each block and that the facilities may be located anywhere on the streets.

Approach One

  • Euclidean Distance

  • Treat all points as those they can be reached by all of the other points

  • Treat the barries as blocks with no incidencies

  • Number of incidences gets multiplied by the N/S weight and E/W weight to get the total distance for traveling to that location.

  • All distances and added up acr0ss the board and the two smallest distances are where the stations should be placed.

  • Below we can see the matrix of incidencies

  • Below that we see the distances. Block (5,3) and (5,4) have the smallest distances, so the stations should be placed there.

Ranchero
      [,1] [,2] [,3] [,4] [,5]
 [1,]    3    1    4    2    5
 [2,]    3    2    3    3    2
 [3,]    2    0    3    3    2
 [4,]    3    0    0    3    1
 [5,]    3    4    3    3    5
 [6,]    2    3    4    4    0
 [7,]    1    2    0    1    3
 [8,]    0    2    0    3    2
 [9,]    3    0    0    0    4
[10,]    3    1    0    4    2
print("/n")
[1] "/n"
raw.distance
       [,1]  [,2]  [,3]  [,4]  [,5]
 [1,] 11040  9820  9200  9260 10360
 [2,]  9885  8665  8045  8105  9205
 [3,]  9120  7900  7280  7340  8440
 [4,]  8655  7435  6815  6875  7975
 [5,]  8400  7180  6560  6620  7720
 [6,]  8685  7465  6845  6905  8005
 [7,]  9360  8140  7520  7580  8680
 [8,] 10245  9025  8405  8465  9565
 [9,] 11340 10120  9500  9560 10660
[10,] 12645 11425 10805 10865 11965

Approach Two

  • Treat the matrix as a transition matrix

  • Use graphical based algorithms instead of euclidean distance, because not all points are reached by other points

  • Obstacles are left outside of the transition matrix, meaning they are not able to be reached

  • Create a digraph of all the points in the grid

  • Use Djikstras shorest path algorithm to find the shorest distances between all of the points

  • Add up all of those distance to find the total sitances between all the possible points in the model

  • The N/S and E/W distances are divided by the number of incidencies at that location. Since there are 0 incidencies on some locations, 1 was added to everything to make the math possible.

  • Transition Matrix

Rio.Rancho
       [,1] [,2] [,3]     [,4]     [,5] [,6]     [,7] [,8] [,9]    [,10]
 [1,]  0.00   10    0 0.000000 0.000000 3.75 0.000000 0.00 0.00 0.000000
 [2,]  5.00    0    4 0.000000 0.000000 0.00 5.000000 0.00 0.00 0.000000
 [3,]  0.00   10    0 5.000000 0.000000 0.00 0.000000 3.75 0.00 0.000000
 [4,]  0.00    0    4 0.000000 3.333333 0.00 0.000000 0.00 3.75 0.000000
 [5,]  0.00    0    0 6.666667 0.000000 0.00 0.000000 0.00 0.00 5.000000
 [6,]  3.75    0    0 0.000000 0.000000 0.00 6.666667 0.00 0.00 0.000000
 [7,]  0.00    5    0 0.000000 0.000000 5.00 0.000000 5.00 0.00 0.000000
 [8,]  0.00    0    3 0.000000 0.000000 0.00 6.666667 0.00 5.00 0.000000
 [9,]  0.00    0    0 5.000000 0.000000 0.00 0.000000 5.00 0.00 6.666667
[10,]  0.00    0    0 0.000000 2.500000 0.00 0.000000 0.00 5.00 0.000000
[11,]  0.00    0    0 0.000000 0.000000 3.75 0.000000 0.00 0.00 0.000000
[12,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[13,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 3.75 0.00 0.000000
[14,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 3.75 0.000000
[15,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 5.000000
[16,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[17,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[18,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[19,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[20,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[21,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[22,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[23,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[24,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[25,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[26,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[27,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[28,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[29,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[30,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[31,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[32,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[33,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[34,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[35,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[36,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[37,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[38,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[39,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[40,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[41,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[42,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[43,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[44,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[45,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[46,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[47,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[48,] 10.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[49,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
[50,]  0.00    0    0 0.000000 0.000000 0.00 0.000000 0.00 0.00 0.000000
         [,11] [,12] [,13] [,14]    [,15] [,16] [,17] [,18] [,19] [,20]
 [1,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
 [2,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
 [3,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
 [4,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
 [5,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
 [6,] 6.666667     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
 [7,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
 [8,] 0.000000     0  3.75  0.00 0.000000  0.00     0     0  0.00   0.0
 [9,] 0.000000     0  0.00  3.75 0.000000  0.00     0     0  0.00   0.0
[10,] 0.000000     0  0.00  0.00 5.000000  0.00     0     0  0.00   0.0
[11,] 0.000000     0  0.00  0.00 0.000000  3.75     0     0  0.00   0.0
[12,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[13,] 0.000000     0  0.00  5.00 0.000000  0.00     0     0  0.00   0.0
[14,] 0.000000     0  5.00  0.00 6.666667  0.00     0     0  3.75   0.0
[15,] 0.000000     0  0.00  5.00 0.000000  0.00     0     0  0.00   7.5
[16,] 5.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[17,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[18,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[19,] 0.000000     0  0.00  0.00 3.750000  0.00     0     0  0.00  10.0
[20,] 0.000000     0  0.00  0.00 5.000000  0.00     0     0  5.00   0.0
[21,] 0.000000     0  0.00  0.00 0.000000  3.75     0     0  0.00   0.0
[22,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[23,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[24,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  3.75   0.0
[25,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   7.5
[26,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[27,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[28,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[29,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[30,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[31,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[32,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[33,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[34,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[35,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[36,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[37,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[38,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[39,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[40,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[41,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[42,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[43,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[44,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[45,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[46,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[47,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[48,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[49,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
[50,] 0.000000     0  0.00  0.00 0.000000  0.00     0     0  0.00   0.0
      [,21] [,22]    [,23] [,24]    [,25]    [,26] [,27] [,28] [,29] [,30]
 [1,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
 [2,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
 [3,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
 [4,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
 [5,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
 [6,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
 [7,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
 [8,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
 [9,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[10,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[11,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[12,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[13,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[14,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[15,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[16,]  3.75     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[17,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[18,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[19,]  0.00     0 0.000000  3.75 0.000000 0.000000  0.00     0     0     0
[20,]  0.00     0 0.000000  0.00 2.500000 0.000000  0.00     0     0     0
[21,]  0.00     4 0.000000  0.00 0.000000 5.000000  0.00     0     0     0
[22,]  0.00     0 6.666667  0.00 0.000000 0.000000  3.75     0     0     0
[23,]  0.00     0 0.000000  5.00 0.000000 0.000000  0.00     3     0     0
[24,]  0.00     0 0.000000  0.00 3.333333 0.000000  0.00     0     3     0
[25,]  0.00     0 0.000000  5.00 0.000000 0.000000  0.00     0     0    15
[26,]  3.75     0 0.000000  0.00 0.000000 0.000000  5.00     0     0     0
[27,]  0.00     3 0.000000  0.00 0.000000 6.666667  0.00     4     0     0
[28,]  0.00     0 3.750000  0.00 0.000000 0.000000  5.00     0     4     0
[29,]  0.00     0 0.000000  3.75 0.000000 0.000000  0.00     4     0    20
[30,]  0.00     0 0.000000  0.00 2.500000 0.000000  0.00     0     4     0
[31,]  0.00     0 0.000000  0.00 0.000000 5.000000  0.00     0     0     0
[32,]  0.00     0 0.000000  0.00 0.000000 0.000000  3.75     0     0     0
[33,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     3     0     0
[34,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     3     0
[35,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0    15
[36,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[37,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[38,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[39,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[40,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[41,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[42,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[43,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[44,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[45,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[46,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[47,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[48,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[49,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
[50,]  0.00     0 0.000000  0.00 0.000000 0.000000  0.00     0     0     0
      [,31]    [,32] [,33] [,34] [,35] [,36]    [,37] [,38] [,39]    [,40]
 [1,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
 [2,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
 [3,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
 [4,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
 [5,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
 [6,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
 [7,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
 [8,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
 [9,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[10,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[11,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[12,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[13,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[14,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[15,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[16,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[17,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[18,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[19,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[20,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[21,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[22,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[23,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[24,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[25,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[26,]   7.5 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[27,]   0.0 5.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[28,]   0.0 0.000000    15   0.0  0.00     0 0.000000     0  0.00 0.000000
[29,]   0.0 0.000000     0   7.5  0.00     0 0.000000     0  0.00 0.000000
[30,]   0.0 0.000000     0   0.0  3.75     0 0.000000     0  0.00 0.000000
[31,]   0.0 6.666667     0   0.0  0.00    15 0.000000     0  0.00 0.000000
[32,]  10.0 0.000000    20   0.0  0.00     0 5.000000     0  0.00 0.000000
[33,]   0.0 6.666667     0  10.0  0.00     0 0.000000     0  0.00 0.000000
[34,]   0.0 0.000000    20   0.0  5.00     0 0.000000     0  3.75 0.000000
[35,]   0.0 0.000000     0  10.0  0.00     0 0.000000     0  0.00 5.000000
[36,]   7.5 0.000000     0   0.0  0.00     0 6.666667     0  0.00 0.000000
[37,]   0.0 5.000000     0   0.0  0.00    20 0.000000     0  0.00 0.000000
[38,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[39,]   0.0 0.000000     0   7.5  0.00     0 0.000000     0  0.00 6.666667
[40,]   0.0 0.000000     0   0.0  3.75     0 0.000000     0  5.00 0.000000
[41,]   0.0 0.000000     0   0.0  0.00    15 0.000000     0  0.00 0.000000
[42,]   0.0 0.000000     0   0.0  0.00     0 5.000000     0  0.00 0.000000
[43,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[44,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  3.75 0.000000
[45,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 5.000000
[46,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[47,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[48,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[49,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
[50,]   0.0 0.000000     0   0.0  0.00     0 0.000000     0  0.00 0.000000
      [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49]    [,50]
 [1,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
 [2,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
 [3,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
 [4,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
 [5,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
 [6,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
 [7,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
 [8,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
 [9,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[10,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[11,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[12,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[13,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[14,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[15,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[16,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[17,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[18,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[19,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[20,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[21,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[22,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[23,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[24,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[25,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[26,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[27,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[28,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[29,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[30,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[31,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[32,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[33,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[34,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[35,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[36,]  3.75     0     0     0     0  0.00   0.0     0     0 0.000000
[37,]  0.00    15     0     0     0  0.00   0.0     0     0 0.000000
[38,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[39,]  0.00     0     0    15     0  0.00   0.0     0     0 0.000000
[40,]  0.00     0     0     0     3  0.00   0.0     0     0 0.000000
[41,]  0.00    20     0     0     0  3.75   0.0     0     0 0.000000
[42,]  5.00     0     0     0     0  0.00   7.5     0     0 0.000000
[43,]  0.00     0     0     0     0  0.00   0.0     0     0 0.000000
[44,]  0.00     0     0     0     4  0.00   0.0     0     3 0.000000
[45,]  0.00     0     0    20     0  0.00   0.0     0     0 5.000000
[46,]  3.75     0     0     0     0  0.00  10.0     0     0 0.000000
[47,]  0.00    15     0     0     0  5.00   0.0    20     0 0.000000
[48,]  0.00     0     0     0     0  0.00   0.0     0     4 0.000000
[49,]  0.00     0     0    15     0  0.00   0.0    20     0 6.666667
[50,]  0.00     0     0     0     4  0.00   0.0     0     4 0.000000

Plot

  • Digraph, showing the connection between all of the nodes

  • The nodes off to the sides are the barriers, they cannot be reached by anything

graph <- graph.adjacency(Rio.Rancho, mode = "directed", weight = TRUE)

plot.igraph(graph, layout = layout.fruchterman.reingold, vertex.label.color = "black", edge.color = "black")

plot of chunk unnamed-chunk-14

Shortest Path

  • Below we see the distance matrix, shows the distacne from the nodes on the rows, to each node across the columns.
  • Below that is the two smallest distances. The translate to the points (5,4) and (6,4)
  • The stations should now be stacked, compared to being next to each other. That could be due to the boundaries getting in the way.
s.path
          [,1]     [,2]      [,3]      [,4]      [,5]     [,6]     [,7]
 [1,]  0.00000  5.00000  9.000000 13.000000 16.333333  3.75000  8.75000
 [2,]  5.00000  0.00000  4.000000  8.000000 11.333333  8.75000  5.00000
 [3,]  9.00000  4.00000  0.000000  4.000000  7.333333 12.75000  8.00000
 [4,] 13.00000  8.00000  4.000000  0.000000  3.333333 16.75000 12.00000
 [5,] 16.33333 11.33333  7.333333  3.333333  0.000000 20.08333 15.33333
 [6,]  3.75000  8.75000 12.750000 16.750000 20.083333  0.00000  5.00000
 [7,]  8.75000  5.00000  8.000000 12.000000 15.333333  5.00000  0.00000
 [8,] 12.00000  7.00000  3.000000  7.000000 10.333333 10.00000  5.00000
 [9,] 16.75000 11.75000  7.750000  3.750000  7.083333 15.00000 10.00000
[10,] 18.83333 13.83333  9.833333  5.833333  2.500000 20.00000 15.00000
[11,]  7.50000 12.50000 16.500000 20.500000 23.833333  3.75000  8.75000
[12,]       NA       NA        NA        NA        NA       NA       NA
[13,] 15.75000 10.75000  6.750000 10.750000 14.083333 13.75000  8.75000
[14,] 20.50000 15.50000 11.500000  7.500000 10.833333 18.75000 13.75000
[15,] 23.83333 18.83333 14.833333 10.833333  7.500000 23.75000 18.75000
[16,] 11.25000 16.25000 20.250000 24.250000 27.583333  7.50000 12.50000
[17,]       NA       NA        NA        NA        NA       NA       NA
[18,]       NA       NA        NA        NA        NA       NA       NA
[19,] 24.25000 19.25000 15.250000 11.250000 11.250000 22.50000 17.50000
[20,] 28.83333 23.83333 19.833333 15.833333 12.500000 27.50000 22.50000
[21,] 15.00000 20.00000 24.000000 28.000000 30.666667 11.25000 16.25000
[22,] 19.00000 24.00000 28.000000 26.666667 26.666667 15.25000 20.25000
[23,] 25.66667 28.00000 24.000000 20.000000 20.000000 21.91667 26.25000
[24,] 28.00000 23.00000 19.000000 15.000000 15.000000 26.25000 21.25000
[25,] 31.33333 26.33333 22.333333 18.333333 15.000000 29.58333 24.58333
[26,] 18.75000 23.75000 27.750000 31.000000 31.000000 15.00000 20.00000
[27,] 22.00000 27.00000 30.000000 26.000000 26.000000 18.25000 23.25000
[28,] 26.00000 30.00000 26.000000 22.000000 22.000000 22.25000 27.25000
[29,] 27.50000 26.00000 22.000000 18.000000 18.000000 26.25000 24.25000
[30,] 31.50000 28.83333 24.833333 20.833333 17.500000 30.25000 27.08333
[31,] 23.75000 28.75000 32.750000 36.000000 36.000000 20.00000 25.00000
[32,] 25.75000 30.75000 33.750000 29.750000 29.750000 22.00000 27.00000
[33,] 29.00000 33.00000 29.000000 25.000000 25.000000 25.25000 30.25000
[34,] 24.50000 29.00000 25.000000 21.000000 21.000000 28.25000 27.25000
[35,] 27.75000 32.58333 28.583333 24.583333 21.250000 31.50000 30.83333
[36,] 31.25000 36.25000 40.250000 41.416667 41.416667 27.50000 32.50000
[37,] 30.75000 35.75000 38.750000 34.750000 34.750000 27.00000 32.00000
[38,]       NA       NA        NA        NA        NA       NA       NA
[39,] 20.75000 25.75000 28.750000 24.750000 24.750000 24.50000 29.50000
[40,] 24.00000 29.00000 32.333333 28.333333 25.000000 27.75000 32.75000
[41,] 35.00000 40.00000 44.000000 44.750000 44.750000 31.25000 36.25000
[42,] 35.75000 40.75000 43.750000 39.750000 39.750000 32.00000 37.00000
[43,]       NA       NA        NA        NA        NA       NA       NA
[44,] 17.00000 22.00000 26.000000 28.500000 28.500000 20.75000 25.75000
[45,] 21.00000 26.00000 30.000000 31.333333 28.000000 24.75000 29.75000
[46,] 35.00000 40.00000 44.000000 48.000000 48.500000 35.00000 40.00000
[47,] 30.00000 35.00000 39.000000 43.000000 46.333333 33.75000 38.75000
[48,] 10.00000 15.00000 19.000000 23.000000 26.333333 13.75000 18.75000
[49,] 14.00000 19.00000 23.000000 27.000000 30.333333 17.75000 22.75000
[50,] 18.00000 23.00000 27.000000 31.000000 32.000000 21.75000 26.75000
          [,8]      [,9]     [,10]    [,11] [,12]    [,13]    [,14]
 [1,] 12.00000 16.750000 18.833333  7.50000    NA 15.75000 20.50000
 [2,]  7.00000 11.750000 13.833333 12.50000    NA 10.75000 15.50000
 [3,]  3.00000  7.750000  9.833333 16.50000    NA  6.75000 11.50000
 [4,]  7.00000  3.750000  5.833333 20.50000    NA 10.75000  7.50000
 [5,] 10.33333  7.083333  2.500000 23.83333    NA 14.08333 10.83333
 [6,] 10.00000 15.000000 20.000000  3.75000    NA 13.75000 18.75000
 [7,]  5.00000 10.000000 15.000000  8.75000    NA  8.75000 13.75000
 [8,]  0.00000  5.000000 10.000000 13.75000    NA  3.75000  8.75000
 [9,]  5.00000  0.000000  5.000000 18.75000    NA  8.75000  3.75000
[10,] 10.00000  5.000000  0.000000 23.75000    NA 13.75000  8.75000
[11,] 13.75000 18.750000 23.750000  0.00000    NA 17.50000 22.50000
[12,]       NA        NA        NA       NA     0       NA       NA
[13,]  3.75000  8.750000 13.750000 17.50000    NA  0.00000  5.00000
[14,]  8.75000  3.750000  8.750000 22.50000    NA  5.00000  0.00000
[15,] 13.75000  8.750000  5.000000 27.50000    NA 10.00000  5.00000
[16,] 17.50000 22.500000 27.500000  3.75000    NA 21.25000 26.25000
[17,]       NA        NA        NA       NA    NA       NA       NA
[18,]       NA        NA        NA       NA    NA       NA       NA
[19,] 12.50000  7.500000  8.750000 26.25000    NA  8.75000  3.75000
[20,] 17.50000 12.500000 10.000000 29.00000    NA 13.75000  8.75000
[21,] 21.25000 26.250000 28.166667  7.50000    NA 25.00000 23.16667
[22,] 25.25000 22.916667 24.166667 11.50000    NA 24.16667 19.16667
[23,] 21.25000 16.250000 17.500000 18.16667    NA 17.50000 12.50000
[24,] 16.25000 11.250000 12.500000 23.16667    NA 12.50000  7.50000
[25,] 19.58333 14.583333 12.500000 26.50000    NA 15.83333 10.83333
[26,] 25.00000 27.250000 28.500000 11.25000    NA 28.50000 23.50000
[27,] 27.25000 22.250000 23.500000 14.50000    NA 23.50000 18.50000
[28,] 23.25000 18.250000 19.500000 18.50000    NA 19.50000 14.50000
[29,] 19.25000 14.250000 15.500000 22.50000    NA 15.50000 10.50000
[30,] 22.08333 17.083333 15.000000 26.50000    NA 18.33333 13.33333
[31,] 30.00000 32.250000 33.500000 16.25000    NA 33.50000 28.50000
[32,] 31.00000 26.000000 27.250000 18.25000    NA 27.25000 22.25000
[33,] 26.25000 21.250000 22.500000 21.50000    NA 22.50000 17.50000
[34,] 22.25000 17.250000 18.500000 25.50000    NA 18.50000 13.50000
[35,] 25.83333 20.833333 18.750000 30.25000    NA 22.08333 17.08333
[36,] 37.50000 37.666667 38.916667 23.75000    NA 38.91667 33.91667
[37,] 36.00000 31.000000 32.250000 23.25000    NA 32.25000 27.25000
[38,]       NA        NA        NA       NA    NA       NA       NA
[39,] 26.00000 21.000000 22.250000 28.25000    NA 22.25000 17.25000
[40,] 29.58333 24.583333 22.500000 31.50000    NA 25.83333 20.83333
[41,] 41.25000 41.000000 42.250000 27.50000    NA 42.25000 37.25000
[42,] 41.00000 36.000000 37.250000 28.25000    NA 37.25000 32.25000
[43,]       NA        NA        NA       NA    NA       NA       NA
[44,] 29.00000 24.750000 26.000000 24.50000    NA 26.00000 21.00000
[45,] 32.58333 27.583333 25.500000 28.50000    NA 28.83333 23.83333
[46,] 45.00000 44.750000 46.000000 31.25000    NA 46.00000 41.00000
[47,] 42.00000 43.500000 44.750000 35.75000    NA 44.75000 39.75000
[48,] 22.00000 26.750000 28.833333 17.50000    NA 25.75000 28.00000
[49,] 26.00000 27.750000 29.000000 21.50000    NA 29.00000 24.00000
[50,] 30.00000 31.583333 29.500000 25.50000    NA 32.83333 27.83333
         [,15]    [,16] [,17] [,18]     [,19]     [,20]    [,21]     [,22]
 [1,] 23.83333 11.25000    NA    NA 24.250000 28.833333 15.00000 19.000000
 [2,] 18.83333 16.25000    NA    NA 19.250000 23.833333 20.00000 24.000000
 [3,] 14.83333 20.25000    NA    NA 15.250000 19.833333 24.00000 28.000000
 [4,] 10.83333 24.25000    NA    NA 11.250000 15.833333 28.00000 26.666667
 [5,]  7.50000 27.58333    NA    NA 11.250000 12.500000 30.66667 26.666667
 [6,] 23.75000  7.50000    NA    NA 22.500000 27.500000 11.25000 15.250000
 [7,] 18.75000 12.50000    NA    NA 17.500000 22.500000 16.25000 20.250000
 [8,] 13.75000 17.50000    NA    NA 12.500000 17.500000 21.25000 25.250000
 [9,]  8.75000 22.50000    NA    NA  7.500000 12.500000 26.25000 22.916667
[10,]  5.00000 27.50000    NA    NA  8.750000 10.000000 28.16667 24.166667
[11,] 27.50000  3.75000    NA    NA 26.250000 29.000000  7.50000 11.500000
[12,]       NA       NA    NA    NA        NA        NA       NA        NA
[13,] 10.00000 21.25000    NA    NA  8.750000 13.750000 25.00000 24.166667
[14,]  5.00000 26.25000    NA    NA  3.750000  8.750000 23.16667 19.166667
[15,]  0.00000 26.91667    NA    NA  3.750000  5.000000 23.16667 19.166667
[16,] 26.91667  0.00000    NA    NA 23.166667 25.250000  3.75000  7.750000
[17,]       NA       NA     0    NA        NA        NA       NA        NA
[18,]       NA       NA    NA     0        NA        NA       NA        NA
[19,]  3.75000 23.16667    NA    NA  0.000000  5.000000 19.41667 15.416667
[20,]  5.00000 25.25000    NA    NA  5.000000  0.000000 21.50000 17.500000
[21,] 23.16667  3.75000    NA    NA 19.416667 21.500000  0.00000  4.000000
[22,] 19.16667  7.75000    NA    NA 15.416667 17.500000  4.00000  0.000000
[23,] 12.50000 14.41667    NA    NA  8.750000 10.833333 10.66667  6.666667
[24,]  7.50000 19.41667    NA    NA  3.750000  5.833333 15.66667 11.666667
[25,]  7.50000 22.75000    NA    NA  7.083333  2.500000 19.00000 15.000000
[26,] 23.50000  7.50000    NA    NA 19.750000 21.833333  3.75000  7.750000
[27,] 18.50000 10.75000    NA    NA 14.750000 16.833333  7.00000  3.000000
[28,] 14.50000 14.75000    NA    NA 10.750000 12.833333 11.00000  7.000000
[29,] 10.50000 18.75000    NA    NA  6.750000  8.833333 15.00000 11.000000
[30,] 10.00000 22.75000    NA    NA  9.583333  5.000000 19.00000 15.000000
[31,] 28.50000 12.50000    NA    NA 24.750000 26.833333  8.75000 12.750000
[32,] 22.25000 14.50000    NA    NA 18.500000 20.583333 10.75000  6.750000
[33,] 17.50000 17.75000    NA    NA 13.750000 15.833333 14.00000 10.000000
[34,] 13.50000 21.75000    NA    NA  9.750000 11.833333 18.00000 14.000000
[35,] 13.75000 26.50000    NA    NA 13.333333  8.750000 22.75000 18.750000
[36,] 33.91667 20.00000    NA    NA 30.166667 32.250000 16.25000 18.416667
[37,] 27.25000 19.50000    NA    NA 23.500000 25.583333 15.75000 11.750000
[38,]       NA       NA    NA    NA        NA        NA       NA        NA
[39,] 17.25000 25.50000    NA    NA 13.500000 15.583333 21.75000 17.750000
[40,] 17.50000 30.25000    NA    NA 17.083333 12.500000 26.50000 22.500000
[41,] 37.25000 23.75000    NA    NA 33.500000 35.583333 20.00000 21.750000
[42,] 32.25000 24.50000    NA    NA 28.500000 30.583333 20.75000 16.750000
[43,]       NA       NA    NA    NA        NA        NA       NA        NA
[44,] 21.00000 28.25000    NA    NA 17.250000 19.333333 25.50000 21.500000
[45,] 20.50000 32.25000    NA    NA 20.083333 15.500000 29.50000 25.500000
[46,] 41.00000 27.50000    NA    NA 37.250000 39.333333 23.75000 25.500000
[47,] 39.75000 32.00000    NA    NA 36.000000 38.083333 28.25000 24.250000
[48,] 28.00000 21.25000    NA    NA 24.250000 26.333333 25.00000 28.500000
[49,] 24.00000 25.25000    NA    NA 20.250000 22.333333 28.50000 24.500000
[50,] 24.50000 29.25000    NA    NA 24.083333 19.500000 32.50000 28.500000
          [,23]     [,24]     [,25]    [,26]    [,27]    [,28]     [,29]
 [1,] 25.666667 28.000000 31.333333 18.75000 22.00000 26.00000 27.500000
 [2,] 28.000000 23.000000 26.333333 23.75000 27.00000 30.00000 26.000000
 [3,] 24.000000 19.000000 22.333333 27.75000 30.00000 26.00000 22.000000
 [4,] 20.000000 15.000000 18.333333 31.00000 26.00000 22.00000 18.000000
 [5,] 20.000000 15.000000 15.000000 31.00000 26.00000 22.00000 18.000000
 [6,] 21.916667 26.250000 29.583333 15.00000 18.25000 22.25000 26.250000
 [7,] 26.250000 21.250000 24.583333 20.00000 23.25000 27.25000 24.250000
 [8,] 21.250000 16.250000 19.583333 25.00000 27.25000 23.25000 19.250000
 [9,] 16.250000 11.250000 14.583333 27.25000 22.25000 18.25000 14.250000
[10,] 17.500000 12.500000 12.500000 28.50000 23.50000 19.50000 15.500000
[11,] 18.166667 23.166667 26.500000 11.25000 14.50000 18.50000 22.500000
[12,]        NA        NA        NA       NA       NA       NA        NA
[13,] 17.500000 12.500000 15.833333 28.50000 23.50000 19.50000 15.500000
[14,] 12.500000  7.500000 10.833333 23.50000 18.50000 14.50000 10.500000
[15,] 12.500000  7.500000  7.500000 23.50000 18.50000 14.50000 10.500000
[16,] 14.416667 19.416667 22.750000  7.50000 10.75000 14.75000 18.750000
[17,]        NA        NA        NA       NA       NA       NA        NA
[18,]        NA        NA        NA       NA       NA       NA        NA
[19,]  8.750000  3.750000  7.083333 19.75000 14.75000 10.75000  6.750000
[20,] 10.833333  5.833333  2.500000 21.83333 16.83333 12.83333  8.833333
[21,] 10.666667 15.666667 19.000000  3.75000  7.00000 11.00000 15.000000
[22,]  6.666667 11.666667 15.000000  7.75000  3.00000  7.00000 11.000000
[23,]  0.000000  5.000000  8.333333 12.00000  7.00000  3.00000  7.000000
[24,]  5.000000  0.000000  3.333333 16.00000 11.00000  7.00000  3.000000
[25,]  8.333333  3.333333  0.000000 19.33333 14.33333 10.33333  6.333333
[26,] 12.000000 16.000000 19.333333  0.00000  5.00000  9.00000 13.000000
[27,]  7.000000 11.000000 14.333333  5.00000  0.00000  4.00000  8.000000
[28,]  3.000000  7.000000 10.333333  9.00000  4.00000  0.00000  4.000000
[29,]  7.000000  3.000000  6.333333 13.00000  8.00000  4.00000  0.000000
[30,] 10.833333  5.833333  2.500000 17.00000 12.00000  8.00000  4.000000
[31,] 17.000000 21.000000 24.333333  5.00000 10.00000 14.00000 18.000000
[32,] 10.750000 14.750000 18.083333  8.75000  3.75000  7.75000 11.750000
[33,]  6.000000 10.000000 13.333333 12.00000  7.00000  3.00000  7.000000
[34,] 10.000000  6.000000  9.333333 16.00000 11.00000  7.00000  3.000000
[35,] 14.583333  9.583333  6.250000 20.75000 15.75000 11.75000  7.750000
[36,] 22.416667 26.416667 29.750000 12.50000 15.41667 19.41667 23.416667
[37,] 15.750000 19.750000 23.083333 13.75000  8.75000 12.75000 16.750000
[38,]        NA        NA        NA       NA       NA       NA        NA
[39,] 13.750000  9.750000 13.083333 19.75000 14.75000 10.75000  6.750000
[40,] 18.333333 13.333333 10.000000 24.50000 19.50000 15.50000 11.500000
[41,] 25.750000 29.750000 33.083333 16.25000 18.75000 22.75000 26.750000
[42,] 20.750000 24.750000 28.083333 18.75000 13.75000 17.75000 21.750000
[43,]        NA        NA        NA       NA       NA       NA        NA
[44,] 17.500000 13.500000 16.833333 23.50000 18.50000 14.50000 10.500000
[45,] 21.333333 16.333333 13.000000 27.50000 22.50000 18.50000 14.500000
[46,] 29.500000 33.500000 36.833333 20.00000 22.50000 26.50000 30.500000
[47,] 28.250000 32.250000 35.583333 25.00000 21.25000 25.25000 29.250000
[48,] 24.500000 20.500000 23.833333 28.75000 25.50000 21.50000 17.500000
[49,] 20.500000 16.500000 19.833333 26.50000 21.50000 17.50000 13.500000
[50,] 24.500000 20.333333 17.000000 30.50000 25.50000 21.50000 17.500000
          [,30]     [,31]     [,32]     [,33]     [,34]     [,35]
 [1,] 31.500000 23.750000 25.750000 29.000000 24.500000 27.750000
 [2,] 28.833333 28.750000 30.750000 33.000000 29.000000 32.583333
 [3,] 24.833333 32.750000 33.750000 29.000000 25.000000 28.583333
 [4,] 20.833333 36.000000 29.750000 25.000000 21.000000 24.583333
 [5,] 17.500000 36.000000 29.750000 25.000000 21.000000 21.250000
 [6,] 30.250000 20.000000 22.000000 25.250000 28.250000 31.500000
 [7,] 27.083333 25.000000 27.000000 30.250000 27.250000 30.833333
 [8,] 22.083333 30.000000 31.000000 26.250000 22.250000 25.833333
 [9,] 17.083333 32.250000 26.000000 21.250000 17.250000 20.833333
[10,] 15.000000 33.500000 27.250000 22.500000 18.500000 18.750000
[11,] 26.500000 16.250000 18.250000 21.500000 25.500000 30.250000
[12,]        NA        NA        NA        NA        NA        NA
[13,] 18.333333 33.500000 27.250000 22.500000 18.500000 22.083333
[14,] 13.333333 28.500000 22.250000 17.500000 13.500000 17.083333
[15,] 10.000000 28.500000 22.250000 17.500000 13.500000 13.750000
[16,] 22.750000 12.500000 14.500000 17.750000 21.750000 26.500000
[17,]        NA        NA        NA        NA        NA        NA
[18,]        NA        NA        NA        NA        NA        NA
[19,]  9.583333 24.750000 18.500000 13.750000  9.750000 13.333333
[20,]  5.000000 26.833333 20.583333 15.833333 11.833333  8.750000
[21,] 19.000000  8.750000 10.750000 14.000000 18.000000 22.750000
[22,] 15.000000 12.750000  6.750000 10.000000 14.000000 18.750000
[23,] 10.833333 17.000000 10.750000  6.000000 10.000000 14.583333
[24,]  5.833333 21.000000 14.750000 10.000000  6.000000  9.583333
[25,]  2.500000 24.333333 18.083333 13.333333  9.333333  6.250000
[26,] 17.000000  5.000000  8.750000 12.000000 16.000000 20.750000
[27,] 12.000000 10.000000  3.750000  7.000000 11.000000 15.750000
[28,]  8.000000 14.000000  7.750000  3.000000  7.000000 11.750000
[29,]  4.000000 18.000000 11.750000  7.000000  3.000000  7.750000
[30,]  0.000000 22.000000 15.750000 11.000000  7.000000  3.750000
[31,] 22.000000  0.000000  6.666667 13.333333 21.000000 25.750000
[32,] 15.750000  6.666667  0.000000  6.666667 14.750000 19.500000
[33,] 11.000000 13.333333  6.666667  0.000000 10.000000 14.750000
[34,]  7.000000 21.000000 14.750000 10.000000  0.000000  5.000000
[35,]  3.750000 25.750000 19.500000 14.750000  5.000000  0.000000
[36,] 27.416667  7.500000 11.666667 18.333333 26.416667 31.166667
[37,] 20.750000 11.666667  5.000000 11.666667 19.750000 24.500000
[38,]        NA        NA        NA        NA        NA        NA
[39,] 10.750000 24.750000 18.500000 13.750000  3.750000  8.750000
[40,]  7.500000 29.500000 23.250000 18.500000  8.750000  3.750000
[41,] 30.750000 11.250000 15.000000 21.666667 29.750000 34.500000
[42,] 25.750000 16.250000 10.000000 16.666667 24.750000 29.500000
[43,]        NA        NA        NA        NA        NA        NA
[44,] 14.500000 28.500000 22.250000 17.500000  7.500000 10.750000
[45,] 10.500000 32.500000 26.250000 21.500000 11.500000  6.750000
[46,] 34.500000 15.000000 18.750000 25.416667 33.500000 38.250000
[47,] 33.250000 20.000000 17.500000 24.166667 32.250000 37.000000
[48,] 21.500000 33.750000 29.250000 24.500000 14.500000 17.750000
[49,] 17.500000 31.500000 25.250000 20.500000 10.500000 13.750000
[50,] 14.500000 35.500000 29.250000 24.500000 14.500000 10.750000
          [,36]     [,37] [,38]    [,39]    [,40]    [,41]    [,42] [,43]
 [1,] 31.250000 30.750000    NA 20.75000 24.00000 35.00000 35.75000    NA
 [2,] 36.250000 35.750000    NA 25.75000 29.00000 40.00000 40.75000    NA
 [3,] 40.250000 38.750000    NA 28.75000 32.33333 44.00000 43.75000    NA
 [4,] 41.416667 34.750000    NA 24.75000 28.33333 44.75000 39.75000    NA
 [5,] 41.416667 34.750000    NA 24.75000 25.00000 44.75000 39.75000    NA
 [6,] 27.500000 27.000000    NA 24.50000 27.75000 31.25000 32.00000    NA
 [7,] 32.500000 32.000000    NA 29.50000 32.75000 36.25000 37.00000    NA
 [8,] 37.500000 36.000000    NA 26.00000 29.58333 41.25000 41.00000    NA
 [9,] 37.666667 31.000000    NA 21.00000 24.58333 41.00000 36.00000    NA
[10,] 38.916667 32.250000    NA 22.25000 22.50000 42.25000 37.25000    NA
[11,] 23.750000 23.250000    NA 28.25000 31.50000 27.50000 28.25000    NA
[12,]        NA        NA    NA       NA       NA       NA       NA    NA
[13,] 38.916667 32.250000    NA 22.25000 25.83333 42.25000 37.25000    NA
[14,] 33.916667 27.250000    NA 17.25000 20.83333 37.25000 32.25000    NA
[15,] 33.916667 27.250000    NA 17.25000 17.50000 37.25000 32.25000    NA
[16,] 20.000000 19.500000    NA 25.50000 30.25000 23.75000 24.50000    NA
[17,]        NA        NA    NA       NA       NA       NA       NA    NA
[18,]        NA        NA    NA       NA       NA       NA       NA    NA
[19,] 30.166667 23.500000    NA 13.50000 17.08333 33.50000 28.50000    NA
[20,] 32.250000 25.583333    NA 15.58333 12.50000 35.58333 30.58333    NA
[21,] 16.250000 15.750000    NA 21.75000 26.50000 20.00000 20.75000    NA
[22,] 18.416667 11.750000    NA 17.75000 22.50000 21.75000 16.75000    NA
[23,] 22.416667 15.750000    NA 13.75000 18.33333 25.75000 20.75000    NA
[24,] 26.416667 19.750000    NA  9.75000 13.33333 29.75000 24.75000    NA
[25,] 29.750000 23.083333    NA 13.08333 10.00000 33.08333 28.08333    NA
[26,] 12.500000 13.750000    NA 19.75000 24.50000 16.25000 18.75000    NA
[27,] 15.416667  8.750000    NA 14.75000 19.50000 18.75000 13.75000    NA
[28,] 19.416667 12.750000    NA 10.75000 15.50000 22.75000 17.75000    NA
[29,] 23.416667 16.750000    NA  6.75000 11.50000 26.75000 21.75000    NA
[30,] 27.416667 20.750000    NA 10.75000  7.50000 30.75000 25.75000    NA
[31,]  7.500000 11.666667    NA 24.75000 29.50000 11.25000 16.25000    NA
[32,] 11.666667  5.000000    NA 18.50000 23.25000 15.00000 10.00000    NA
[33,] 18.333333 11.666667    NA 13.75000 18.50000 21.66667 16.66667    NA
[34,] 26.416667 19.750000    NA  3.75000  8.75000 29.75000 24.75000    NA
[35,] 31.166667 24.500000    NA  8.75000  3.75000 34.50000 29.50000    NA
[36,]  0.000000  6.666667    NA 30.16667 34.91667  3.75000  8.75000    NA
[37,]  6.666667  0.000000    NA 23.50000 28.25000 10.00000  5.00000    NA
[38,]        NA        NA     0       NA       NA       NA       NA    NA
[39,] 30.166667 23.500000    NA  0.00000  5.00000 33.50000 28.50000    NA
[40,] 34.916667 28.250000    NA  5.00000  0.00000 38.25000 33.25000    NA
[41,]  3.750000 10.000000    NA 33.50000 38.25000  0.00000  5.00000    NA
[42,]  8.750000  5.000000    NA 28.50000 33.25000  5.00000  0.00000    NA
[43,]        NA        NA    NA       NA       NA       NA       NA     0
[44,] 33.916667 27.250000    NA  3.75000  7.00000 35.75000 32.25000    NA
[45,] 37.916667 31.250000    NA  7.75000  3.00000 39.75000 36.25000    NA
[46,]  7.500000 13.750000    NA 35.75000 39.00000  3.75000  8.75000    NA
[47,] 12.500000 12.500000    NA 30.75000 34.00000  8.75000  7.50000    NA
[48,] 32.500000 32.500000    NA 10.75000 14.00000 28.75000 27.50000    NA
[49,] 36.500000 30.250000    NA  6.75000 10.00000 32.75000 31.50000    NA
[50,] 40.500000 34.250000    NA 10.75000  7.00000 36.75000 35.50000    NA
         [,44]    [,45]    [,46]    [,47]    [,48]    [,49]    [,50]
 [1,] 17.00000 21.00000 35.00000 30.00000 10.00000 14.00000 18.00000
 [2,] 22.00000 26.00000 40.00000 35.00000 15.00000 19.00000 23.00000
 [3,] 26.00000 30.00000 44.00000 39.00000 19.00000 23.00000 27.00000
 [4,] 28.50000 31.33333 48.00000 43.00000 23.00000 27.00000 31.00000
 [5,] 28.50000 28.00000 48.50000 46.33333 26.33333 30.33333 32.00000
 [6,] 20.75000 24.75000 35.00000 33.75000 13.75000 17.75000 21.75000
 [7,] 25.75000 29.75000 40.00000 38.75000 18.75000 22.75000 26.75000
 [8,] 29.00000 32.58333 45.00000 42.00000 22.00000 26.00000 30.00000
 [9,] 24.75000 27.58333 44.75000 43.50000 26.75000 27.75000 31.58333
[10,] 26.00000 25.50000 46.00000 44.75000 28.83333 29.00000 29.50000
[11,] 24.50000 28.50000 31.25000 35.75000 17.50000 21.50000 25.50000
[12,]       NA       NA       NA       NA       NA       NA       NA
[13,] 26.00000 28.83333 46.00000 44.75000 25.75000 29.00000 32.83333
[14,] 21.00000 23.83333 41.00000 39.75000 28.00000 24.00000 27.83333
[15,] 21.00000 20.50000 41.00000 39.75000 28.00000 24.00000 24.50000
[16,] 28.25000 32.25000 27.50000 32.00000 21.25000 25.25000 29.25000
[17,]       NA       NA       NA       NA       NA       NA       NA
[18,]       NA       NA       NA       NA       NA       NA       NA
[19,] 17.25000 20.08333 37.25000 36.00000 24.25000 20.25000 24.08333
[20,] 19.33333 15.50000 39.33333 38.08333 26.33333 22.33333 19.50000
[21,] 25.50000 29.50000 23.75000 28.25000 25.00000 28.50000 32.50000
[22,] 21.50000 25.50000 25.50000 24.25000 28.50000 24.50000 28.50000
[23,] 17.50000 21.33333 29.50000 28.25000 24.50000 20.50000 24.50000
[24,] 13.50000 16.33333 33.50000 32.25000 20.50000 16.50000 20.33333
[25,] 16.83333 13.00000 36.83333 35.58333 23.83333 19.83333 17.00000
[26,] 23.50000 27.50000 20.00000 25.00000 28.75000 26.50000 30.50000
[27,] 18.50000 22.50000 22.50000 21.25000 25.50000 21.50000 25.50000
[28,] 14.50000 18.50000 26.50000 25.25000 21.50000 17.50000 21.50000
[29,] 10.50000 14.50000 30.50000 29.25000 17.50000 13.50000 17.50000
[30,] 14.50000 10.50000 34.50000 33.25000 21.50000 17.50000 14.50000
[31,] 28.50000 32.50000 15.00000 20.00000 33.75000 31.50000 35.50000
[32,] 22.25000 26.25000 18.75000 17.50000 29.25000 25.25000 29.25000
[33,] 17.50000 21.50000 25.41667 24.16667 24.50000 20.50000 24.50000
[34,]  7.50000 11.50000 33.50000 32.25000 14.50000 10.50000 14.50000
[35,] 10.75000  6.75000 38.25000 37.00000 17.75000 13.75000 10.75000
[36,] 33.91667 37.91667  7.50000 12.50000 32.50000 36.50000 40.50000
[37,] 27.25000 31.25000 13.75000 12.50000 32.50000 30.25000 34.25000
[38,]       NA       NA       NA       NA       NA       NA       NA
[39,]  3.75000  7.75000 35.75000 30.75000 10.75000  6.75000 10.75000
[40,]  7.00000  3.00000 39.00000 34.00000 14.00000 10.00000  7.00000
[41,] 35.75000 39.75000  3.75000  8.75000 28.75000 32.75000 36.75000
[42,] 32.25000 36.25000  8.75000  7.50000 27.50000 31.50000 35.50000
[43,]       NA       NA       NA       NA       NA       NA       NA
[44,]  0.00000  4.00000 32.00000 27.00000  7.00000  3.00000  7.00000
[45,]  4.00000  0.00000 36.00000 31.00000 11.00000  7.00000  4.00000
[46,] 32.00000 36.00000  0.00000  5.00000 25.00000 29.00000 33.00000
[47,] 27.00000 31.00000  5.00000  0.00000 20.00000 24.00000 28.00000
[48,]  7.00000 11.00000 25.00000 20.00000  0.00000  4.00000  8.00000
[49,]  3.00000  7.00000 29.00000 24.00000  4.00000  0.00000  4.00000
[50,]  7.00000  4.00000 33.00000 28.00000  8.00000  4.00000  0.00000
print("/n")
[1] "/n"
mins
[1] 663.3333 683.6667