forked from ChalkboardSonata/J2SummerResearch
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path.Rhistory
50 lines (50 loc) · 1.5 KB
/
.Rhistory
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
knitr::opts_chunk$set(echo = TRUE)
library(tidyverse)
library(knitr)
library(readr)
library(dict)
devtools::install_github("mkuhn/dict")
knitr::opts_chunk$set(echo = TRUE)
library(tidyverse)
library(knitr)
library(readr)
library(dict)
devtools::install_github("mkuhn/dict")
knitr::opts_chunk$set(echo = TRUE)
library(tidyverse)
library(knitr)
library(readr)
library(dict)
knitr::opts_chunk$set(echo = TRUE)
library(tidyverse)
install.packages("Rcpp")
knitr::opts_chunk$set(echo = TRUE)
library(tidyverse)
library(knitr)
library(readr)
library(dict)
devtools::install_github("mkuhn/dict")
1
knitr::opts_chunk$set(echo = TRUE)
library(tidyverse)
library(knitr)
library(readr)
library(dict)
library(fs)
# Imprt csv's as one data frame
mydir = "Form_Data/"
files<-dir_ls(mydir)
forms<-files %>% map_dfr(read_csv, col_names = F, col_types = cols(.default = "c"))
#Rename columns
forms <- as_tibble(forms)
forms<-rename(forms, form = X1, invertible = X2, inverse = X3, dinverse = X4, rank = X5,rank_diff = X6, PF = X7, stable = X8)
#Treat rank and pfaffian as numbers
forms<- forms %>% mutate(rank = as.numeric(rank), PF = map_dbl(PF, ~eval(parse(text = .x))))
forms_invert<-forms %>% count(rank, invertible) %>% group_by(rank) %>% mutate(prop = n/ sum(n))
forms_invert %>% filter(invertible =="True") %>%
ggplot(aes( x = rank, y = prop)) +
geom_point()
forms_stable<-forms %>% count(rank, stable) %>% group_by(rank) %>% mutate(prop = n/ sum(n))
forms_stable %>% filter(stable =="True") %>%
ggplot(aes( x = rank, y = prop)) +
geom_point()