add final model plot to readme

Author: mpadge

README.Rmd

@@ -28,6 +28,8 @@ City](https://www1.nyc.gov/html/dot/html/about/datafeeds.shtml#Pedestrians).
 
 Final result looks like this:
 ```{r load, echo = FALSE}
+library (ggplot2)
+theme_set (theme_minimal ())
 devtools::load_all (".", export_all = FALSE)
 data_dir <- "/data/data/moveability/nyc"
 ```
@@ -48,3 +50,15 @@ The final flow model looks like this:
 
 ![](./man/figures/flowmap.png)
 
+And a final statistical relationship between modelled and observed pedestrian
+counts that looks like this:
+```{r model-plot}
+data_dir <- "/data/data/moveability/nyc"
+dat <- readRDS (file.path (data_dir, "ped-model-final.Rds"))
+mod <- lm (dat$p ~ dat$flowvars)
+res <- data.frame (predicted = predict (mod),
+                   actual = dat$p)
+ggplot (res, aes (x = predicted, y = actual)) +
+    geom_point () +
+    geom_smooth (method = "lm")
+```

README.md

@@ -5,8 +5,7 @@ Status](https://travis-ci.org/moveability/calibration.svg)](https://travis-ci.or
 [![Project Status: Concept - Minimal or no implementation has been done
 yet.](http://www.repostatus.org/badges/0.1.0/concept.svg)](http://www.repostatus.org/#concept)
 
-calibration
-===========
+# calibration
 
 <!-- badges: start -->
 
@@ -21,7 +20,7 @@ City](https://www1.nyc.gov/html/dot/html/about/datafeeds.shtml#Pedestrians).
 
 Final result looks like this:
 
-    #> Loading calibration
+    #> Loading nyped
 
 ``` r
 x <- calc_layer_scales (data_dir, model = "final")
@@ -31,22 +30,22 @@ x <- calc_layer_scales (data_dir, model = "final")
 knitr::kable (x)
 ```
 
-|     | layer   |      estimates|   full\_rel|    ped\_rel|     full\_abs|      ped\_abs|          t|          p|
-|-----|:--------|--------------:|-----------:|-----------:|-------------:|-------------:|----------:|----------:|
-| 13  | hea-dis |  105658.191105|   0.5719363|   0.3459949|  60429.755766|  36557.196583|   9.868915|  0.0000000|
-| 1   | sub-dis |      22.566344|   0.5219142|   0.5682277|     11.777695|     12.822822|   8.986039|  0.0000000|
-| 12  | sub-cen |     -10.087818|  -0.2347841|  -0.3304276|      2.368459|      3.333293|  -6.994975|  0.0000000|
-| 6   | sub-hea |       8.378134|   0.4338598|   0.4335501|      3.634935|      3.632341|   6.658314|  0.0000000|
-| 7   | sus-edu |   -5923.748211|  -0.6416411|  -0.5697332|   3800.920378|   3374.956231|  -6.056855|  0.0000000|
-| 8   | edu-hea |  -24921.480025|  -0.4563520|  -0.3789464|  11372.966534|   9443.905594|  -5.607032|  0.0000002|
-| 4   | edu-tra |   23977.299359|   0.3225789|   0.2620421|   7734.569685|   6283.060961|   5.347316|  0.0000006|
-| 2   | sub-tra |       6.444493|   0.2392262|   0.2839541|      1.541692|      1.829940|   5.080766|  0.0000018|
-| 9   | sus-res |    6258.235319|   0.3904308|   0.4362252|   2443.407704|   2729.999677|   5.080036|  0.0000018|
-| 11  | sus-sub |   -1336.884757|  -0.1903741|  -0.2170844|    254.508256|    290.216796|  -4.039428|  0.0001066|
-| 5   | sus-ent |    1445.795365|   0.2070682|   0.1704797|    299.378267|    246.478763|   4.002238|  0.0001220|
-| 3   | edu-dis |  -78056.981732|  -0.5872607|  -0.3261469|  45839.794420|  25458.042278|  -3.183290|  0.0019518|
-| 14  | ent-tra |   38179.006598|   0.1893925|   0.1540936|   7230.818735|   5883.141715|   3.176518|  0.0019936|
-| 10  | edu-sus |   16903.889029|   0.2340051|   0.1677712|   3955.596308|   2835.985535|   3.033825|  0.0030914|
+|    | layer   |      estimates |   full\_rel |    ped\_rel |    full\_abs |     ped\_abs |          t |         p |
+| -- | :------ | -------------: | ----------: | ----------: | -----------: | -----------: | ---------: | --------: |
+| 13 | hea-dis |  105658.191105 |   0.5719363 |   0.3459949 | 60429.755766 | 36557.196583 |   9.868915 | 0.0000000 |
+| 1  | sub-dis |      22.566344 |   0.5219142 |   0.5682277 |    11.777695 |    12.822822 |   8.986039 | 0.0000000 |
+| 12 | sub-cen |    \-10.087818 | \-0.2347841 | \-0.3304276 |     2.368459 |     3.333293 | \-6.994975 | 0.0000000 |
+| 6  | sub-hea |       8.378134 |   0.4338598 |   0.4335501 |     3.634935 |     3.632341 |   6.658314 | 0.0000000 |
+| 7  | sus-edu |  \-5923.748211 | \-0.6416411 | \-0.5697332 |  3800.920378 |  3374.956231 | \-6.056855 | 0.0000000 |
+| 8  | edu-hea | \-24921.480025 | \-0.4563520 | \-0.3789464 | 11372.966534 |  9443.905594 | \-5.607032 | 0.0000002 |
+| 4  | edu-tra |   23977.299359 |   0.3225789 |   0.2620421 |  7734.569685 |  6283.060961 |   5.347316 | 0.0000006 |
+| 2  | sub-tra |       6.444493 |   0.2392262 |   0.2839541 |     1.541692 |     1.829940 |   5.080766 | 0.0000018 |
+| 9  | sus-res |    6258.235319 |   0.3904308 |   0.4362252 |  2443.407704 |  2729.999677 |   5.080036 | 0.0000018 |
+| 11 | sus-sub |  \-1336.884757 | \-0.1903741 | \-0.2170844 |   254.508256 |   290.216796 | \-4.039428 | 0.0001066 |
+| 5  | sus-ent |    1445.795365 |   0.2070682 |   0.1704797 |   299.378267 |   246.478763 |   4.002238 | 0.0001220 |
+| 3  | edu-dis | \-78056.981732 | \-0.5872607 | \-0.3261469 | 45839.794420 | 25458.042278 | \-3.183290 | 0.0019518 |
+| 14 | ent-tra |   38179.006598 |   0.1893925 |   0.1540936 |  7230.818735 |  5883.141715 |   3.176518 | 0.0019936 |
+| 10 | edu-sus |   16903.889029 |   0.2340051 |   0.1677712 |  3955.596308 |  2835.985535 |   3.033825 | 0.0030914 |
 
 With the following model statistics:
 
@@ -90,3 +89,19 @@ summary (lm (dat$p ~ f_))
 The final flow model looks like this:
 
 ![](./man/figures/flowmap.png)
+
+And a final statistical relationship between modelled and observed
+pedestrian counts that looks like this:
+
+``` r
+data_dir <- "/data/data/moveability/nyc"
+dat <- readRDS (file.path (data_dir, "ped-model-final.Rds"))
+mod <- lm (dat$p ~ dat$flowvars)
+res <- data.frame (predicted = predict (mod),
+                   actual = dat$p)
+ggplot (res, aes (x = predicted, y = actual)) +
+    geom_point () +
+    geom_smooth (method = "lm")
+```
+
+<img src="man/figures/README-model-plot-1.png" width="100%" />