Explore
For Enterprise
Join for Free
Log In

艺术与人文
商务
计算机科学
数据科学
信息技术
健康
数学和逻辑
个人发展
物理科学与工程
社会科学
语言学习
學位
證書

瀏覽 Coursera 的全部課程

瀏覽
搜索
企業版
登錄
免費加入

Exploratory Analysis in R Part III (7:26)

Loading...

基因组数据科学所需的统计学

约翰霍普金斯大学

4.1（182 個評分） | 16K 名學生已註冊

課程 7（共 8 門，基因组数据科学專項課程）

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.

查看授課大綱

您將學習的技能

Statistics, Data Analysis, R Programming, Biostatistics

審閱

4.1（182 個評分）

5 stars
93 ratings
4 stars
52 ratings
3 stars
18 ratings
2 stars
7 ratings
1 star
12 ratings

LR

May 23, 2016

I have really enjoyed the course and I have learnt different concepts relevant for my current study.\n\nYurany

TK

Mar 25, 2019

Very good course and useful understanding statistical aspects of data.

從本節課中

Module 1

This course is structured to hit the key conceptual ideas of normalization, exploratory analysis, linear modeling, testing, and multiple testing that arise over and over in genomic studies.

Module 1 Overview (1:07)1:07

Reproducible Research (3:42)3:42

Achieving Reproducible Research (5:02)5:02

R Markdown (6:26)6:26

The Three Tables in Genomics (2:10)2:10

The Three Tables in Genomics (in R) (3:46)3:46

Experimental Design: Variability, Replication, and Power (14:17)14:17

Experimental Design: Confounding and Randomization (9:26)9:26

Exploratory Analysis (9:21)9:21

Exploratory Analysis in R Part I (7:22)7:22

Exploratory Analysis in R Part II (10:07)10:07

Exploratory Analysis in R Part III (7:26)7:26

Data Transforms (7:31)7:31

Clustering (8:43)8:43

Clustering in R (9:09)9:09

教學方

Jeff Leek, PhD
Associate Professor, Biostatistics

腳本

This is the third lecture on exploratory data analysis so again,

I'm going to set things up a little bit quickly.

You can watch the first lecture on exploratory data analysis if you want to

see why I'm setting up the parameters to be so

that the plots are pretty and then loading all the basic libraries.

And then, at least the library this lecture depends on, and

then I'm going to load the data set in from the web.

So now that I'm up to speed, I'm going to do in the sort of the third

part of the exploration of a data set, I usually check for consistency.

And so you usually want to be able to check for

consistency using some kind of external data, if you have it available.

As an example of that here, I think what I'm going to do is I'm going to show you

how I compare a data set that I have to some meta data or some annotation that

I have from somewhere else to make sure they appear to be doing the right thing.

So here, what I'm going to do is I'm going to get the IDs for

the features and so here I've got these IDs, that are ensemble IDs.

Then, I'm going to extract the chromosome information from those ensemble IDs.

So if I look at that information, I've got chromosome labels now, for

each of the ensemble IDs in the data set.

Once I've done that,

hopefully that has the same dimension as the data set, but it doesn't.

So okay, obviously, we're already kind of stuck in a problem.

And one reason why is you can see that some of the chromosome,

some of the samples are duplicate labeled.

So I'm going to remove all of the duplicated

values from the chromosome data set.

So once I remove those duplicated values, I can see

that the chromosome data set has the same dimension as the expression data set.

But we want to go one step further and make sure that the chromosome data

matches the information in the expression data, so

what I'm going to do is I'm going to say, do all.

I'm going to need the all command to basically check do all of the ensemble IDs

that are in this chromosome data set match the real names of the expression data set.

That's true, so that's good.

Basically, again first, I basically check to make sure that this ensemble ID

matches the ensemble ID in the expression data, so that the IDs all match up.

Now, what I'm going to do, is I'm going to filter to just the chromosome Y samples.

To do that, I'm going to have to convert the data frame here.

When you're applying this filtering command,

you're going to have to convert it to a data frame.

And then, you should be able to apply the filtering command in it, so now,

I've got a smaller data set that just consists of

all the genes measured on the Y chromosome.

