pollution <- read.csv2("../R Scripts/avgpm25.csv")
colClasses = c("numeric", "character", "factor", "numeric", "numeric")
summary(pollution$pm25)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 3.383 8.549 10.047 9.836 11.356 18.441
pollution
## ï..X pm25 fips region longitude latitude
## 1 NA 9.771185 1003 east -87.74826 30.59278
## 2 NA 9.993817 1027 east -85.84286 33.26581
## 3 NA 10.688618 1033 east -87.72596 34.73148
## 4 NA 11.337424 1049 east -85.79892 34.45913
## 5 NA 12.119764 1055 east -86.03212 34.01860
## 6 NA 10.827805 1069 east -85.35039 31.18973
## 7 NA 11.583928 1073 east -86.82805 33.52787
## 8 NA 11.261996 1089 east -86.58823 34.73079
## 9 NA 9.414423 1097 east -88.13967 30.72226
## 10 NA 11.391494 1103 east -86.91892 34.50702
## 11 NA 12.384795 1113 east -85.10110 32.37600
## 12 NA 10.649500 1117 east -86.69867 33.26912
## 13 NA 11.333821 1121 east -86.17828 33.36850
## 14 NA 12.302436 1125 east -87.51169 33.23560
## 15 NA 11.024508 1127 east -87.28541 33.81989
## 16 NA 10.835515 5001 east -91.42941 34.36000
## 17 NA 10.436235 5003 east -91.78535 33.18542
## 18 NA 11.114743 5035 east -90.27280 35.19770
## 19 NA 10.823157 5045 east -92.37994 35.11948
## 20 NA 10.745523 5051 east -93.07589 34.54531
## 21 NA 10.450351 5067 east -91.22758 35.61326
## 22 NA 10.723959 5107 east -90.76197 34.46837
## 23 NA 10.441737 5113 east -94.25363 34.48549
## 24 NA 11.097645 5115 east -93.08445 35.33048
## 25 NA 11.840086 5119 east -92.29447 34.76654
## 26 NA 11.105003 5131 east -94.35177 35.29392
## 27 NA 10.887024 5139 east -92.61279 33.20718
## 28 NA 10.696524 5143 east -94.19868 36.05138
## 29 NA 10.641453 5145 east -91.73056 35.25011
## 30 NA 9.821778 9001 east -73.35482 41.21290
## 31 NA 8.894819 9003 east -72.71545 41.76997
## 32 NA 7.029871 9005 east -73.21269 41.75426
## 33 NA 10.001609 9009 east -72.94032 41.39076
## 34 NA 8.821875 9011 east -72.09295 41.44611
## 35 NA 10.412954 10001 east -75.56027 39.09271
## 36 NA 11.551756 10003 east -75.61171 39.69561
## 37 NA 10.837718 10005 east -75.34227 38.65142
## 38 NA 11.023823 11001 east -77.01322 38.91361
## 39 NA 7.503461 12001 east -82.37995 29.67644
## 40 NA 6.877201 12009 east -80.68519 28.22008
## 41 NA 6.977536 12011 east -80.20610 26.13479
## 42 NA 7.375229 12017 east -82.46916 28.89464
## 43 NA 8.378378 12031 east -81.63887 30.32087
## 44 NA 8.961842 12033 east -87.27736 30.50204
## 45 NA 8.341545 12057 east -82.41208 27.96635
## 46 NA 6.744966 12071 east -81.84626 26.60254
## 47 NA 9.795967 12073 east -84.26089 30.46293
## 48 NA 5.000000 12086 east -80.30226 25.75643
## 49 NA 7.285535 12095 east -81.40416 28.54713
## 50 NA 6.348540 12099 east -80.17639 26.63774
## 51 NA 7.738196 12103 east -82.72777 27.88965
## 52 NA 7.679519 12105 east -81.75830 27.99244
## 53 NA 7.584909 12111 east -80.36087 27.34825
## 54 NA 6.800815 12115 east -82.39888 27.16493
## 55 NA 7.500355 12117 east -81.31083 28.70991
## 56 NA 7.213773 12127 east -81.13895 29.04904
## 57 NA 11.959208 13021 east -83.66521 32.82969
## 58 NA 10.824470 13051 east -81.10716 32.04756
## 59 NA 11.292279 13059 east -83.38011 33.95475
## 60 NA 12.913903 13063 east -84.36560 33.56618
## 61 NA 12.018416 13067 east -84.55928 33.93554
## 62 NA 12.145414 13089 east -84.25768 33.79094
## 63 NA 11.746439 13095 east -84.16321 31.57287
## 64 NA 11.817821 13121 east -84.41933 33.77960
## 65 NA 9.793812 13127 east -81.49629 31.18763
## 66 NA 12.145278 13135 east -84.05845 33.96706
## 67 NA 11.188334 13139 east -83.84296 34.29680
## 68 NA 11.422325 13153 east -83.66323 32.55806
## 69 NA 9.709880 13185 east -83.27718 30.84853
## 70 NA 12.200660 13215 east -84.94218 32.49006
## 71 NA 10.949059 13223 east -84.83112 33.90278
## 72 NA 12.392875 13245 east -82.02840 33.42949
## 73 NA 10.672707 13295 east -85.29183 34.81869
## 74 NA 11.358464 13303 east -82.78261 32.94936
## 75 NA 12.121381 13319 east -83.22660 32.83589
## 76 NA 9.353740 17001 east -91.26805 39.96545
## 77 NA 11.059005 17019 east -88.21675 40.13682
## 78 NA 11.878008 17031 east -87.76782 41.83765
## 79 NA 10.901309 17043 east -88.06359 41.85603
## 80 NA 11.269838 17065 east -88.53624 38.09678
## 81 NA 10.336829 17083 east -90.34214 39.08361
## 82 NA 10.581457 17089 east -88.33844 41.91823
## 83 NA 9.256495 17097 east -87.98380 42.31332
## 84 NA 10.275219 17099 east -88.90429 41.33216
## 85 NA 10.003862 17111 east -88.35377 42.29640
## 86 NA 10.469588 17113 east -88.88358 40.49525
## 87 NA 11.631926 17115 east -88.95323 39.85314
## 88 NA 12.702123 17119 east -90.01884 38.81251
## 89 NA 11.116070 17143 east -89.66506 40.75249
## 90 NA 10.104734 17157 east -89.79688 38.07921
## 91 NA 9.655263 17161 east -90.50753 41.49143
## 92 NA 12.331060 17163 east -90.00207 38.53268
## 93 NA 11.039065 17167 east -89.64650 39.76385
## 94 NA 10.505179 17197 east -88.02984 41.52743
## 95 NA 10.034725 17201 east -89.08163 42.31112
## 96 NA 13.093119 18003 east -85.10007 41.09028
## 97 NA 13.144659 18019 east -85.72746 38.38648
## 98 NA 11.855302 18035 east -85.39608 40.21189
## 99 NA 13.023923 18037 east -86.90055 38.35009
## 100 NA 12.046176 18039 east -85.89940 41.62609
## 101 NA 12.367141 18043 east -85.86329 38.31693
## 102 NA 11.630029 18051 east -87.55836 38.31770
## 103 NA 11.443252 18065 east -85.40599 39.92353
## 104 NA 11.692452 18067 east -86.12641 40.47544
## 105 NA 11.811781 18083 east -87.41894 38.69742
## 106 NA 12.753287 18089 east -87.39525 41.52960
## 107 NA 10.646946 18091 east -86.77926 41.59922
## 108 NA 12.089511 18095 east -85.71920 40.14594
## 109 NA 13.883217 18097 east -86.14181 39.79279
## 110 NA 12.197354 18127 east -87.08665 41.52351
## 111 NA 13.435254 18141 east -86.24193 41.66478
## 112 NA 12.224040 18147 east -87.01320 38.02437
## 113 NA 12.195055 18157 east -86.88088 40.40948
## 114 NA 12.828878 18163 east -87.55696 37.98989
## 115 NA 12.260916 18167 east -87.39127 39.45689
## 116 NA 10.473696 19013 east -92.34473 42.48759
## 117 NA 11.699045 19045 east -90.42740 41.87117
## 118 NA 10.917263 19103 east -91.56887 41.67491
## 119 NA 11.123664 19111 east -91.41661 40.57683
## 120 NA 10.407355 19113 east -91.63056 42.02999
## 121 NA 9.288892 19137 east -95.15925 41.01310
## 122 NA 12.679211 19139 east -91.08480 41.47246
## 123 NA 8.828504 19147 east -94.68872 43.09030
## 124 NA 9.402015 19153 east -93.63094 41.62794
## 125 NA 10.831346 19155 east -95.68040 41.30853
## 126 NA 11.960760 19163 east -90.57886 41.57213
## 127 NA 9.281105 19177 east -91.94881 40.74022
## 128 NA 8.938962 20091 east -94.75377 38.93218
## 129 NA 9.337790 20107 east -94.80532 38.22383
## 130 NA 9.057115 20173 east -97.36347 37.68133
## 131 NA 9.322162 20177 east -95.71026 39.03915
## 132 NA 8.681833 20191 east -97.43256 37.26461
## 133 NA 9.741182 20209 east -94.69276 39.10371
## 134 NA 11.281826 21013 east -83.70225 36.69139
## 135 NA 11.405537 21019 east -82.66458 38.42536
## 136 NA 12.665325 21029 east -85.67205 38.00518
## 137 NA 11.801188 21037 east -84.43330 39.03652
## 138 NA 9.685132 21043 east -83.04932 38.30785
## 139 NA 11.322465 21047 east -87.48199 36.85444
## 140 NA 12.051753 21059 east -87.10510 37.75065
## 141 NA 12.395794 21067 east -84.49464 38.02963
## 142 NA 11.103714 21073 east -84.86839 38.20365
## 143 NA 11.765029 21093 east -85.91493 37.73892
## 144 NA 11.859195 21101 east -87.56808 37.80262
## 145 NA 13.058814 21111 east -85.70402 38.20924
## 146 NA 11.903645 21117 east -84.53935 39.02674
## 147 NA 11.670030 21145 east -88.65111 37.05971
## 148 NA 10.343433 21151 east -84.27224 37.69968
## 149 NA 11.554662 21183 east -86.87936 37.45651
## 150 NA 10.090150 21195 east -82.41111 37.44779
## 151 NA 11.511463 21227 east -86.42455 36.98691
## 152 NA 10.851174 22017 east -93.81391 32.52756
## 153 NA 8.874310 22019 east -93.29176 30.23744
## 154 NA 9.970397 22033 east -91.11965 30.48012
## 155 NA 9.761818 22047 east -91.28719 30.27517
## 156 NA 9.079453 22051 east -90.15140 29.93524
## 157 NA 9.799970 22055 east -92.03888 30.20739
## 158 NA 9.703961 22073 east -92.11586 32.50938
## 159 NA 8.898718 22079 east -92.47659 31.26208
## 160 NA 12.178354 22087 east -89.89026 29.92209
## 161 NA 9.134486 22105 east -90.46295 30.60740
## 162 NA 8.528276 22109 east -90.73592 29.54117
## 163 NA 11.802678 22121 east -91.26305 30.45209
## 164 NA 7.794286 23001 east -70.21766 44.13507
## 165 NA 6.162443 23003 east -68.26166 46.63736
## 166 NA 8.182713 23005 east -70.32798 43.78448
## 167 NA 4.504539 23009 east -68.41571 44.51274
## 168 NA 7.334426 23011 east -69.78455 44.39996
## 169 NA 8.259664 23017 east -70.70517 44.32452
## 170 NA 7.384147 23019 east -68.72247 45.17193
## 171 NA 5.186250 23021 east -69.33369 45.56518
## 172 NA 11.439402 24003 east -76.57689 39.05786
## 173 NA 11.023192 24005 east -76.61263 39.37221
## 174 NA 11.665486 24015 east -75.95146 39.59263
## 175 NA 9.529465 24025 east -76.29024 39.53467
## 176 NA 10.669508 24031 east -77.12439 39.07261
## 177 NA 10.926625 24033 east -76.88218 38.90633
## 178 NA 10.919743 24043 east -77.77444 39.61470
## 179 NA 11.370975 24510 east -76.61702 39.30796
## 180 NA 9.194519 25003 east -73.20989 42.39613
## 181 NA 8.282754 25005 east -71.10213 41.77833
## 182 NA 7.949325 25009 east -70.97083 42.63547
## 183 NA 9.244576 25013 east -72.57131 42.12756
## 184 NA 8.031137 25017 east -71.27557 42.45909
## 185 NA 8.550007 25023 east -70.81896 41.97888
## 186 NA 9.788672 25025 east -71.07349 42.33495
## 187 NA 8.834979 25027 east -71.84021 42.32964
## 188 NA 8.935905 26005 east -85.90387 42.59273
## 189 NA 8.404685 26017 east -83.91741 43.63666
## 190 NA 9.208338 26021 east -86.42682 41.97124
## 191 NA 14.623500 26033 east -84.49799 46.32227
## 192 NA 9.140998 26049 east -83.70050 43.01188
## 193 NA 9.253624 26065 east -84.46983 42.66775
## 194 NA 10.071916 26077 east -85.56132 42.25964
## 195 NA 10.229916 26081 east -85.61270 42.98864
## 196 NA 9.535325 26091 east -84.06593 41.91763
## 197 NA 9.669258 26099 east -82.95176 42.59114
## 198 NA 6.745091 26101 east -86.13540 44.32379
## 199 NA 6.153301 26113 east -85.14557 44.32747
## 200 NA 10.370516 26115 east -83.48413 41.91905
## 201 NA 9.090500 26121 east -86.21617 43.25468
## 202 NA 10.020057 26125 east -83.30560 42.58823
## 203 NA 9.993553 26139 east -86.04975 42.92997
## 204 NA 9.888886 26147 east -82.59785 42.91695
## 205 NA 10.104034 26161 east -83.75947 42.26033
## 206 NA 10.880669 26163 east -83.19795 42.33724
## 207 NA 9.568676 27037 east -93.12013 44.75880
## 208 NA 9.649893 27053 east -93.36236 44.97133
## 209 NA 9.594459 27109 east -92.44255 44.00748
## 210 NA 10.205015 27123 east -93.10659 44.99092
## 211 NA 6.263498 27137 east -92.40455 47.36915
## 212 NA 9.422119 27139 east -93.51083 44.69634
## 213 NA 8.559982 27145 east -94.49589 45.54086
## 214 NA 10.287751 28001 east -91.35435 31.53554
## 215 NA 10.758182 28011 east -90.83014 33.80328
## 216 NA 10.645153 28033 east -89.99207 34.90434
## 217 NA 11.998575 28035 east -89.27621 31.26302
## 218 NA 9.768253 28043 east -89.80345 33.77504
## 219 NA 8.993276 28047 east -89.07954 30.42790
## 220 NA 11.447412 28049 east -90.30537 32.29433
## 221 NA 9.705737 28059 east -88.63314 30.45571
## 222 NA 12.157404 28067 east -89.17378 31.65227
## 223 NA 11.497955 28075 east -88.67819 32.40008
## 224 NA 11.304253 28081 east -88.68251 34.27690
## 225 NA 11.375990 29021 east -94.82701 39.72254
## 226 NA 9.823971 29037 east -94.37666 38.68937
## 227 NA 10.199264 29039 east -93.88540 37.73050
## 228 NA 9.609904 29047 east -94.47403 39.25367
## 229 NA 10.135240 29077 east -93.31609 37.21630
## 230 NA 11.173592 29095 east -94.46324 39.02401
## 231 NA 10.935242 29099 east -90.49927 38.28584
## 232 NA 12.523549 29510 east -90.24281 38.62772
## 233 NA 9.376956 31055 east -96.05094 41.26403
## 234 NA 7.367368 31079 east -98.41773 40.89416
## 235 NA 8.503896 31109 east -96.68139 40.79670
## 236 NA 10.006475 31153 east -96.03110 41.14287
## 237 NA 8.759718 31177 east -96.19888 41.51628
## 238 NA 5.711976 33001 east -71.43608 43.51520
## 239 NA 10.208299 33005 east -72.24208 42.91180
## 240 NA 6.834998 33009 east -71.89463 43.89607
## 241 NA 8.349729 33011 east -71.58274 42.89558
## 242 NA 8.515929 33013 east -71.63725 43.27770
## 243 NA 7.598258 33015 east -71.07762 42.96287
## 244 NA 8.827018 34001 east -74.62157 39.44772
## 245 NA 9.035905 34003 east -74.05904 40.93046
## 246 NA 9.876226 34007 east -75.02926 39.86921
## 247 NA 9.117745 34015 east -75.13620 39.75278
## 248 NA 10.656258 34017 east -74.06900 40.74035
## 249 NA 9.024159 34021 east -74.71165 40.25291
## 250 NA 8.632513 34023 east -74.38392 40.49404
## 251 NA 8.073528 34027 east -74.50044 40.86717
## 252 NA 8.677116 34029 east -74.21090 39.91882
## 253 NA 9.538970 34031 east -74.21683 40.94649
## 254 NA 10.432357 34039 east -74.30165 40.65775
## 255 NA 9.551347 34041 east -75.03070 40.80302
## 256 NA 8.233488 36001 east -73.84908 42.66210
## 257 NA 11.367129 36005 east -73.87321 40.84676
## 258 NA 7.771279 36013 east -79.34572 42.24070
## 259 NA 10.034822 36029 east -78.79630 42.85991
## 260 NA 4.625279 36031 east -73.73759 44.16219
## 261 NA 10.884358 36047 east -73.95225 40.64981
## 262 NA 8.251590 36055 east -77.63032 43.15996
## 263 NA 9.704766 36059 east -73.60254 40.71595
## 264 NA 11.784900 36061 east -73.97353 40.76966
## 265 NA 8.904244 36063 east -78.84810 43.14609
## 266 NA 7.926522 36067 east -76.17912 43.04843
## 267 NA 8.591969 36071 east -74.27948 41.40612
## 268 NA 10.114641 36081 east -73.82038 40.71429
## 269 NA 9.936354 36085 east -74.14092 40.58220
## 270 NA 7.500272 36101 east -77.36315 42.28362
## 271 NA 8.924974 36103 east -73.02623 40.83331
## 272 NA 9.611338 36119 east -73.79401 41.05717
## 273 NA 11.272163 37001 east -79.41392 36.07604
## 274 NA 9.319477 37021 east -82.53146 35.59526
## 275 NA 10.383642 37033 east -79.32001 36.40641
## 276 NA 11.806276 37035 east -81.24757 35.69561
## 277 NA 9.707310 37037 east -79.28286 35.70231
## 278 NA 11.376835 37051 east -78.90725 35.05621
## 279 NA 12.433921 37057 east -80.19781 35.81681
## 280 NA 9.750804 37061 east -77.98801 34.92412
## 281 NA 10.954336 37063 east -78.90036 35.99834
## 282 NA 10.062751 37065 east -77.64031 35.91624
## 283 NA 10.971491 37067 east -80.24460 36.11420
## 284 NA 11.308123 37071 east -81.16645 35.28455
## 285 NA 11.022960 37081 east -79.84698 36.05264
## 286 NA 11.045066 37087 east -82.96053 35.52338
## 287 NA 9.445623 37099 east -83.18386 35.31347
## 288 NA 9.809391 37107 east -77.63059 35.25056
## 289 NA 11.033742 37111 east -82.03840 35.68694
## 290 NA 9.444095 37117 east -77.12383 35.83833
## 291 NA 11.315136 37119 east -80.83086 35.23179
## 292 NA 10.194477 37121 east -82.14284 35.97769
## 293 NA 10.506809 37123 east -79.90331 35.34010
## 294 NA 9.971608 37147 east -77.39108 35.58991
## 295 NA 10.710838 37155 east -79.10653 34.65837
## 296 NA 11.681061 37159 east -80.51550 35.61734
## 297 NA 10.150782 37173 east -83.44653 35.42182
## 298 NA 10.747831 37183 east -78.66575 35.79651
## 299 NA 9.434592 37189 east -81.70011 36.21793
## 300 NA 11.059944 37191 east -77.99534 35.36471
## 301 NA 8.598159 38017 east -97.08444 46.89958
## 302 NA 9.719606 39009 east -82.07508 39.35691
## 303 NA 13.424673 39017 east -84.51319 39.42786
## 304 NA 12.776898 39023 east -83.82636 39.92131
## 305 NA 11.575359 39025 east -84.17742 39.05957
## 306 NA 12.593666 39035 east -81.66385 41.46030
## 307 NA 11.909242 39049 east -82.99653 39.99051
## 308 NA 11.875767 39057 east -83.94940 39.71256
## 309 NA 13.770987 39061 east -84.49379 39.17953
## 310 NA 12.807829 39081 east -80.72585 40.36070
## 311 NA 12.228452 39087 east -82.54751 38.53040
## 312 NA 10.415077 39093 east -82.13582 41.37027
## 313 NA 11.556389 39095 east -83.60847 41.64927
## 314 NA 12.434767 39099 east -80.70650 41.05256
## 315 NA 11.149050 39103 east -81.87205 41.11258
## 316 NA 13.217910 39113 east -84.21859 39.75026
## 317 NA 11.494418 39133 east -81.25372 41.17219
## 318 NA 11.685097 39135 east -84.63719 39.74386
## 319 NA 11.594706 39145 east -82.94307 38.78893
## 320 NA 13.224641 39151 east -81.38360 40.81941
## 321 NA 13.082987 39153 east -81.52064 41.09560
## 322 NA 11.864689 39165 east -84.22176 39.44015
## 323 NA 12.232856 40001 east -94.64069 35.89642
## 324 NA 8.950644 40015 east -98.34474 35.13654
## 325 NA 10.716408 40097 east -95.21124 36.31482
## 326 NA 10.549134 40101 east -95.36713 35.68659
## 327 NA 9.630629 40109 east -97.49653 35.50967
## 328 NA 9.467113 40115 east -94.83983 36.86137
## 329 NA 10.217874 40121 east -95.71430 34.93971
## 330 NA 10.245121 40135 east -94.78890 35.48435
## 331 NA 11.147216 40143 east -95.93900 36.12476
## 332 NA 11.355195 42001 east -77.18254 39.85438
## 333 NA 13.136062 42003 east -79.95970 40.44057
## 334 NA 12.897254 42007 east -80.31267 40.69541
## 335 NA 11.132102 42011 east -75.91424 40.37774
## 336 NA 11.342894 42017 east -75.02784 40.25494
## 337 NA 13.195630 42021 east -78.78923 40.43547
## 338 NA 9.923461 42027 east -77.83446 40.88728
## 339 NA 13.654429 42029 east -75.68330 39.98881
## 340 NA 13.919460 42041 east -77.13368 40.20120
## 341 NA 12.420925 42043 east -76.80466 40.32449
## 342 NA 13.308193 42045 east -75.35064 39.91298
## 343 NA 10.511073 42049 east -80.05698 42.05203
## 344 NA 9.677219 42069 east -75.63221 41.43348
## 345 NA 12.702044 42071 east -76.28122 40.06767
## 346 NA 11.020534 42085 east -80.31313 41.29443
## 347 NA 10.774001 42091 east -75.31743 40.17060
## 348 NA 12.709182 42095 east -75.32207 40.70286
## 349 NA 11.500995 42101 east -75.14479 39.99801
## 350 NA 11.838172 42125 east -80.13855 40.18261
## 351 NA 13.281318 42129 east -79.56187 40.32611
## 352 NA 12.459742 42133 east -76.74218 39.92730
## 353 NA 6.615918 44003 east -71.47953 41.70412
## 354 NA 8.690118 44007 east -71.45616 41.85285
## 355 NA 9.201633 45019 east -79.99138 32.80964
## 356 NA 10.297114 45025 east -80.15629 34.65417
## 357 NA 10.740024 45037 east -81.90446 33.77124
## 358 NA 10.563050 45041 east -79.74030 34.06507
## 359 NA 11.558594 45045 east -82.35311 34.86127
## 360 NA 11.422803 45063 east -81.23232 33.94922
## 361 NA 8.600291 45073 east -83.02142 34.71504
## 362 NA 11.203930 45079 east -80.97433 34.03540
## 363 NA 8.319356 46011 east -96.79849 44.34542
## 364 NA 8.136039 46013 east -98.38199 45.53867
## 365 NA 8.631982 46029 east -97.15738 44.93933
## 366 NA 9.101323 46099 east -96.74220 43.60144
## 367 NA 11.782392 47009 east -83.96924 35.75428
## 368 NA 10.815339 47037 east -86.76576 36.15660
## 369 NA 10.104499 47045 east -89.35898 36.06046
## 370 NA 11.480995 47065 east -85.23010 35.09552
## 371 NA 12.155269 47093 east -83.95593 35.97311
## 372 NA 9.043446 47099 east -87.37581 35.20526
## 373 NA 12.445652 47105 east -84.29551 35.74490
## 374 NA 11.721120 47107 east -84.59591 35.41638
## 375 NA 9.721455 47113 east -88.82738 35.62865
## 376 NA 9.755862 47119 east -87.06484 35.61640
## 377 NA 10.607539 47125 east -87.37043 36.51062
## 378 NA 10.459144 47141 east -85.49528 36.15783
## 379 NA 12.035383 47145 east -84.52673 35.87976
## 380 NA 10.658054 47157 east -89.92690 35.13338
## 381 NA 11.053976 47163 east -82.38360 36.53036
## 382 NA 10.431050 47165 east -86.48751 36.43644
## 383 NA 8.911126 48029 east -98.50894 29.45323
## 384 NA 11.277835 48037 east -94.25020 33.43807
## 385 NA 10.711332 48061 east -97.55744 26.08299
## 386 NA 10.102954 48113 east -96.78883 32.80358
## 387 NA 9.747652 48139 east -96.79844 32.39143
## 388 NA 11.727956 48201 east -95.36306 29.78778
## 389 NA 9.873713 48203 east -94.35308 32.54868
## 390 NA 10.204761 48215 east -98.15980 26.24229
## 391 NA 9.866045 48245 east -94.06177 29.98423
## 392 NA 9.722218 48355 east -97.45894 27.74715
## 393 NA 9.835390 48361 east -93.84206 30.11460
## 394 NA 10.304592 48439 east -97.26583 32.76383
## 395 NA 9.214648 48453 east -97.76981 30.32181
## 396 NA 6.902222 50003 east -73.10796 43.02942
## 397 NA 6.296837 50007 east -73.12727 44.47568
## 398 NA 9.782239 50021 east -73.05306 43.59508
## 399 NA 9.386446 51003 east -78.53528 38.02418
## 400 NA 10.881701 51013 east -77.10826 38.88034
## 401 NA 9.487579 51036 east -77.07161 37.34487
## 402 NA 10.236205 51041 east -77.54590 37.41994
## 403 NA 10.353459 51059 east -77.24957 38.84085
## 404 NA 11.075066 51069 east -78.24462 39.19275
## 405 NA 9.767818 51087 east -77.47699 37.58475
## 406 NA 10.190762 51107 east -77.56882 39.06463
## 407 NA 9.844875 51139 east -78.48847 38.61346
## 408 NA 10.801909 51165 east -78.84438 38.46651
## 409 NA 10.221530 51520 east -82.17619 36.61112
## 410 NA 9.327619 51680 east -79.17020 37.40367
## 411 NA 10.257793 51710 east -76.25990 36.88575
## 412 NA 10.314598 51770 east -79.95571 37.27690
## 413 NA 9.928187 51775 east -80.05584 37.28690
## 414 NA 9.802806 51810 east -76.08718 36.83450
## 415 NA 12.876411 54003 east -77.98393 39.46299
## 416 NA 13.339803 54009 east -80.59096 40.29642
## 417 NA 13.043723 54011 east -82.32324 38.41292
## 418 NA 12.405182 54029 east -80.57548 40.50472
## 419 NA 11.783240 54033 east -80.34195 39.29704
## 420 NA 12.420487 54039 east -81.60569 38.34893
## 421 NA 12.843694 54049 east -80.19794 39.49985
## 422 NA 13.131861 54051 east -80.71289 39.90745
## 423 NA 11.556837 54061 east -79.99478 39.63479
## 424 NA 12.423321 54069 east -80.67612 40.07877
## 425 NA 10.114421 54081 east -81.22827 37.77134
## 426 NA 12.840365 54107 east -81.53170 39.25943
## 427 NA 5.715239 55003 east -90.72189 46.33752
## 428 NA 11.366288 55009 east -88.03277 44.49349
## 429 NA 11.098744 55025 east -89.39367 43.07150
## 430 NA 6.308198 55041 east -88.77048 45.60437
## 431 NA 11.219191 55043 east -90.66734 42.87649
## 432 NA 10.440923 55059 east -87.95243 42.56867
## 433 NA 10.058572 55063 east -91.18913 43.85941
## 434 NA 9.417947 55071 east -87.75858 44.10232
## 435 NA 11.575934 55079 east -87.96071 43.03411
## 436 NA 10.248115 55087 east -88.41002 44.33653
## 437 NA 10.030949 55089 east -87.94472 43.32925
## 438 NA 9.749491 55109 east -92.50697 45.02847
## 439 NA 9.208060 55111 east -89.90087 43.44451
## 440 NA 8.483095 55119 east -90.49127 45.19183
## 441 NA 6.111594 55125 east -89.53184 46.02458
## 442 NA 11.975980 55133 east -88.25589 43.03635
## 443 NA 6.058860 2020 west -149.76210 61.19190
## 444 NA 11.101467 2090 west -147.56838 64.81859
## 445 NA 7.308113 2110 west -134.51158 58.35142
## 446 NA 7.147626 2170 west -149.48109 61.76274
## 447 NA 6.929844 4003 west -109.90432 31.75027
## 448 NA 6.132351 4005 west -111.51106 35.77144
## 449 NA 8.228339 4013 west -112.08791 33.49451
## 450 NA 5.328475 4019 west -111.08862 32.17841
## 451 NA 10.502862 4021 west -111.49811 32.96567
## 452 NA 11.326499 4023 west -110.90573 31.48190
## 453 NA 5.165413 4025 west -112.41471 34.65033
## 454 NA 9.465902 6001 west -122.09615 37.71666
## 455 NA 9.842297 6007 west -121.64631 39.64835
## 456 NA 6.768743 6009 west -120.55724 38.16757
## 457 NA 6.944759 6011 west -122.15239 39.18749
## 458 NA 8.545173 6013 west -122.06089 37.94242
## 459 NA 16.194524 6019 west -119.90347 36.63837
## 460 NA 7.196966 6023 west -123.98826 40.74740
## 461 NA 9.417947 6025 west -115.48759 32.96310
## 462 NA 6.431331 6027 west -117.69333 36.70904
## 463 NA 15.803780 6029 west -118.68334 35.29602
## 464 NA 18.440731 6031 west -119.81127 36.15514
## 465 NA 3.494351 6033 west -122.75179 39.02326
## 466 NA 16.661801 6037 west -118.23422 34.08851
## 467 NA 8.151318 6045 west -123.43022 39.36310
## 468 NA 15.015733 6047 west -120.67412 37.24578
## 469 NA 6.455655 6053 west -121.52895 36.45073
## 470 NA 5.280999 6057 west -120.81427 39.27078
## 471 NA 14.506272 6059 west -117.86166 33.73152
## 472 NA 7.969183 6061 west -120.94762 38.97158
## 473 NA 9.678393 6063 west -120.85120 40.02015
## 474 NA 17.429049 6065 west -116.80361 33.78331
## 475 NA 9.679342 6067 west -121.38927 38.56111
## 476 NA 5.995774 6069 west -121.28884 36.74575
## 477 NA 9.843627 6073 west -117.06299 32.90462
## 478 NA 10.741087 6075 west -122.43739 37.75988
## 479 NA 11.560022 6077 west -121.28131 37.94596
## 480 NA 7.308466 6079 west -120.58746 35.37391
## 481 NA 8.999636 6081 west -122.33193 37.53104
## 482 NA 8.415964 6083 west -120.07143 34.64938
## 483 NA 9.841446 6085 west -121.91282 37.30651
## 484 NA 6.196792 6087 west -121.97376 37.00235
## 485 NA 5.560546 6089 west -122.12120 40.68870
## 486 NA 5.702414 6093 west -122.50282 41.61494
## 487 NA 10.099460 6095 west -122.05183 38.23593
## 488 NA 6.925403 6097 west -122.76862 38.44009
## 489 NA 16.251904 6099 west -120.95880 37.61380
## 490 NA 7.309657 6101 west -121.65816 39.07619
## 491 NA 16.183583 6107 west -119.16609 36.23465
## 492 NA 9.617136 6111 west -119.04533 34.27705
## 493 NA 7.705003 6113 west -121.77571 38.63415
## 494 NA 8.469519 8001 west -104.86995 39.86672
## 495 NA 6.654664 8005 west -104.84145 39.64190
## 496 NA 6.771765 8013 west -105.19515 40.05809
## 497 NA 7.733877 8031 west -104.96549 39.72629
## 498 NA 5.732005 8035 west -104.88872 39.46117
## 499 NA 4.186090 8039 west -104.29216 39.29364
## 500 NA 5.971900 8041 west -104.74883 38.86544
## 501 NA 6.614878 8069 west -105.21301 40.53410
## 502 NA 9.328776 8077 west -108.50910 39.09523
## 503 NA 6.275940 8083 west -108.57584 37.35369
## 504 NA 7.680745 8123 west -104.72821 40.34668
## 505 NA 6.824419 15001 west -155.39316 19.68389
## 506 NA 5.783519 15003 west -158.03586 21.42613
## 507 NA 5.693945 15009 west -156.61584 20.90103
## 508 NA 10.329405 16001 west -116.28007 43.56030
## 509 NA 4.917140 16005 west -112.34627 42.77779
## 510 NA 10.792448 16009 west -116.63764 47.24004
## 511 NA 8.846641 16027 west -116.67363 43.61814
## 512 NA 10.376327 16041 west -111.85559 42.15321
## 513 NA 12.437898 16059 west -113.77824 44.94882
## 514 NA 13.346910 16079 west -116.04049 47.41788
## 515 NA 8.185202 30029 west -114.32128 48.26788
## 516 NA 6.875453 30031 west -111.18679 45.66387
## 517 NA 8.030033 30049 west -112.15393 46.81387
## 518 NA 11.464039 30053 west -115.35935 48.55862
## 519 NA 8.294273 30063 west -113.96861 46.93997
## 520 NA 7.187243 30081 west -114.11724 46.21017
## 521 NA 5.135916 30083 west -104.45825 47.77934
## 522 NA 6.170566 30089 west -115.08446 47.63311
## 523 NA 8.858588 30093 west -112.56254 45.96093
## 524 NA 6.672426 32003 west -115.11790 36.13922
## 525 NA 7.203337 32031 west -119.76363 39.70417
## 526 NA 5.425800 35001 west -106.61247 35.10105
## 527 NA 5.689097 35005 west -104.44499 33.35589
## 528 NA 7.838288 35013 west -106.78167 32.27412
## 529 NA 4.793644 35017 west -108.24675 32.70857
## 530 NA 5.678738 35025 west -103.32776 32.72844
## 531 NA 4.601408 35045 west -108.30406 36.57760
## 532 NA 4.195688 35049 west -106.00713 35.58883
## 533 NA 4.460193 38007 west -103.35967 46.99011
## 534 NA 7.057676 38015 west -100.64104 46.89068
## 535 NA 6.121487 38057 west -101.68637 47.28470
## 536 NA 4.978397 41017 west -121.33720 43.99606
## 537 NA 11.422585 41025 west -119.03972 43.46143
## 538 NA 6.544022 41029 west -122.78566 42.38983
## 539 NA 7.857464 41033 west -123.43236 42.36216
## 540 NA 11.201131 41035 west -121.65163 42.66203
## 541 NA 9.606104 41037 west -120.57698 42.71338
## 542 NA 7.981166 41039 west -123.08081 43.98925
## 543 NA 9.559339 41043 west -122.73143 44.54483
## 544 NA 7.138978 41051 west -122.61323 45.52247
## 545 NA 7.512537 41059 west -118.81731 45.63730
## 546 NA 6.722891 41061 west -118.03131 45.32501
## 547 NA 8.617499 41067 west -122.91514 45.49680
## 548 NA 4.324736 46033 west -103.51523 43.69311
## 549 NA 4.175901 46071 west -101.66929 43.76851
## 550 NA 6.105285 46103 west -103.10550 44.04519
## 551 NA 5.349604 48043 west -103.40222 29.99413
## 552 NA 7.412242 48135 west -102.43482 31.86721
## 553 NA 10.305085 48141 west -106.35181 31.77016
## 554 NA 7.740238 48303 west -101.85207 33.57008
## 555 NA 5.998857 48375 west -101.84221 35.25705
## 556 NA 7.974543 49003 west -112.48032 41.59579
## 557 NA 9.757741 49005 west -111.84212 41.75108
## 558 NA 9.631552 49011 west -111.94075 41.00435
## 559 NA 9.632232 49035 west -111.91557 40.65580
## 560 NA 6.624555 49045 west -112.69003 40.41542
## 561 NA 8.866858 49049 west -111.73671 40.23660
## 562 NA 9.195895 49057 west -111.97287 41.23218
## 563 NA 7.930779 53011 west -122.54783 45.69219
## 564 NA 7.939444 53033 west -122.19685 47.54874
## 565 NA 8.843770 53053 west -122.38357 47.17157
## 566 NA 7.821492 53061 west -122.09418 47.97684
## 567 NA 9.462253 53077 west -120.51294 46.48245
## 568 NA 5.272252 56005 west -105.49616 44.08580
## 569 NA 3.382626 56009 west -105.50272 42.85738
## 570 NA 8.505596 56013 west -108.68112 43.09821
## 571 NA 4.132739 56021 west -104.76208 41.19620
## 572 NA 4.955570 56029 west -108.99902 44.56896
## 573 NA 6.549239 56033 west -106.96813 44.78363
## 574 NA 5.632587 56035 west -109.99115 42.75095
## 575 NA 6.349710 56037 west -109.16820 41.64071
## 576 NA 4.565808 56039 west -110.67353 43.61212
boxplot(pollution$pm25, col = "blue")
abline(h=12)
hist(pollution$pm25, col="green")
abline(v = 12, lwd =2)
abline(v = median(pollution$pm25), col = "magenta", lwd = 4)
abline(v = mean(pollution$pm25), col = "blue", lwd = 4)
hist(pollution$pm25, col="green", breaks = 100)
rug(pollution$pm25)
barplot(table(pollution$region), col = "wheat", main = "Number of Counties in Each Region")
boxplot(pm25 ~ region, data=pollution, col = "red")
par(mfrow = c(2,1), mar = c(4,4,2,1))
hist(subset(pollution, region=="east")$pm25, col="green")
hist(subset(pollution, region=="west")$pm25, col="green")
with(pollution, plot(latitude, pm25))
abline(h=12, lwd=4,lty=5)
par(mfrow = c(1,2), mar = c(5,4,2,1))
with(subset(pollution, region =="west"),
plot(latitude, pm25, main = "Oeste"))
with(subset(pollution, region =="east"),
plot(latitude, pm25, main = "Este"))
library(datasets)
data("cars")
with(cars, plot(speed, dist))
library(lattice)
state <- data.frame(state.x77, region=state.region)
xyplot(Life.Exp~Income | region, data = state, layout=c(4,1))
library(ggplot2)
data(mpg)
qplot(displ, hwy, data=mpg)
library(datasets)
hist(airquality$Ozone)
library(datasets)
with(airquality, plot(Wind, Ozone))
library(datasets)
airquality <- transform(airquality, Month=factor(Month))
boxplot(Ozone~Month, airquality, xlab="Month", ylab="Ozone (ppb)")
par("lty")
## [1] "solid"
par("col")
## [1] "black"
par("pch")
## [1] 1
par("bg")
## [1] "white"
par("mar")
## [1] 5.1 4.1 4.1 2.1
par("mfrow")
## [1] 1 1
library(datasets)
with(airquality, plot(Wind, Ozone))
title(main = "Ozone and Wind in New York City")
library(datasets)
with(airquality, plot(Wind, Ozone, main="Ozone and Wind in New York City"))
with(subset(airquality, Month == 5), points(Wind, Ozone, col = "blue"))
library(datasets)
with(airquality, plot(Wind, Ozone, main="Ozono y viento en la ciudad de Nueva York", type="n"))
with(subset(airquality, Month == 5), points(Wind, Ozone, col = "blue"))
with(subset(airquality, Month != 5), points(Wind, Ozone, col = "red"))
legend("topright", pch=1, col=c("blue","red"), legend=c("May", "Other Months"))
library(datasets)
with(airquality, plot(Wind, Ozone, main="Ozono y viento en la ciudad de Nueva York", pch=20))
model <- lm(Ozone ~ Wind, airquality)
abline(model, lwd=2)
summary(model)
##
## Call:
## lm(formula = Ozone ~ Wind, data = airquality)
##
## Residuals:
## Min 1Q Median 3Q Max
## -51.572 -18.854 -4.868 15.234 90.000
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 96.8729 7.2387 13.38 < 2e-16 ***
## Wind -5.5509 0.6904 -8.04 9.27e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 26.47 on 114 degrees of freedom
## (37 observations deleted due to missingness)
## Multiple R-squared: 0.3619, Adjusted R-squared: 0.3563
## F-statistic: 64.64 on 1 and 114 DF, p-value: 9.272e-13
par(mfrow=c(1,2))
with(airquality, {
plot(Wind, Ozone, main="Ozone and Wind")
plot(Solar.R, Ozone, main="Ozone and Solar Radiation")
})
par(mfrow=c(1,3), mar=c(4,4,2,1), oma=c(0,0,2,0))
with(airquality,{
plot(Wind,Ozone,main="Ozone and Wind")
plot(Solar.R, Ozone, main = "Ozone and Solar Radiation")
plot(Temp, Ozone, main = "Ozone and Temperature")
mtext("Ozone and Weather in New York City", outer=TRUE)
})
Histograma aleatorio
x <- rnorm(100)
hist(x)
boxplot(x)
summary(x)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -2.22742 -0.63403 -0.01828 0.01784 0.61539 2.09843
y <- rnorm(100)
plot(x,y)
z <- rnorm(100)
plot(x, z)
example(points)
##
## points> require(stats) # for rnorm
##
## points> plot(-4:4, -4:4, type = "n") # setting up coord. system
##
## points> points(rnorm(200), rnorm(200), col = "red")
##
## points> points(rnorm(100)/2, rnorm(100)/2, col = "blue", cex = 1.5)
##
## points> op <- par(bg = "light blue")
##
## points> x <- seq(0, 2*pi, len = 51)
##
## points> ## something "between type='b' and type='o'":
## points> plot(x, sin(x), type = "o", pch = 21, bg = par("bg"), col = "blue", cex = .6,
## points+ main = 'plot(..., type="o", pch=21, bg=par("bg"))')
##
## points> par(op)
##
## points> ## Not run:
## points> ##D ## The figure was produced by calls like
## points> ##D png("pch.png", height = 0.7, width = 7, res = 100, units = "in")
## points> ##D par(mar = rep(0,4))
## points> ##D plot(c(-1, 26), 0:1, type = "n", axes = FALSE)
## points> ##D text(0:25, 0.6, 0:25, cex = 0.5)
## points> ##D points(0:25, rep(0.3, 26), pch = 0:25, bg = "grey")
## points> ## End(Not run)
## points>
## points> ##-------- Showing all the extra & some char graphics symbols ---------
## points> pchShow <-
## points+ function(extras = c("*",".", "o","O","0","+","-","|","%","#"),
## points+ cex = 3, ## good for both .Device=="postscript" and "x11"
## points+ col = "red3", bg = "gold", coltext = "brown", cextext = 1.2,
## points+ main = paste("plot symbols : points (... pch = *, cex =",
## points+ cex,")"))
## points+ {
## points+ nex <- length(extras)
## points+ np <- 26 + nex
## points+ ipch <- 0:(np-1)
## points+ k <- floor(sqrt(np))
## points+ dd <- c(-1,1)/2
## points+ rx <- dd + range(ix <- ipch %/% k)
## points+ ry <- dd + range(iy <- 3 + (k-1)- ipch %% k)
## points+ pch <- as.list(ipch) # list with integers & strings
## points+ if(nex > 0) pch[26+ 1:nex] <- as.list(extras)
## points+ plot(rx, ry, type = "n", axes = FALSE, xlab = "", ylab = "", main = main)
## points+ abline(v = ix, h = iy, col = "lightgray", lty = "dotted")
## points+ for(i in 1:np) {
## points+ pc <- pch[[i]]
## points+ ## 'col' symbols with a 'bg'-colored interior (where available) :
## points+ points(ix[i], iy[i], pch = pc, col = col, bg = bg, cex = cex)
## points+ if(cextext > 0)
## points+ text(ix[i] - 0.3, iy[i], pc, col = coltext, cex = cextext)
## points+ }
## points+ }
##
## points> pchShow()
##
## points> pchShow(c("o","O","0"), cex = 2.5)
##
## points> pchShow(NULL, cex = 4, cextext = 0, main = NULL)
##
## points> ## No test:
## points> ##D ## ------------ test code for various pch specifications -------------
## points> ##D # Try this in various font families (including Hershey)
## points> ##D # and locales. Use sign = -1 asserts we want Latin-1.
## points> ##D # Standard cases in a MBCS locale will not plot the top half.
## points> ##D TestChars <- function(sign = 1, font = 1, ...)
## points> ##D {
## points> ##D MB <- l10n_info()$MBCS
## points> ##D r <- if(font == 5) { sign <- 1; c(32:126, 160:254)
## points> ##D } else if(MB) 32:126 else 32:255
## points> ##D if (sign == -1) r <- c(32:126, 160:255)
## points> ##D par(pty = "s")
## points> ##D plot(c(-1,16), c(-1,16), type = "n", xlab = "", ylab = "",
## points> ##D xaxs = "i", yaxs = "i",
## points> ##D main = sprintf("sign = %d, font = %d", sign, font))
## points> ##D grid(17, 17, lty = 1) ; mtext(paste("MBCS:", MB))
## points> ##D for(i in r) try(points(i%%16, i%/%16, pch = sign*i, font = font,...))
## points> ##D }
## points> ##D TestChars()
## points> ##D try(TestChars(sign = -1))
## points> ##D TestChars(font = 5) # Euro might be at 160 (0+10*16).
## points> ##D # macOS has apple at 240 (0+15*16).
## points> ##D try(TestChars(-1, font = 2)) # bold
## points> ## End(No test)
## points>
## points>