Session 4 R practice

Logistic models

To run your code, click run. It will let you know if your code is correct or not, and will offer hints on what to fix if needed. If you are stuck, you can see the solution by clicking Solution.

Exercise 1: Using mutate()

Reminder:

\[ yoy\_growth = \frac{level_t - level_{t-4}}{level_{t-4}} \]

eyJsYW5ndWFnZSI6InIiLCJwcmVfZXhlcmNpc2VfY29kZSI6ImRmIDwtIGRhdGEuZnJhbWUodGljPXJlcChcIldNVFwiLDE5NCksXG4gICAgICAgICAgICAgICAgIGZ5ZWFycT1jKDE5NzAsMTk3MCwxOTcwLDE5NzAsMTk3MSwxOTcxLDE5NzEsMTk3MSwxOTcyLDE5NzIsMTk3MiwxOTcyLDE5NzMsMTk3MywxOTczLDE5NzMsMTk3NCwxOTc0LDE5NzQsMTk3NCwxOTc1LDE5NzUsMTk3NSwxOTc1LDE5NzYsMTk3NiwxOTc2LDE5NzYsMTk3NywxOTc3LDE5NzcsMTk3NywxOTc4LDE5NzgsMTk3OCwxOTc4LDE5NzksMTk3OSwxOTc5LDE5NzksMTk4MCwxOTgwLDE5ODAsMTk4MCwxOTgxLDE5ODEsMTk4MSwxOTgxLDE5ODIsMTk4MiwxOTgyLDE5ODIsMTk4MywxOTgzLDE5ODMsMTk4MywxOTg0LDE5ODQsMTk4NCwxOTg0LDE5ODUsMTk4NSwxOTg1LDE5ODUsMTk4NiwxOTg2LDE5ODYsMTk4NiwxOTg3LDE5ODcsMTk4NywxOTg3LDE5ODgsMTk4OCwxOTg4LDE5ODgsMTk4OSwxOTg5LDE5ODksMTk4OSwxOTkwLDE5OTAsMTk5MCwxOTkwLDE5OTEsMTk5MSwxOTkxLDE5OTEsMTk5MiwxOTkyLDE5OTIsMTk5MiwxOTkzLDE5OTMsMTk5MywxOTkzLDE5OTQsMTk5NCwxOTk0LDE5OTQsMTk5NSwxOTk1LDE5OTUsMTk5NSwxOTk2LDE5OTYsMTk5NiwxOTk2LDE5OTcsMTk5NywxOTk3LDE5OTcsMTk5OCwxOTk4LDE5OTgsMTk5OCwxOTk5LDE5OTksMTk5OSwxOTk5LDIwMDAsMjAwMCwyMDAwLDIwMDAsMjAwMSwyMDAxLDIwMDEsMjAwMSwyMDAyLDIwMDIsMjAwMiwyMDAyLDIwMDMsMjAwMywyMDAzLDIwMDMsMjAwNCwyMDA0LDIwMDQsMjAwNCwyMDA1LDIwMDUsMjAwNSwyMDA1LDIwMDYsMjAwNiwyMDA2LDIwMDYsMjAwNywyMDA3LDIwMDcsMjAwNywyMDA4LDIwMDgsMjAwOCwyMDA4LDIwMDksMjAwOSwyMDA5LDIwMDksMjAxMCwyMDEwLDIwMTAsMjAxMCwyMDExLDIwMTEsMjAxMSwyMDExLDIwMTIsMjAxMiwyMDEyLDIwMTIsMjAxMywyMDEzLDIwMTMsMjAxMywyMDE0LDIwMTQsMjAxNCwyMDE0LDIwMTUsMjAxNSwyMDE1LDIwMTUsMjAxNiwyMDE2LDIwMTYsMjAxNiwyMDE3LDIwMTcsMjAxNywyMDE3LDIwMTgsMjAxOCksXG4gICAgICAgICAgICAgICAgIGZxdHI9YygxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIpLFxuICAgICAgICAgICAgICAgICByZXZ0cT1jKDcuMjY3LDkuOTM5LDExLjQ3LDE1LjYxLDEyLjc1OCwxOC4yMDcsMTkuMjA1LDI3Ljg0NCwyMS4wMTEsMjkuOTE4LDMyLjY4Myw0MS4yNzcsMzIuNzM3LDM3LjY1NSw0Mi4wOTgsNTUuMDcxLDQ2LjM1NSw1Ni4xNzcsNTkuNDYyLDc0LjIxNSw2NC42NSw3OC44NTIsODYuMTA4LDExMC43MjEsOTMuOTcxLDExNS4xODYsMTIxLjg3NSwxNDcuNzc1LDEyMi42NTUsMTUyLjM4MSwxNzIuNjM3LDIzMC43ODMsMTY3LjI3NCwyMDQuMDg1LDIyMS45ODIsMzA2Ljk1NywyMzguNTY5LDI5MS42ODUsMzAwLjc0Nyw0MTcuMTc1LDMxNC41NzUsMzcyLjU5OCw0MDEuODQyLDU1NC4xODQsNDM0LjI1MSw1MzMuNjUzLDYzNC40NzYsODQyLjYxNyw2NjUuMjI5LDc4OC43MDUsODI3Ljc2NywxMDk0LjU1MSw4NTQuNzM0LDEwOTguOTQyLDExNjYuNTk0LDE1NDYuNjM5LDEyMzQuNzY4LDE1MDguNTMzLDE1ODMuNTczLDIwNzMuOTgzLDE2NTUuNjYxLDE5MzQuMzU4LDIwODcuNTMzLDI3NzMuOTM2LDIzNDQuMjUzLDI3NzAuMjc2LDI5NDkuMDA4LDM4NDUuNTU3LDMyMTEuNzYxLDM3MzIuMjM3LDQwMTUuMDE3LDUwMDAuMjM4LDQyODIuMDEyLDQ4NDUuMTMzLDQ5OTEuMTEzLDY1MzAuNzM0LDUzNzMuMjU4LDYwNDYuNDEsNjI4My40OTYsODEwNy40ODQsNjc2OC4xOTEsNzU0My41MDgsNzkzMC45NDksMTAzNTguODk4LDkyODAuNTY2LDEwMzM5Ljg5OCwxMDYyNy41LDEzNjM4Ljc5NywxMTY0OS4zOTgsMTMwMjguMzk4LDEzNjgzLjc5NywxNzEyMiwxMzkyMC4zOTgsMTYyMzYuNSwxNjgyNi44OTgsMjAzNjAuNSwxNzY4Ni4wOTgsMTk5NDIuMjk3LDIwNDE3LjY5OSwyNDQ0Ny43OTcsMjA0NDAsMjI3MjMsMjI5MTQsMjc1NTAsMjI3NzIsMjU1ODcsMjU2NDQsMzA4NTYsMjU0MDksMjgzODYsMjg3NzcsMzUzODYsMjk4MTksMzM1MjEsMzM1MDksNDA3ODUsMzUxMjksMzg5MTMsNDA4OTksNTE4NjgsNDM0NDcsNDY1ODgsNDYxODEsNTcwNzksNDg1NjUsNTMxODcsNTMxODUsNjQ3MzUsNTIxMjYsNTY3ODEsNTU3NjUsNjY5MDUsNTcyMjQsNjMyMzEsNjMwMzYsNzUxOTAsNjU0MzQsNzA1MTYsNjkyODIsODI5MDAsNzA3NTUsNzY2OTcsNzUzOTcsODkyNTIsNzk2NzUsODU0MzAsODQ0NjcsOTkwNzgsODYzNzgsOTI4MjcsOTE3OTQsMTA3Mjg5LDk0OTQwLDEwMjM0Miw5ODM0NSwxMDg3NDcsOTQyMTQsMTAwODc2LDk5MzczLDExMzYyMiw5OTgxMSwxMDM3MjYsMTAxOTUyLDExNjM2MCwxMDQxODksMTA5MzY2LDExMDIyNiwxMjMxNjksMTEzMDEwLDExNDA2NSwxMTM4MDAsMTI3Nzc2LDExNDA3MCwxMTY4MzAsMTE1Njg4LDEyOTcwNiwxMTQ5NjAsMTIwMTI1LDExOTAwMSwxMzE1NjUsMTE0ODI2LDEyMDIyOSwxMTc0MDgsMTI5NjY3LDExNTkwNCwxMjAzMTksMTE4MTc5LDEzMDkzNiwxMTc1NDIsMTIyOTY4LDEyMzE3OSwxMzYyNjcsMTIyNjkwLDEyODAyOCksXG4gICAgICAgICAgICAgICAgIHJlY2Vzc2lvbj1jKDEsIDEsIDEsIDEsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDEsIDEsIDEsIDEsIDEsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDEsIDEsIDAsIDAsIDAsIDAsIDEsIDEsIDEsIDEsIDEsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDEsIDEsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDEsIDEsIDEsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDEsIDEsIDEsIDEsIDEsIDEsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDApKSIsInNhbXBsZSI6IiMgVGhlIGRhdGFmcmFtZSBiZWxvdyBjb250YWlucyBoaXN0b3JpY2FsIHJldmVudWUgZm9yIFdhbG1hcnQgc2luY2UgMTk3MFxuIyBBbGwgZGF0YSBpcyBxdWFydGVybHlcbiMgUmVjZXNzaW9uIGRhdGEgY29tZXMgZnJvbSB0aGUgVVMgRmVkZXJhbCBSZXNlcnZlXG4jIFRoZSBmaWxlIHVzZXMgQ29tcHVzdGF0J3Mgbm9tZW5jbGF0dXJlOlxuIyAgICAgLSByZXZlbnVlIGlzIHJldnRxXG4jICAgICAtIHRpY2tlciBpcyB0aWNcbiMgICAgIC0gZnllYXJxIGlzIHllYXJcbiMgICAgIC0gZnF0ciBpcyBxdWFydGVyXG4jICAgICAtIHJlY2Vzc2lvbiBpcyBhbiBpbmRpY2F0b3IgZm9yIFVTIHJlY2Vzc2lvbnNcbmRmWzEsNV1cblxubGlicmFyeShkcGx5cikgIyBsb2FkcyBkcGx5ciBmb3IgZGF0YSBtYW5pcHVsYXRpb25cblxuIyBDcmVhdGUgYSB2YXJpYWJsZSBmb3IgeWVhciBvdmVyIHllYXIgcmV2ZW51ZSBncm93dGggYnkgcXVhcnRlclxuIyBTdG9yZSB0aGUgcmVzdWx0IGluIGEgZGF0YSBmcmFtZSBjb2x1bW4gbmFtZWQgcmV2dHFfeW95XG5cblxuIyBDcmVhdGUgYW4gaW5kaWNhdG9yIGVxdWFsIHRvIDEgaWYgcmV2dHFfeW95IGlzIGxlc3MgdGhhbiAwLCBhbmQgMCBvdGhlcndpc2VcbiMgVXNpbmcgaWZlbHNlIG1ha2VzIHRoaXMgZWFzeSB0byBkb1xuIyBTdG9yZSB0aGUgcmVzdWx0IGluIGEgZGF0YSBmcmFtZSBjb2x1bW4gbmFtZWQgcmV2dHFfZG93blxuXG5cbiMgQ3JlYXRlIGEgbGFnIG9mIHJlY2Vzc2lvbiwgY2FsbCBpdCByZWNlc3Npb25fbGFnXG5cblxuI1ZpZXcgdGhlIG91dHB1dFxuaGVhZChkZilcblxuI0VORCIsInNvbHV0aW9uIjoiIyBUaGUgZGF0YWZyYW1lIGJlbG93IGNvbnRhaW5zIGhpc3RvcmljYWwgcmV2ZW51ZSBmb3IgV2FsbWFydCBzaW5jZSAxOTcwXG4jIEFsbCBkYXRhIGlzIHF1YXJ0ZXJseVxuIyBSZWNlc3Npb24gZGF0YSBjb21lcyBmcm9tIHRoZSBVUyBGZWRlcmFsIFJlc2VydmVcbiMgVGhlIGZpbGUgdXNlcyBDb21wdXN0YXQncyBub21lbmNsYXR1cmU6XG4jICAgICAtIHJldmVudWUgaXMgcmV2dHFcbiMgICAgIC0gdGlja2VyIGlzIHRpY1xuIyAgICAgLSBmeWVhcnEgaXMgeWVhclxuIyAgICAgLSBmcXRyIGlzIHF1YXJ0ZXJcbiMgICAgIC0gcmVjZXNzaW9uIGlzIGFuIGluZGljYXRvciBmb3IgVVMgcmVjZXNzaW9uc1xuZGZbMSw1XVxuXG5saWJyYXJ5KGRwbHlyKSAjIGxvYWRzIGRwbHlyIGZvciBkYXRhIG1hbmlwdWxhdGlvblxuXG4jIENyZWF0ZSBhIHZhcmlhYmxlIGZvciB5ZWFyIG92ZXIgeWVhciByZXZlbnVlIGdyb3d0aCBieSBxdWFydGVyXG4jIFN0b3JlIHRoZSByZXN1bHQgaW4gYSBkYXRhIGZyYW1lIGNvbHVtbiBuYW1lZCByZXZ0cV95b3lcbmRmIDwtIGRmICU+JSBtdXRhdGUocmV2dHFfeW95ID0gKHJldnRxIC0gbGFnKHJldnRxLCA0KSkgLyBsYWcocmV2dHEsIDQpKVxuXG4jIENyZWF0ZSBhbiBpbmRpY2F0b3IgZXF1YWwgdG8gMSBpZiByZXZ0cV95b3kgaXMgbGVzcyB0aGFuIDAsIGFuZCAwIG90aGVyd2lzZVxuIyBTdG9yZSB0aGUgcmVzdWx0IGluIGEgZGF0YSBmcmFtZSBjb2x1bW4gbmFtZWQgcmV2dHFfZG93blxuZGYgPC0gZGYgJT4lIG11dGF0ZShyZXZ0cV9kb3duID0gaWZlbHNlKHJldnRxX3lveTwwLDEsMCkpXG5cbiMgQ3JlYXRlIGEgMSBxdWFydGVyIGxhZyBvZiByZWNlc3Npb24sIGNhbGwgaXQgcmVjZXNzaW9uX2xhZ1xuZGYgPC0gZGYgJT4lIG11dGF0ZShyZWNlc3Npb25fbGFnID0gbGFnKHJlY2Vzc2lvbikpXG5cbiNWaWV3IHRoZSBvdXRwdXRcbmhlYWQoZGYpXG5cbiNFTkQiLCJzY3QiOiIjIFRlbXBsYXRlIGJhc2VkIG9uIGh0dHBzOi8vd3d3LnJkb2N1bWVudGF0aW9uLm9yZy9wYWNrYWdlcy90ZXN0d2hhdC92ZXJzaW9ucy80LjEuMVxuIyBDaGVjayBpZiBzb21ldGhpbmcgaXMgZXhwbGljaXRseSB0eXBlZFxuXG50ZXN0X3N0dWRlbnRfdHlwZWQoXCI8LVwiLCBub3RfdHlwZWRfbXNnPSdEaWQgeW91IHJlbWVtYmVyIHRvIHN0b3JlIHRoZSByZXN1bHRzIGJhY2sgaW50byB0aGUgZGF0YSBmcmFtZT8nKVxudGVzdF9leHByZXNzaW9uX291dHB1dChcIidyZXZ0cV95b3knICVpbiUgY29sbmFtZXMoZGYpXCIsIGluY29ycmVjdF9tc2c9XCJEaWQgeW91IHJlbWVtYmVyIHRvIHNhdmUgdGhlIG91dHB1dCB0byB0aGUgZGF0YWZyYW1lLCBhbmQgdG8gbmFtZSB0aGUgY29sdW1uIGByZXZ0cV95b3lgP1wiKVxudGVzdF9leHByZXNzaW9uX291dHB1dChcImRmJHJldnRxX3lveVwiLCBpbmNvcnJlY3RfbXNnPVwiSXQgbG9va3MgbGlrZSBlaXRoZXIgdGhlIGNvbHVtbiBuYW1lIGlzIHdyb25nIG9yIHRoZSBjYWxjdWxhdGlvbiBpcyB3cm9uZyBmb3IgYHJldnRxX3lveWAuXCIpXG5cbnRlc3RfZXhwcmVzc2lvbl9vdXRwdXQoXCIncmV2dHFfZG93bicgJWluJSBjb2xuYW1lcyhkZilcIiwgaW5jb3JyZWN0X21zZz1cIkRpZCB5b3UgcmVtZW1iZXIgdG8gc2F2ZSB0aGUgb3V0cHV0IHRvIHRoZSBkYXRhZnJhbWUsIGFuZCB0byBuYW1lIHRoZSBjb2x1bW4gYHJldnRxX2Rvd25gP1wiKVxudGVzdF9leHByZXNzaW9uX291dHB1dChcImRmJHJldnRxX2Rvd25cIiwgaW5jb3JyZWN0X21zZz1cIkl0IGxvb2tzIGxpa2UgZWl0aGVyIHRoZSBjb2x1bW4gbmFtZSBpcyB3cm9uZyBvciB0aGUgY2FsY3VsYXRpb24gaXMgd3JvbmcgZm9yIGByZXZ0cV9kb3duYC5cIilcblxudGVzdF9leHByZXNzaW9uX291dHB1dChcIidyZWNlc3Npb25fbGFnJyAlaW4lIGNvbG5hbWVzKGRmKVwiLCBpbmNvcnJlY3RfbXNnPVwiRGlkIHlvdSByZW1lbWJlciB0byBzYXZlIHRoZSBvdXRwdXQgdG8gdGhlIGRhdGFmcmFtZSwgYW5kIHRvIG5hbWUgdGhlIGNvbHVtbiBgcmVjZXNzaW9uX2xhZ2A/XCIpXG50ZXN0X2V4cHJlc3Npb25fb3V0cHV0KFwiZGYkcmVjZXNzaW9uX2xhZ1wiLCBpbmNvcnJlY3RfbXNnPVwiSXQgbG9va3MgbGlrZSBlaXRoZXIgdGhlIGNvbHVtbiBuYW1lIGlzIHdyb25nIG9yIHRoZSBjYWxjdWxhdGlvbiBpcyB3cm9uZyBmb3IgYHJlY2Vzc2lvbl9sYWdgLlwiKVxudGVzdF9zdHVkZW50X3R5cGVkKCdoZWFkKGRmKScsIG5vdF90eXBlZF9tc2c9J0RvblxcJ3QgcmVtb3ZlIHRoZSBjb2RlIGBoZWFkKGRmKWAuJylcblxuIyB0ZXN0X3N0dWRlbnRfdHlwZWQoJ3ggPC0gMicsIG5vdF90eXBlZF9tc2c9JycpXG5cbiMgQ2hlY2sgaWYgZnVuY3Rpb24gd2FzIHVzZWQgaW4gaW5wdXQgY29kZVxuIyB0ZXN0X2Z1bmN0aW9uKCdjJyxpbmNvcnJlY3RfbXNnPScnKSAgXG5cbiMgUmVxdWlyZXMgYW4gb2JqZWN0IGB4YCB0byBoYXZlIHRoZSBzYW1lIHZhbHVlIGFzIHRoZSBzb2x1dGlvblxuIyB0ZXN0X29iamVjdChcInhcIixpbmNvcnJlY3RfbXNnID0gXCJcIix1bmRlZmluZWRfbXNnID0gXCJcIikgIFxuXG4jIFJlcXVpcmVzIGFuIG9uamVjdCB3aXRoIHRoZSBzYW1lIHZhbHVlIG9mIGB4YCBpbiB0aGUgc29sdXRpb25cbiMgdGVzdF9hbl9vYmplY3QoXCJ4XCIsdW5kZWZpbmVkX21zZz1cIlwiKVxuXG4jIENoZWNrcyBpZiBvdXRwdXQgb2Ygc3R1ZGVudCdzIGNvZGUgY29udGFpbnMgZ2l2ZW4gZXZhbHVhdGVkIGV4cHJlc3Npb25cbiMgdGVzdF9vdXRwdXRfY29udGFpbnMoXCJ4XCIsaW5jb3JyZWN0X21zZyA9IFwiXCIpXG5cbiMgQ2hlY2sgaWYgYSB2ZWN0b3Igb2YgcHJlZGVmaW5lZCBvYmplY3RzIGFyZSB1bmNoYW5nZWRcbiMgdGVzdF9wcmVkZWZpbmVkX29iamVjdHMoYygneCcsJ3knKSxpbmNvcnJlY3RfbXNnPVwiRG9uJ3Qgb252ZXJ3cml0ZSB0aGUgcHJlZGVmaW5lZCB2YXJpYWJsZXNcIilcblxuIyBDaGVja3MgZm9yIGEgcmVnZXggcGF0dGVybiBpbiB0cmhlIG91dHB1dFxuIyB0ZXN0X291dHB1dF9yZWdleChwYXR0ZXJuLGZpeGVkPUYsIHRpbWVzPTEsIGluY29ycmVjdF9tc2c9JycpXG5cbiMgQ2FuIGNoZWNrIGFuIGFyYml0cmFyeSBleHByZXNzaW9uIGFjcm9zcyBib3RoIHNvbHV0aW9uIGFuZCBzdHVkZW50IGNvZGVcbiN0ZXN0X2V4cHJlc3Npb25fb3V0cHV0KFwidHlwZW9mKGNvbXBhbnlfbmFtZSlcIiwgaW5jb3JyZWN0X21zZz1cIkRpZCB5b3Ugc3RvcmUgdGV4dHVhbCBkYXRhIGluIGBjb21wYW55X25hbWVgP1wiKVxuXG50ZXN0X2Vycm9yKClcbnN1Y2Nlc3NfbXNnKFwiQXdlc29tZSFcIilcblxuIyBPdGhlciBmdW5jdGlvbnMgdG8gbm90ZTpcbiMgICAgIC0gdGVzdF9vcihhLGIpIC0tIGNoZWNrcyBpZiBlaXRoZXIgdGVzdCBhIG9yIHRlc3QgYiBwYXNzXG4jICAgICAtIHRlc3RfZ2dwbG90KCkgLS0gY2FuIGNoZWNrIGlmIHBsb3RzIGFyZSBjb3JyZWN0XG4jICAgICAtIHRlc3RfZnVuY3Rpb24oKSAtLSBjYW4gYWxzbyBjaGVjayBpbmNsdWRlZCBwYXJhbWV0ZXJzXG4jICAgICAtIHRlc3RfbG9vcCgpIC0tIGNoZWNraW5nIGZvciBhbmQgd2hpbGUgbG9vcHNcbiMgICAgIC0gdGVzdF9saWJyYXJ5X2Z1bmN0aW9uKCdwYWNrYWdlJywgbm90X2NhbGxlZF9tc2c9JycsaW5jb3JyZWN0X21zZz0nJylcbiMgICAgIC0gdGVzdF9pZl9lbHNlKCkgLS0gY2hlY2tpbmcgaWYgc3RhdGVtZW50c1xuIyAgICAgLSB0ZXN0X2V4cHJlc3Npb25fZXJyb3IoKSAtLSBjYW4gY2hlY2sgaWYgZnVuY3Rpb25zIGFyZSBwcm9wZXJseSBkZWZpbmVkXG4jICAgICAtIHRlc3Rfb3BlcmF0b3IoJ29wZXJhdG9yJywpLCBub3RfY2FsbGVkX21zZz0nJyxpbmNvcnJlY3RfbXNnPScnKVxuIyAgICAgLSB0ZXN0X2Z1bmN0aW9uX2RlZmluaXRpb24oKSAtLSByaWdvcm91c2x5IGNoZWNrIGRlZmluZWQgZnVuY3Rpb25cbiMgICAgIC0gdGVzdF9kYXRhX2ZyYW1lKCkgLS0gY2hlY2sgaWYgZGF0YWZyYW1lIFtjb2x1bW5zXSBhcmUgZXF1aXZhbGVudFxuIyAgICAgLSB0ZXN0X2Z1bmN0aW9uX3Jlc3VsdCwgdGVzdF9leHByZXNzaW9uX3Jlc3VsdCJ9

