Introduction to Bayesian Statistics and an Application and an - - PowerPoint PPT Presentation

introduction to bayesian statistics
SMART_READER_LITE
LIVE PREVIEW

Introduction to Bayesian Statistics and an Application and an - - PowerPoint PPT Presentation

Introduction to Bayesian Statistics Introduction to Bayesian Statistics and an Application and an Application Timothy M. Bahr Unconfounding the Confounded: Separating Introduction Treatment and Batch Effects in Confounded Definitions


slide-1
SLIDE 1

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Introduction to Bayesian Statistics and an Application

Unconfounding the Confounded: Separating Treatment and Batch Effects in Confounded Microarray Experiments Timothy M. Bahr

Department of Statistics Brigham Young University

March 16, 2009

slide-2
SLIDE 2

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Introduction Who am I?

Tim Bahr, Undergrad...

slide-3
SLIDE 3

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Introduction Who am I?

Tim Bahr, Undergrad...

◮ 22, B.S. in Statistics,

emphasis: Biostat

slide-4
SLIDE 4

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Introduction Who am I?

Tim Bahr, Undergrad...

◮ 22, B.S. in Statistics,

emphasis: Biostat

◮ My first intro to Statistics

in High School

slide-5
SLIDE 5

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Introduction Who am I?

Tim Bahr, Undergrad...

◮ 22, B.S. in Statistics,

emphasis: Biostat

◮ My first intro to Statistics

in High School

◮ Fascination with the

Numerical Patterns in Science

slide-6
SLIDE 6

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Introduction Who am I?

Tim Bahr, Undergrad...

◮ 22, B.S. in Statistics,

emphasis: Biostat

◮ My first intro to Statistics

in High School

◮ Fascination with the

Numerical Patterns in Science

◮ Future Goals

slide-7
SLIDE 7

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Introduction Who are you?

Bioinformatics

◮ Majors? ◮ Math/Stat Background? ◮ Microarrays? ◮ Research? ◮ Why Bioinformatics? ◮ Can I tell you what I think

about Bioinformatics?

slide-8
SLIDE 8

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

slide-9
SLIDE 9

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

◮ Bayesian Statistics >>> statistical inferences on

experimental data + prior knowledge.

slide-10
SLIDE 10

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

◮ Bayesian Statistics >>> statistical inferences on

experimental data + prior knowledge.

◮ Classical (Frequentist) Statistics >>> data from

  • bservations or experiments only.
slide-11
SLIDE 11

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

◮ Bayesian Statistics >>> statistical inferences on

experimental data + prior knowledge.

◮ Classical (Frequentist) Statistics >>> data from

  • bservations or experiments only.

◮ Prior Distribution: The distribution we assume our

parameters come from.

slide-12
SLIDE 12

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

◮ Bayesian Statistics >>> statistical inferences on

experimental data + prior knowledge.

◮ Classical (Frequentist) Statistics >>> data from

  • bservations or experiments only.

◮ Prior Distribution: The distribution we assume our

parameters come from.

◮ Gibbs Sampling (simplification): An algorithm that

allows us to give interatively infer point estimates for “random” parameters.

slide-13
SLIDE 13

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

slide-14
SLIDE 14

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

◮ Biostatistics: The application of statistics to a wide

range of topics in biology.

slide-15
SLIDE 15

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

◮ Biostatistics: The application of statistics to a wide

range of topics in biology.

◮ Gene Expression Microarray: A high-throughput

technology in molecular biology used to detect gene expression levels in a cellular sample.

slide-16
SLIDE 16

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Definitions

◮ Biostatistics: The application of statistics to a wide

range of topics in biology.

◮ Gene Expression Microarray: A high-throughput

technology in molecular biology used to detect gene expression levels in a cellular sample.

◮ Confounded Experiment: when two or more variables

vary together so that it is impossible to separate their unique effects.

slide-17
SLIDE 17

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Probabilistic inference that computes the distribution of the model parameters and gives prediction for previously unseen input values probabilistically. Freqentist

slide-18
SLIDE 18

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Probabilistic inference that computes the distribution of the model parameters and gives prediction for previously unseen input values probabilistically. Freqentist

◮ θ, parameters, are fixed

and unknown

slide-19
SLIDE 19

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Probabilistic inference that computes the distribution of the model parameters and gives prediction for previously unseen input values probabilistically. Freqentist

◮ θ, parameters, are fixed

and unknown

◮ X, random variables

(data), are random Bayesian

slide-20
SLIDE 20

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Probabilistic inference that computes the distribution of the model parameters and gives prediction for previously unseen input values probabilistically. Freqentist

◮ θ, parameters, are fixed

and unknown

◮ X, random variables

(data), are random Bayesian

◮ θ, parameters, are random

and unknown

