########################## library(plyr) library(reshape) library(tidyr) library(dplyr) library(multcomp) ## Need to merge in necropsy data for litter weight setwd("C:/Users/blakebe/Documents/CD1 PFOA GenX/R CSV files") df<-read.csv("Combined necropsy data E11.5 and E17.5.csv") df<-subset(df,pregnant=="Yes") df.11<-subset(df, timepoint=="E11.5") df.17<-subset(df, timepoint=="E17.5") m11<-lm(d.wt.p.11~group+litter.n,data=df.11) summary(m11) confint(m11) coef(m11) summary(glht(lm(d.wt.p.11~group+litter.n,data=df.11), linfct=mcp(group="Tukey"))) m17<-lm(d.wt.pt.17~group+litter.n,data=df.17) summary(m17) confint(m17) coef(m17) summary(glht(lm(d.wt.pt.17~group+litter.n,data=df.17), linfct=mcp(group="Tukey")))