Tutorial 4, Probability, Statistics and Modelling 2, BSc Security and Crime Science, UCL
Outcomes of this tutorial
In this tutorial, you will:
- apply NHST and Bayesian hypothesis testing
Task 1-3 are to be done in R. Task 4-6 in JASP.
Task 1: Evaluation of driving campaign
You are tasked to evaluate a safe driving campaign. On 150 road crossing, a local police force implemented a safe-driving cmapign through ads that promote better driving practices. As a control, they chose 150 comparable crossings where no such campaign was run. At each locaiton, police officers evaluated the driving behavior by assigning a score from 0-100. Higher scores, indicate better driving practices.
You are now asked to decide whether or not the funding for the campaign should be extended.
Provide an answer in the NHST framework and in the Bayesian framework.
The data are located in ./data/tutorial4/driving_data.csv
.
#your code here
What do you advice? Should the funding be continued?
#type your answer here
Task 2: Revisiting fraud data
Load the dataframe damage_data
from ./data/tutorial4/analyst_data.RData
. In that dataframe, you have 1000 observations each for junior and senior analysts and the damage they caused for their employers in fine that needed to be paid for compliance issues. You are in charge of leading a new, better statistical insights team and your employer asks you whether junior or senior analysts are more problematic.
Provide an answer from both the frequentist and the Bayesian perspective.
#your code here
Task 3: Original data
You will use data from a paper published in Applied Cognitive Psychology. The paper (OA at https://onlinelibrary.wiley.com/doi/full/10.1002/acp.3439) replicated an experiment that found that people who are deceptive about their future plans struggle to be as detailed as truth-tellers. Here, the outcome variable of interest was the number of specific time references in the participants’ answers to interview questions (e.g. “9am”, “half past seven”) as counted by instructed human coders.
Load the data from ./data/tutorial4/replication_data.csv
and test:
- whether the findings replicated using the NHST framework
- what Bayesian statistics concludes.
#your code here
What is your conclusion?
#write your answer here
Task 4: “Evaluation of driving campaign” in JASP
You are tasked to evaluate a safe driving campaign. On 150 road crossing, a local police force implemented a safe-driving cmapign through ads that promote better driving practices. As a control, they chose 150 comparable crossings where no such campaign was run. At each locaiton, police officers evaluated the driving behavior by assigning a score from 0-100. Higher scores, indicate better driving practices.
You are now asked to decide whether or not the funding for the campaign should be extended.
Provide an answer in the NHST framework and in the Bayesian framework.
Conduct these analyses in JASP
Task 5: “Revisiting fraud data” in JASP
Load the dataframe damage_data
from ./data/tutorial4/analyst_data.RData
. In that dataframe, you have 1000 observations each for junior and senior analysts and the damage they caused for their employers in fine that needed to be paid for compliance issues. You are in charge of leading a new, better statistical insights team and your employer asks you whether junior or senior analysts are more problematic.
Provide an answer from both the frequentist and the Bayesian perspective.
Conduct these analyses in JASP
Task 6: “Original data” in JASP
You will use data from a paper published in Applied Cognitive Psychology. The paper (OA at https://onlinelibrary.wiley.com/doi/full/10.1002/acp.3439) replicated an experiment that found that people who are deceptive about their future plans struggle to be as detailed as truth-tellers. Here, the outcome variable of interest was the number of specific time references in the participants’ answers to interview questions (e.g. “9am”, “half past seven”) as counted by instructed human coders.
Load the data from ./data/tutorial4/replication_data.csv
and test:
- whether the findings replicated using the NHST framework
- what Bayesian statistics concludes.
Conduct these analyses in JASP
END
LS0tCnRpdGxlOiBCYXllc2lhbiBoeXBvdGhlc2lzIHRlc3RpbmcKYXV0aG9yOiBCIEtsZWluYmVyZwpkYXRlOiA1IE1hcmNoIDIwMTkKc3VidGl0bGU6IERlcHQgb2YgU2VjdXJpdHkgYW5kIENyaW1lIFNjaWVuY2UsIFVDTApvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgotLS0KClR1dG9yaWFsIDQsIFByb2JhYmlsaXR5LCBTdGF0aXN0aWNzIGFuZCBNb2RlbGxpbmcgMiwgQlNjIFNlY3VyaXR5IGFuZCBDcmltZSBTY2llbmNlLCBVQ0wKCi0tLQoKIyMgT3V0Y29tZXMgb2YgdGhpcyB0dXRvcmlhbAoKSW4gdGhpcyB0dXRvcmlhbCwgeW91IHdpbGw6CgotIGFwcGx5IE5IU1QgYW5kIEJheWVzaWFuIGh5cG90aGVzaXMgdGVzdGluZwogICAgLSB1c2luZyBSCiAgICAtIHVzaW5nIEpBU1AKCi0tLQoKKipUYXNrIDEtMyBhcmUgdG8gYmUgZG9uZSBpbiBSLiBUYXNrIDQtNiBpbiBKQVNQLioqCgojIyBUYXNrIDE6IEV2YWx1YXRpb24gb2YgZHJpdmluZyBjYW1wYWlnbgoKWW91IGFyZSB0YXNrZWQgdG8gZXZhbHVhdGUgYSBzYWZlIGRyaXZpbmcgY2FtcGFpZ24uIE9uIDE1MCByb2FkIGNyb3NzaW5nLCBhIGxvY2FsIHBvbGljZSBmb3JjZSBpbXBsZW1lbnRlZCBhIHNhZmUtZHJpdmluZyBjbWFwaWduIHRocm91Z2ggYWRzIHRoYXQgcHJvbW90ZSBiZXR0ZXIgZHJpdmluZyBwcmFjdGljZXMuIEFzIGEgY29udHJvbCwgdGhleSBjaG9zZSAxNTAgY29tcGFyYWJsZSBjcm9zc2luZ3Mgd2hlcmUgbm8gc3VjaCBjYW1wYWlnbiB3YXMgcnVuLiBBdCBlYWNoIGxvY2FpdG9uLCBwb2xpY2Ugb2ZmaWNlcnMgZXZhbHVhdGVkIHRoZSBkcml2aW5nIGJlaGF2aW9yIGJ5IGFzc2lnbmluZyBhIHNjb3JlIGZyb20gMC0xMDAuIEhpZ2hlciBzY29yZXMsIGluZGljYXRlIGJldHRlciBkcml2aW5nIHByYWN0aWNlcy4KCllvdSBhcmUgbm93IGFza2VkIHRvIGRlY2lkZSB3aGV0aGVyIG9yIG5vdCB0aGUgZnVuZGluZyBmb3IgdGhlIGNhbXBhaWduIHNob3VsZCBiZSBleHRlbmRlZC4KClByb3ZpZGUgYW4gYW5zd2VyIGluIHRoZSBOSFNUIGZyYW1ld29yayBhbmQgaW4gdGhlIEJheWVzaWFuIGZyYW1ld29yay4KClRoZSBkYXRhIGFyZSBsb2NhdGVkIGluIGAuL2RhdGEvdHV0b3JpYWw0L2RyaXZpbmdfZGF0YS5jc3ZgLgoKYGBge3J9CiN5b3VyIGNvZGUgaGVyZQoKCmBgYAoKV2hhdCBkbyB5b3UgYWR2aWNlPyBTaG91bGQgdGhlIGZ1bmRpbmcgYmUgY29udGludWVkPwoKYGBge3J9CiN0eXBlIHlvdXIgYW5zd2VyIGhlcmUKCgpgYGAKCgojIyBUYXNrIDI6IFJldmlzaXRpbmcgZnJhdWQgZGF0YQoKTG9hZCB0aGUgZGF0YWZyYW1lIGBkYW1hZ2VfZGF0YWAgZnJvbSBgLi9kYXRhL3R1dG9yaWFsNC9hbmFseXN0X2RhdGEuUkRhdGFgLiBJbiB0aGF0IGRhdGFmcmFtZSwgeW91IGhhdmUgMTAwMCBvYnNlcnZhdGlvbnMgZWFjaCBmb3IganVuaW9yIGFuZCBzZW5pb3IgYW5hbHlzdHMgYW5kIHRoZSBkYW1hZ2UgdGhleSBjYXVzZWQgZm9yIHRoZWlyIGVtcGxveWVycyBpbiBmaW5lIHRoYXQgbmVlZGVkIHRvIGJlIHBhaWQgZm9yIGNvbXBsaWFuY2UgaXNzdWVzLiBZb3UgYXJlIGluIGNoYXJnZSBvZiBsZWFkaW5nIGEgbmV3LCBiZXR0ZXIgc3RhdGlzdGljYWwgaW5zaWdodHMgdGVhbSBhbmQgeW91ciBlbXBsb3llciBhc2tzIHlvdSB3aGV0aGVyIGp1bmlvciBvciBzZW5pb3IgYW5hbHlzdHMgYXJlIG1vcmUgcHJvYmxlbWF0aWMuCgpQcm92aWRlIGFuIGFuc3dlciBmcm9tIGJvdGggdGhlIGZyZXF1ZW50aXN0IGFuZCB0aGUgQmF5ZXNpYW4gcGVyc3BlY3RpdmUuCgpgYGB7cn0KI3lvdXIgY29kZSBoZXJlCgoKYGBgCgojIyBUYXNrIDM6IE9yaWdpbmFsIGRhdGEKCllvdSB3aWxsIHVzZSBkYXRhIGZyb20gYSBwYXBlciBwdWJsaXNoZWQgaW4gQXBwbGllZCBDb2duaXRpdmUgUHN5Y2hvbG9neS4gVGhlIHBhcGVyIFsoT0EgYXQgaHR0cHM6Ly9vbmxpbmVsaWJyYXJ5LndpbGV5LmNvbS9kb2kvZnVsbC8xMC4xMDAyL2FjcC4zNDM5KV0oaHR0cHM6Ly9vbmxpbmVsaWJyYXJ5LndpbGV5LmNvbS9kb2kvZnVsbC8xMC4xMDAyL2FjcC4zNDM5KSByZXBsaWNhdGVkIGFuIGV4cGVyaW1lbnQgdGhhdCBmb3VuZCB0aGF0IHBlb3BsZSB3aG8gYXJlIGRlY2VwdGl2ZSBhYm91dCB0aGVpciBmdXR1cmUgcGxhbnMgc3RydWdnbGUgdG8gYmUgYXMgZGV0YWlsZWQgYXMgdHJ1dGgtdGVsbGVycy4gSGVyZSwgdGhlIG91dGNvbWUgdmFyaWFibGUgb2YgaW50ZXJlc3Qgd2FzIHRoZSBudW1iZXIgb2Ygc3BlY2lmaWMgdGltZSByZWZlcmVuY2VzIGluIHRoZSBwYXJ0aWNpcGFudHMnIGFuc3dlcnMgdG8gaW50ZXJ2aWV3IHF1ZXN0aW9ucyAoZS5nLiAiOWFtIiwgImhhbGYgcGFzdCBzZXZlbiIpIGFzIGNvdW50ZWQgYnkgaW5zdHJ1Y3RlZCBodW1hbiBjb2RlcnMuCgpMb2FkIHRoZSBkYXRhIGZyb20gYC4vZGF0YS90dXRvcmlhbDQvcmVwbGljYXRpb25fZGF0YS5jc3ZgIGFuZCB0ZXN0OgoKMS4gd2hldGhlciB0aGUgZmluZGluZ3MgcmVwbGljYXRlZCB1c2luZyB0aGUgTkhTVCBmcmFtZXdvcmsKMi4gd2hhdCBCYXllc2lhbiBzdGF0aXN0aWNzIGNvbmNsdWRlcy4KCmBgYHtyfQojeW91ciBjb2RlIGhlcmUKCmBgYAoKV2hhdCBpcyB5b3VyIGNvbmNsdXNpb24/CgpgYGB7cn0KI3dyaXRlIHlvdXIgYW5zd2VyIGhlcmUKCgpgYGAKCgotLS0KCiMjIFRhc2sgNDogIkV2YWx1YXRpb24gb2YgZHJpdmluZyBjYW1wYWlnbiIgaW4gSkFTUAoKWW91IGFyZSB0YXNrZWQgdG8gZXZhbHVhdGUgYSBzYWZlIGRyaXZpbmcgY2FtcGFpZ24uIE9uIDE1MCByb2FkIGNyb3NzaW5nLCBhIGxvY2FsIHBvbGljZSBmb3JjZSBpbXBsZW1lbnRlZCBhIHNhZmUtZHJpdmluZyBjbWFwaWduIHRocm91Z2ggYWRzIHRoYXQgcHJvbW90ZSBiZXR0ZXIgZHJpdmluZyBwcmFjdGljZXMuIEFzIGEgY29udHJvbCwgdGhleSBjaG9zZSAxNTAgY29tcGFyYWJsZSBjcm9zc2luZ3Mgd2hlcmUgbm8gc3VjaCBjYW1wYWlnbiB3YXMgcnVuLiBBdCBlYWNoIGxvY2FpdG9uLCBwb2xpY2Ugb2ZmaWNlcnMgZXZhbHVhdGVkIHRoZSBkcml2aW5nIGJlaGF2aW9yIGJ5IGFzc2lnbmluZyBhIHNjb3JlIGZyb20gMC0xMDAuIEhpZ2hlciBzY29yZXMsIGluZGljYXRlIGJldHRlciBkcml2aW5nIHByYWN0aWNlcy4KCllvdSBhcmUgbm93IGFza2VkIHRvIGRlY2lkZSB3aGV0aGVyIG9yIG5vdCB0aGUgZnVuZGluZyBmb3IgdGhlIGNhbXBhaWduIHNob3VsZCBiZSBleHRlbmRlZC4KClByb3ZpZGUgYW4gYW5zd2VyIGluIHRoZSBOSFNUIGZyYW1ld29yayBhbmQgaW4gdGhlIEJheWVzaWFuIGZyYW1ld29yay4KCioqQ29uZHVjdCB0aGVzZSBhbmFseXNlcyBpbiBKQVNQKioKCiMjIFRhc2sgNTogIlJldmlzaXRpbmcgZnJhdWQgZGF0YSIgaW4gSkFTUAoKTG9hZCB0aGUgZGF0YWZyYW1lIGBkYW1hZ2VfZGF0YWAgZnJvbSBgLi9kYXRhL3R1dG9yaWFsNC9hbmFseXN0X2RhdGEuUkRhdGFgLiBJbiB0aGF0IGRhdGFmcmFtZSwgeW91IGhhdmUgMTAwMCBvYnNlcnZhdGlvbnMgZWFjaCBmb3IganVuaW9yIGFuZCBzZW5pb3IgYW5hbHlzdHMgYW5kIHRoZSBkYW1hZ2UgdGhleSBjYXVzZWQgZm9yIHRoZWlyIGVtcGxveWVycyBpbiBmaW5lIHRoYXQgbmVlZGVkIHRvIGJlIHBhaWQgZm9yIGNvbXBsaWFuY2UgaXNzdWVzLiBZb3UgYXJlIGluIGNoYXJnZSBvZiBsZWFkaW5nIGEgbmV3LCBiZXR0ZXIgc3RhdGlzdGljYWwgaW5zaWdodHMgdGVhbSBhbmQgeW91ciBlbXBsb3llciBhc2tzIHlvdSB3aGV0aGVyIGp1bmlvciBvciBzZW5pb3IgYW5hbHlzdHMgYXJlIG1vcmUgcHJvYmxlbWF0aWMuCgpQcm92aWRlIGFuIGFuc3dlciBmcm9tIGJvdGggdGhlIGZyZXF1ZW50aXN0IGFuZCB0aGUgQmF5ZXNpYW4gcGVyc3BlY3RpdmUuCgoKKipDb25kdWN0IHRoZXNlIGFuYWx5c2VzIGluIEpBU1AqKgoKIyMgVGFzayA2OiAiT3JpZ2luYWwgZGF0YSIgaW4gSkFTUAoKWW91IHdpbGwgdXNlIGRhdGEgZnJvbSBhIHBhcGVyIHB1Ymxpc2hlZCBpbiBBcHBsaWVkIENvZ25pdGl2ZSBQc3ljaG9sb2d5LiBUaGUgcGFwZXIgWyhPQSBhdCBodHRwczovL29ubGluZWxpYnJhcnkud2lsZXkuY29tL2RvaS9mdWxsLzEwLjEwMDIvYWNwLjM0MzkpXShodHRwczovL29ubGluZWxpYnJhcnkud2lsZXkuY29tL2RvaS9mdWxsLzEwLjEwMDIvYWNwLjM0MzkpIHJlcGxpY2F0ZWQgYW4gZXhwZXJpbWVudCB0aGF0IGZvdW5kIHRoYXQgcGVvcGxlIHdobyBhcmUgZGVjZXB0aXZlIGFib3V0IHRoZWlyIGZ1dHVyZSBwbGFucyBzdHJ1Z2dsZSB0byBiZSBhcyBkZXRhaWxlZCBhcyB0cnV0aC10ZWxsZXJzLiBIZXJlLCB0aGUgb3V0Y29tZSB2YXJpYWJsZSBvZiBpbnRlcmVzdCB3YXMgdGhlIG51bWJlciBvZiBzcGVjaWZpYyB0aW1lIHJlZmVyZW5jZXMgaW4gdGhlIHBhcnRpY2lwYW50cycgYW5zd2VycyB0byBpbnRlcnZpZXcgcXVlc3Rpb25zIChlLmcuICI5YW0iLCAiaGFsZiBwYXN0IHNldmVuIikgYXMgY291bnRlZCBieSBpbnN0cnVjdGVkIGh1bWFuIGNvZGVycy4KCkxvYWQgdGhlIGRhdGEgZnJvbSBgLi9kYXRhL3R1dG9yaWFsNC9yZXBsaWNhdGlvbl9kYXRhLmNzdmAgYW5kIHRlc3Q6CgoxLiB3aGV0aGVyIHRoZSBmaW5kaW5ncyByZXBsaWNhdGVkIHVzaW5nIHRoZSBOSFNUIGZyYW1ld29yawoyLiB3aGF0IEJheWVzaWFuIHN0YXRpc3RpY3MgY29uY2x1ZGVzLgoKCioqQ29uZHVjdCB0aGVzZSBhbmFseXNlcyBpbiBKQVNQKioKCiMjIEVORAoKLS0tCgo=