slide-21
SLIDE 21

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Probabilistic inference that computes the distribution of the model parameters and gives prediction for previously unseen input values probabilistically. Freqentist

◮ θ, parameters, are fixed

and unknown

◮ X, random variables

(data), are random Bayesian

◮ θ, parameters, are random

and unknown

◮ X, random variables

(data), are random “If you want to work on really interesting problems [Bayesian Inference] is where those problems lie”

  • Don Rubin, Ph.D., Dept. Chair, Harvard Statistics
slide-22
SLIDE 22

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference The idea of a prior

slide-23
SLIDE 23

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference The idea of a prior

◮ Frequentists assume a parameter is fixed:

◮ For example X ∼ N(µ, σ2) ◮ µ is a fixed unknown value

slide-24
SLIDE 24

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference The idea of a prior

◮ Frequentists assume a parameter is fixed:

◮ For example X ∼ N(µ, σ2) ◮ µ is a fixed unknown value

◮ What if µ is not fixed? What if it too can assume a

distribution with variation

slide-25
SLIDE 25

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference The idea of a prior

◮ Frequentists assume a parameter is fixed:

◮ For example X ∼ N(µ, σ2) ◮ µ is a fixed unknown value

◮ What if µ is not fixed? What if it too can assume a

distribution with variation

◮ We assume a prior on µ. i.e. µ ∼ N(mµ, s2 µ)

slide-26
SLIDE 26

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Bayes’ Theorem, based on basic theories of probability: π(θ|x) = f(x|θ)π(θ)

  • f(x|θ)π(θ)dθ

(1)

slide-27
SLIDE 27

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Bayes’ Theorem, based on basic theories of probability: π(θ|x) = f(x|θ)π(θ)

  • f(x|θ)π(θ)dθ

(1)

◮ π(θ|x) is the posterior distribution of our parameters, θ.

slide-28
SLIDE 28

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Bayes’ Theorem, based on basic theories of probability: π(θ|x) = f(x|θ)π(θ)

  • f(x|θ)π(θ)dθ

(1)

◮ π(θ|x) is the posterior distribution of our parameters, θ. ◮ f(x|θ) is the likelihood of the data

slide-29
SLIDE 29

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

Bayes’ Theorem, based on basic theories of probability: π(θ|x) = f(x|θ)π(θ)

  • f(x|θ)π(θ)dθ

(1)

◮ π(θ|x) is the posterior distribution of our parameters, θ. ◮ f(x|θ) is the likelihood of the data ◮ π(θ) is the prior distribution assumed on our parameters,

θ.

slide-30
SLIDE 30

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Bayesian Inference

In the End: Estimate Parameters

◮ We solve for the posterior of the parameters ◮ Use different methods to estimate an “optimum” value of

  • ur parameters.

◮ Take the Expected Value of a Parameter ◮ Gibbs Sampling ◮ Metropolis-Hastings

slide-31
SLIDE 31

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Microarrays What is a Microarray?

slide-32
SLIDE 32

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Microarrays What is a Microarray?

◮ We use microarrays to detect gene expression levels for a

given cellular sample.

slide-33
SLIDE 33

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Microarrays

slide-34
SLIDE 34

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Microarrays

slide-35
SLIDE 35

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Microarrays

slide-36
SLIDE 36

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments

◮ Consider a fertilizer experiment with corn:

First, An “unconfounded” experiment.

◮ 1 plot of corn; left half- control (no fertilizer), right half-

treatment (Fertilizer)

◮ Differences in corn quality can be attributed to the

treatment effect.

slide-37
SLIDE 37

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn:

slide-38
SLIDE 38

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn: ◮ Plot 1 (Batch 1)

slide-39
SLIDE 39

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn: ◮ Plot 1 (Batch 1) ◮ Control (no fertilizer)

slide-40
SLIDE 40

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn: ◮ Plot 1 (Batch 1) ◮ Control (no fertilizer) ◮ Plot 2 (Batch 2) - 1 mi. away

slide-41
SLIDE 41

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn: ◮ Plot 1 (Batch 1) ◮ Control (no fertilizer) ◮ Plot 2 (Batch 2) - 1 mi. away ◮ Treatment (New Fertilizer)

slide-42
SLIDE 42

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn:

If we observe a significant difference between the corn quality

  • f the two plots (batches), can we attribute this difference to

the fertilizer?

slide-43
SLIDE 43

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn:

  • No. The difference may be due to the treatment effect, the

plot (batch effect), or a combination of the two.

slide-44
SLIDE 44

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn:

The Treatment Effect is confounded with the Plot or Batch Effect.

slide-45
SLIDE 45

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments What is a Confounded Experiment?

◮ Consider a fertilizer experiment with corn:

The same principle applies to microarray experiments.

slide-46
SLIDE 46

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

slide-47
SLIDE 47

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

