Residual Disease Paper
========================================================
Comparing Ovary and Omentum Expression
-------------------------------------------------------

by Shelley Herbrich

*May 13, 2013*


```{r options, echo=FALSE}

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

```


## 1 Executive Summary

### 1.1 Introduction

In this report, we compare the PCR expression of our target genes between ovary and omental samples. 

We also compare FABP4 exression between ovary and omental samples from patients who had both available.

### 1.2 Data & Methods

In comparing PCR expression of our target genes between ovary and omental samples, we use PCR quantification summaries of samples run on the 17 plates containing tissue from both ovary and omentum. 

### 1.3 Results

We produce density plots of ADH1B and FABP4 expression by tissue. 

We see that for both genes, omental levels are markedly higher than ovary.

We produce a scatterplot of paired FABP4 values for 4 cases with data from both ovary and omentum available.

The measurement from omentum is higher than that from ovary in all 4 cases, with percent differences of 2%, 85%, 110% and 131%.

## 2 Loading PCR Data

```{r libraries}
library(RColorBrewer)
```
```{r loadPCRFiles}
load(file.path("RDataObjects","OMResults.RData"))
load(file.path("RDataObjects","PCRResults.RData"))
```

```{r matchPlates}
#match plates
plate <- unique(OMResults$Plate)
ovFABP4 <- PCRResults$FABP4[which(PCRResults$Plate %in% plate)]
ovADH1B <- PCRResults$ADH1B[which(PCRResults$Plate %in% plate)]
```

## 3 Analyses

### 3.1 Plotting Densities of ADH1B and FABP4 Expression

```{r plotDensities}
#pdf(file="OVOMExpression.pdf",paper="USr")
plot(density(ovADH1B), col="red", xlim=c(-30,-5), main="ADH1B Expression")
lines(density(OMResults$ADH1B),col="grey")
legend("topright",c("Ovary","Omentum"),col=c("red","grey"),lty=1,bty="n")

plot(density(ovFABP4), col="red", xlim=c(-30,-5), main="FABP4 Expression")
lines(density(OMResults$FABP4),col="grey")
legend("topright",c("Ovary","Omentum"),col=c("red","grey"),lty=1,bty="n")
#dev.off()
```

### 3.2 Scatterplot of FABP4 and ADH1B from omentum versus primary

The following scatterplot shows FABP4 values in primary tumor and omentum for the 4 patients in the validation cohort who had both types of tumor tissue assayed.

Solid symbols show cases with RD. The dashed line indicates equality.

```{r plotOmentumPrimary}

pairsOV <- c("M20","M59","M65","M34")
pairsOM <- c("OM1","OM20","OM40","OM3")
sym <- c(16,16,1,16)

ovPaired <- c()
omPaired <- c()

for (i in 1:4)  {
  ovPaired <- c(ovPaired,PCRResults$FABP4[which(PCRResults$Sample.Name == pairsOV[i])])
  omPaired <- c(omPaired,OMResults$FABP4[which(OMResults$Sample.Name == pairsOM[i])])
}

plot(ovPaired, omPaired, xlab="FABP4 Expression from Primary Tumor", ylab="FABP4 Expression from Omental Tumor", xlim=c(-30,-10), ylim=c(-30,-10),pch=sym)

abline(0,1,col="red",lty=2)

ovPaired / omPaired

```

## Appendix

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

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