Residual Disease Paper
========================================================
Plotting RD Status
--------------------------------------------------------

by Shelley Herbrich

*May 8, 2013*


```{r options, echo=FALSE}

opts_chunk$set(tidy=TRUE, message=FALSE, warning=FALSE)

```


## 1 Executive Summary

### 1.1 Introduction

In this report, we present the script used to display RD by ADH1B and FABP4.

### 1.2 Data & Methods

We use the "PCRResults.RData" file  which contains the sample-specific source, plate, randomized sample identifier, ADH1B and FABP4 quantifications, RD call, and true RD status. 

### 1.3 Results

We generate a plot of the ADH1B versus FABP4 PCR values indicating samples with RD by red points. The solid vertical line corresponds to the final threshold used to identify samples with high FABP4 (top 25%). The dashed line shows the partitioning of samples using both FABP4 and ADH1B.

## 2 Loading Data

First, we load the deidentified PCR results.

```{r loadRawPCRData}
load(file.path("RDataObjects","PCRResults.RData"))
```

## 3 Plotting PCR Data

We plot the PCR measurements for ADH1B against FABP4 indicating RD in red.

```{r plotPCR}
pdf(file="plottingRD.pdf",paper="USr")
rd <- factor(PCRResults$RDStatus)
plot(PCRResults$FABP4,PCRResults$ADH1B,pch=21,bg=c("grey","red")[rd],xlab="FABP4",ylab="ADH1B")
legend("topleft",c("RD","no RD"),pch=19,col=c("red","grey"),bty="n",cex=0.8)
abline(v=-20.05)
abline(a=-39.5,b=-1,lty=2)
dev.off()
```

## 4 Appendix

```{r getLocation}
getwd()
```

```{r sessionInfo}
sessionInfo()
```