git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/MineICA@83834 bc3139a8-67e5-0310-9ffc-ced21a209358
... | ... |
@@ -562,7 +562,7 @@ To obtain a biological interpretation of the component, it can be useful to stud |
562 | 562 |
In order to identify the gene sets which are enriched in the list of selected (contributing) genes, the function \Robject{runEnrich} uses \verb$R$ \verb$GOstats$ package \cite{Falcon2007Using} which makes use of a hypergeometric distribution to test the over-representation of a gene set in a given list of genes. |
563 | 563 |
<<runEnrich, echo=TRUE, eval=FALSE>>= |
564 | 564 |
## run enrichment analysis on the first three components of icaSetMainz, |
565 |
-## using gene sets from KEGG and ontology 'Biological Process' (BP) of Gene Ontology (GO) |
|
565 |
+## using gene sets from the ontology 'Biological Process' (BP) of Gene Ontology (GO) |
|
566 | 566 |
resEnrich <- runEnrich(params=params,icaSet=icaSetMainz[,,1:3], |
567 | 567 |
dbs=c("GO"), ontos="BP") |
568 | 568 |
@ |
... | ... |
@@ -694,13 +694,13 @@ When a variable is quantitative, its association with a component can be studied |
694 | 694 |
## for correlations exceeding this threshold. |
695 | 695 |
resQuant <- quantVarAnalysis(params=params, icaSet=icaSetMainz, keepVar="age", |
696 | 696 |
typeCor="pearson", cutoffOn="cor", |
697 |
- cutoff=0.3, adjustBy="none", |
|
697 |
+ cutoff=0.3, adjustBy="none", |
|
698 | 698 |
path="quantVarAnalysis/", filename="quantVar") |
699 | 699 |
@ |
700 | 700 |
|
701 |
-The absolute correlation between age and sample contributions exceeds 0.3 only for the third component. |
|
701 |
+The absolute correlation between age and sample contributions exceeds 0.3 only for the second component. |
|
702 | 702 |
<<quantVarAncor, echo=TRUE, print=FALSE>>= |
703 |
-resQuant$cor[3] |
|
703 |
+resQuant$cor[2] |
|
704 | 704 |
@ |
705 | 705 |
|
706 | 706 |
The corresponding scatter plot is available in Figure~\ref{fig:exscatter}. A tendency of the women whose tumors are located at the positive end of the component to be younger indeed appears. |
... | ... |
@@ -709,7 +709,7 @@ The function creates a HTML file "quantVar.htm" containing correlations values, |
709 | 709 |
|
710 | 710 |
\begin{figure}[htbp] |
711 | 711 |
\centering |
712 |
-\includegraphics[width=0.8\linewidth]{mainz/quantVarAnalysis/plots/3_age.png} |
|
712 |
+\includegraphics[width=0.8\linewidth]{mainz/quantVarAnalysis/plots/2_age.png} |
|
713 | 713 |
\caption[Scatter plot of age vs sample contributions.]{Scatter plot of AGE vs sample contributions. The gene witness is \textit{KRT16}. At the bottom of the plot, each sample is represented by a square point whose colour denotes the expression value of the \textit{KRT16} gene. The scale of these colors is denoted by a legend at the upper right of the graph. Note that the gene expression profiles were centered to have mean zero. } |
714 | 714 |
\label{fig:exscatter} |
715 | 715 |
\end{figure} |