An introduction to the statistics behind the most popular genomic data science projects. This is the sixth course in the Genomic Big Data Science Specialization from Johns Hopkins University.

Loading...

From the course by Johns Hopkins University

Statistics for Genomic Data Science

110 ratings

An introduction to the statistics behind the most popular genomic data science projects. This is the sixth course in the Genomic Big Data Science Specialization from Johns Hopkins University.

From the lesson

Module 3

This week we will cover modeling non-continuous outcomes (like binary or count data), hypothesis testing, and multiple hypothesis testing.

- Jeff Leek, PhDAssociate Professor, Biostatistics

Bloomberg School of Public Health

For the most common kind of statistical inference that's performed in genomics and

Â indeed across a large range of genomic and biological sciences hypothesis testing.

Â Hypothesis testing gets a bad rap amongst statisticians because sometimes it can be

Â a little bit difficult to interpret.

Â And in particular, it's important that you get a null and

Â alternative hypotheses right when you're performing an analysis.

Â So what do I mean by a null and alternative hypothesis?

Â So, first, you're going to set up a hypothesis that you'd like to reject.

Â So in this case, you would say, maybe if you're doing a gene expression experiment,

Â and you're looking for genes associated with age,

Â the null hypothesis might be that there's no relationship at all between age and

Â the expression for a particular gene.

Â The alternative hypothesis would be that there is a relationship between age and

Â gene expression so that the relationship is non zero.

Â So when you set up these hypotheses, basically you're trying to set up

Â the thing that you're trying to reject if you collect enough data and

Â the thing that you're trying to see if that might be true.

Â So often you actually define the null hypothesis alone in that and

Â then you don't define the alternative hypothesis.

Â Some people don't like that in particular set up but that's often how it's done.

Â In other words, almost always the hypothesis that's being

Â tested in genomic experiments is that the relationship is exactly zero

Â versus the alternative that it's not equal to zero.

Â So that's not really defining a specific alternative you're just

Â trying to reject the case that the null hypothesis makes sure that there's exactly

Â a zero relationship.

Â So if we look back at our linear regression model so

Â suppose we're modeling gene expression as a linear function of age.

Â Then the way that we define these hypotheses in terms of mathematical

Â notation is that the null hypothesis is that there's no relationship.

Â Between age and expression.

Â In other words, this b coefficient is exactly zero.

Â And the alternative is that it's anything other than zero, so

Â it could be positive or negative.

Â We don't have a specific alternative in this case.

Â So the goal is to basically set up a scenario where we can quantify,

Â is b1 equal to zero and then say something about sort of

Â the basically try to reject that hypothesis and say something

Â about what is the distribution of the data look like under that hypothesis?

Â And are we seeing something that's very extreme with respect to that distribution.

Â So in general it's not possible to accept the null.

Â You're going to try to reject it, but it's not possible to accept it.

Â And it's definitely not possible to accept the alternative.

Â You can often reject the null, and you're going to say,

Â we're claiming in favor of the alternative.

Â But it doesn't mean that you've necessarily identified exactly what

Â the alternative is.

Â Ideally, when you're building your statistic to measure this relationship,

Â it's set up so that it's monotone.

Â In other words, as that statistic gets really big, or it gets really small,

Â then it's sort of less likely to be null.

Â In other words, it's more of an extreme observation from the null distribution.

Â If you change the variables that you're adjusting for

Â in your regression model, actually the null hypothesis change

Â because the null hypothesis depends on all the regression coefficients you're saying.

Â In the previous model we only had one coefficient that we were fitting and so

Â the null and the alternative hypothesis were very easy to define.

Â But if you have a large number of adjustment variables

Â then you're saying that the effect is equal to zero.

Â Once you adjust for all of the other variables variables in your model.

Â In general,

Â you have to be very careful to make sure that the null makes intuitive sense.

Â And you can really twist yourself into knots if you're not careful about

Â defining, very clearly in advance, what would be the no effect scenario.

Â In this particular modeling strategy.

Â And it's very important to get this step right.

Â This is the reason why hypothesis testing is often very highly

Â criticized by a number of statisticians is because it's very easy

Â to get the null hypothesis wrong.

Â Or at least to argue about what the null hypothesis is.

Â This is actually a case where this was actually discussed in very great depth.

Â On a blog, on Lior Pachter's blog,

Â he talked about a particular null hypotheses from a particular paper

Â that he disagreed with the way they had defined the null hypotheses, and

Â then discussed the how you would define that null hypothesis in a case.

Â I'm pointing you to this article because there was this very long discussion

Â in the comments section on this blog post about how to define the null hypothesis.

Â And I think it gives you some insight into looking at

Â how difficult it is to get the null hypothesis right,

Â particularly when dealing with complex and high dimensional genomic data.

Â But this is a point that's worth paying a lot of attention to.

Â Again I'd point you to this inference class if you want to learn a lot

Â more about hypothesis testing.

Â And the Statistics and R for the Life Sciences course also has a lot more about

Â inference and hypothesis testing in particular if you care about that.

Â But again the key point is to remember to get your null and

Â alternative hypothesis right.

Â And the null hypothesis in general, usually for genomics,

Â is defining what would this model look like if there was no effect?

Â Coursera provides universal access to the worldâ€™s best education,
partnering with top universities and organizations to offer courses online.