-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathSec-BiocCase.Rnw
45 lines (40 loc) · 1.22 KB
/
Sec-BiocCase.Rnw
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
\section{A Bioc script}
\begin{frame}
\begin{block}{}
A typical Bioconductor script for microarray data analysis.
\begin{itemize}
\item Getting data
\item Import data into \R using dedicated infrastructure
\item Analyse
\item Save script, plots and objects
\end{itemize}
\end{block}
\end{frame}
\begin{frame}
\begin{block}{Using a subset of the \textit{tg-gates} data}
\begin{itemize}
\item \texttt{E-MTAB-800}: transcription profiling by array of rat liver and kidney after exposure to approximately 130 chemicals collected from repeat dosing studies\footnote{\url{http://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-800/}}
\item Downloaded and unzipped \texttt{E-MTAB-800.raw.1.zip}
\item Using only a subset of files below.
\end{itemize}
\end{block}
\end{frame}
\begin{frame}[fragile]{\exercise}
\begin{block}{Loading libraries}
<<biocenv, eval = FALSE>>=
library("Biobase")
library("affy")
@
\end{block}
\begin{block}{Reading data}
<<readcel, eval = FALSE>>=
rawdata <- ReadAffy(filenames = flnms)
@
\end{block}
\begin{block}{Normalisation}
<<norm, eval = FALSE>>=
eset <- rma(rawdata)
@
\end{block}
See \texttt{Exercise-06.R}
\end{frame}