diff --git a/isleOfWight/beisIoWData.Rmd b/isleOfWight/beisIoWData.Rmd index 39a2a807c770406f3c04f4f59f023f76f7581cff..c18feb11754569793ffde847abb807aec6db2d78 100644 --- a/isleOfWight/beisIoWData.Rmd +++ b/isleOfWight/beisIoWData.Rmd @@ -199,10 +199,15 @@ nrow(dtMsoaSpineIoW[elecRatio > 0.80]) # Ellis' question... -1) join income deprivation score to LSOA shapefiles. If you could produce a map that would be great. -2) Count dwellings in LSOAs with have EPC < C (D-G) and main fuel type not mains gas. If you could map all EPCs and then the EPCs that match the criteria. -3) Join tenure to LSOA shapefiles. Multiply count by LSOA level census figure for % of private dwellings. -4) Multiply by 1/0.687 to correct for dwellings without an EPC. +For the Isle of Wight, can we estimate the number of private (owned or rented) addresses in high income deprivation areas which do not have gas? + +Our method: + + 1. Join LSOA data to EPC data to + join IMD income deprivation score to LSOA. + 2. Count dwellings in LSOAs with have EPC < C (D-G) and main fuel type not mains gas. If you could map all EPCs and then the EPCs that match the criteria. + 3. Join tenure to LSOA shapefiles. Multiply count by LSOA level census figure for % of private dwellings. + 4. Multiply by 1/0.687 to correct for dwellings without an EPC. ```{r loadData} # this will break here if you don't have the data obvs @@ -427,21 +432,26 @@ Table \@ref(tab:depr10pc) shows the same results but for the 10% most income dep Next we assume that the EPC addresses of interest are as likely to be social/private rented as owner-occupied and multiply the EPC address counts by the Census-derived household ownership tenure rate for each LSOA. Without knowing the tensure status of the EPC we cannot do anything else. ```{r select20pc} -lsoaDT[, nHHsSelected := selectedEPCAddresses*`Owned (%)`] +# multiply by tenure needed +lsoaDT[, nEPCAddressessSelectedOwned := selectedEPCAddresses*`Owned (%)`] # just owned +lsoaDT[, nEPCAddressessSelectedPrivate := selectedEPCAddresses*(`Owned (%)` + `Private rented (%)`)] # just owned -t <- lsoaDT[, .(nHHsSelected = sum(nHHsSelected), +t <- lsoaDT[, .(nEPCAddressessSelectedOwned = sum(nEPCAddressessSelectedOwned), + nEPCAddressessSelectedPrivate = sum(nEPCAddressessSelectedPrivate), nHHsCensus = sum(nHHsCensus), nLSOAs = .N), keyby = .(deprived20pc)] kableExtra::kable(t, digits = 2, caption = "Summary of EPC households and all census households by LSOA deprivation type (20% most income deprived)") %>% kable_styling() -res20pc <- round(t[deprived20pc == "highDep", nHHsSelected]) +res20pcOwned <- round(t[deprived20pc == "highDep", nEPCAddressessSelectedOwned]) +res20pcPrivate <- round(t[deprived20pc == "highDep", nEPCAddressessSelectedPrivate]) ``` ```{r select10pc} -t <- lsoaDT[, .(nHHsSelected = sum(nHHsSelected), +t <- lsoaDT[, .(nEPCAddressessSelectedOwned = sum(nEPCAddressessSelectedOwned), + nEPCAddressessSelectedPrivate = sum(nEPCAddressessSelectedPrivate), nHHsCensus = sum(nHHsCensus), nLSOAs = .N), keyby = .(deprived10pc)] @@ -449,21 +459,26 @@ kableExtra::kable(t, digits = 2, caption = "Summary of EPC households and all census households by LSOA deprivation type (10% most income deprived)") %>% kable_styling() -res10pc <- round(t[deprived10pc == "highDep", nHHsSelected]) - +res10pcOwned <- round(t[deprived10pc == "highDep", nEPCAddressessSelectedOwned]) +res10pcPrivate <- round(t[deprived10pc == "highDep", nEPCAddressessSelectedPrivate]) epcMissingPC <- 100*(1-(1/epcMissingFactor)) ``` -> In summary using the EPC data we have estimated the number of addresses in low income areas which are owned and do not have gas to be: +> In summary using the EPC data we have estimated the number of addresses in high income deprivation areas which do not have gas to be: - * 10% most deprived: `r res10pc` or, if we uprate to allow for the `r epcMissingPC` % missing EPCs, `r round(res10pc * epcMissingFactor)`; - * 20% most deprived: `r res20pc` or, if we uprate to allow for the `r epcMissingPC` % missing EPCs, `r round(res20pc * epcMissingFactor)`; + * 10% most deprived: + * Owned: `r res10pcOwned` or, if we uprate to allow for the `r epcMissingPC` % missing EPCs, `r round(res10pcOwned * epcMissingFactor)` + * Private (owned & private rent): `r res10pcPrivate` or, if we uprate to allow for the `r epcMissingPC` % missing EPCs, `r round(res10pcPrivate * epcMissingFactor)` + * 20% most deprived: + * Owned: `r res20pcOwned` or, if we uprate to allow for the `r epcMissingPC` % missing EPCs, `r round(res20pcOwned * epcMissingFactor)` + * Private (owned & private rent): `r res20pcPrivate` or, if we uprate to allow for the `r epcMissingPC` % missing EPCs, `r round(res20pcPrivate * epcMissingFactor)` Table \@ref(tab:lsoaTable) shows the LSOA level results. ```{r lsoaTable} -t<- lsoaDT[, .(lsoaCode, nAddressesEPC = nHHsEPC, notSelectedEPCAddresses, selectedEPCAddresses, pcOwned=pcOwned*100, deprived20pc, deprived10pc,nHHsSelected)] +t<- lsoaDT[, .(lsoaCode, nAddressesEPC = nHHsEPC, notSelectedEPCAddresses, selectedEPCAddresses, pcOwned=pcOwned*100, deprived20pc, deprived10pc, + nEPCAddressessSelectedOwned,nEPCAddressessSelectedPrivate)] setkey(t, deprived20pc, deprived10pc) kableExtra::kable(t, caption = "Final results at LSOA level", @@ -476,20 +491,33 @@ data.table::fwrite(t, file = "lsoaTableForEllis.csv") ```{r resultsChecks} -totSelectedHH <- sum(lsoaDT[!is.na(lsoaCode)]$selectedEPCAddresses) +totPrivSelected <- sum(lsoaDT[!is.na(lsoaCode)]$nEPCAddressessSelectedPrivate) +totOwnedSelected <- sum(lsoaDT[!is.na(lsoaCode)]$nEPCAddressessSelectedOwned) + +meanPrivate <- mean(lsoaDT$pcOwned + lsoaDT$`Private rented (%)`) + +res10own <- totOwnedSelected * meanPrivate * 0.1 +res20own <- totOwnedSelected * meanPrivate * 0.2 +res10priv <- totPrivSelected * meanPrivate * 0.1 +res20priv <- totPrivSelected * meanPrivate * 0.2 -meanOwnership <- mean(lsoaDT$pcOwned) +res10priv +res10own -res10 <- totSelectedHH * meanOwnership * 0.1 -res20 <- totSelectedHH * meanOwnership * 0.2 +res20priv +res20own ``` -As a rough check, if we multiply the total number of EPC addresses which fit the criteria by the mean ownership rate we get: +As a rough check, if we multiply the total number of EPC addresses which fit the criteria by the mean private ownership rate (owned + private rental = `r round(100*meanPrivate)`%) we get: - * 10% of households = `r round(res10)` - * 20% of households = `r round(res20)` + * 10% of households: + * Private: `r round(res10priv)` + * Owned: `r round(res10own)` + * 20% of households: + * Private: `r round(res20priv)` + * Owned: `r round(res20own)` -Note that this is only a rough guide since it assumes that LSOAs have exactly the same number of addresses/households in them. +Note that these values have not been uprated to correct for missing EPCs and that they are only a rough guide since they assume that LSOAs have exactly the same number of addresses/households in them. # Annexes diff --git a/isleOfWight/beisIoWData.html b/isleOfWight/beisIoWData.html index ed249ad606a8368e56a68f858300495d2317a066..531dd8164bcaf10da80657a7533363578e6506e3 100644 --- a/isleOfWight/beisIoWData.html +++ b/isleOfWight/beisIoWData.html @@ -1592,7 +1592,7 @@ div.tocify { <h1 class="title toc-ignore">Exploring Isle of Wight Data</h1> <h3 class="subtitle">Code and notes</h3> <h4 class="author">Ben Anderson (<a href="mailto:b.anderson@soton.ac.uk" class="email">b.anderson@soton.ac.uk</a>), <a href="http://www.energy.soton.ac.uk/">SERG, Energy & Climate Change</a>, University of Southampton</h4> -<h4 class="date">Last run at: 2020-04-23 23:06:26</h4> +<h4 class="date">Last run at: 2020-04-24 12:09:42</h4> </div> @@ -2362,8 +2362,11 @@ nrow(dtMsoaSpineIoW[elecRatio > 0.80])</code></pre> </div> <div id="ellis-question" class="section level1"> <h1><span class="header-section-number">3</span> Ellis’ question…</h1> +<p>For the Isle of Wight, can we estimate the number of private (owned or rented) addresses in high income deprivation areas which do not have gas?</p> +<p>Our method:</p> <ol style="list-style-type: decimal"> -<li>join income deprivation score to LSOA shapefiles. If you could produce a map that would be great.</li> +<li>Join LSOA data to EPC data to +join IMD income deprivation score to LSOA.</li> <li>Count dwellings in LSOAs with have EPC < C (D-G) and main fuel type not mains gas. If you could map all EPCs and then the EPCs that match the criteria.</li> <li>Join tenure to LSOA shapefiles. Multiply count by LSOA level census figure for % of private dwellings.</li> <li>Multiply by 1/0.687 to correct for dwellings without an EPC.</li> @@ -3025,9 +3028,12 @@ notHighDep </table> <p>Table <a href="#tab:depr10pc">3.3</a> shows the same results but for the 10% most income deprived LSOAs. As we can see there are 9 LSOAs classed as highly income deprived using this threshold.</p> <p>Next we assume that the EPC addresses of interest are as likely to be social/private rented as owner-occupied and multiply the EPC address counts by the Census-derived household ownership tenure rate for each LSOA. Without knowing the tensure status of the EPC we cannot do anything else.</p> -<pre class="r"><code>lsoaDT[, nHHsSelected := selectedEPCAddresses*`Owned (%)`] +<pre class="r"><code># multiply by tenure needed +lsoaDT[, nEPCAddressessSelectedOwned := selectedEPCAddresses*`Owned (%)`] # just owned +lsoaDT[, nEPCAddressessSelectedPrivate := selectedEPCAddresses*(`Owned (%)` + `Private rented (%)`)] # just owned -t <- lsoaDT[, .(nHHsSelected = sum(nHHsSelected), +t <- lsoaDT[, .(nEPCAddressessSelectedOwned = sum(nEPCAddressessSelectedOwned), + nEPCAddressessSelectedPrivate = sum(nEPCAddressessSelectedPrivate), nHHsCensus = sum(nHHsCensus), nLSOAs = .N), keyby = .(deprived20pc)] @@ -3044,7 +3050,10 @@ kableExtra::kable(t, digits = 2, deprived20pc </th> <th style="text-align:right;"> -nHHsSelected +nEPCAddressessSelectedOwned +</th> +<th style="text-align:right;"> +nEPCAddressessSelectedPrivate </th> <th style="text-align:right;"> nHHsCensus @@ -3063,6 +3072,9 @@ NA 158.58 </td> <td style="text-align:right;"> +177.13 +</td> +<td style="text-align:right;"> 2288 </td> <td style="text-align:right;"> @@ -3077,6 +3089,9 @@ highDep 821.68 </td> <td style="text-align:right;"> +1340.95 +</td> +<td style="text-align:right;"> 11475 </td> <td style="text-align:right;"> @@ -3091,6 +3106,9 @@ notHighDep 5297.50 </td> <td style="text-align:right;"> +6456.72 +</td> +<td style="text-align:right;"> 46216 </td> <td style="text-align:right;"> @@ -3099,8 +3117,10 @@ notHighDep </tr> </tbody> </table> -<pre class="r"><code>res20pc <- round(t[deprived20pc == "highDep", nHHsSelected])</code></pre> -<pre class="r"><code>t <- lsoaDT[, .(nHHsSelected = sum(nHHsSelected), +<pre class="r"><code>res20pcOwned <- round(t[deprived20pc == "highDep", nEPCAddressessSelectedOwned]) +res20pcPrivate <- round(t[deprived20pc == "highDep", nEPCAddressessSelectedPrivate])</code></pre> +<pre class="r"><code>t <- lsoaDT[, .(nEPCAddressessSelectedOwned = sum(nEPCAddressessSelectedOwned), + nEPCAddressessSelectedPrivate = sum(nEPCAddressessSelectedPrivate), nHHsCensus = sum(nHHsCensus), nLSOAs = .N), keyby = .(deprived10pc)] @@ -3117,7 +3137,10 @@ kableExtra::kable(t, digits = 2, deprived10pc </th> <th style="text-align:right;"> -nHHsSelected +nEPCAddressessSelectedOwned +</th> +<th style="text-align:right;"> +nEPCAddressessSelectedPrivate </th> <th style="text-align:right;"> nHHsCensus @@ -3136,6 +3159,9 @@ NA 158.58 </td> <td style="text-align:right;"> +177.13 +</td> +<td style="text-align:right;"> 2288 </td> <td style="text-align:right;"> @@ -3150,6 +3176,9 @@ highDep 496.50 </td> <td style="text-align:right;"> +847.32 +</td> +<td style="text-align:right;"> 6066 </td> <td style="text-align:right;"> @@ -3164,6 +3193,9 @@ notHighDep 5622.68 </td> <td style="text-align:right;"> +6950.34 +</td> +<td style="text-align:right;"> 51625 </td> <td style="text-align:right;"> @@ -3172,18 +3204,27 @@ notHighDep </tr> </tbody> </table> -<pre class="r"><code>res10pc <- round(t[deprived10pc == "highDep", nHHsSelected]) - +<pre class="r"><code>res10pcOwned <- round(t[deprived10pc == "highDep", nEPCAddressessSelectedOwned]) +res10pcPrivate <- round(t[deprived10pc == "highDep", nEPCAddressessSelectedPrivate]) epcMissingPC <- 100*(1-(1/epcMissingFactor))</code></pre> <blockquote> -<p>In summary using the EPC data we have estimated the number of addresses in low income areas which are owned and do not have gas to be:</p> +<p>In summary using the EPC data we have estimated the number of addresses in high income deprivation areas which do not have gas to be:</p> </blockquote> <ul> -<li>10% most deprived: 496 or, if we uprate to allow for the 31.3 % missing EPCs, 722;</li> -<li>20% most deprived: 822 or, if we uprate to allow for the 31.3 % missing EPCs, 1197;</li> +<li>10% most deprived: +<ul> +<li>Owned: 496 or, if we uprate to allow for the 31.3 % missing EPCs, 722</li> +<li>Private (owned & private rent): 847 or, if we uprate to allow for the 31.3 % missing EPCs, 1233</li> +</ul></li> +<li>20% most deprived: +<ul> +<li>Owned: 822 or, if we uprate to allow for the 31.3 % missing EPCs, 1197</li> +<li>Private (owned & private rent): 1341 or, if we uprate to allow for the 31.3 % missing EPCs, 1952</li> +</ul></li> </ul> <p>Table <a href="#tab:lsoaTable">3.6</a> shows the LSOA level results.</p> -<pre class="r"><code>t<- lsoaDT[, .(lsoaCode, nAddressesEPC = nHHsEPC, notSelectedEPCAddresses, selectedEPCAddresses, pcOwned=pcOwned*100, deprived20pc, deprived10pc,nHHsSelected)] +<pre class="r"><code>t<- lsoaDT[, .(lsoaCode, nAddressesEPC = nHHsEPC, notSelectedEPCAddresses, selectedEPCAddresses, pcOwned=pcOwned*100, deprived20pc, deprived10pc, + nEPCAddressessSelectedOwned,nEPCAddressessSelectedPrivate)] setkey(t, deprived20pc, deprived10pc) kableExtra::kable(t, caption = "Final results at LSOA level", @@ -3217,7 +3258,10 @@ deprived20pc deprived10pc </th> <th style="text-align:right;"> -nHHsSelected +nEPCAddressessSelectedOwned +</th> +<th style="text-align:right;"> +nEPCAddressessSelectedPrivate </th> </tr> </thead> @@ -3247,6 +3291,9 @@ NA <td style="text-align:right;"> 32.56 </td> +<td style="text-align:right;"> +37.24 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3273,6 +3320,9 @@ NA <td style="text-align:right;"> 13.74 </td> +<td style="text-align:right;"> +16.62 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3299,6 +3349,9 @@ NA <td style="text-align:right;"> 112.28 </td> +<td style="text-align:right;"> +123.28 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3325,6 +3378,9 @@ highDep <td style="text-align:right;"> 38.58 </td> +<td style="text-align:right;"> +52.68 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3351,6 +3407,9 @@ highDep <td style="text-align:right;"> 66.76 </td> +<td style="text-align:right;"> +120.89 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3377,6 +3436,9 @@ highDep <td style="text-align:right;"> 33.23 </td> +<td style="text-align:right;"> +39.93 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3403,6 +3465,9 @@ highDep <td style="text-align:right;"> 10.48 </td> +<td style="text-align:right;"> +12.28 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3429,6 +3494,9 @@ highDep <td style="text-align:right;"> 109.37 </td> +<td style="text-align:right;"> +266.59 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3455,6 +3523,9 @@ highDep <td style="text-align:right;"> 51.96 </td> +<td style="text-align:right;"> +66.87 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3481,6 +3552,9 @@ highDep <td style="text-align:right;"> 44.99 </td> +<td style="text-align:right;"> +69.14 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3507,6 +3581,9 @@ highDep <td style="text-align:right;"> 110.55 </td> +<td style="text-align:right;"> +181.22 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3533,6 +3610,9 @@ highDep <td style="text-align:right;"> 30.58 </td> +<td style="text-align:right;"> +37.72 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3559,6 +3639,9 @@ notHighDep <td style="text-align:right;"> 26.01 </td> +<td style="text-align:right;"> +33.89 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3585,6 +3668,9 @@ notHighDep <td style="text-align:right;"> 8.30 </td> +<td style="text-align:right;"> +10.34 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3611,6 +3697,9 @@ notHighDep <td style="text-align:right;"> 76.87 </td> +<td style="text-align:right;"> +93.16 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3637,6 +3726,9 @@ notHighDep <td style="text-align:right;"> 70.18 </td> +<td style="text-align:right;"> +142.73 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3663,6 +3755,9 @@ notHighDep <td style="text-align:right;"> 52.61 </td> +<td style="text-align:right;"> +79.78 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3689,6 +3784,9 @@ notHighDep <td style="text-align:right;"> 16.96 </td> +<td style="text-align:right;"> +19.90 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3715,6 +3813,9 @@ notHighDep <td style="text-align:right;"> 14.54 </td> +<td style="text-align:right;"> +18.12 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3741,6 +3842,9 @@ notHighDep <td style="text-align:right;"> 59.71 </td> +<td style="text-align:right;"> +95.69 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3767,6 +3871,9 @@ notHighDep <td style="text-align:right;"> 20.55 </td> +<td style="text-align:right;"> +24.51 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3793,6 +3900,9 @@ notHighDep <td style="text-align:right;"> 88.79 </td> +<td style="text-align:right;"> +104.03 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3819,6 +3929,9 @@ notHighDep <td style="text-align:right;"> 46.20 </td> +<td style="text-align:right;"> +53.30 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3845,6 +3958,9 @@ notHighDep <td style="text-align:right;"> 70.03 </td> +<td style="text-align:right;"> +78.56 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3871,6 +3987,9 @@ notHighDep <td style="text-align:right;"> 36.30 </td> +<td style="text-align:right;"> +41.31 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3897,6 +4016,9 @@ notHighDep <td style="text-align:right;"> 14.41 </td> +<td style="text-align:right;"> +15.87 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3923,6 +4045,9 @@ notHighDep <td style="text-align:right;"> 45.26 </td> +<td style="text-align:right;"> +54.31 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3949,6 +4074,9 @@ notHighDep <td style="text-align:right;"> 60.03 </td> +<td style="text-align:right;"> +68.22 +</td> </tr> <tr> <td style="text-align:left;"> @@ -3975,6 +4103,9 @@ notHighDep <td style="text-align:right;"> 276.76 </td> +<td style="text-align:right;"> +316.30 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4001,6 +4132,9 @@ notHighDep <td style="text-align:right;"> 300.36 </td> +<td style="text-align:right;"> +322.61 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4027,6 +4161,9 @@ notHighDep <td style="text-align:right;"> 11.41 </td> +<td style="text-align:right;"> +13.91 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4053,6 +4190,9 @@ notHighDep <td style="text-align:right;"> 6.00 </td> +<td style="text-align:right;"> +6.76 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4079,6 +4219,9 @@ notHighDep <td style="text-align:right;"> 78.31 </td> +<td style="text-align:right;"> +91.38 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4105,6 +4248,9 @@ notHighDep <td style="text-align:right;"> 170.82 </td> +<td style="text-align:right;"> +191.83 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4131,6 +4277,9 @@ notHighDep <td style="text-align:right;"> 255.00 </td> +<td style="text-align:right;"> +290.04 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4157,6 +4306,9 @@ notHighDep <td style="text-align:right;"> 327.59 </td> +<td style="text-align:right;"> +364.62 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4183,6 +4335,9 @@ notHighDep <td style="text-align:right;"> 255.23 </td> +<td style="text-align:right;"> +298.22 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4209,6 +4364,9 @@ notHighDep <td style="text-align:right;"> 157.15 </td> +<td style="text-align:right;"> +230.89 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4235,6 +4393,9 @@ notHighDep <td style="text-align:right;"> 33.56 </td> +<td style="text-align:right;"> +36.77 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4261,6 +4422,9 @@ notHighDep <td style="text-align:right;"> 11.91 </td> +<td style="text-align:right;"> +13.36 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4287,6 +4451,9 @@ notHighDep <td style="text-align:right;"> 28.31 </td> +<td style="text-align:right;"> +38.97 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4313,6 +4480,9 @@ notHighDep <td style="text-align:right;"> 19.52 </td> +<td style="text-align:right;"> +22.86 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4339,6 +4509,9 @@ notHighDep <td style="text-align:right;"> 30.09 </td> +<td style="text-align:right;"> +48.15 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4365,6 +4538,9 @@ notHighDep <td style="text-align:right;"> 27.59 </td> +<td style="text-align:right;"> +34.91 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4391,6 +4567,9 @@ notHighDep <td style="text-align:right;"> 63.44 </td> +<td style="text-align:right;"> +75.44 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4417,6 +4596,9 @@ notHighDep <td style="text-align:right;"> 38.20 </td> +<td style="text-align:right;"> +48.69 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4443,6 +4625,9 @@ notHighDep <td style="text-align:right;"> 35.85 </td> +<td style="text-align:right;"> +50.14 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4469,6 +4654,9 @@ notHighDep <td style="text-align:right;"> 38.47 </td> +<td style="text-align:right;"> +42.08 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4495,6 +4683,9 @@ notHighDep <td style="text-align:right;"> 35.13 </td> +<td style="text-align:right;"> +46.03 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4521,6 +4712,9 @@ notHighDep <td style="text-align:right;"> 52.46 </td> +<td style="text-align:right;"> +69.04 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4547,6 +4741,9 @@ notHighDep <td style="text-align:right;"> 44.90 </td> +<td style="text-align:right;"> +50.63 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4573,6 +4770,9 @@ notHighDep <td style="text-align:right;"> 36.18 </td> +<td style="text-align:right;"> +42.52 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4599,6 +4799,9 @@ notHighDep <td style="text-align:right;"> 17.68 </td> +<td style="text-align:right;"> +20.69 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4625,6 +4828,9 @@ notHighDep <td style="text-align:right;"> 20.50 </td> +<td style="text-align:right;"> +22.21 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4651,6 +4857,9 @@ notHighDep <td style="text-align:right;"> 97.16 </td> +<td style="text-align:right;"> +107.49 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4677,6 +4886,9 @@ notHighDep <td style="text-align:right;"> 39.98 </td> +<td style="text-align:right;"> +42.51 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4703,6 +4915,9 @@ notHighDep <td style="text-align:right;"> 40.79 </td> +<td style="text-align:right;"> +46.05 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4729,6 +4944,9 @@ notHighDep <td style="text-align:right;"> 51.64 </td> +<td style="text-align:right;"> +69.61 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4755,6 +4973,9 @@ notHighDep <td style="text-align:right;"> 33.09 </td> +<td style="text-align:right;"> +40.36 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4781,6 +5002,9 @@ notHighDep <td style="text-align:right;"> 110.38 </td> +<td style="text-align:right;"> +119.88 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4807,6 +5031,9 @@ notHighDep <td style="text-align:right;"> 20.89 </td> +<td style="text-align:right;"> +25.47 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4833,6 +5060,9 @@ notHighDep <td style="text-align:right;"> 38.79 </td> +<td style="text-align:right;"> +44.19 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4859,6 +5089,9 @@ notHighDep <td style="text-align:right;"> 74.11 </td> +<td style="text-align:right;"> +122.50 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4885,6 +5118,9 @@ notHighDep <td style="text-align:right;"> 126.37 </td> +<td style="text-align:right;"> +234.64 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4911,6 +5147,9 @@ notHighDep <td style="text-align:right;"> 31.74 </td> +<td style="text-align:right;"> +37.72 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4937,6 +5176,9 @@ notHighDep <td style="text-align:right;"> 46.08 </td> +<td style="text-align:right;"> +57.72 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4963,6 +5205,9 @@ notHighDep <td style="text-align:right;"> 35.44 </td> +<td style="text-align:right;"> +45.03 +</td> </tr> <tr> <td style="text-align:left;"> @@ -4989,6 +5234,9 @@ notHighDep <td style="text-align:right;"> 66.13 </td> +<td style="text-align:right;"> +77.56 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5015,6 +5263,9 @@ notHighDep <td style="text-align:right;"> 63.21 </td> +<td style="text-align:right;"> +104.68 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5041,6 +5292,9 @@ notHighDep <td style="text-align:right;"> 46.16 </td> +<td style="text-align:right;"> +56.11 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5067,6 +5321,9 @@ notHighDep <td style="text-align:right;"> 42.64 </td> +<td style="text-align:right;"> +63.83 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5093,6 +5350,9 @@ notHighDep <td style="text-align:right;"> 74.90 </td> +<td style="text-align:right;"> +117.99 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5119,6 +5379,9 @@ notHighDep <td style="text-align:right;"> 152.88 </td> +<td style="text-align:right;"> +181.38 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5145,6 +5408,9 @@ notHighDep <td style="text-align:right;"> 15.87 </td> +<td style="text-align:right;"> +16.88 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5171,6 +5437,9 @@ notHighDep <td style="text-align:right;"> 186.31 </td> +<td style="text-align:right;"> +196.41 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5197,6 +5466,9 @@ notHighDep <td style="text-align:right;"> 84.04 </td> +<td style="text-align:right;"> +106.64 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5223,6 +5495,9 @@ notHighDep <td style="text-align:right;"> 72.95 </td> +<td style="text-align:right;"> +83.74 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5249,6 +5524,9 @@ notHighDep <td style="text-align:right;"> 44.57 </td> +<td style="text-align:right;"> +53.74 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5275,6 +5553,9 @@ notHighDep <td style="text-align:right;"> 31.97 </td> +<td style="text-align:right;"> +34.13 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5301,6 +5582,9 @@ notHighDep <td style="text-align:right;"> 75.12 </td> +<td style="text-align:right;"> +94.98 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5327,6 +5611,9 @@ notHighDep <td style="text-align:right;"> 110.49 </td> +<td style="text-align:right;"> +188.60 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5353,6 +5640,9 @@ notHighDep <td style="text-align:right;"> 82.01 </td> +<td style="text-align:right;"> +92.98 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5379,6 +5669,9 @@ notHighDep <td style="text-align:right;"> 59.91 </td> +<td style="text-align:right;"> +78.46 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5405,6 +5698,9 @@ notHighDep <td style="text-align:right;"> 95.10 </td> +<td style="text-align:right;"> +119.06 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5431,6 +5727,9 @@ notHighDep <td style="text-align:right;"> 208.24 </td> +<td style="text-align:right;"> +243.07 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5457,6 +5756,9 @@ notHighDep <td style="text-align:right;"> 90.32 </td> +<td style="text-align:right;"> +101.52 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5483,6 +5785,9 @@ notHighDep <td style="text-align:right;"> 30.10 </td> +<td style="text-align:right;"> +36.81 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5509,6 +5814,9 @@ notHighDep <td style="text-align:right;"> 111.74 </td> +<td style="text-align:right;"> +128.84 +</td> </tr> <tr> <td style="text-align:left;"> @@ -5535,21 +5843,45 @@ notHighDep <td style="text-align:right;"> 52.40 </td> +<td style="text-align:right;"> +56.97 +</td> </tr> </tbody> </table> -<pre class="r"><code>totSelectedHH <- sum(lsoaDT[!is.na(lsoaCode)]$selectedEPCAddresses) - -meanOwnership <- mean(lsoaDT$pcOwned) - -res10 <- totSelectedHH * meanOwnership * 0.1 -res20 <- totSelectedHH * meanOwnership * 0.2</code></pre> -<p>As a rough check, if we multiply the total number of EPC addresses which fit the criteria by the mean ownership rate we get:</p> +<pre class="r"><code>data.table::fwrite(t, file = "lsoaTableForEllis.csv")</code></pre> +<pre class="r"><code>totPrivSelected <- sum(lsoaDT[!is.na(lsoaCode)]$nEPCAddressessSelectedPrivate) +totOwnedSelected <- sum(lsoaDT[!is.na(lsoaCode)]$nEPCAddressessSelectedOwned) + +meanPrivate <- mean(lsoaDT$pcOwned + lsoaDT$`Private rented (%)`) + +res10own <- totOwnedSelected * meanPrivate * 0.1 +res20own <- totOwnedSelected * meanPrivate * 0.2 +res10priv <- totPrivSelected * meanPrivate * 0.1 +res20priv <- totPrivSelected * meanPrivate * 0.2 + +res10priv</code></pre> +<pre><code>## [1] 708.1836</code></pre> +<pre class="r"><code>res10own</code></pre> +<pre><code>## [1] 557.4818</code></pre> +<pre class="r"><code>res20priv</code></pre> +<pre><code>## [1] 1416.367</code></pre> +<pre class="r"><code>res20own</code></pre> +<pre><code>## [1] 1114.964</code></pre> +<p>As a rough check, if we multiply the total number of EPC addresses which fit the criteria by the mean private ownership rate (owned + private rental = 89%) we get:</p> +<ul> +<li>10% of households: <ul> -<li>10% of households = 635</li> -<li>20% of households = 1269</li> +<li>Private: 708</li> +<li>Owned: 557</li> +</ul></li> +<li>20% of households: +<ul> +<li>Private: 1416</li> +<li>Owned: 1115</li> +</ul></li> </ul> -<p>Note that this is only a rough guide since it assumes that LSOAs have exactly the same number of addresses/households in them.</p> +<p>Note that these values have not been uprated to correct for missing EPCs and that they are only a rough guide since they assume that LSOAs have exactly the same number of addresses/households in them.</p> </div> <div id="annexes" class="section level1"> <h1><span class="header-section-number">4</span> Annexes</h1> @@ -8927,7 +9259,7 @@ kwhNDMedianGas </div> <div id="runtime" class="section level1"> <h1><span class="header-section-number">5</span> Runtime</h1> -<p>Analysis completed in 54.74 seconds ( 0.91 minutes) using <a href="https://cran.r-project.org/package=knitr">knitr</a> in <a href="http://www.rstudio.com">RStudio</a> with R version 3.6.3 (2020-02-29) running on x86_64-apple-darwin15.6.0.</p> +<p>Analysis completed in 51.23 seconds ( 0.85 minutes) using <a href="https://cran.r-project.org/package=knitr">knitr</a> in <a href="http://www.rstudio.com">RStudio</a> with R version 3.6.3 (2020-02-29) running on x86_64-apple-darwin15.6.0.</p> </div> <div id="r-environment" class="section level1"> <h1><span class="header-section-number">6</span> R environment</h1> diff --git a/isleOfWight/lsoaTableForEllis.csv b/isleOfWight/lsoaTableForEllis.csv index ac977780baa0136f426a81d41899c9519870d298..3a2e52607623a34adec20efe97c081dddd65e4df 100644 --- a/isleOfWight/lsoaTableForEllis.csv +++ b/isleOfWight/lsoaTableForEllis.csv @@ -1,90 +1,90 @@ -lsoaCode,nAddressesEPC,notSelectedEPCAddresses,selectedEPCAddresses,pcOwned,deprived20pc,deprived10pc,nHHsSelected -E01033229,615,576,39,83.4890966,,,32.560747674 -E01033237,566,531,35,39.2523364,,,13.73831774 -E01033239,904,778,126,89.1089109,,,112.277227734 -E01017322,558,473,85,45.3936348,highDep,highDep,38.58458958 -E01017328,597,458,139,48.030303,highDep,highDep,66.76212117 -E01017331,589,531,58,57.2948328,highDep,highDep,33.231003024 -E01017332,710,688,22,47.6190476,highDep,highDep,10.476190472 -E01017336,905,622,283,38.647343,highDep,highDep,109.37198069 -E01017343,437,348,89,58.377425,highDep,highDep,51.95590825 -E01017345,753,636,117,38.4520885,highDep,highDep,44.988943545 -E01017363,887,667,220,50.25,highDep,highDep,110.55 -E01017366,455,406,49,62.3987034,highDep,highDep,30.575364666 -E01017295,356,315,41,63.4328358,highDep,notHighDep,26.007462678 -E01017318,320,306,14,59.2592593,highDep,notHighDep,8.296296302 -E01017319,533,421,112,68.6363636,highDep,notHighDep,76.872727232 -E01017326,815,636,179,39.2086331,highDep,notHighDep,70.183453249 -E01017340,736,645,91,57.8144854,highDep,notHighDep,52.611181714 -E01017341,411,384,27,62.8304821,highDep,notHighDep,16.964230167 -E01017350,559,536,23,63.2198953,highDep,notHighDep,14.540575919 -E01017356,672,565,107,55.8035714,highDep,notHighDep,59.709821398 -E01017282,458,430,28,73.3788396,notHighDep,notHighDep,20.546075088 -E01017283,506,397,109,81.4589666,notHighDep,notHighDep,88.790273594 -E01017284,714,660,54,85.546875,notHighDep,notHighDep,46.1953125 -E01017285,664,577,87,80.490524,notHighDep,notHighDep,70.02675588 -E01017286,457,414,43,84.4250364,notHighDep,notHighDep,36.302765652 -E01017287,439,423,16,90.0815217,notHighDep,notHighDep,14.413043472 -E01017288,489,432,57,79.4066318,notHighDep,notHighDep,45.261780126 -E01017289,650,571,79,75.9864713,notHighDep,notHighDep,60.029312327 -E01017290,366,34,332,83.3625219,notHighDep,notHighDep,276.763572708 -E01017291,412,69,343,87.5675676,notHighDep,notHighDep,300.356756868 -E01017292,549,534,15,76.0432767,notHighDep,notHighDep,11.406491505 -E01017293,471,464,7,85.7142857,notHighDep,notHighDep,5.999999999 -E01017294,414,320,94,83.3075734,notHighDep,notHighDep,78.309118996 -E01017296,330,130,200,85.4085603,notHighDep,notHighDep,170.8171206 -E01017297,485,161,324,78.7021631,notHighDep,notHighDep,254.995008444 -E01017298,427,47,380,86.2068966,notHighDep,notHighDep,327.58620708 -E01017299,423,65,358,71.2945591,notHighDep,notHighDep,255.234521578 -E01017300,1324,1079,245,64.144454,notHighDep,notHighDep,157.1539123 -E01017301,277,240,37,90.7024793,notHighDep,notHighDep,33.559917341 -E01017302,462,448,14,85.0609756,notHighDep,notHighDep,11.908536584 -E01017303,521,481,40,70.7720588,notHighDep,notHighDep,28.30882352 -E01017304,388,362,26,75.0902527,notHighDep,notHighDep,19.523465702 -E01017305,717,663,54,55.7275542,notHighDep,notHighDep,30.092879268 -E01017306,360,318,42,65.6880734,notHighDep,notHighDep,27.588990828 -E01017307,815,737,78,81.3348416,notHighDep,notHighDep,63.441176448 -E01017308,598,538,60,63.6612022,notHighDep,notHighDep,38.19672132 -E01017309,543,489,54,66.3805436,notHighDep,notHighDep,35.845493544 -E01017310,349,306,43,89.4736842,notHighDep,notHighDep,38.473684206 -E01017311,439,387,52,67.5637394,notHighDep,notHighDep,35.133144488 -E01017312,755,674,81,64.7606383,notHighDep,notHighDep,52.456117023 -E01017315,443,373,70,64.1369048,notHighDep,notHighDep,44.89583336 -E01017316,747,703,44,82.2360248,notHighDep,notHighDep,36.183850912 -E01017317,428,405,23,76.8768769,notHighDep,notHighDep,17.681681687 -E01017320,344,321,23,89.1304348,notHighDep,notHighDep,20.500000004 -E01017321,416,301,115,84.4897959,notHighDep,notHighDep,97.163265285 -E01017323,358,314,44,90.8548708,notHighDep,notHighDep,39.976143152 -E01017324,413,364,49,83.2520325,notHighDep,notHighDep,40.793495925 -E01017325,483,391,92,56.1258278,notHighDep,notHighDep,51.635761576 -E01017327,443,398,45,73.5426009,notHighDep,notHighDep,33.094170405 -E01017329,569,442,127,86.9158879,notHighDep,notHighDep,110.383177633 -E01017333,391,364,27,77.3700306,notHighDep,notHighDep,20.889908262 -E01017334,319,268,51,76.059322,notHighDep,notHighDep,38.79025422 -E01017335,627,492,135,54.8941799,notHighDep,notHighDep,74.107142865 -E01017337,740,496,244,51.7902813,notHighDep,notHighDep,126.368286372 -E01017338,475,435,40,79.3495935,notHighDep,notHighDep,31.7398374 -E01017339,497,435,62,74.3243243,notHighDep,notHighDep,46.081081066 -E01017342,560,509,51,69.4986072,notHighDep,notHighDep,35.444289672 -E01017344,464,384,80,82.6629681,notHighDep,notHighDep,66.13037448 -E01017346,623,510,113,55.9345157,notHighDep,notHighDep,63.206002741 -E01017347,495,431,64,72.1288515,notHighDep,notHighDep,46.16246496 -E01017348,486,417,69,61.8042226,notHighDep,notHighDep,42.644913594 -E01017349,783,663,120,62.4161074,notHighDep,notHighDep,74.89932888 -E01017351,708,520,188,81.3169985,notHighDep,notHighDep,152.87595718 -E01017352,434,417,17,93.3566434,notHighDep,notHighDep,15.870629378 -E01017353,353,138,215,86.6541353,notHighDep,notHighDep,186.306390895 -E01017354,546,426,120,70.0342466,notHighDep,notHighDep,84.04109592 -E01017355,504,415,89,81.969697,notHighDep,notHighDep,72.95303033 -E01017357,472,402,70,63.6650869,notHighDep,notHighDep,44.56556083 -E01017358,425,387,38,84.137931,notHighDep,notHighDep,31.97241378 -E01017359,648,550,98,76.6519824,notHighDep,notHighDep,75.118942752 -E01017360,892,695,197,56.0853199,notHighDep,notHighDep,110.488080203 -E01017361,545,447,98,83.68,notHighDep,notHighDep,82.0064 -E01017362,686,599,87,68.8650307,notHighDep,notHighDep,59.912576709 -E01017364,632,508,124,76.6961652,notHighDep,notHighDep,95.103244848 -E01017365,731,474,257,81.0271041,notHighDep,notHighDep,208.239657537 -E01017367,585,483,102,88.5478158,notHighDep,notHighDep,90.318772116 -E01017368,475,432,43,70,notHighDep,notHighDep,30.1 -E01017369,439,294,145,77.0597738,notHighDep,notHighDep,111.73667201 -E01017370,486,418,68,77.0595691,notHighDep,notHighDep,52.400506988 +lsoaCode,nAddressesEPC,notSelectedEPCAddresses,selectedEPCAddresses,pcOwned,deprived20pc,deprived10pc,nEPCAddressessSelectedOwned,nEPCAddressessSelectedPrivate +E01033229,615,576,39,83.4890966,,,32.560747674,37.238317779 +E01033237,566,531,35,39.2523364,,,13.73831774,16.616822425 +E01033239,904,778,126,89.1089109,,,112.277227734,123.27812784 +E01017322,558,473,85,45.3936348,highDep,highDep,38.58458958,52.68006699 +E01017328,597,458,139,48.030303,highDep,highDep,66.76212117,120.887878691 +E01017331,589,531,58,57.2948328,highDep,highDep,33.231003024,39.930091184 +E01017332,710,688,22,47.6190476,highDep,highDep,10.476190472,12.278095236 +E01017336,905,622,283,38.647343,highDep,highDep,109.37198069,266.594203038 +E01017343,437,348,89,58.377425,highDep,highDep,51.95590825,66.867724839 +E01017345,753,636,117,38.4520885,highDep,highDep,44.988943545,69.136363647 +E01017363,887,667,220,50.25,highDep,highDep,110.55,181.225 +E01017366,455,406,49,62.3987034,highDep,highDep,30.575364666,37.722852517 +E01017295,356,315,41,63.4328358,highDep,notHighDep,26.007462678,33.886194017 +E01017318,320,306,14,59.2592593,highDep,notHighDep,8.296296302,10.343079922 +E01017319,533,421,112,68.6363636,highDep,notHighDep,76.872727232,93.163636272 +E01017326,815,636,179,39.2086331,highDep,notHighDep,70.183453249,142.72781769 +E01017340,736,645,91,57.8144854,highDep,notHighDep,52.611181714,79.783989831 +E01017341,411,384,27,62.8304821,highDep,notHighDep,16.964230167,19.903576986 +E01017350,559,536,23,63.2198953,highDep,notHighDep,14.540575919,18.123036645 +E01017356,672,565,107,55.8035714,highDep,notHighDep,59.709821398,95.694940466 +E01017282,458,430,28,73.3788396,notHighDep,notHighDep,20.546075088,24.511945388 +E01017283,506,397,109,81.4589666,notHighDep,notHighDep,88.790273594,104.030395155 +E01017284,714,660,54,85.546875,notHighDep,notHighDep,46.1953125,53.296875018 +E01017285,664,577,87,80.490524,notHighDep,notHighDep,70.02675588,78.561872958 +E01017286,457,414,43,84.4250364,notHighDep,notHighDep,36.302765652,41.31004367 +E01017287,439,423,16,90.0815217,notHighDep,notHighDep,14.413043472,15.869565216 +E01017288,489,432,57,79.4066318,notHighDep,notHighDep,45.261780126,54.314136174 +E01017289,650,571,79,75.9864713,notHighDep,notHighDep,60.029312327,68.223224401 +E01017290,366,34,332,83.3625219,notHighDep,notHighDep,276.763572708,316.301225952 +E01017291,412,69,343,87.5675676,notHighDep,notHighDep,300.356756868,322.605405563 +E01017292,549,534,15,76.0432767,notHighDep,notHighDep,11.406491505,13.910355495 +E01017293,471,464,7,85.7142857,notHighDep,notHighDep,5.999999999,6.761467888 +E01017294,414,320,94,83.3075734,notHighDep,notHighDep,78.309118996,91.384853166 +E01017296,330,130,200,85.4085603,notHighDep,notHighDep,170.8171206,191.8287938 +E01017297,485,161,324,78.7021631,notHighDep,notHighDep,254.995008444,290.036605716 +E01017298,427,47,380,86.2068966,notHighDep,notHighDep,327.58620708,364.6176912 +E01017299,423,65,358,71.2945591,notHighDep,notHighDep,255.234521578,298.221388404 +E01017300,1324,1079,245,64.144454,notHighDep,notHighDep,157.1539123,230.88564065 +E01017301,277,240,37,90.7024793,notHighDep,notHighDep,33.559917341,36.770661161 +E01017302,462,448,14,85.0609756,notHighDep,notHighDep,11.908536584,13.359756102 +E01017303,521,481,40,70.7720588,notHighDep,notHighDep,28.30882352,38.97058824 +E01017304,388,362,26,75.0902527,notHighDep,notHighDep,19.523465702,22.855595672 +E01017305,717,663,54,55.7275542,notHighDep,notHighDep,30.092879268,48.148606818 +E01017306,360,318,42,65.6880734,notHighDep,notHighDep,27.588990828,34.910091762 +E01017307,815,737,78,81.3348416,notHighDep,notHighDep,63.441176448,75.44117646 +E01017308,598,538,60,63.6612022,notHighDep,notHighDep,38.19672132,48.6885246 +E01017309,543,489,54,66.3805436,notHighDep,notHighDep,35.845493544,50.137339014 +E01017310,349,306,43,89.4736842,notHighDep,notHighDep,38.473684206,42.077972715 +E01017311,439,387,52,67.5637394,notHighDep,notHighDep,35.133144488,46.033994344 +E01017312,755,674,81,64.7606383,notHighDep,notHighDep,52.456117023,69.043882977 +E01017315,443,373,70,64.1369048,notHighDep,notHighDep,44.89583336,50.62500002 +E01017316,747,703,44,82.2360248,notHighDep,notHighDep,36.183850912,42.524223588 +E01017317,428,405,23,76.8768769,notHighDep,notHighDep,17.681681687,20.6861862 +E01017320,344,321,23,89.1304348,notHighDep,notHighDep,20.500000004,22.20833333 +E01017321,416,301,115,84.4897959,notHighDep,notHighDep,97.163265285,107.489795855 +E01017323,358,314,44,90.8548708,notHighDep,notHighDep,39.976143152,42.512922496 +E01017324,413,364,49,83.2520325,notHighDep,notHighDep,40.793495925,46.052032502 +E01017325,483,391,92,56.1258278,notHighDep,notHighDep,51.635761576,69.609271472 +E01017327,443,398,45,73.5426009,notHighDep,notHighDep,33.094170405,40.358744415 +E01017329,569,442,127,86.9158879,notHighDep,notHighDep,110.383177633,119.878504718 +E01017333,391,364,27,77.3700306,notHighDep,notHighDep,20.889908262,25.472477079 +E01017334,319,268,51,76.059322,notHighDep,notHighDep,38.79025422,44.192796573 +E01017335,627,492,135,54.8941799,notHighDep,notHighDep,74.107142865,122.499999945 +E01017337,740,496,244,51.7902813,notHighDep,notHighDep,126.368286372,234.639386032 +E01017338,475,435,40,79.3495935,notHighDep,notHighDep,31.7398374,37.72357724 +E01017339,497,435,62,74.3243243,notHighDep,notHighDep,46.081081066,57.717717722 +E01017342,560,509,51,69.4986072,notHighDep,notHighDep,35.444289672,45.033426156 +E01017344,464,384,80,82.6629681,notHighDep,notHighDep,66.13037448,77.55894592 +E01017346,623,510,113,55.9345157,notHighDep,notHighDep,63.206002741,104.67530695 +E01017347,495,431,64,72.1288515,notHighDep,notHighDep,46.16246496,56.1120448 +E01017348,486,417,69,61.8042226,notHighDep,notHighDep,42.644913594,63.834932757 +E01017349,783,663,120,62.4161074,notHighDep,notHighDep,74.89932888,117.98657724 +E01017351,708,520,188,81.3169985,notHighDep,notHighDep,152.87595718,181.378254224 +E01017352,434,417,17,93.3566434,notHighDep,notHighDep,15.870629378,16.881118881 +E01017353,353,138,215,86.6541353,notHighDep,notHighDep,186.306390895,196.40977431 +E01017354,546,426,120,70.0342466,notHighDep,notHighDep,84.04109592,106.6438356 +E01017355,504,415,89,81.969697,notHighDep,notHighDep,72.95303033,83.740909099 +E01017357,472,402,70,63.6650869,notHighDep,notHighDep,44.56556083,53.74407584 +E01017358,425,387,38,84.137931,notHighDep,notHighDep,31.97241378,34.134482756 +E01017359,648,550,98,76.6519824,notHighDep,notHighDep,75.118942752,94.977973608 +E01017360,892,695,197,56.0853199,notHighDep,notHighDep,110.488080203,188.595984925 +E01017361,545,447,98,83.68,notHighDep,notHighDep,82.0064,92.9824 +E01017362,686,599,87,68.8650307,notHighDep,notHighDep,59.912576709,78.460122717 +E01017364,632,508,124,76.6961652,notHighDep,notHighDep,95.103244848,119.061946884 +E01017365,731,474,257,81.0271041,notHighDep,notHighDep,208.239657537,243.068473582 +E01017367,585,483,102,88.5478158,notHighDep,notHighDep,90.318772116,101.5182999 +E01017368,475,432,43,70,notHighDep,notHighDep,30.1,36.810606058 +E01017369,439,294,145,77.0597738,notHighDep,notHighDep,111.73667201,128.836833615 +E01017370,486,418,68,77.0595691,notHighDep,notHighDep,52.400506988,56.968314372