This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.

Try executing this chunk by clicking the Run button within the chunk or by placing your cursor inside it and pressing Ctrl+Shift+Enter.

Load & install packages

library(curl)

package 㤼㸱curl㤼㸲 was built under R version 3.3.3

library(devtools)
install_github("weatherData", "Ram-N")

Username parameter is deprecated. Please use Ram-N/weatherDataSkipping install of 'weatherData' from a github remote, the SHA1 (76401386) has not changed since last install.
  Use `force = TRUE` to force installation

library(weatherData)
#install.packages("weatherData") - the personal weather station data is only available via github version
#package info at https://ram-n.github.io/weatherData/
library(ggplot2)

Stackoverflow is a great place to get help: http://stackoverflow.com/tags/ggplot2.

library(psych)


Attaching package: 㤼㸱psych㤼㸲

The following objects are masked from 㤼㸱package:ggplot2㤼㸲:

    %+%, alpha

library(reshape2)
library(knitr)
library(lattice) #just to illustrate another histogram function

package 㤼㸱lattice㤼㸲 was built under R version 3.3.3

Add a new chunk by clicking the Insert Chunk button on the toolbar or by pressing Ctrl+Alt+I.

When you save the notebook, an HTML file containing the code and output will be saved alongside it (click the Preview button or press Ctrl+Shift+K to preview the HTML file).

checkDataAvailability("SYD", "2017-02-15", station_type = "airportCode")

Getting data from:
 https://www.wunderground.com/history/airport/SYD/2017/2/15/DailyHistory.html?format=1

[1] 72

The following columns are available for:2017-02-15

 [1] "TimeAEDT"              "TemperatureC"          "Dew_PointC"           
 [4] "Humidity"              "Sea_Level_PressurehPa" "VisibilityKm"         
 [7] "Wind_Direction"        "Wind_SpeedKm_h"        "Gust_SpeedKm_h"       
[10] "Precipitationmm"       "Events"                "Conditions"           
[13] "WindDirDegrees"        "DateUTC"              
[1] 36 14
        V1 V2 V3 V4   V5 V6   V7   V8 V9 V10  V11                V12 V13                 V14
1 12:00 AM 19 17 86 1018 15 West  7.4        Rain Light Rain Showers 260 2017-02-14 13:00:00
2 12:30 AM 19 17 88 1017 10  SSW 14.8  - N/A Rain Light Rain Showers 200 2017-02-14 13:30:00
3  1:00 AM 19 16 83 1017 10  WSW 11.1  - N/A           Mostly Cloudy 240 2017-02-14 14:00:00
4  2:00 AM 19 16 81 1017 NA West 11.1                                260 2017-02-14 15:00:00
5  2:30 AM 19 16 83 1016 10   NW 11.1  - N/A Rain Light Rain Showers 310 2017-02-14 15:30:00
6  3:00 AM 18 15 78 1017  7  NNW 11.1        Rain Light Rain Showers 340 2017-02-14 16:00:00

Checking Data Availability For SYD
Found Records for 2017-02-15
Data is Available

[1] 1

getWeatherForDate("SYD", "2017-02-15", opt_detailed=TRUE)

Getting data from:
 https://www.wunderground.com/history/airport/SYD/2017/2/15/DailyHistory.html?format=1

[1] 72

The following columns are available for:2017-02-15

 [1] "TimeAEDT"              "TemperatureC"          "Dew_PointC"           
 [4] "Humidity"              "Sea_Level_PressurehPa" "VisibilityKm"         
 [7] "Wind_Direction"        "Wind_SpeedKm_h"        "Gust_SpeedKm_h"       
[10] "Precipitationmm"       "Events"                "Conditions"           
[13] "WindDirDegrees"        "DateUTC"              
[1] 36 14
        V1 V2 V3 V4   V5 V6   V7   V8 V9 V10  V11                V12 V13                 V14
1 12:00 AM 19 17 86 1018 15 West  7.4        Rain Light Rain Showers 260 2017-02-14 13:00:00
2 12:30 AM 19 17 88 1017 10  SSW 14.8  - N/A Rain Light Rain Showers 200 2017-02-14 13:30:00
3  1:00 AM 19 16 83 1017 10  WSW 11.1  - N/A           Mostly Cloudy 240 2017-02-14 14:00:00
4  2:00 AM 19 16 81 1017 NA West 11.1                                260 2017-02-14 15:00:00
5  2:30 AM 19 16 83 1016 10   NW 11.1  - N/A Rain Light Rain Showers 310 2017-02-14 15:30:00
6  3:00 AM 18 15 78 1017  7  NNW 11.1        Rain Light Rain Showers 340 2017-02-14 16:00:00

Checking Data Availability For SYD
Found Records for 2017-02-15
Data is Available
Will be fetching these Columns:

[1] "Time"         "TemperatureC"

Begin getting Daily Data for SYD
SYD 1 2017-02-15 : Fetching 71 Rows with 2 Column(s)

getWeatherForDate("SYD", "2017-02-10", opt_detailed=TRUE)

Getting data from:
 https://www.wunderground.com/history/airport/SYD/2017/2/10/DailyHistory.html?format=1

[1] 75

The following columns are available for:2017-02-10

 [1] "TimeAEDT"              "TemperatureC"          "Dew_PointC"           
 [4] "Humidity"              "Sea_Level_PressurehPa" "VisibilityKm"         
 [7] "Wind_Direction"        "Wind_SpeedKm_h"        "Gust_SpeedKm_h"       
[10] "Precipitationmm"       "Events"                "Conditions"           
[13] "WindDirDegrees"        "DateUTC"              
[1] 37 14
        V1 V2 V3 V4   V5    V6  V7   V8 V9 V10 V11   V12 V13                 V14
1 12:00 AM 25 21 72 1011    25 NNE 13.0         NA  Haze  30 2017-02-09 13:00:00
2 12:30 AM 25 21 78 1010 -9999  NE  9.3  - N/A  NA Clear  40 2017-02-09 13:30:00
3  1:00 AM 25 21 78 1010 -9999  NE 11.1  - N/A  NA Clear  50 2017-02-09 14:00:00
4  2:00 AM 24 22 83 1010    NA  NE 16.7         NA        40 2017-02-09 15:00:00
5  2:30 AM 24 22 89 1010 -9999  NE 13.0  - N/A  NA Clear  50 2017-02-09 15:30:00
6  3:00 AM 24 21 83 1010 -9999  NE  9.3  - N/A  NA Clear  40 2017-02-09 16:00:00

Checking Data Availability For SYD
Found Records for 2017-02-10
Data is Available
Will be fetching these Columns:

[1] "Time"         "TemperatureC"

Begin getting Daily Data for SYD
SYD 1 2017-02-10 : Fetching 74 Rows with 2 Column(s)

