## CD-1 PFOA/GenX gestational exposure and latent health outcomes study
## Week 6 necropsy data
## Body weight, liver weight, relative liver weight
## HAC 01/22/2020 (w/ QC data)

## Load packages (needs to be done every time R is opened)
library(dplyr)
library(lme4)
library(ggplot2)
library(multcomp)
library(tidyverse)

# Set your working directory (replace this example with your file folder destination)
setwd("/Users/harliecope/Desktop/GenX:PFOA Study/GenX PFOA Mouse Study/GenXPFOA Study")

# Read in your files
x<- read_csv("~/Desktop/GenX:PFOA Study/GenX PFOA Mouse Study/QC Data Files/csv for R/QCWeek18NecropsyWeights.csv")
x$ID <-as.factor(x$ID)
x$Dam.ID<-as.factor(x$Dam.ID)
x$Group<-as.factor(x$Group)

str(x)

w18necrop<-x


names(w18necrop)

w18necrop.summary.stats<-w18necrop%>%
  group_by(Group,Sex,Diet)%>%
  summarise(
    bw.mean=mean(BW, na.rm=T),
    bw.sd=sd(BW, na.rm=T),
    liv.mean=mean(LW, na.rm=T),
    liv.sd=sd(LW, na.rm=T),
    rel.liv.mean=mean(RelLivWt, na.rm=T),
    rel.liv.sd=sd(RelLivWt, na.rm=T))
w18necrop.summary.stats


write.csv(w18necrop.summary.stats, "Week18Necropsy Weights Summary Statistics.csv",row.names = F)


### Statistical Analyses

## Statistical analysis
## Males

w18m<-subset(w18necrop, Sex=="M")
w18f<-subset(w18necrop, Sex=="F")

w18mhf<-subset(w18m, Diet=="H")
w18mlf<-subset(w18m, Diet=="L")
w18fhf<-subset(w18f, Diet=="H")
w18flf<-subset(w18f, Diet=="L")

w18fveh<-subset(w18f, Group=="Control")
w18mveh<-subset(w18m, Group=="Control")
## Male, HFD
summary(glht(aov(BW~Group,data=w18mhf), linfct=mcp(Group="Dunnett")))
## No differences

summary(glht(aov(LW~Group,data=w18mhf), linfct=mcp(Group="Dunnett")))
## No differences

summary(glht(aov(RelLivWt~Group,data=w18mhf), linfct=mcp(Group="Dunnett")))
## No differences


## Male, LFD
summary(glht(aov(BW~Group,data=w18mlf), linfct=mcp(Group="Dunnett")))
## Difference in 2.0 mg/kg GenX and control (p=0.0261) and 0.1 mg/kg PFOA and control (p=0.0456) (1.0 mg/kg PFOA 0.0788)

summary(glht(aov(LW~Group,data=w18mlf), linfct=mcp(Group="Dunnett")))
## Difference in 2.0 mg/kg GenX compared to control (p=0.0160) and 1.0 mg/kg PFOA and control (p=0.0189)

summary(glht(aov(RelLivWt~Group,data=w18mlf), linfct=mcp(Group="Dunnett")))
## No differences


## Female, HFD
summary(glht(aov(BW~Group,data=w18fhf), linfct=mcp(Group="Dunnett")))
## No differences

summary(glht(aov(LW~Group,data=w18fhf), linfct=mcp(Group="Dunnett")))
## No differences (1.0 mg/kg PFOA compared to control p=0.0567)

summary(glht(aov(RelLivWt~Group,data=w18fhf), linfct=mcp(Group="Dunnett")))
## Difference in 1.0 mg/kg PFOA compared to control, (p=0.0111)


## Female, LFD
summary(glht(aov(BW~Group,data=w18flf), linfct=mcp(Group="Dunnett")))
## No differences

summary(glht(aov(LW~Group,data=w18flf), linfct=mcp(Group="Dunnett")))
## No differences

summary(glht(aov(RelLivWt~Group,data=w18flf), linfct=mcp(Group="Dunnett")))
## No differences

summary(glht(aov(BW~Diet, data=w18fveh)))
summary(glht(aov(BW~Diet, data=w18mveh)))
summary(glht(aov(LW~Diet, data=w18fveh)))
summary(glht(aov(LW~Diet, data=w18mveh)))
summary(glht(aov(RelLivWt~Diet, data=w18fveh)))
summary(glht(aov(RelLivWt~Diet, data=w18mveh)))