Exercise 2: Running a logit regression

eyJsYW5ndWFnZSI6InIiLCJwcmVfZXhlcmNpc2VfY29kZSI6ImxpYnJhcnkoZHBseXIpICMgbG9hZHMgZHBseXIgZm9yIGRhdGEgbWFuaXB1bGF0aW9uXG5kZiA8LSBkYXRhLmZyYW1lKHRpYz1yZXAoXCJXTVRcIiwxOTQpLFxuICAgICAgICAgICAgICAgICBmeWVhcnE9YygxOTcwLDE5NzAsMTk3MCwxOTcwLDE5NzEsMTk3MSwxOTcxLDE5NzEsMTk3MiwxOTcyLDE5NzIsMTk3MiwxOTczLDE5NzMsMTk3MywxOTczLDE5NzQsMTk3NCwxOTc0LDE5NzQsMTk3NSwxOTc1LDE5NzUsMTk3NSwxOTc2LDE5NzYsMTk3NiwxOTc2LDE5NzcsMTk3NywxOTc3LDE5NzcsMTk3OCwxOTc4LDE5NzgsMTk3OCwxOTc5LDE5NzksMTk3OSwxOTc5LDE5ODAsMTk4MCwxOTgwLDE5ODAsMTk4MSwxOTgxLDE5ODEsMTk4MSwxOTgyLDE5ODIsMTk4MiwxOTgyLDE5ODMsMTk4MywxOTgzLDE5ODMsMTk4NCwxOTg0LDE5ODQsMTk4NCwxOTg1LDE5ODUsMTk4NSwxOTg1LDE5ODYsMTk4NiwxOTg2LDE5ODYsMTk4NywxOTg3LDE5ODcsMTk4NywxOTg4LDE5ODgsMTk4OCwxOTg4LDE5ODksMTk4OSwxOTg5LDE5ODksMTk5MCwxOTkwLDE5OTAsMTk5MCwxOTkxLDE5OTEsMTk5MSwxOTkxLDE5OTIsMTk5MiwxOTkyLDE5OTIsMTk5MywxOTkzLDE5OTMsMTk5MywxOTk0LDE5OTQsMTk5NCwxOTk0LDE5OTUsMTk5NSwxOTk1LDE5OTUsMTk5NiwxOTk2LDE5OTYsMTk5NiwxOTk3LDE5OTcsMTk5NywxOTk3LDE5OTgsMTk5OCwxOTk4LDE5OTgsMTk5OSwxOTk5LDE5OTksMTk5OSwyMDAwLDIwMDAsMjAwMCwyMDAwLDIwMDEsMjAwMSwyMDAxLDIwMDEsMjAwMiwyMDAyLDIwMDIsMjAwMiwyMDAzLDIwMDMsMjAwMywyMDAzLDIwMDQsMjAwNCwyMDA0LDIwMDQsMjAwNSwyMDA1LDIwMDUsMjAwNSwyMDA2LDIwMDYsMjAwNiwyMDA2LDIwMDcsMjAwNywyMDA3LDIwMDcsMjAwOCwyMDA4LDIwMDgsMjAwOCwyMDA5LDIwMDksMjAwOSwyMDA5LDIwMTAsMjAxMCwyMDEwLDIwMTAsMjAxMSwyMDExLDIwMTEsMjAxMSwyMDEyLDIwMTIsMjAxMiwyMDEyLDIwMTMsMjAxMywyMDEzLDIwMTMsMjAxNCwyMDE0LDIwMTQsMjAxNCwyMDE1LDIwMTUsMjAxNSwyMDE1LDIwMTYsMjAxNiwyMDE2LDIwMTYsMjAxNywyMDE3LDIwMTcsMjAxNywyMDE4LDIwMTgpLFxuICAgICAgICAgICAgICAgICBmcXRyPWMoMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyLDMsNCwxLDIsMyw0LDEsMiwzLDQsMSwyKSxcbiAgICAgICAgICAgICAgICAgcmV2dHE9Yyg3LjI2Nyw5LjkzOSwxMS40NywxNS42MSwxMi43NTgsMTguMjA3LDE5LjIwNSwyNy44NDQsMjEuMDExLDI5LjkxOCwzMi42ODMsNDEuMjc3LDMyLjczNywzNy42NTUsNDIuMDk4LDU1LjA3MSw0Ni4zNTUsNTYuMTc3LDU5LjQ2Miw3NC4yMTUsNjQuNjUsNzguODUyLDg2LjEwOCwxMTAuNzIxLDkzLjk3MSwxMTUuMTg2LDEyMS44NzUsMTQ3Ljc3NSwxMjIuNjU1LDE1Mi4zODEsMTcyLjYzNywyMzAuNzgzLDE2Ny4yNzQsMjA0LjA4NSwyMjEuOTgyLDMwNi45NTcsMjM4LjU2OSwyOTEuNjg1LDMwMC43NDcsNDE3LjE3NSwzMTQuNTc1LDM3Mi41OTgsNDAxLjg0Miw1NTQuMTg0LDQzNC4yNTEsNTMzLjY1Myw2MzQuNDc2LDg0Mi42MTcsNjY1LjIyOSw3ODguNzA1LDgyNy43NjcsMTA5NC41NTEsODU0LjczNCwxMDk4Ljk0MiwxMTY2LjU5NCwxNTQ2LjYzOSwxMjM0Ljc2OCwxNTA4LjUzMywxNTgzLjU3MywyMDczLjk4MywxNjU1LjY2MSwxOTM0LjM1OCwyMDg3LjUzMywyNzczLjkzNiwyMzQ0LjI1MywyNzcwLjI3NiwyOTQ5LjAwOCwzODQ1LjU1NywzMjExLjc2MSwzNzMyLjIzNyw0MDE1LjAxNyw1MDAwLjIzOCw0MjgyLjAxMiw0ODQ1LjEzMyw0OTkxLjExMyw2NTMwLjczNCw1MzczLjI1OCw2MDQ2LjQxLDYyODMuNDk2LDgxMDcuNDg0LDY3NjguMTkxLDc1NDMuNTA4LDc5MzAuOTQ5LDEwMzU4Ljg5OCw5MjgwLjU2NiwxMDMzOS44OTgsMTA2MjcuNSwxMzYzOC43OTcsMTE2NDkuMzk4LDEzMDI4LjM5OCwxMzY4My43OTcsMTcxMjIsMTM5MjAuMzk4LDE2MjM2LjUsMTY4MjYuODk4LDIwMzYwLjUsMTc2ODYuMDk4LDE5OTQyLjI5NywyMDQxNy42OTksMjQ0NDcuNzk3LDIwNDQwLDIyNzIzLDIyOTE0LDI3NTUwLDIyNzcyLDI1NTg3LDI1NjQ0LDMwODU2LDI1NDA5LDI4Mzg2LDI4Nzc3LDM1Mzg2LDI5ODE5LDMzNTIxLDMzNTA5LDQwNzg1LDM1MTI5LDM4OTEzLDQwODk5LDUxODY4LDQzNDQ3LDQ2NTg4LDQ2MTgxLDU3MDc5LDQ4NTY1LDUzMTg3LDUzMTg1LDY0NzM1LDUyMTI2LDU2NzgxLDU1NzY1LDY2OTA1LDU3MjI0LDYzMjMxLDYzMDM2LDc1MTkwLDY1NDM0LDcwNTE2LDY5MjgyLDgyOTAwLDcwNzU1LDc2Njk3LDc1Mzk3LDg5MjUyLDc5Njc1LDg1NDMwLDg0NDY3LDk5MDc4LDg2Mzc4LDkyODI3LDkxNzk0LDEwNzI4OSw5NDk0MCwxMDIzNDIsOTgzNDUsMTA4NzQ3LDk0MjE0LDEwMDg3Niw5OTM3MywxMTM2MjIsOTk4MTEsMTAzNzI2LDEwMTk1MiwxMTYzNjAsMTA0MTg5LDEwOTM2NiwxMTAyMjYsMTIzMTY5LDExMzAxMCwxMTQwNjUsMTEzODAwLDEyNzc3NiwxMTQwNzAsMTE2ODMwLDExNTY4OCwxMjk3MDYsMTE0OTYwLDEyMDEyNSwxMTkwMDEsMTMxNTY1LDExNDgyNiwxMjAyMjksMTE3NDA4LDEyOTY2NywxMTU5MDQsMTIwMzE5LDExODE3OSwxMzA5MzYsMTE3NTQyLDEyMjk2OCwxMjMxNzksMTM2MjY3LDEyMjY5MCwxMjgwMjgpLFxuICAgICAgICAgICAgICAgICByZWNlc3Npb249YygxLCAxLCAxLCAxLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAxLCAxLCAxLCAxLCAxLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAxLCAxLCAwLCAwLCAwLCAwLCAxLCAxLCAxLCAxLCAxLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAxLCAxLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAxLCAxLCAxLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAxLCAxLCAxLCAxLCAxLCAxLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwKSlcbmRmIDwtIGRmICU+JSBtdXRhdGUocmV2dHFfeW95ID0gKHJldnRxIC0gbGFnKHJldnRxLCA0KSkgLyBsYWcocmV2dHEsIDQpKVxuZGYgPC0gZGYgJT4lIG11dGF0ZShyZWNlc3Npb25fbGFnID0gbGFnKHJlY2Vzc2lvbikpXG5kZiA8LSBkZiAlPiUgbXV0YXRlKHJldnRxX2Rvd24gPSBpZmVsc2UocmV2dHFfeW95PDAsMSwwKSkiLCJzYW1wbGUiOiIjIFRoZSBkYXRhZnJhbWUgYmVsb3cgY29udGFpbnMgaGlzdG9yaWNhbCByZXZlbnVlIGZvciBXYWxtYXJ0IHNpbmNlIDE5NzBcbiMgQWxsIGRhdGEgaXMgcXVhcnRlcmx5XG4jIFJlY2Vzc2lvbiBkYXRhIGNvbWVzIGZyb20gdGhlIFVTIEZlZGVyYWwgUmVzZXJ2ZVxuIyBUaGUgZmlsZSB1c2VzIENvbXB1c3RhdCdzIG5vbWVuY2xhdHVyZTpcbiMgICAgIC0gcmV2ZW51ZSBpcyByZXZ0cVxuIyAgICAgLSB0aWNrZXIgaXMgdGljXG4jICAgICAtIGZ5ZWFycSBpcyB5ZWFyXG4jICAgICAtIGZxdHIgaXMgcXVhcnRlclxuIyAgICAgLSByZWNlc3Npb24gaXMgYW4gaW5kaWNhdG9yIGZvciBVUyByZWNlc3Npb25zXG4jIEl0IGFsc28gY29udGFpbnMgcmV2dHFfeW95LCByZXZ0cV9kb3duLCBhbmQgcmVjZXNzaW9uX2xhZyBmcm9tIHRoZSBmaXJzdCBleGVyY2lzZVxuZGZbMTo1LF1cblxuIyBDYWxjdWxhdGUgYSBsb2dpc3RpYyBtb2RlbCBvZiBhIGRlY3JlYXNlIGluIHllYXIgb3ZlciB5ZWFyIHF1YXJ0ZXJseSByZXZlbnVlXG4jIG9uIHRoZSBwcmVzZW5jZSBvZiBhIHJlY2Vzc2lvbiBsYXN0IHF1YXJ0ZXIsIHN0b3JlIGl0IGluIG1vZDFcbm1vZDEgPC0gZ2xtKCB+ICwgZGF0YT0sIGZhbWlseT0pXG5cbiMgUHJpbnQgb3V0IGEgc3VtbWFyeSBvZiB0aGUgbW9kZWwgdXNpbmcgc3VtbWFyeVxuXG5cbiNFTkQiLCJzb2x1dGlvbiI6IiMgVGhlIGRhdGFmcmFtZSBiZWxvdyBjb250YWlucyBoaXN0b3JpY2FsIHJldmVudWUgZm9yIFdhbG1hcnQgc2luY2UgMTk3MFxuIyBBbGwgZGF0YSBpcyBxdWFydGVybHlcbiMgUmVjZXNzaW9uIGRhdGEgY29tZXMgZnJvbSB0aGUgVVMgRmVkZXJhbCBSZXNlcnZlXG4jIFRoZSBmaWxlIHVzZXMgQ29tcHVzdGF0J3Mgbm9tZW5jbGF0dXJlOlxuIyAgICAgLSByZXZlbnVlIGlzIHJldnRxXG4jICAgICAtIHRpY2tlciBpcyB0aWNcbiMgICAgIC0gZnllYXJxIGlzIHllYXJcbiMgICAgIC0gZnF0ciBpcyBxdWFydGVyXG4jICAgICAtIHJlY2Vzc2lvbiBpcyBhbiBpbmRpY2F0b3IgZm9yIFVTIHJlY2Vzc2lvbnNcbiMgSXQgYWxzbyBjb250YWlucyByZXZ0cV95b3ksIHJldnRxX2Rvd24sIGFuZCByZWNlc3Npb25fbGFnIGZyb20gdGhlIGZpcnN0IGV4ZXJjaXNlXG5kZlsxOjUsXVxuXG4jIENhbGN1bGF0ZSBhIGxvZ2lzdGljIG1vZGVsIG9mIGEgZGVjcmVhc2UgaW4geWVhciBvdmVyIHllYXIgcXVhcnRlcmx5IHJldmVudWVcbiMgb24gdGhlIHByZXNlbmNlIG9mIGEgcmVjZXNzaW9uIGxhc3QgcXVhcnRlciwgc3RvcmUgaXQgaW4gbW9kMVxubW9kMSA8LSBnbG0ocmV2dHFfZG93biB+IHJlY2Vzc2lvbl9sYWcsIGRhdGE9ZGYsIGZhbWlseT1iaW5vbWlhbClcblxuIyBQcmludCBvdXQgYSBzdW1tYXJ5IG9mIHRoZSBtb2RlbCB1c2luZyBzdW1tYXJ5XG5zdW1tYXJ5KG1vZDEpXG5cbiNFTkQiLCJzY3QiOiIjIFRlbXBsYXRlIGJhc2VkIG9uIGh0dHBzOi8vd3d3LnJkb2N1bWVudGF0aW9uLm9yZy9wYWNrYWdlcy90ZXN0d2hhdC92ZXJzaW9ucy80LjEuMVxuIyBDaGVjayBpZiBzb21ldGhpbmcgaXMgZXhwbGljaXRseSB0eXBlZFxuXG50ZXN0X3N0dWRlbnRfdHlwZWQoJ3JldnRxX2Rvd24gfiByZWNlc3Npb25fbGFnJyxub3RfdHlwZWRfbXNnPSdNYWtlIHN1cmUgdG8gcmVncmVzcyBgcmV2dHFfZG93bmAgb24gYHJlY2Vzc2lvbl9sYWdgLicpXG50ZXN0X3N0dWRlbnRfdHlwZWQoJ2RhdGE9ZGYnLG5vdF90eXBlZF9tc2c9J01ha2Ugc3VyZSB0byB1c2luZyB0aGUgYGRmYCBkYXRhIGZyYW1lIGZvciB0aGUgdGVzdC4nKVxudGVzdF9zdHVkZW50X3R5cGVkKCdmYW1pbHk9Ymlub21pYWwnLG5vdF90eXBlZF9tc2c9J01ha2Ugc3VyZSB0byBzZXQgZmFtaWx5IHRvIGJlIGBiaW5vbWlhbGAsIHNvIGFzIHRvIHJ1biBhIGxvZ2lzdGljIHJlZ3Jlc3Npb24uJylcbnRlc3RfZXhwcmVzc2lvbl9vdXRwdXQoXCJtb2QxXCIsIGluY29ycmVjdF9tc2c9XCJZb3VyIHJlZ3Jlc3Npb24gbW9kZWwgaXNuJ3QgcXVpdGUgcmlnaHQuXCIpXG50ZXN0X3N0dWRlbnRfdHlwZWQoJ3N1bW1hcnkobW9kMSknLG5vdF90eXBlZF9tc2c9J1lvdXIgbW9kZWwgaXMgY29ycmVjdCwgYnV0IGRpZCB5b3UgcmVtZW1iZXIgdG8gc3VtbWFyaXplIGl0IHdpdGggYHN1bW1hcnkoKWA/JylcblxuIyB0ZXN0X3N0dWRlbnRfdHlwZWQoJ3ggPC0gMicsIG5vdF90eXBlZF9tc2c9JycpXG5cbiMgQ2hlY2sgaWYgZnVuY3Rpb24gd2FzIHVzZWQgaW4gaW5wdXQgY29kZVxuIyB0ZXN0X2Z1bmN0aW9uKCdjJyxpbmNvcnJlY3RfbXNnPScnKSAgXG5cbiMgUmVxdWlyZXMgYW4gb2JqZWN0IGB4YCB0byBoYXZlIHRoZSBzYW1lIHZhbHVlIGFzIHRoZSBzb2x1dGlvblxuIyB0ZXN0X29iamVjdChcInhcIixpbmNvcnJlY3RfbXNnID0gXCJcIix1bmRlZmluZWRfbXNnID0gXCJcIikgIFxuXG4jIFJlcXVpcmVzIGFuIG9uamVjdCB3aXRoIHRoZSBzYW1lIHZhbHVlIG9mIGB4YCBpbiB0aGUgc29sdXRpb25cbiMgdGVzdF9hbl9vYmplY3QoXCJ4XCIsdW5kZWZpbmVkX21zZz1cIlwiKVxuXG4jIENoZWNrcyBpZiBvdXRwdXQgb2Ygc3R1ZGVudCdzIGNvZGUgY29udGFpbnMgZ2l2ZW4gZXZhbHVhdGVkIGV4cHJlc3Npb25cbiMgdGVzdF9vdXRwdXRfY29udGFpbnMoXCJ4XCIsaW5jb3JyZWN0X21zZyA9IFwiXCIpXG5cbiMgQ2hlY2sgaWYgYSB2ZWN0b3Igb2YgcHJlZGVmaW5lZCBvYmplY3RzIGFyZSB1bmNoYW5nZWRcbiMgdGVzdF9wcmVkZWZpbmVkX29iamVjdHMoYygneCcsJ3knKSxpbmNvcnJlY3RfbXNnPVwiRG9uJ3Qgb252ZXJ3cml0ZSB0aGUgcHJlZGVmaW5lZCB2YXJpYWJsZXNcIilcblxuIyBDaGVja3MgZm9yIGEgcmVnZXggcGF0dGVybiBpbiB0cmhlIG91dHB1dFxuIyB0ZXN0X291dHB1dF9yZWdleChwYXR0ZXJuLGZpeGVkPUYsIHRpbWVzPTEsIGluY29ycmVjdF9tc2c9JycpXG5cbiMgQ2FuIGNoZWNrIGFuIGFyYml0cmFyeSBleHByZXNzaW9uIGFjcm9zcyBib3RoIHNvbHV0aW9uIGFuZCBzdHVkZW50IGNvZGVcbiN0ZXN0X2V4cHJlc3Npb25fb3V0cHV0KFwidHlwZW9mKGNvbXBhbnlfbmFtZSlcIiwgaW5jb3JyZWN0X21zZz1cIkRpZCB5b3Ugc3RvcmUgdGV4dHVhbCBkYXRhIGluIGBjb21wYW55X25hbWVgP1wiKVxuXG50ZXN0X2Vycm9yKClcbnN1Y2Nlc3NfbXNnKFwiQXdlc29tZSFcIilcblxuIyBPdGhlciBmdW5jdGlvbnMgdG8gbm90ZTpcbiMgICAgIC0gdGVzdF9vcihhLGIpIC0tIGNoZWNrcyBpZiBlaXRoZXIgdGVzdCBhIG9yIHRlc3QgYiBwYXNzXG4jICAgICAtIHRlc3RfZ2dwbG90KCkgLS0gY2FuIGNoZWNrIGlmIHBsb3RzIGFyZSBjb3JyZWN0XG4jICAgICAtIHRlc3RfZnVuY3Rpb24oKSAtLSBjYW4gYWxzbyBjaGVjayBpbmNsdWRlZCBwYXJhbWV0ZXJzXG4jICAgICAtIHRlc3RfbG9vcCgpIC0tIGNoZWNraW5nIGZvciBhbmQgd2hpbGUgbG9vcHNcbiMgICAgIC0gdGVzdF9saWJyYXJ5X2Z1bmN0aW9uKCdwYWNrYWdlJywgbm90X2NhbGxlZF9tc2c9JycsaW5jb3JyZWN0X21zZz0nJylcbiMgICAgIC0gdGVzdF9pZl9lbHNlKCkgLS0gY2hlY2tpbmcgaWYgc3RhdGVtZW50c1xuIyAgICAgLSB0ZXN0X2V4cHJlc3Npb25fZXJyb3IoKSAtLSBjYW4gY2hlY2sgaWYgZnVuY3Rpb25zIGFyZSBwcm9wZXJseSBkZWZpbmVkXG4jICAgICAtIHRlc3Rfb3BlcmF0b3IoJ29wZXJhdG9yJywpLCBub3RfY2FsbGVkX21zZz0nJyxpbmNvcnJlY3RfbXNnPScnKVxuIyAgICAgLSB0ZXN0X2Z1bmN0aW9uX2RlZmluaXRpb24oKSAtLSByaWdvcm91c2x5IGNoZWNrIGRlZmluZWQgZnVuY3Rpb25cbiMgICAgIC0gdGVzdF9kYXRhX2ZyYW1lKCkgLS0gY2hlY2sgaWYgZGF0YWZyYW1lIFtjb2x1bW5zXSBhcmUgZXF1aXZhbGVudFxuIyAgICAgLSB0ZXN0X2Z1bmN0aW9uX3Jlc3VsdCwgdGVzdF9leHByZXNzaW9uX3Jlc3VsdCJ9