getSummarizedWeather("SYD", "2017-02-10")

Grab some data

SYD_weather <- getWeatherForDate("SYD", "2017-02-15", end_date = "2017-04-04", opt_detailed=TRUE, opt_all_columns = T)

Getting data from:
 https://www.wunderground.com/history/airport/SYD/2017/2/15/DailyHistory.html?format=1

[1] 72

The following columns are available for:2017-02-15

 [1] "TimeAEDT"              "TemperatureC"          "Dew_PointC"           
 [4] "Humidity"              "Sea_Level_PressurehPa" "VisibilityKm"         
 [7] "Wind_Direction"        "Wind_SpeedKm_h"        "Gust_SpeedKm_h"       
[10] "Precipitationmm"       "Events"                "Conditions"           
[13] "WindDirDegrees"        "DateUTC"              
[1] 36 14
        V1 V2 V3 V4   V5 V6   V7   V8 V9 V10  V11                V12 V13                 V14
1 12:00 AM 19 17 86 1018 15 West  7.4        Rain Light Rain Showers 260 2017-02-14 13:00:00
2 12:30 AM 19 17 88 1017 10  SSW 14.8  - N/A Rain Light Rain Showers 200 2017-02-14 13:30:00
3  1:00 AM 19 16 83 1017 10  WSW 11.1  - N/A           Mostly Cloudy 240 2017-02-14 14:00:00
4  2:00 AM 19 16 81 1017 NA West 11.1                                260 2017-02-14 15:00:00
5  2:30 AM 19 16 83 1016 10   NW 11.1  - N/A Rain Light Rain Showers 310 2017-02-14 15:30:00
6  3:00 AM 18 15 78 1017  7  NNW 11.1        Rain Light Rain Showers 340 2017-02-14 16:00:00

Getting data from:
 https://www.wunderground.com/history/airport/SYD/2017/4/4/DailyHistory.html?format=1

[1] 77

The following columns are available for:2017-04-04

 [1] "TimeAEST"              "TemperatureC"          "Dew_PointC"           
 [4] "Humidity"              "Sea_Level_PressurehPa" "VisibilityKm"         
 [7] "Wind_Direction"        "Wind_SpeedKm_h"        "Gust_SpeedKm_h"       
[10] "Precipitationmm"       "Events"                "Conditions"           
[13] "WindDirDegrees"        "DateUTC"              
[1] 38 14
        V1 V2 V3 V4   V5 V6  V7   V8   V9 V10  V11                V12 V13                 V14
1 12:00 AM 17 13 67 1028 15  SE 35.2          Rain Light Rain Showers 130 2017-04-03 14:00:00
2 12:30 AM 18 14 77 1027 10 ESE 29.6    - N/A        Scattered Clouds 120 2017-04-03 14:30:00
3  1:00 AM 18 14 77 1027 10  SE 25.9    - N/A        Scattered Clouds 130 2017-04-03 15:00:00
4  2:00 AM 17 14 78 1027 NA ESE 27.8                                  110 2017-04-03 16:00:00
5  2:30 AM 17 13 77 1026  8 ESE 35.2 53.7 N/A Rain Light Rain Showers 120 2017-04-03 16:30:00
6  3:00 AM 18 15 83 1026 10 SSE 25.9    - N/A           Mostly Cloudy 150 2017-04-03 17:00:00

Checking Data Availability For SYD
Found 36 records for 2017-02-15
Found 38 records for 2017-04-04

Data is Available for the interval.

Will be fetching these Columns:

 [1] "Time"                  "TimeAEDT"              "TemperatureC"         
 [4] "Dew_PointC"            "Humidity"              "Sea_Level_PressurehPa"
 [7] "VisibilityKm"          "Wind_Direction"        "Wind_SpeedKm_h"       
[10] "Gust_SpeedKm_h"        "Precipitationmm"       "Events"               
[13] "Conditions"            "WindDirDegrees"        "DateUTC"

Begin getting Daily Data for SYD
SYD 1 2017-02-15 : Fetching 71 Rows with 15 Column(s)
SYD 2 2017-02-16 : Fetching 69 Rows with 15 Column(s)
SYD 3 2017-02-17 : Fetching 77 Rows with 15 Column(s)
SYD 4 2017-02-18 : Fetching 71 Rows with 15 Column(s)
SYD 5 2017-02-19 : Fetching 74 Rows with 15 Column(s)
SYD 6 2017-02-20 : Fetching 70 Rows with 15 Column(s)
SYD 7 2017-02-21 : Fetching 66 Rows with 15 Column(s)
SYD 8 2017-02-22 : Fetching 67 Rows with 15 Column(s)
SYD 9 2017-02-23 : Fetching 67 Rows with 15 Column(s)
SYD 10 2017-02-24 : Fetching 69 Rows with 15 Column(s)
SYD 11 2017-02-25 : Fetching 90 Rows with 15 Column(s)
SYD 12 2017-02-26 : Fetching 82 Rows with 15 Column(s)
SYD 13 2017-02-27 : Fetching 79 Rows with 15 Column(s)
SYD 14 2017-02-28 : Fetching 69 Rows with 15 Column(s)
SYD 15 2017-03-01 : Fetching 78 Rows with 15 Column(s)
SYD 16 2017-03-02 : Fetching 72 Rows with 15 Column(s)
SYD 17 2017-03-03 : Fetching 78 Rows with 15 Column(s)
SYD 18 2017-03-04 : Fetching 76 Rows with 15 Column(s)
SYD 19 2017-03-05 : Fetching 77 Rows with 15 Column(s)
SYD 20 2017-03-06 : Fetching 69 Rows with 15 Column(s)
SYD 21 2017-03-07 : Fetching 75 Rows with 15 Column(s)
SYD 22 2017-03-08 : Fetching 69 Rows with 15 Column(s)
SYD 23 2017-03-09 : Fetching 66 Rows with 15 Column(s)
SYD 24 2017-03-10 : Fetching 70 Rows with 15 Column(s)
SYD 25 2017-03-11 : Fetching 68 Rows with 15 Column(s)
SYD 26 2017-03-12 : Fetching 72 Rows with 15 Column(s)
SYD 27 2017-03-13 : Fetching 69 Rows with 15 Column(s)
SYD 28 2017-03-14 : Fetching 72 Rows with 15 Column(s)
SYD 29 2017-03-15 : Fetching 72 Rows with 15 Column(s)
SYD 30 2017-03-16 : Fetching 78 Rows with 15 Column(s)
SYD 31 2017-03-17 : Fetching 88 Rows with 15 Column(s)
SYD 32 2017-03-18 : Fetching 80 Rows with 15 Column(s)
SYD 33 2017-03-19 : Fetching 81 Rows with 15 Column(s)
SYD 34 2017-03-20 : Fetching 67 Rows with 15 Column(s)
SYD 35 2017-03-21 : Fetching 72 Rows with 15 Column(s)
SYD 36 2017-03-22 : Fetching 76 Rows with 15 Column(s)
SYD 37 2017-03-23 : Fetching 72 Rows with 15 Column(s)
SYD 38 2017-03-24 : Fetching 67 Rows with 15 Column(s)
SYD 39 2017-03-25 : Fetching 67 Rows with 15 Column(s)
SYD 40 2017-03-26 : Fetching 64 Rows with 15 Column(s)
SYD 41 2017-03-27 : Fetching 71 Rows with 15 Column(s)
SYD 42 2017-03-28 : Fetching 70 Rows with 15 Column(s)
SYD 43 2017-03-29 : Fetching 67 Rows with 15 Column(s)
SYD 44 2017-03-30 : Fetching 81 Rows with 15 Column(s)
SYD 45 2017-03-31 : Fetching 70 Rows with 15 Column(s)
SYD 46 2017-04-01 : Fetching 68 Rows with 15 Column(s)
SYD 47 2017-04-02 : Fetching 76 Rows with 15 Column(s)
SYD 48 2017-04-03 : Fetching 73 Rows with 15 Column(s)
SYD 49 2017-04-04 : Fetching 76 Rows with 15 Column(s)

