-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathindex.Rmd
More file actions
45 lines (30 loc) · 1.42 KB
/
Copy pathindex.Rmd
File metadata and controls
45 lines (30 loc) · 1.42 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
---
title: "Front Door"
author: "Edward Lopatto"
date: "Last updated July 24, 2018"
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```
## Welcome
Welcome to Edward Lopatto's tutorial site. This site was created by Ed Lopatto to document bioinformatic skills he preformed during a MAP in the Summer of 2018. The main purpose of this site is to provide documentation for analysis for both Ed (when he forgets what he did) and Prof. Shannon Hinsa's lab.
This site is a work in progress as Ed slowly crawls on his way to learning github and Rmarkdown.
Any questions can be directed to lopattoe@grinnell.edu
```{r, echo=FALSE, warning=FALSE, message=FALSE}
library(ggplot2)
library(ggtree)
library(treeio)
library(dplyr)
tree <- read.tree("~/desktop/fungene\ str/tree_file")
meta <- read.csv("~/desktop/fungene\ str/strtreedata.csv")
#tree
p1 <- ggtree(tree, layout="circular")
treemeta<- meta[c("ID","Cluster", "Primer", "Number")]
p1 %<+% treemeta+
geom_tiplab(aes(label=Cluster, color=Primer), size=2, offset = 0.02, align=TRUE) +
geom_point(aes(size=Number, color=Primer))+
scale_colour_manual(labels = c("None","strA","strB",""), values = c("black","darkgreen","blue"))+
geom_cladelabel(node=120, label="strB", align=TRUE, color="blue", offset = 0.3, barsize = 0) +
geom_cladelabel(node=152, label="strA", align=TRUE, color = "darkgreen", offset = 0.3, barsize=0)+
theme(legend.position = "right")
```