git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/branches/RELEASE_2_13/madman/Rpacks/MineICA@84468 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 |
@ |
... | ... |
@@ -696,13 +696,13 @@ When a variable is quantitative, its association with a component can be studied |
696 | 696 |
## for correlations exceeding this threshold. |
697 | 697 |
resQuant <- quantVarAnalysis(params=params, icaSet=icaSetMainz, keepVar="age", |
698 | 698 |
typeCor="pearson", cutoffOn="cor", |
699 |
- cutoff=0.3, adjustBy="none", |
|
699 |
+ cutoff=0.3, adjustBy="none", |
|
700 | 700 |
path="quantVarAnalysis/", filename="quantVar") |
701 | 701 |
@ |
702 | 702 |
|
703 |
-The absolute correlation between age and sample contributions exceeds 0.3 only for the third component. |
|
703 |
+The absolute correlation between age and sample contributions exceeds 0.3 only for the second component. |
|
704 | 704 |
<<quantVarAncor, echo=TRUE, print=FALSE>>= |
705 |
-resQuant$cor[3] |
|
705 |
+resQuant$cor[2] |
|
706 | 706 |
@ |
707 | 707 |
|
708 | 708 |
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. |
... | ... |
@@ -711,7 +711,7 @@ The function creates a HTML file "quantVar.htm" containing correlations values, |
711 | 711 |
|
712 | 712 |
\begin{figure}[htbp] |
713 | 713 |
\centering |
714 |
-\includegraphics[width=0.8\linewidth]{mainz/quantVarAnalysis/plots/3_age.png} |
|
714 |
+\includegraphics[width=0.8\linewidth]{mainz/quantVarAnalysis/plots/2_age.png} |
|
715 | 715 |
\caption[Scatter plot of age vs sample contributions.]{Scatter plot of AGE vs sample contributions. The witness gene 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. } |
716 | 716 |
\label{fig:exscatter} |
717 | 717 |
\end{figure} |