Biostatistics Case Studies 2010 - PowerPoint PPT Presentation

Biostatistics case studies 2010
1 / 24

  • Uploaded on
  • Presentation posted in: General

Biostatistics Case Studies 2010. Session 5: Microarray Statistics. Peter D. Christenson Biostatistician http://gcrc. /Biostat. Case Study. A compound found in red grapes improves the health and lifespan of mice on a high calorie diet. Treatment Groups.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

Download Presentation

Biostatistics Case Studies 2010

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript

Biostatistics case studies 2010

Biostatistics Case Studies 2010

Session 5: Microarray Statistics

Peter D. Christenson


Case study

Case Study

Biostatistics case studies 2010

A compound found in red grapes improves the health and lifespan of mice on a high calorie diet.

Treatment groups

Treatment Groups

  • Male middle-aged mice (11 months) were randomized to:

  • Standard Diet (SD): known as AIN-93G

    • N=60; N=5 for gene expression

  • High Calorie (HC): SD + coconut oil → 60% fat.

  • N=55; N=5 for gene expression

  • Resveritrol (HCR): HC + 22.4 mg resveritrol/ kg/day.

  • N=55; N=4 for gene expression (+1 w/ degraded sample)

  • Outcome mortality

    Outcome: Mortality

    0.42 died

    0.58 died

    • Methods: Survival Analysis.

    • 114 week mortality ratio for HCR/HC = 0.42/0.58=0.72, which is a 28% reduction at 114 weeks.

    • Paper reports a “hazard ratio”, which is similar at 0.69, with a p-value of 0.02. How does it differ from the mortality ratio?

    Outcome agility

    Outcome: Agility

    What statistical analysis was done here?

    Outcome clinical markers

    Outcome: Clinical Markers

    What statistical analyses?

    Other outcomes

    Other Outcomes

    Gene expression in liver at age 18 months

    Gene Expression in Liver at age 18 months

    Fourteen Microarray “Experiments”: each of 5+5+4 mice had a separate array run for ~40,000 genes.

    Gene expression data 536 872 numbers

    Gene Expression Data: 536,872 Numbers

    First 2 SD mice. 12 others →

    38,348 rows: each a gene

    Gene expression results

    Gene Expression Results

    HCR over-expressed, compared to HC.

    How were results for (a) and (b) calculated?

    HCR under-expressed, compared to HC.

    Microarray analysis

    Microarray Analysis

    • How can we analyze these data?

    • What are “experimental units”: mice or genes?

    • Consider each gene independently?

    • If so, Ns of 4 and 5 seem small to say much - low power.

    • So, maybe combine genes for larger Ns?

    • Pair up HCR and HC mice, find ratio, and average?

    • Ratio of mean for N=4 HCR and mean for N=5 HC?

    • If p<0.05 is used for each gene, expect many false positives among 38,348 genes.

    • SD among only 5 mice could be large just due to differences from array to array, not biologic diff, and thus miss finding important genes.

    Detectable effects with n 5 per group

    Detectable Effects with N=5 per Group

    Suppose we compare the mean of 5 appropriately scaled #s for a gene’s expression with the mean of 5 in another group, using a t-test.


    So, we need ~ 2SD difference in gene expression to be fairly sure (80%) of detecting this gene with only N=5+5.

    This is a large effect – see next slide.

    Detectable effects with n 5 per group1

    Detectable Effects with N=5 per Group


    Normal Range

    2SD Effect corresponds to 50th→ 97th percentile, about 2/5 of normal range

    Detectable effects with n 5 per group2

    Detectable Effects with N=5 per Group

    So, how can we try to avoid missing genes that are important, but are not detected with p<0.05?

    Recall that p<0.05 corresponds to approximately:

    |t| =|effect/SE(effect)| = |Δ/SE(Δ)| = |signal/noise| >2

    where noise is a function of ~ SD/ sqrt(N).

    • Thus, if ↑ N is not possible to reduce noise, we can:

    • Try to reduce SD, or

    • Ignore SD and base the decision for gene selection on the signal, i.e., effect, i.e., mean differential expression, only.

    Microarray analysis 1 try to reduce sd

    Microarray Analysis: 1. Try to reduce SD

    Here, SD is the SD among the expressions for 5 mice in a group.

    How can we “reduce SD”? Isn’t it natural subject-to-subject heterogeneity, a characteristic of the population?

    This SD is among measured expression, which includes both array-to-array error and subject-to-subject heterogeneity. (Confounded-there is no internal control.)

    We try to statistically remove some of the inherent array-to-array error through normalization.

    Side point on microarray design

    Side Point on Microarray Design

    • Single Channel Chip: 1 sample, many probes.

      • No replicated measures. This study.

    • Two Channel Chip: 2 samples, possibly fewer probes that are common for both samples.

      • One sample may be an internal control.

      • The two samples may be matched, e.g.,

        • Two conditions, times, etc, for the same subject.

        • Twins, littermates, etc, treated differently.



    • There are many ways to normalize. They exploit the assumption that most of 1000s of genes will be the same in many subjects. Two common methods:

    • Global: All genes in an array are multiplied by the ratio of the (global) mean over all genes for all arrays to the mean over all genes for this array. E.g., array1 has mean 1000 and fourteen arrays have mean 900, multiply by 0.90.

    • Z-score: Replace expression x by z=(x-mean)/SD, where mean and SD are over genes for this array. Expression becomes # of SDs deviant from gene mean.

    Microarray analysis 2 ignore sd

    Microarray Analysis: 2. Ignore SD

    • Here, SD is the SD among the expressions for 5 mice in a group.

    • Use an effect measure for each gene, such as the ratio of mean of 4 HCR to the mean of 5 HC, usually standardized to a “normal range” as with z-scores.

    • Usually select genes with either:

      • Ratio>c or <1/c, some c such as 1.5 or 2.

      • A specified number or percent of genes with largest or smallest ratios.

    Microarray analysis this study

    Microarray Analysis: This study

    • Genes selected with both:

      • “Z-Ratio” >1.5 or <-1.5.

      • The p-value from a z-test for comparing the mean z-score of 4 HCR mice to the mean of 5 HC mice is <0.05.

    Raw expression is normalized within each array by z-scores on log(expression).

    The Z-Ratio is the difference between the mean z-score of 4 HCR mice to the mean of 5 HC mice (which is the numerator for the z-test), divided by the SD of these differences over different genes.

    Microarray analysis gene hsd3b5

    Microarray Analysis: Gene Hsd3b5

    Use raw data to generate results for the most up-regulated gene.

    Microarray analysis gene hsd3b51

    Microarray Analysis: Gene Hsd3b5

    Microarray analysis gene hsd3b52

    Microarray Analysis: Gene Hsd3b5

    Two Sample T-Test for HCR vs. HC on Gene Hsd3b5

    N Mean SD SE

    HCR 4 1.136 0.634 0.32

    HC 5 -0.555 0.107 0.048

    Diff 1.691

    95% CI for Diff: ( 1.02, 2.362)

    T-Test T = 5.96 P = 0.0006

    Antilog(1.691) =~ 5.42 fold greater HCR expression

    “Z-Ratio” = Diff of logs/SD = 1.691/0.14 = 11.99

    Here, SD=0.14 is among these diffs over genes.

    Expected identified genes among 38 348 genes using p values

    Expected Identified Genes among 38,348 Genes using p-values

    Suppose the decision rule is to declare a particular gene important if its mean expression in HCR mice differs enough from that for HC mice so that p<0.05:

    Significantly less → down-regulated.

    Significantly more → up-regulated.

    Then the expected number of identified genes among, say, 38,000 that are not affected (false positives) is:

    0.05*38,000 = 1900

    Thus, confirmatory analyses such as PCR are needed.

  • Login