Visualization with maps

Exercise 3: Using Plotly for maps

NOTE: Unforetunately, though DataCamp light has plotly installed, it doesn’t work at present (plots don’t open). As such, we can’t run this here and see any output. Instead, you can download this file, load the file into R with read.csv(), and then run the code given below. Alternatively, you can download this Rmd file containing this exercise’s code. The output of the code is available here.

You will need plotly installed to run the code. If you do not have it installed, run install.packages("plotly").

eyJsYW5ndWFnZSI6InIiLCJzYW1wbGUiOiIjIFRoZSBmaWxlIG1nLmNzdiBjb250YWlucyB0aGUgZm9sbG93aW5nIGluZm9ybWF0aW9uIGZvciBUeXBob29uIE1hbmdraHV0OlxuIyAgICAgLSBsYXQ6IGxhdGl0dWRlXG4jICAgICAtIGxvbjogbG9uZ2l0dWRlIChmaXhlZClcbiMgICAgIC0gZHQ6IERhdGUgdGltZSBpbiBmb3JtYXQgXCIlWSViJWQgJUglTSVTXCJcbm1nIDwtIHJlYWQuY3N2KFwibWcuY3N2XCIpXG5cbiMgU2FtcGxlIG9mIHRoZSBkYXRhOlxubWdbMTo1LF1cblxuIyBMb2FkIHBsb3RseSBmb3IgcGxvdHRpbmdcbmxpYnJhcnkocGxvdGx5KVxuXG4jIFNldHRpbmdzIGZvciB0aGUgcGxvdCAtLSBjb21wbGV0ZSwgYnV0IHNvbWUgb3RoZXIgb3B0aW9ucyBhcmUgbGlzdGVkIGluIHRoZSBjb21tZW50c1xuZ2VvIDwtIGxpc3QoXG4gIHNob3dsYW5kID0gVFJVRSxcbiAgc2hvd2xha2VzID0gVFJVRSxcbiAgc2hvd2NvdW50cmllcyA9IFRSVUUsXG4gIHNob3dvY2VhbiA9IFRSVUUsXG4gIGNvdW50cnl3aWR0aCA9IDAuNSxcbiAgbGFuZGNvbG9yID0gdG9SR0IoXCJncmV5OTBcIiksICAjIFJ1biBjb2xvcnMoKSBhZnRlciBsb2FkaW5nIHBsb3RseSB0byBzZWUgYXZhaWxhYmxlIGNvbG9yc1xuICBsYWtlY29sb3IgPSB0b1JHQihcImFsaWNlYmx1ZVwiKSxcbiAgb2NlYW5jb2xvciA9IHRvUkdCKFwiYWxpY2VibHVlXCIpLFxuICBwcm9qZWN0aW9uID0gbGlzdChcbiAgICB0eXBlID0gJ29ydGhvZ3JhcGhpYycsICAjIENvbnNpZGVyIHRyeWluZyAnZXF1aXJlY3Rhbmd1bGFyJyBvciAnc3RlcmlvZ3JhcGhpYydcbiAgICByb3RhdGlvbiA9IGxpc3QoXG4gICAgICBsb24gPSAxMDAsICAjIFRoaXMgc2V0cyB0aGUgY2VudGVyIGxvbmdpdHVkZVxuICAgICAgbGF0ID0gMSwgICAgIyBUaGlzIHNldHMgdGhlIGNlbnRlciBsYXRpdHVkZVxuICAgICAgcm9sbCA9IDBcbiAgICApXG4gICksXG4gIGxvbmF4aXMgPSBsaXN0KFxuICAgIHNob3dncmlkID0gVFJVRSxcbiAgICBncmlkY29sb3IgPSB0b1JHQihcImdyYXk0MFwiKSxcbiAgICBncmlkd2lkdGggPSAwLjVcbiAgKSxcbiAgbGF0YXhpcyA9IGxpc3QoXG4gICAgc2hvd2dyaWQgPSBUUlVFLFxuICAgIGdyaWRjb2xvciA9IHRvUkdCKFwiZ3JheTQwXCIpLFxuICAgIGdyaWR3aWR0aCA9IDAuNVxuICApXG4pXG5cblxuIyBQbG90XG5wbG90X2dlbygpICU+JVxuICBhZGRfbWFya2VycyhkYXRhPW1nLFxuICAgICAgICAgICAgICB4ID0gfmxvbixcbiAgICAgICAgICAgICAgeSA9IH5sYXQsXG4gICAgICAgICAgICAgIHRleHQ9fnBhc3RlKFwiVGltZTogXCIsIGR0KSkgJT4lXG4gIGxheW91dChzaG93bGVnZW5kID0gVFJVRSwgZ2VvID0gZ2VvLFxuICAgIHRpdGxlID0gJ1R5cGhvb24gTWFuZ2todXQgKFNlcHQgMjAxOCksIGZ1bGwgcGF0aCcpIn0=

Exercise 4: Plotting with leaflet…

NOTE: Unforetunately, DataCamp light doesn’t have leaflet installed. As such, we can’t run this here. If you would like to see a nice graph of the typhoon path in leaflet, you can download this file (same csv as for Exercise 3) and this Rmd file containing this exercise’s code. The output of the code is available here.

You will need leaflet installed to run the code. If you do not have it installed, run install.packages("leaflet")

The final map also requires leaflet.extras installed. If you do not have it installed, run install.packages("leaflet.extras")