% \VignetteEngine{knitr::knitr}
% \VignetteIndexEntry{TPP_introduction_NPARC}
% \VignettePackage{TPP}
\documentclass[10pt,a4paper,twoside]{article}

\usepackage[utf8]{inputenc}
\usepackage{graphicx}
\usepackage{forloop}

%\pagestyle{empty}
<<style-knitr, eval=TRUE, echo=FALSE, results='asis'>>=
BiocStyle::latex()
@

\bioctitle[TPP]{Analyzing Thermal Proteome Profiling experiments by NPARC (Non-Parametric Analysis of Response Curves)}
\author{Dorothee Childs, Nils Kurzawa\\
European Molecular Biology Laboratory (EMBL),
\\ Heidelberg, Germany\\
\texttt{dorothee.childs@embl.de}}

\date{\Rpackage{TPP} version \Sexpr{packageDescription("TPP")$Version} (Last revision 2016-08-01)}

\begin{document}
<<opts, include=FALSE, echo=FALSE>>=
knitr::opts_chunk$set(concordance=TRUE, 
                      eval = TRUE,
                      cache = TRUE,
                      resize.width="0.45\\textwidth",
                      fig.align='center',
                      tidy = FALSE,
                      message=FALSE)
@

\maketitle

\begin{abstract}
This document demonstrates how to analyze TPP-TR (temperature range) experiments by the NPARC approach. NPARC is a recent extension to the \Rpackage{TPP} package. It offers a novel methodology to model the temperature dependent melting behaviour of each protein, and to detect significant changes in this behavior due to changes in experimental conditions like drug treatment \cite{savitski2014tracking}.

In brief, the melting curve of each protein is represented by natural splines, either once for all conditions (null model), or separately for different experimental conditions (alternative model). Condition specific effects are then detected by testing for significant improvements in goodnes-of-fit of the alternative model relative to the null model. % to do: cite our manuscript
\end{abstract}

\tableofcontents

\section{Installation and input data}
This document describes the analysis of TPP-TR data by the NPARC approach.
General information regarding package installation and format of the input data is provided in the complementary \href{https://www.bioconductor.org/packages/release/bioc/vignettes/TPP/inst/doc/TPP_introduction.pdf}{introductory vignette}.

\section{Getting started}
First, we load the \Rpackage{TPP} package by
<<package>>=
library("TPP")
@
Note that this command only works if the package has been correctly installed beforehand, as described in the \href{https://www.bioconductor.org/packages/release/bioc/vignettes/TPP/inst/doc/TPP_introduction.pdf}{introductory vignette}.

Next, we load the example data for the analysis (treatment of K562 cells using the clinical HDAC
inhibitor panobinostat\cite{franken2015}):
<<load_tr_data>>=
data("hdacTR_smallExample")
@
This command loads two objects: 
\begin{enumerate}
\item \texttt{hdacTR\_data}: a list of data frames that contain the measurements to be analyzed,
\item \texttt{hdacTR\_config}: a configuration table with details about each experiment.
\end{enumerate}

Now we define the desired output path for the results:
<<result_path_TR>>=
resultPath = file.path(getwd(), 'NPARC_Vignette_Example')
@

We start the workflow by typing
<<analyzeTR, eval = TRUE>>=
TRresults <- analyzeTPPTR(configTable = hdacTR_config, data = hdacTR_data, 
                          method = "splinefit", nCores = 2, resultPath = resultPath, 
                          plotCurves = FALSE) 
@
Here, the new methodology is activated by the argument \texttt{method = "splinefit"} (available since package version 3.0.0)
This performs the spline fitting procedure in parallel on a maximum of two CPUs (requirement for package vignettes). Without specifying the \texttt{nCores} argument, fitting is performed by default on the maximum number of CPUs on your device.



\bibliography{TPP_references}
\end{document}