MEL_weather <- getWeatherForDate("MEL", "2017-02-15", end_date = "2017-04-04", opt_detailed=TRUE, opt_all_columns = T)

Getting data from:
 https://www.wunderground.com/history/airport/MEL/2017/2/15/DailyHistory.html?format=1

[1] 5

The following columns are available for:2017-02-15

 [1] "Time_0545"             "TemperatureC"          "Dew_PointC"           
 [4] "Humidity"              "Sea_Level_PressurehPa" "VisibilityKm"         
 [7] "Wind_Direction"        "Wind_SpeedKm_h"        "Gust_SpeedKm_h"       
[10] "Precipitationmm"       "Events"                "Conditions"           
[13] "WindDirDegrees"        "DateUTC"              
[1]  2 14
       V1 V2 V3 V4 V5 V6    V7   V8 V9 V10 V11 V12 V13                 V14
1 8:45 AM  3 -4 50 NA  8  Calm Calm NA  NA  NA  NA  NA 2017-02-15 03:00:00
2 2:45 PM 21 -7  7 NA 12 South 13.0 NA  NA  NA  NA 190 2017-02-15 09:00:00

Getting data from:
 https://www.wunderground.com/history/airport/MEL/2017/4/4/DailyHistory.html?format=1

[1] 6

The following columns are available for:2017-04-04

 [1] "Time_0545"             "TemperatureC"          "Dew_PointC"           
 [4] "Humidity"              "Sea_Level_PressurehPa" "VisibilityKm"         
 [7] "Wind_Direction"        "Wind_SpeedKm_h"        "Gust_SpeedKm_h"       
[10] "Precipitationmm"       "Events"                "Conditions"           
[13] "WindDirDegrees"        "DateUTC"              
[1]  3 14
        V1 V2 V3 V4 V5 V6    V7  V8 V9 V10 V11              V12 V13                 V14
1  5:45 AM  6  4 81 NA  8 South 3.7 NA  NA  NA                  190 2017-04-04 00:00:00
2 11:45 AM 22  5 23 NA 10    SW 1.9 NA  NA  NA Scattered Clouds 220 2017-04-04 06:00:00
3  5:45 PM 19  2 21 NA  8   WSW 3.7 NA  NA  NA    Partly Cloudy 250 2017-04-04 12:00:00

Checking Data Availability For MEL
Found 2 records for 2017-02-15
Found 3 records for 2017-04-04

Data is Available for the interval.

Will be fetching these Columns:

 [1] "Time"                  "Time_0545"             "TemperatureC"         
 [4] "Dew_PointC"            "Humidity"              "Sea_Level_PressurehPa"
 [7] "VisibilityKm"          "Wind_Direction"        "Wind_SpeedKm_h"       
[10] "Gust_SpeedKm_h"        "Precipitationmm"       "Events"               
[13] "Conditions"            "WindDirDegrees"        "DateUTC"

Begin getting Daily Data for MEL
MEL 1 2017-02-15 : Fetching 4 Rows with 15 Column(s)
MEL 2 2017-02-16 : Fetching 5 Rows with 15 Column(s)
MEL 3 2017-02-17 : Fetching 5 Rows with 15 Column(s)
MEL 4 2017-02-18 : Fetching 5 Rows with 15 Column(s)
MEL 5 2017-02-19 : Fetching 5 Rows with 15 Column(s)
MEL 6 2017-02-20 : Fetching 5 Rows with 15 Column(s)
MEL 7 2017-02-21 : Fetching 5 Rows with 15 Column(s)
MEL 8 2017-02-22 : Fetching 5 Rows with 15 Column(s)
MEL 9 2017-02-23 : Fetching 4 Rows with 15 Column(s)
MEL 10 2017-02-24 : Fetching 5 Rows with 15 Column(s)
MEL 11 2017-02-25 : Fetching 4 Rows with 15 Column(s)
MEL 12 2017-02-26 : Fetching 5 Rows with 15 Column(s)
MEL 13 2017-02-27 : Fetching 5 Rows with 15 Column(s)
MEL 14 2017-02-28 : Fetching 5 Rows with 15 Column(s)
MEL 15 2017-03-01 : Fetching 5 Rows with 15 Column(s)
MEL 16 2017-03-02 : Fetching 4 Rows with 15 Column(s)
MEL 17 2017-03-03 : Fetching 4 Rows with 15 Column(s)
MEL 18 2017-03-04 : Fetching 5 Rows with 15 Column(s)
MEL 19 2017-03-05 : Fetching 4 Rows with 15 Column(s)
MEL 20 2017-03-06 : Fetching 5 Rows with 15 Column(s)
MEL 21 2017-03-07 : Fetching 4 Rows with 15 Column(s)
MEL 22 2017-03-08 : Fetching 4 Rows with 15 Column(s)
MEL 23 2017-03-09 : Fetching 5 Rows with 15 Column(s)
MEL 24 2017-03-10 : Fetching 3 Rows with 15 Column(s)
MEL 25 2017-03-11 : Fetching 4 Rows with 15 Column(s)
MEL 26 2017-03-12 : Fetching 4 Rows with 15 Column(s)
MEL 27 2017-03-13 : Fetching 5 Rows with 15 Column(s)
MEL 28 2017-03-14 : Fetching 3 Rows with 15 Column(s)
MEL 29 2017-03-15 : Fetching 4 Rows with 15 Column(s)
MEL 30 2017-03-16 : Fetching 4 Rows with 15 Column(s)
MEL 31 2017-03-17 : Fetching 5 Rows with 15 Column(s)
MEL 32 2017-03-18 : Fetching 5 Rows with 15 Column(s)
MEL 33 2017-03-19 : Fetching 4 Rows with 15 Column(s)
MEL 34 2017-03-20 : Fetching 5 Rows with 15 Column(s)
MEL 35 2017-03-21 : Fetching 4 Rows with 15 Column(s)
MEL 36 2017-03-22 : Fetching 4 Rows with 15 Column(s)
MEL 37 2017-03-23 : Fetching 5 Rows with 15 Column(s)
MEL 38 2017-03-24 : Fetching 5 Rows with 15 Column(s)
MEL 39 2017-03-25 : Fetching 4 Rows with 15 Column(s)
MEL 40 2017-03-26 : Fetching 5 Rows with 15 Column(s)
MEL 41 2017-03-27 : Fetching 5 Rows with 15 Column(s)
MEL 42 2017-03-28 : Fetching 4 Rows with 15 Column(s)
MEL 43 2017-03-29 : Fetching 5 Rows with 15 Column(s)
MEL 44 2017-03-30 : Fetching 4 Rows with 15 Column(s)
MEL 45 2017-03-31 : Fetching 5 Rows with 15 Column(s)
MEL 46 2017-04-01 : Fetching 4 Rows with 15 Column(s)
MEL 47 2017-04-02 : Fetching 5 Rows with 15 Column(s)
MEL 48 2017-04-03 : Fetching 5 Rows with 15 Column(s)
MEL 49 2017-04-04 : Fetching 5 Rows with 15 Column(s)

