-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathqPCR_Analysis.Rmd
More file actions
55 lines (40 loc) · 1.45 KB
/
qPCR_Analysis.Rmd
File metadata and controls
55 lines (40 loc) · 1.45 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
---
title: "qPCR_Analysis"
author: "Nathan Brandt"
date: "2/6/2018"
output: html_document
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```
******
##Step 1: Load relevant libraries
```{r}
# This is a function that just makes sure you have a package, or installs it for you without prompting
requireInstall <- function(packageName,isBioconductor=F) {
if ( !try(require(packageName,character.only=T)) ) {
print(paste0("You don't have ",packageName," accessible, ",
"I'm gonna install it"))
if (isBioconductor) {
source("http://bioconductor.org/biocLite.R")
biocLite(packageName)
} else {
install.packages("packageName", repos = "http://cran.us.r-project.org")
}
}
return(1)
}
#Load libraries
requireInstall("tidyverse")
requireInstall("googlesheets")
```
###Step 2: Read in Data
```{r}
path = "/Users/Brandt/Google Drive/GAT1 evolution paper/Additional_Exp/qPCR_Data/"
filename = "MEP2_CNV_Expression_2nd_02122018.txt"
sample_filename = "MEP2_CNV_Expression_2nd_02122018_SampleSheet"
data <- left_join(gs_read(ss = gs_title(sample_filename)), select(read_tsv(paste(path, filename, sep = ""), skip = 1, col_names = TRUE), c("Name","Cp")), by = c("Sample_Number" = "Name"))
#Figure out a way to read files from a folder to get name_list????
#Or might need another master list to control flow in
#Create qPCR Sample Sheet Template
```