And so, the first thing that I want to do is then check to see, so

I'm going to take the sum in each column.

So take out the genes from the Y chromosome.

I'm going to take the sums for

each sample, that's the total count for the Y chromosome genes for each sample.

I'm going to plot that versus gender.

When I do that, I can see that the males

have way more counts than the females, so that's kind of a good independent check.

And then, I can overlay the data points if I want to, actually.

If I make that boxplot without coloring it and then I overlay the data points,

so here, points like lines from the previous lecture,

basically just overlays points on top of the current plot.

I'm going to again, plot the column sums of the total counts for the Y chromosome.

Here, I'm using a jitter command,

you'll see that again in a lot of plots that I make like this, with box plots.

So jitter basically adds random noise to this value for gender, so it basically is

going to make the points so they don't all land on top of each other.

And then I'm going to color it by the gender, and so you can see that

the females have counts mostly of zero, and the males have some positive counts.

So that makes us think that maybe there is

sort of an external validation that the plot is looking okay.

The next thing that I might want to do,

is look at some kind of multi-variate plots of the data.

The way that people often do that is with the heatmap command.

To do that, the first thing I'm going to do is I'm going to filter the data set

down to just the genes that have a very high roamings, so

I'm going to do that because it'll make it quicker to make this heatmap.

But you can obviously make it with much larger number of values,

if you have time to wait for the plot to make.

So then, I'm going to apply the heatmap command to that matrix.

When I do that, I can see here, I've got sample by sample.

And then in each of the columns and then in each of the rows I've got the genes.

What it's done already is it's clustered these data and

we're going to talk about clustering in a bit.

But one thing that you might want to do first of all,

is maybe make the colors be different than the heat colors.

You can do that by defining a new set of colors, a new color ramp.

You can do that with the color ramp palette function.

So here, I'm going to say blend the colors from the third color in the tropical

palette that I defined, to white to the second color.

And I'm going to say make a nine-color color ramp so you can see.

Now this has defined the colors in hex decimal format.

The next thing that we can do is we can make a heat map, but I can tell it,

instead of using the colors you were using before, use these new colors.

So now you can see, you got this kind of, nice, white, blue, pink colors format.

The other thing you can do is if you don't want to cluster them,

you actually want to see the samples in the order that they originally were in,

you can again, use the heatmap command.

And we're going to use the pretty colors I just defined.

Then, if you say, Rowv=NA,

Colv=NA, it will remove those sort of clustering diagrams.

So now, these are actually the samples in the order that they are in the data

set you have.

So one thing you can already kind of detect from this is that you can see,

right here, that these samples,

here, have a little bit higher values than the other samples.

That's helpful to be able to immediately check that out.

The other thing you can do in the G+ package is use the heatmap_2

function to add, say, for example, a scale to this plot,

because sometimes I said that these values were higher.

But how do you know what's higher or lower in terms of the colors?

And so I'm going to again, use basically,

almost the exact same syntax as you would have with the heatmap.

But now, I'm going to tell it not to cluster by using dendogram=none.

And then, I'm going to tell it to scale each row.

Now, what it's going to do is it's going to make this plot same as it did

before, but now it's added this color key histogram, here, for you.

It says most of the values are near zero in the matrix and

that blue is high values and pink is low values.

And so then, you can kind of see it right there, in the data set.

So that's the sort of multivariate way I go about looking at the data.

After I've done all of those checks, and I've sort of done all that exploration,

I start thinking about plots that I'm going to make related to model it.

探索我們的目錄

免費加入並獲得個性化推薦、更新和優惠。

Coursera 致力於普及全世界最好的教育，它與全球一流大學和機構合作提供在線課程。

© 2019 Coursera Inc. 保留所有權利。

通過 App Store 下載

通過 Google Play 獲取

Coursera
關於
管理團隊
工作機會
目錄
證書
MasterTrack™ 證書
學位
企業版
政府版

社區
Learners
Partners
Developers
Beta Testers
Translators

連接
博客
Facebook
領英
Twitter
Instagram
YouTube
技術博客

更多
條款
隱私
幫助
內容訪問
媒體
聯繫我們
目錄
附屬公司