#getWeatherForDate(station_id = "ISYDNEY143", start_date = "2017-02-10", end_date = "2017-04-04", opt_detailed=TRUE, opt_all_columns = T, station_type = "ID")

How might we explore some of this data?

kable(rbind(describe(SYD_weather$TemperatureC),describe(MEL_weather$TemperatureC)))

	vars	n	mean	sd	median	trimmed	mad	min	max	range	skew	kurtosis	se
X1	1	3515	-356.99943	1912.3546	22	21.89264	2.9652	-9999	32	10031	-4.842284	21.45387	32.25565
X11	1	221	-78.54751	950.2152	13	11.91525	8.8956	-9999	30	10029	-10.297442	104.51675	63.91839

Let’s have a basic look at the distribution

hist(MEL_weather$TemperatureC)

hist(SYD_weather$TemperatureC)

Ah, there’s something wrong! Let’s fix that and get the melbourne and sydney data into the same frame

SYD_temp <- as.data.frame(as.numeric(unlist(subset(SYD_weather, TemperatureC >-300, select=c("TemperatureC"))))) #you could also replace these wiht NA, but here we're just going to exclude the missing data
colnames(SYD_temp)[1] <- "temp"
SYD_temp$loc <- "SYD"
MEL_temp <- as.data.frame(as.numeric(unlist(subset(MEL_weather, TemperatureC >-300, select=c("TemperatureC")))))
colnames(MEL_temp)[1] <- "temp"
MEL_temp$loc <- "MEL"
temps <- rbind(SYD_temp, MEL_temp)
temps$temp <- as.numeric(temps$temp)

Let’s look at it again

ggplot(temps, aes(x = temp, fill = loc)) + geom_histogram(alpha = .5, position = 'identity') #ah, what's the problem here?

Try that again

ggplot(temps, aes(x = temp, fill = loc)) + geom_histogram(alpha = .5, aes(y = ..density..), position = 'identity') #note use of 'density' because we have unequal temperature counts in each dataset. Alpha is the transparency level.

histogram(~ temp | loc, data=temps)

What’s wrong with this?

ggplot(temps) + 
  geom_bar(aes(x = loc, y = temp, fill = loc),
           position = "dodge", stat = "summary", fun.y = "mean")

More informative?

ggplot(temps, aes(x=loc, y=temp, fill=loc)) + geom_boxplot() +
    guides(fill=FALSE)+
    stat_summary(fun.y=mean, geom="point", shape=5, size=4)

Couple of useful things - let’s pull the date out to its own value, and this time we’ll replace missing values (-9999) with NA

ggplot(SYD_weather, aes(x=TemperatureC, y=Dew_PointC)) +
    geom_point(shape=1)      # Use hollow circles

SYD_weather[SYD_weather == -9999] <- NA
SYD_weather$date <- as.Date(SYD_weather$DateUTC)
MEL_weather[MEL_weather == -9999] <- NA
MEL_weather$date <- as.Date(MEL_weather$DateUTC)
ggplot(SYD_weather, aes(x=TemperatureC, y=Dew_PointC)) +
    geom_point(shape=1)      # Use hollow circles

What if we want to explore the relationship between Dew_PointC and other features https://support.office.com/en-us/article/Present-your-data-in-a-scatter-chart-or-a-line-chart-4570a80f-599a-4d6b-a155-104a9018b86e

One way you might be tempted to do this…

bad_example <- subset(SYD_weather, !is.na(Humidity), select=c("Humidity", "Dew_PointC","date"))
bad_example[c("Humidity","Dew_PointC")] <- lapply(bad_example[c("Humidity","Dew_PointC")],as.numeric)

NAs introduced by coercion

bad_example <- aggregate(. ~ date, bad_example, FUN=mean)
#convert to long
bad_example <- melt(bad_example, id.vars = c("date"))
ggplot(data=bad_example, aes(x=date, y=value, group=variable, colour=variable)) +
    geom_line() +
    geom_point()

#Is date an important variable in this analysis? Does the scaling of the data gives us the best available insight into relationships of paired values? Is the use of a line to join datapoints appropriate given missing data?

A better way?

ggplot(SYD_weather, aes(x=Humidity, y=Dew_PointC)) +
    geom_point(shape=1)      # Use hollow circles

cor.test(as.numeric(SYD_weather$Humidity),as.numeric(SYD_weather$Dew_PointC))

NAs introduced by coercion


    Pearson's product-moment correlation

data:  as.numeric(SYD_weather$Humidity) and as.numeric(SYD_weather$Dew_PointC)
t = 45.252, df = 3380, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 0.5927923 0.6347879
sample estimates:
      cor 
0.6142248

Ok what if we want to look at how weather varies over time and place?

SYD_weather$loc <- "Sydney"
MEL_weather$loc <- "Melbourne"

weather <- rbind(SYD_weather[c("TemperatureC","Dew_PointC","Humidity","Wind_SpeedKm_h","Precipitationmm","Events","Conditions","DateUTC","loc")],MEL_weather[c("TemperatureC","Dew_PointC","Humidity","Wind_SpeedKm_h","Precipitationmm","Events","Conditions","DateUTC","loc")])
weather$month <- format(as.Date(weather$DateUTC), "%m")
ggplot(weather, aes(x=month, y=TemperatureC, fill=loc)) + geom_boxplot() +
    guides(fill=FALSE) +
    stat_summary(fun.y=mean, geom="point", shape=5, size=4) +
    facet_wrap(~loc)

