I had a few quick questions about the limma models. The first is when writing these to a file, the fold changes would be the "Coef" column, not the A column, which is the average across all of the arrays? The second is when I write my bayes model to a file, it only attaches the probe ID, I wanted to see if there was a way to attach the remaining gene annotation but the limma documentation isn't clear on how to accomplish that.

Here is how I am currently performing the fitting:
    contmat   <- makeContrasts(d-s,
                               levels = design)
    fit1      <- lmFit(x[[1]]$avg$E, design)
    fit2      <- eBayes(contrasts.fit(fit1, contmat))
    diff      <- topTable(fit2,
                          coef     = 1,
                          adjust   = "BH",
                          genelist = x[[1]]$avg$genes )
    res       <- decideTests(fit2)
    comp      <- vennDiagram(res)
    write.fit(fit1,
              results  = res,
              sep      = sep,
              file     = paste(c(x[[1]]$bg, x[[1]]$norm, "lmfit.csv"),
              collapse = ""))

    write.fit(fit2,
              results  = res,
              sep      = sep,
              file     = paste(c(x[[1]]$bg, x[[1]]$norm, "be.csv"),
              collapse = ""))

where x[[1]]$avg is the corrected, normalized, replicate averaged array.

Here is some example output from write.fit(fit2...):

A

Coef

t

p.value

F

F.p.value

d - s

10.99

5.352

6.82

6.00E-04

46.55

6.00E-04

0

12.48

5.347

14.13

1.00E-05

199.78

1.00E-05

1

12.84

4.775

11.99

3.00E-05

143.65

3.00E-05

1

13

4.646

10.98

5.00E-05

120.49

5.00E-05

1



Joseph Cornish
Post-Bac IRTA
NIAID EVPS


	[[alternative HTML version deleted]]

