12.07.2015 Views

Analysis of microarray data - VSN International

Analysis of microarray data - VSN International

Analysis of microarray data - VSN International

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

The final example in this section examines the F values inthe column FRatio[1] from the file “Hyb-ANOVA.gwb”,which contains the treatment F-ratio for the results fromthe Single Channel ANOVA for the 9 Affymetrix slidesconsidered earlier. These are compared with the expectedvalues from an F distribution with 2 and 6 degrees <strong>of</strong>freedom. As F-ratios are very skewed, this comparison isbest plotted on a probability scale. The resulting graph isshown to the right, and this shows the same result asfound in the False Discovery Rate menu for the same<strong>data</strong>, i.e. many more <strong>of</strong> the probes shown significantdifferential expression than expected under the nullhypothesis <strong>of</strong> no differential expression.<strong>Analysis</strong> <strong>of</strong> Microarray Data in GenStat 49Display volcano plotThis menu shows jointly the level <strong>of</strong> differentialexpression and the significance on a single graph. It iscalled a volcano plot as the points typically take up a vshape looking like the ash spewing out <strong>of</strong> a volcano. Thisis because it is rare to get points with very small levels <strong>of</strong>differential expression, but a high level <strong>of</strong> significance,and vice versa. In addition, as the significance is onlypositive, the negative and positive expressions are plottedagainst positive values, generally giving a symmetric plot(unless there is a big difference in the number <strong>of</strong>positively and negatively expressed log-ratios). Openingthe menu gives the window to the right. In this case, wewill look at the CEst[1], the contrast between DM andControl from the file “Estimates13-6-9.gsh”, andCPVal[1], the corresponding probability under a twosidednull hypothesis <strong>of</strong> no differential expression. Todemonstrate the ability to colour the points on anothervariate, we will use the mean intensity <strong>of</strong> the spots for aprobe over the slides. To add this variate to the Estimatesspreadsheet, open the full <strong>data</strong> set “Data13-6-9.gwb” anduse the Spread | Calculate | Summary Stats menu toobtain the means <strong>of</strong> intensity for each probe. Thecompleted dialog that does this is shown to the right.Clicking OK will create a spreadsheet containing themeans (shown on the next page).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!