Or at how weather events vary by place

unique(weather$Events)

[1] "Rain"              ""                  NA                  "Thunderstorm"     
[5] "Rain-Thunderstorm" "Snow"              "Hail-Thunderstorm"

unique(weather$Conditions)

 [1] "Light Rain Showers"           ""                            
 [3] "Mostly Cloudy"                "Rain Showers"                
 [5] "Partly Cloudy"                "Scattered Clouds"            
 [7] "Clear"                        "Unknown"                     
 [9] "Haze"                         "Light Thunderstorms and Rain"
[11] "Thunderstorms and Rain"       "Light Rain"                  
[13] "Heavy Rain Showers"           "Heavy Thunderstorms and Rain"
[15] "Smoke"                        "Light Drizzle"               
[17] "Drizzle"                      "Thunderstorm"                
[19] "Heavy Rain"                   "Rain"                        
[21] "Overcast"                     "Heavy Drizzle"               
[23] NA                             "Heavy Snow"                  
[25] "Light Snow"                   "Thunderstorms with Hail"

table(weather$Conditions,weather$loc)

                              
                               Melbourne Sydney
                                      49    762
  Clear                                0    394
  Drizzle                              2      6
  Haze                                 1     42
  Heavy Drizzle                        0      1
  Heavy Rain                           0      2
  Heavy Rain Showers                   0     25
  Heavy Snow                           2      0
  Heavy Thunderstorms and Rain         0      5
  Light Drizzle                        4     66
  Light Rain                           1     56
  Light Rain Showers                   0    413
  Light Snow                           2      0
  Light Thunderstorms and Rain         0     15
  Mostly Cloudy                        5    977
  Overcast                             4      9
  Partly Cloudy                       33    428
  Rain                                 0     13
  Rain Showers                         0    108
  Scattered Clouds                    77    226
  Smoke                                0      1
  Thunderstorm                         0      4
  Thunderstorms and Rain               0      4
  Thunderstorms with Hail              1      0
  Unknown                              0     11

weather_con <- unique(subset(weather,select=c("Conditions","DateUTC","loc")))
ggplot(data=weather_con, aes(x=Conditions, fill = loc)) +
    geom_bar(position=position_dodge()) +
    theme(axis.text.x = element_text(angle = 90, vjust = .5, hjust = 1))

weather_event <- unique(subset(weather,select=c("Events","DateUTC","loc")))
ggplot(data=weather_event, aes(x=Events, fill = loc)) +
    geom_bar(position=position_dodge()) +
    #scale_y_continuous(labels=scales::percent) +
    theme(axis.text.x = element_text(angle = 90, vjust = .5, hjust = 1))

We’ve often seen students refer to ‘average mood’. Sometimes this might make sense, but this is an analogous example…

#let's take the weather event data, and code it from best ('no event') to worst ('snow')
weather_event$Event[weather_event$Event==""] <- 6
weather_event$Event[weather_event$Event=="Rain"] <- 5
weather_event$Event[weather_event$Event=="Thunderstorm"] <- 4
weather_event$Event[weather_event$Event=="Rain-Thunderstorm"] <- 3
weather_event$Event[weather_event$Event=="Hail-Thunderstorm"] <- 2
weather_event$Event[weather_event$Event=="Snow"] <- 1
weather_event$Event <- as.numeric(weather_event$Event)
ggplot(weather_event, aes(x = loc, y = Event, fill=loc)) + geom_boxplot() +
    stat_summary(fun.y=mean, geom="point", shape=5, size=4)

LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpUaGlzIGlzIGFuIFtSIE1hcmtkb3duXShodHRwOi8vcm1hcmtkb3duLnJzdHVkaW8uY29tKSBOb3RlYm9vay4gV2hlbiB5b3UgZXhlY3V0ZSBjb2RlIHdpdGhpbiB0aGUgbm90ZWJvb2ssIHRoZSByZXN1bHRzIGFwcGVhciBiZW5lYXRoIHRoZSBjb2RlLiANCg0KVHJ5IGV4ZWN1dGluZyB0aGlzIGNodW5rIGJ5IGNsaWNraW5nIHRoZSAqUnVuKiBidXR0b24gd2l0aGluIHRoZSBjaHVuayBvciBieSBwbGFjaW5nIHlvdXIgY3Vyc29yIGluc2lkZSBpdCBhbmQgcHJlc3NpbmcgKkN0cmwrU2hpZnQrRW50ZXIqLiANCg0KTG9hZCAmIGluc3RhbGwgcGFja2FnZXMNCj09PT09PT09PT09PT09PT09PT09PT09PT09PQ0KYGBge3J9DQpsaWJyYXJ5KGN1cmwpDQpsaWJyYXJ5KGRldnRvb2xzKQ0KaW5zdGFsbF9naXRodWIoIndlYXRoZXJEYXRhIiwgIlJhbS1OIikNCmxpYnJhcnkod2VhdGhlckRhdGEpDQojaW5zdGFsbC5wYWNrYWdlcygid2VhdGhlckRhdGEiKSAtIHRoZSBwZXJzb25hbCB3ZWF0aGVyIHN0YXRpb24gZGF0YSBpcyBvbmx5IGF2YWlsYWJsZSB2aWEgZ2l0aHViIHZlcnNpb24NCiNwYWNrYWdlIGluZm8gYXQgaHR0cHM6Ly9yYW0tbi5naXRodWIuaW8vd2VhdGhlckRhdGEvDQoNCmxpYnJhcnkoZ2dwbG90MikNCmxpYnJhcnkocHN5Y2gpDQpsaWJyYXJ5KHJlc2hhcGUyKQ0KbGlicmFyeShrbml0cikNCmxpYnJhcnkobGF0dGljZSkgI2p1c3QgdG8gaWxsdXN0cmF0ZSBhbm90aGVyIGhpc3RvZ3JhbSBmdW5jdGlvbg0KDQpgYGANCg0KQWRkIGEgbmV3IGNodW5rIGJ5IGNsaWNraW5nIHRoZSAqSW5zZXJ0IENodW5rKiBidXR0b24gb24gdGhlIHRvb2xiYXIgb3IgYnkgcHJlc3NpbmcgKkN0cmwrQWx0K0kqLg0KDQpXaGVuIHlvdSBzYXZlIHRoZSBub3RlYm9vaywgYW4gSFRNTCBmaWxlIGNvbnRhaW5pbmcgdGhlIGNvZGUgYW5kIG91dHB1dCB3aWxsIGJlIHNhdmVkIGFsb25nc2lkZSBpdCAoY2xpY2sgdGhlICpQcmV2aWV3KiBidXR0b24gb3IgcHJlc3MgKkN0cmwrU2hpZnQrSyogdG8gcHJldmlldyB0aGUgSFRNTCBmaWxlKS4NCg0KYGBge3J9DQpjaGVja0RhdGFBdmFpbGFiaWxpdHkoIlNZRCIsICIyMDE3LTAyLTE1Iiwgc3RhdGlvbl90eXBlID0gImFpcnBvcnRDb2RlIikgDQoNCmdldFdlYXRoZXJGb3JEYXRlKCJTWUQiLCAiMjAxNy0wMi0xNSIsIG9wdF9kZXRhaWxlZD1UUlVFKQ0KZ2V0V2VhdGhlckZvckRhdGUoIlNZRCIsICIyMDE3LTAyLTEwIiwgb3B0X2RldGFpbGVkPVRSVUUpDQoNCmdldFN1bW1hcml6ZWRXZWF0aGVyKCJTWUQiLCAiMjAxNy0wMi0xMCIpDQpgYGANCg0KR3JhYiBzb21lIGRhdGENCmBgYHtyfQ0KU1lEX3dlYXRoZXIgPC0gZ2V0V2VhdGhlckZvckRhdGUoIlNZRCIsICIyMDE3LTAyLTE1IiwgZW5kX2RhdGUgPSAiMjAxNy0wNC0wNCIsIG9wdF9kZXRhaWxlZD1UUlVFLCBvcHRfYWxsX2NvbHVtbnMgPSBUKQ0KTUVMX3dlYXRoZXIgPC0gZ2V0V2VhdGhlckZvckRhdGUoIk1FTCIsICIyMDE3LTAyLTE1IiwgZW5kX2RhdGUgPSAiMjAxNy0wNC0wNCIsIG9wdF9kZXRhaWxlZD1UUlVFLCBvcHRfYWxsX2NvbHVtbnMgPSBUKQ0KDQojZ2V0V2VhdGhlckZvckRhdGUoc3RhdGlvbl9pZCA9ICJJU1lETkVZMTQzIiwgc3RhcnRfZGF0ZSA9ICIyMDE3LTAyLTEwIiwgZW5kX2RhdGUgPSAiMjAxNy0wNC0wNCIsIG9wdF9kZXRhaWxlZD1UUlVFLCBvcHRfYWxsX2NvbHVtbnMgPSBULCBzdGF0aW9uX3R5cGUgPSAiSUQiKQ0KYGBgDQoNCg0KSG93IG1pZ2h0IHdlIGV4cGxvcmUgc29tZSBvZiB0aGlzIGRhdGE/DQpgYGB7cn0NCmthYmxlKHJiaW5kKGRlc2NyaWJlKFNZRF93ZWF0aGVyJFRlbXBlcmF0dXJlQyksZGVzY3JpYmUoTUVMX3dlYXRoZXIkVGVtcGVyYXR1cmVDKSkpDQpgYGANCg0KTGV0J3MgaGF2ZSBhIGJhc2ljIGxvb2sgYXQgdGhlIGRpc3RyaWJ1dGlvbg0KYGBge3J9DQpoaXN0KE1FTF93ZWF0aGVyJFRlbXBlcmF0dXJlQykNCmhpc3QoU1lEX3dlYXRoZXIkVGVtcGVyYXR1cmVDKQ0KYGBgDQoNCkFoLCB0aGVyZSdzIHNvbWV0aGluZyB3cm9uZyEgIExldCdzIGZpeCB0aGF0IGFuZCBnZXQgdGhlIG1lbGJvdXJuZSBhbmQgc3lkbmV5IGRhdGEgaW50byB0aGUgc2FtZSBmcmFtZQ0KYGBge3J9DQpTWURfdGVtcCA8LSBhcy5kYXRhLmZyYW1lKGFzLm51bWVyaWModW5saXN0KHN1YnNldChTWURfd2VhdGhlciwgVGVtcGVyYXR1cmVDID4tMzAwLCBzZWxlY3Q9YygiVGVtcGVyYXR1cmVDIikpKSkpICN5b3UgY291bGQgYWxzbyByZXBsYWNlIHRoZXNlIHdpaHQgTkEsIGJ1dCBoZXJlIHdlJ3JlIGp1c3QgZ29pbmcgdG8gZXhjbHVkZSB0aGUgbWlzc2luZyBkYXRhDQpjb2xuYW1lcyhTWURfdGVtcClbMV0gPC0gInRlbXAiDQpTWURfdGVtcCRsb2MgPC0gIlNZRCINCg0KTUVMX3RlbXAgPC0gYXMuZGF0YS5mcmFtZShhcy5udW1lcmljKHVubGlzdChzdWJzZXQoTUVMX3dlYXRoZXIsIFRlbXBlcmF0dXJlQyA+LTMwMCwgc2VsZWN0PWMoIlRlbXBlcmF0dXJlQyIpKSkpKQ0KY29sbmFtZXMoTUVMX3RlbXApWzFdIDwtICJ0ZW1wIg0KTUVMX3RlbXAkbG9jIDwtICJNRUwiDQoNCnRlbXBzIDwtIHJiaW5kKFNZRF90ZW1wLCBNRUxfdGVtcCkNCnRlbXBzJHRlbXAgPC0gYXMubnVtZXJpYyh0ZW1wcyR0ZW1wKQ0KYGBgDQoNCkxldCdzIGxvb2sgYXQgaXQgYWdhaW4NCmBgYHtyfQ0KZ2dwbG90KHRlbXBzLCBhZXMoeCA9IHRlbXAsIGZpbGwgPSBsb2MpKSArIGdlb21faGlzdG9ncmFtKGFscGhhID0gLjUsIHBvc2l0aW9uID0gJ2lkZW50aXR5JykgI2FoLCB3aGF0J3MgdGhlIHByb2JsZW0gaGVyZT8NCmBgYA0KDQpUcnkgdGhhdCBhZ2Fpbg0KYGBge3J9DQpnZ3Bsb3QodGVtcHMsIGFlcyh4ID0gdGVtcCwgZmlsbCA9IGxvYykpICsgZ2VvbV9oaXN0b2dyYW0oYWxwaGEgPSAuNSwgYWVzKHkgPSAuLmRlbnNpdHkuLiksIHBvc2l0aW9uID0gJ2lkZW50aXR5JykgI25vdGUgdXNlIG9mICdkZW5zaXR5JyBiZWNhdXNlIHdlIGhhdmUgdW5lcXVhbCB0ZW1wZXJhdHVyZSBjb3VudHMgaW4gZWFjaCBkYXRhc2V0LiBBbHBoYSBpcyB0aGUgdHJhbnNwYXJlbmN5IGxldmVsLg0KDQpoaXN0b2dyYW0ofiB0ZW1wIHwgbG9jLCBkYXRhPXRlbXBzKQ0KDQpgYGANCg0KV2hhdCdzIHdyb25nIHdpdGggdGhpcz8NCmBgYHtyfQ0KZ2dwbG90KHRlbXBzKSArIA0KICBnZW9tX2JhcihhZXMoeCA9IGxvYywgeSA9IHRlbXAsIGZpbGwgPSBsb2MpLA0KICAgICAgICAgICBwb3NpdGlvbiA9ICJkb2RnZSIsIHN0YXQgPSAic3VtbWFyeSIsIGZ1bi55ID0gIm1lYW4iKQ0KYGBgDQoNCk1vcmUgaW5mb3JtYXRpdmU/DQpgYGB7cn0NCmdncGxvdCh0ZW1wcywgYWVzKHg9bG9jLCB5PXRlbXAsIGZpbGw9bG9jKSkgKyBnZW9tX2JveHBsb3QoKSArDQogICAgZ3VpZGVzKGZpbGw9RkFMU0UpKw0KICAgIHN0YXRfc3VtbWFyeShmdW4ueT1tZWFuLCBnZW9tPSJwb2ludCIsIHNoYXBlPTUsIHNpemU9NCkNCmBgYA0KDQpDb3VwbGUgb2YgdXNlZnVsIHRoaW5ncyAtIGxldCdzIHB1bGwgdGhlIGRhdGUgb3V0IHRvIGl0cyBvd24gdmFsdWUsIGFuZCB0aGlzIHRpbWUgd2UnbGwgcmVwbGFjZSBtaXNzaW5nIHZhbHVlcyAoLTk5OTkpIHdpdGggTkENCmBgYHtyfQ0KZ2dwbG90KFNZRF93ZWF0aGVyLCBhZXMoeD1UZW1wZXJhdHVyZUMsIHk9RGV3X1BvaW50QykpICsNCiAgICBnZW9tX3BvaW50KHNoYXBlPTEpICAgICAgIyBVc2UgaG9sbG93IGNpcmNsZXMNCg0KU1lEX3dlYXRoZXJbU1lEX3dlYXRoZXIgPT0gLTk5OTldIDwtIE5BDQpTWURfd2VhdGhlciRkYXRlIDwtIGFzLkRhdGUoU1lEX3dlYXRoZXIkRGF0ZVVUQykNCg0KTUVMX3dlYXRoZXJbTUVMX3dlYXRoZXIgPT0gLTk5OTldIDwtIE5BDQpNRUxfd2VhdGhlciRkYXRlIDwtIGFzLkRhdGUoTUVMX3dlYXRoZXIkRGF0ZVVUQykNCg0KZ2dwbG90KFNZRF93ZWF0aGVyLCBhZXMoeD1UZW1wZXJhdHVyZUMsIHk9RGV3X1BvaW50QykpICsNCiAgICBnZW9tX3BvaW50KHNoYXBlPTEpICAgICAgIyBVc2UgaG9sbG93IGNpcmNsZXMNCmBgYA0KDQpXaGF0IGlmIHdlIHdhbnQgdG8gZXhwbG9yZSB0aGUgcmVsYXRpb25zaGlwIGJldHdlZW4gRGV3X1BvaW50QyBhbmQgb3RoZXIgZmVhdHVyZXMNCmh0dHBzOi8vc3VwcG9ydC5vZmZpY2UuY29tL2VuLXVzL2FydGljbGUvUHJlc2VudC15b3VyLWRhdGEtaW4tYS1zY2F0dGVyLWNoYXJ0LW9yLWEtbGluZS1jaGFydC00NTcwYTgwZi01OTlhLTRkNmItYTE1NS0xMDRhOTAxOGI4NmUNCg0KT25lIHdheSB5b3UgbWlnaHQgYmUgdGVtcHRlZCB0byBkbyB0aGlzLi4uDQpgYGB7cn0NCmJhZF9leGFtcGxlIDwtIHN1YnNldChTWURfd2VhdGhlciwgIWlzLm5hKEh1bWlkaXR5KSwgc2VsZWN0PWMoIkh1bWlkaXR5IiwgIkRld19Qb2ludEMiLCJkYXRlIikpDQpiYWRfZXhhbXBsZVtjKCJIdW1pZGl0eSIsIkRld19Qb2ludEMiKV0gPC0gbGFwcGx5KGJhZF9leGFtcGxlW2MoIkh1bWlkaXR5IiwiRGV3X1BvaW50QyIpXSxhcy5udW1lcmljKQ0KDQpiYWRfZXhhbXBsZSA8LSBhZ2dyZWdhdGUoLiB+IGRhdGUsIGJhZF9leGFtcGxlLCBGVU49bWVhbikNCg0KI2NvbnZlcnQgdG8gbG9uZw0KYmFkX2V4YW1wbGUgPC0gbWVsdChiYWRfZXhhbXBsZSwgaWQudmFycyA9IGMoImRhdGUiKSkNCg0KZ2dwbG90KGRhdGE9YmFkX2V4YW1wbGUsIGFlcyh4PWRhdGUsIHk9dmFsdWUsIGdyb3VwPXZhcmlhYmxlLCBjb2xvdXI9dmFyaWFibGUpKSArDQogICAgZ2VvbV9saW5lKCkgKw0KICAgIGdlb21fcG9pbnQoKQ0KI0lzIGRhdGUgYW4gaW1wb3J0YW50IHZhcmlhYmxlIGluIHRoaXMgYW5hbHlzaXM/IERvZXMgdGhlIHNjYWxpbmcgb2YgdGhlIGRhdGEgZ2l2ZXMgdXMgdGhlIGJlc3QgYXZhaWxhYmxlIGluc2lnaHQgaW50byByZWxhdGlvbnNoaXBzIG9mIHBhaXJlZCB2YWx1ZXM/IElzIHRoZSB1c2Ugb2YgYSBsaW5lIHRvIGpvaW4gZGF0YXBvaW50cyBhcHByb3ByaWF0ZSBnaXZlbiBtaXNzaW5nIGRhdGE/DQpgYGANCg0KQSBiZXR0ZXIgd2F5Pw0KYGBge3J9DQpnZ3Bsb3QoU1lEX3dlYXRoZXIsIGFlcyh4PUh1bWlkaXR5LCB5PURld19Qb2ludEMpKSArDQogICAgZ2VvbV9wb2ludChzaGFwZT0xKSAgICAgICMgVXNlIGhvbGxvdyBjaXJjbGVzDQoNCmNvci50ZXN0KGFzLm51bWVyaWMoU1lEX3dlYXRoZXIkSHVtaWRpdHkpLGFzLm51bWVyaWMoU1lEX3dlYXRoZXIkRGV3X1BvaW50QykpDQoNCmBgYA0KDQpPayB3aGF0IGlmIHdlIHdhbnQgdG8gbG9vayBhdCBob3cgd2VhdGhlciB2YXJpZXMgb3ZlciB0aW1lIGFuZCBwbGFjZT8NCmBgYHtyfQ0KU1lEX3dlYXRoZXIkbG9jIDwtICJTeWRuZXkiDQpNRUxfd2VhdGhlciRsb2MgPC0gIk1lbGJvdXJuZSINCmBgYA0KDQpgYGB7cn0NCndlYXRoZXIgPC0gcmJpbmQoU1lEX3dlYXRoZXJbYygiVGVtcGVyYXR1cmVDIiwiRGV3X1BvaW50QyIsIkh1bWlkaXR5IiwiV2luZF9TcGVlZEttX2giLCJQcmVjaXBpdGF0aW9ubW0iLCJFdmVudHMiLCJDb25kaXRpb25zIiwiRGF0ZVVUQyIsImxvYyIpXSxNRUxfd2VhdGhlcltjKCJUZW1wZXJhdHVyZUMiLCJEZXdfUG9pbnRDIiwiSHVtaWRpdHkiLCJXaW5kX1NwZWVkS21faCIsIlByZWNpcGl0YXRpb25tbSIsIkV2ZW50cyIsIkNvbmRpdGlvbnMiLCJEYXRlVVRDIiwibG9jIildKQ0Kd2VhdGhlciRtb250aCA8LSBmb3JtYXQoYXMuRGF0ZSh3ZWF0aGVyJERhdGVVVEMpLCAiJW0iKQ0KDQpnZ3Bsb3Qod2VhdGhlciwgYWVzKHg9bW9udGgsIHk9VGVtcGVyYXR1cmVDLCBmaWxsPWxvYykpICsgZ2VvbV9ib3hwbG90KCkgKw0KICAgIGd1aWRlcyhmaWxsPUZBTFNFKSArDQogICAgc3RhdF9zdW1tYXJ5KGZ1bi55PW1lYW4sIGdlb209InBvaW50Iiwgc2hhcGU9NSwgc2l6ZT00KSArDQogICAgZmFjZXRfd3JhcCh+bG9jKQ0KDQpgYGANCg0KT3IgYXQgaG93IHdlYXRoZXIgZXZlbnRzIHZhcnkgYnkgcGxhY2UNCmBgYHtyfQ0KdW5pcXVlKHdlYXRoZXIkRXZlbnRzKQ0KdW5pcXVlKHdlYXRoZXIkQ29uZGl0aW9ucykNCg0KdGFibGUod2VhdGhlciRDb25kaXRpb25zLHdlYXRoZXIkbG9jKQ0KDQp3ZWF0aGVyX2NvbiA8LSB1bmlxdWUoc3Vic2V0KHdlYXRoZXIsc2VsZWN0PWMoIkNvbmRpdGlvbnMiLCJEYXRlVVRDIiwibG9jIikpKQ0KDQpnZ3Bsb3QoZGF0YT13ZWF0aGVyX2NvbiwgYWVzKHg9Q29uZGl0aW9ucywgZmlsbCA9IGxvYykpICsNCiAgICBnZW9tX2Jhcihwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSgpKSArDQogICAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCwgdmp1c3QgPSAuNSwgaGp1c3QgPSAxKSkNCg0Kd2VhdGhlcl9ldmVudCA8LSB1bmlxdWUoc3Vic2V0KHdlYXRoZXIsc2VsZWN0PWMoIkV2ZW50cyIsIkRhdGVVVEMiLCJsb2MiKSkpDQpnZ3Bsb3QoZGF0YT13ZWF0aGVyX2V2ZW50LCBhZXMoeD1FdmVudHMsIGZpbGwgPSBsb2MpKSArDQogICAgZ2VvbV9iYXIocG9zaXRpb249cG9zaXRpb25fZG9kZ2UoKSkgKw0KICAgICNzY2FsZV95X2NvbnRpbnVvdXMobGFiZWxzPXNjYWxlczo6cGVyY2VudCkgKw0KICAgIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTAsIHZqdXN0ID0gLjUsIGhqdXN0ID0gMSkpDQoNCmBgYA0KDQpXZSd2ZSBvZnRlbiBzZWVuIHN0dWRlbnRzIHJlZmVyIHRvICdhdmVyYWdlIG1vb2QnLiBTb21ldGltZXMgdGhpcyBtaWdodCBtYWtlIHNlbnNlLCBidXQgdGhpcyBpcyBhbiBhbmFsb2dvdXMgZXhhbXBsZS4uLg0KYGBge3J9DQojbGV0J3MgdGFrZSB0aGUgd2VhdGhlciBldmVudCBkYXRhLCBhbmQgY29kZSBpdCBmcm9tIGJlc3QgKCdubyBldmVudCcpIHRvIHdvcnN0ICgnc25vdycpDQp3ZWF0aGVyX2V2ZW50JEV2ZW50W3dlYXRoZXJfZXZlbnQkRXZlbnQ9PSIiXSA8LSA2DQp3ZWF0aGVyX2V2ZW50JEV2ZW50W3dlYXRoZXJfZXZlbnQkRXZlbnQ9PSJSYWluIl0gPC0gNQ0Kd2VhdGhlcl9ldmVudCRFdmVudFt3ZWF0aGVyX2V2ZW50JEV2ZW50PT0iVGh1bmRlcnN0b3JtIl0gPC0gNA0Kd2VhdGhlcl9ldmVudCRFdmVudFt3ZWF0aGVyX2V2ZW50JEV2ZW50PT0iUmFpbi1UaHVuZGVyc3Rvcm0iXSA8LSAzDQp3ZWF0aGVyX2V2ZW50JEV2ZW50W3dlYXRoZXJfZXZlbnQkRXZlbnQ9PSJIYWlsLVRodW5kZXJzdG9ybSJdIDwtIDINCndlYXRoZXJfZXZlbnQkRXZlbnRbd2VhdGhlcl9ldmVudCRFdmVudD09IlNub3ciXSA8LSAxDQoNCndlYXRoZXJfZXZlbnQkRXZlbnQgPC0gYXMubnVtZXJpYyh3ZWF0aGVyX2V2ZW50JEV2ZW50KQ0KDQpnZ3Bsb3Qod2VhdGhlcl9ldmVudCwgYWVzKHggPSBsb2MsIHkgPSBFdmVudCwgZmlsbD1sb2MpKSArIGdlb21fYm94cGxvdCgpICsNCiAgICBzdGF0X3N1bW1hcnkoZnVuLnk9bWVhbiwgZ2VvbT0icG9pbnQiLCBzaGFwZT01LCBzaXplPTQpDQpgYGANCg==

R Notebook

Load & install packages