◮ Microarrays prepared at different times, in different places,

by different people etc. ... are often confounded by batch effects.

slide-48
SLIDE 48

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

◮ Microarrays prepared at different times, in different places,

by different people etc. ... are often confounded by batch effects.

◮ We are not interested in the the batch effect. We want to

subtract it out.

slide-49
SLIDE 49

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

◮ Microarrays prepared at different times, in different places,

by different people etc. ... are often confounded by batch effects.

◮ We are not interested in the the batch effect. We want to

subtract it out.

◮ Our algorithm uses statistical methods to adjust for the

Batch effect in confounded microarray experiments.

slide-50
SLIDE 50

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

Why?

slide-51
SLIDE 51

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

Why?

◮ Often times biologists can save money by using data that

was obtained in previous experiments.

slide-52
SLIDE 52

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

Why?

◮ Often times biologists can save money by using data that

was obtained in previous experiments.

◮ Inter-lab collaboration becomes much more reliable when

batch effects are accounted for.

slide-53
SLIDE 53

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

Our Solution

Our method allows precise estimation of the batch effect and the treatment effect.

slide-54
SLIDE 54

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

Our Solution

Our method allows precise estimation of the batch effect and the treatment effect.

◮ A dynamic linear model

slide-55
SLIDE 55

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

Our Solution

Our method allows precise estimation of the batch effect and the treatment effect.

◮ A dynamic linear model ◮ Novel yet Appropriate Assumptions

slide-56
SLIDE 56

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Confounded Experiments: Microarrays

Our Solution

Our method allows precise estimation of the batch effect and the treatment effect.

◮ A dynamic linear model ◮ Novel yet Appropriate Assumptions ◮ Bayesian Statistical Methods

slide-57
SLIDE 57

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Our Model The Model: yig = µg + Xiαg + Ziτg + εig (2)

slide-58
SLIDE 58

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Our Model The Model: yig = µg + Xiαg + Ziτg + εig (2)

◮ yig - the “expression level” for a sample i from gene g

slide-59
SLIDE 59

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Our Model The Model: yig = µg + Xiαg + Ziτg + εig (2)

◮ yig - the “expression level” for a sample i from gene g ◮ µg - an overall average for gene g

slide-60
SLIDE 60

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Our Model The Model: yig = µg + Xiαg + Ziτg + εig (2)

◮ yig - the “expression level” for a sample i from gene g ◮ µg - an overall average for gene g ◮ αg - the Treatment Effect for gene g

slide-61
SLIDE 61

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Our Model The Model: yig = µg + Xiαg + Ziτg + εig (2)

◮ yig - the “expression level” for a sample i from gene g ◮ µg - an overall average for gene g ◮ αg - the Treatment Effect for gene g ◮ τg - the Batch Effect for gene g

slide-62
SLIDE 62

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Our Model The Model: yig = µg + Xiαg + Ziτg + εig (2)

◮ yig - the “expression level” for a sample i from gene g ◮ µg - an overall average for gene g ◮ αg - the Treatment Effect for gene g ◮ τg - the Batch Effect for gene g ◮ εig - error for sample i from gene g

slide-63
SLIDE 63

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Solution Formulation and Assumptions

First, an “unconfounded” formulation.

Difference between treatment and control can be attributed to “treatment effect.”

slide-64
SLIDE 64

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Solution Formulation and Assumptions

We can’t differentiate the values of αg and τg.

slide-65
SLIDE 65

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Solution Formulation and Assumptions

We assume treatment, αg, has no effect on group 2 genes

slide-66
SLIDE 66

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Solution Formulation and Assumptions

Determine which genes in each group >>> estimate αg and τg.

slide-67
SLIDE 67

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

How do we estimate αg and τg?

Gibbs Sampling

◮ A Bayesian Method ◮ Gives us the power to estimate which genes are in each

group

◮ Iteratively estimates values until sequence converges

slide-68
SLIDE 68

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Estimating αg

5 10 15 20 25 30 35 2 4 6 8 10 12

Estimating Alpha

iterations alpha actual

slide-69
SLIDE 69

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Estimating τg

5 10 15 20 25 30 35 2 3 4 5

Estimating Tau

iterations tau actual

slide-70
SLIDE 70

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Application

Possible Applications

◮ Microarrays in Cancer Research ◮ Clinical use of microarrays for diagnosis ◮ Possible applications in non-array experiments

slide-71
SLIDE 71

Introduction to Bayesian Statistics and an Application Timothy M. Bahr Introduction Definitions Bayesian Statistics Microarrays Confounded Experiments Model Gibbs Sampling Application

Acknowledgments

◮ W. Evan Johnson, mentor ◮ Nathaniel Gustafson, programmer ◮ BYU Dept. of Statistics ◮ Johnson Lab

(old picture)