P29.R

# making table data sets
library(dplyr)
library(tidyr)
library(MorpheusData)

#############benchmark 29
attributes <- read.table(text=
"Species	Attribute
Callvulg        MI
Empenigr        MI
Rhodtome        PI
", header=T)

dat <- read.table(text=
"ID  Species Value
1 Callvulg  0.55
2 Callvulg  0.67
3 Callvulg  0.10
1 Empenigr 11.13
2 Empenigr  0.17
3 Empenigr  1.55
1 Rhodtome  0.17
2 Rhodtome  1.55
3 Rhodtome  3.00",
header=T)


write.csv(dat, "data-raw/p29_input1.csv", row.names=FALSE)
write.csv(attributes, "data-raw/p29_input2.csv", row.names=FALSE)

df_out = dat %>% 
  inner_join(attributes) %>% 
  filter(Attribute == "MI") %>% 
  group_by(ID) %>% 
  summarise(Total = sum(Value))

write.csv(df_out, "data-raw/p29_output1.csv", row.names=FALSE)

p29_output1 <- read.csv("data-raw/p29_output1.csv", check.names = FALSE)
fctr.cols <- sapply(p29_output1, is.factor)
int.cols <- sapply(p29_output1, is.integer)

p29_output1[, fctr.cols] <- sapply(p29_output1[, fctr.cols], as.character)
p29_output1[, int.cols] <- sapply(p29_output1[, int.cols], as.numeric)
save(p29_output1, file = "data/p29_output1.rdata")

p29_input1 <- read.csv("data-raw/p29_input1.csv", check.names = FALSE)
fctr.cols <- sapply(p29_input1, is.factor)
int.cols <- sapply(p29_input1, is.integer)

p29_input1[, fctr.cols] <- sapply(p29_input1[, fctr.cols], as.character)
p29_input1[, int.cols] <- sapply(p29_input1[, int.cols], as.numeric)
save(p29_input1, file = "data/p29_input1.rdata")

p29_input2 <- read.csv("data-raw/p29_input2.csv", check.names = FALSE)
fctr.cols <- sapply(p29_input2, is.factor)
int.cols <- sapply(p29_input2, is.integer)

p29_input2[, fctr.cols] <- sapply(p29_input2[, fctr.cols], as.character)
p29_input2[, int.cols] <- sapply(p29_input2[, int.cols], as.numeric)
save(p29_input2, file = "data/p29_input2.rdata")
fredfeng/MorpheusData documentation built on May 16, 2019, 2:42 p.m.