Expected Preparations:

  [FND-STA]
Probability_distribution 		 
  The units listed above are part of this course and contain important preparatory material.  
Keywords: Probability and p-values; significance as a threshold of p-values; deriving probability distributions from simulation and interpreting in terms of significance.
Objectives:

  • Introduce the difference between p-values and event probability;

  • Show how we interpret p-values in terms of “significance”;

  • Illustrate this with an example;

  • Present a permutation example, a strategy that can be used as an alternative to the integration of probability density functions.

  • Discuss a common error that is made when establishing the significance of an observation in the biomedical sciences.

 Outcomes:

After working through this unit you should…

  • Be able to define a “p-value”;

  • Be able to set up a permutation test or a sampling simulation to estimate a probability density;

  • Be able to interpret the frequency of values in that probability density in terms of a p-value;

  • Be able to critically assess whether an observation should be considered “significant” in that context;

  • Be able to identify a common error that is made in the literature when two effects are compared.

Deliverables:

Time management: Before you begin, estimate how long it will take you to complete this unit. Then, record in your course journal: the number of hours you estimated, the number of hours you worked on the unit, and the amount of time that passed between start and completion of this unit.

Journal: Document your progress in your Course Journal. Some tasks may ask you to include specific items in your journal. Don’t overlook these.

Insights: If you find something particularly noteworthy about this unit, make a note in your insights! page.

Evaluation:

NA: This unit is not evaluated for course marks.

Contents

The probability of an event is the chance of it occurring, but how do we relate that to the question whether an observation is significant? In this context we talk about p-values and the meaning of a p-value is not the same as the probability of an observation. The p-value of an observation is the probability that - assuming a null hypothesis is true - an event as extreme or more extreme is observed. This unit contains R code to study this concept.

“Significance” concepts in practice

Here we discuss the idea of a p-value, in particular how to compute “empirical p-values”. These are very easy to compute and simulate in R. There is just one thing to be aware of: while we normally approximate a p-value from observed events r divided by the number of observations N as r / N, if we use this approach to evaluate significance, i.e. we are asking whetehr our r events are taken from the same distribution as the N observations, we need to apply a correction factor: (r + 1) / (N + 1)1.

 

Empirical p-value: (r + 1) / (N + 1)
for r events of interest in N observations.

 

Task…

  • Open RStudio and load the ABC-units R project. If you have loaded it before, choose File ▸ Recent projects ▸ ABC-Units. If you have not loaded it before, follow the instructions in the RPR-Introduction unit.
  • Choose Tools ▸ Version Control ▸ Pull Branches to fetch the most recent version of the project from its GitHub repository with all changes and bug fixes included. This ensures that your data and code remain up to date when we update, or fix bugs.
  • Type init() if requested.
  • Open the file FND-STA-Significance.R and follow the instructions.

     

    Note: take care that you understand all of the code in the script. Evaluation in this course is cumulative and you may be asked to explain any part of code.

 

Controversies

Task…

Examine the papers below that introduce difficulties with P-values and statistical significance. Rephrase the issues in your own words to make sure that you understand what the discussion is about.

Baker, Monya. (2016). “Statisticians issue warning over misuse of P values”. Nature 531(7593):151 .
[PMID: 26961635] [DOI: 10.1038/nature.2016.19503]

Nieuwenhuis, Sander, Birte U Forstmann, and Eric-Jan Wagenmakers. (2011). “Erroneous analyses of interactions in neuroscience: a problem of significance”. Nature Neuroscience 14(9):1105–7 .
[PMID: 21878926] [DOI: 10.1038/nn.2886]

In theory, a comparison of two experimental effects requires a statistical test on their difference. In practice, this comparison is often based on an incorrect procedure involving two separate tests in which researchers conclude that effects differ when one effect is significant (P < 0.05) but the other is not (P > 0.05). We reviewed 513 behavioral, systems and cognitive neuroscience articles in five top-ranking journals (Science, Nature, Nature Neuroscience, Neuron and The Journal of Neuroscience) and found that 78 used the correct procedure and 79 used the incorrect procedure. An additional analysis suggests that incorrect analyses of interactions are even more common in cellular and molecular neuroscience. We discuss scenarios in which the erroneous procedure is particularly beguiling.

Further Reading

J, Murdoch, D, Tsai, Yu-Ling, and Adcock, James. (2008). “P-Values are Random Variables”. The American Statistician Taylor & Francis, Alexandria, VA. 62(3):242–245 .
[DOI: 10.1198/000313008X332421]

P-values are taught in introductory statistics classes in a way that confuses many of the students, leading to common misconceptions about their meaning. In this article, we argue that p-values should be taught through simulation, emphasizing that p-values are random variables. By means of elementary examples we illustrate how to teach students valid interpretations of p-values and give them a deeper understanding of hypothesis testing.

Questions, comments

If in doubt, ask! If anything about this contents is not clear to you, do not proceed but ask for clarification. If you have ideas about how to make this material better, let’s hear them. We are aiming to compile a list of FAQs for all learning units, and your contributions will count towards your participation marks.

Improve this page! If you have questions or comments, please post them on the Quercus Discussion board with a subject line that includes the name of the unit.

References

Page ID: FND-STA-Significance

Author:
Boris Steipe ( <boris.steipe@utoronto.ca> )
Created:
2017-08-05
Last modified:
2022-09-14
Version:
1.1.1
Version History:
–  1.1 2020 Maintenance
–  1.1 Added definition of empirical p-value
–  1.0 First live
–  0.1 First stub
Tagged with:
–  Unit
–  Live
–  Links to R course project
–  Has further reading

 

[END]

  1. ↩︎