Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- mu1_mu2_violin_onset <- age[which(age$Status=="2" & (age$Diet=="wt" | age_onset$Diet=="mu2" | age_onset$Diet=="mu1")), ]
- onset_mu1_mu2_final=ggplot(data=mu1_mu2_violin_onset, aes(x=Diet, y=Age, fill=Diet)) +
- geom_violin() + stat_summary(fun.data=data_summary) + theme_bw() + theme_classic() +
- scale_color_brewer(palette="Dark2") +scale_fill_brewer(palette="Blues")+ theme(axis.text=element_text(family="sans", size=18),axis.line.x = element_line(color="black"), axis.line.y= element_line(color="black"), legend.position="none", axis.title=element_text(family="sans", size=18), panel.grid.major = element_blank(), panel.grid.minor = element_blank()) +
- scale_y_continuous(name="Age") + scale_x_discrete(name="",limits=c("mu2","mu1","control"),labels=expression(italic(mu2),italic(mu1),WT))
- spag= onset_mu1_mu2_final + coord_flip()
- onset_mu1_mu2=spag + facet_grid(Sex~ ., switch="y") +
- theme(strip.text.y = element_text(size=18, angle=180),
- strip.background = element_rect(colour="white"))
- # T-test
- violin_onset_F_mu1 <- mu1_mu2_violin_onset[which(mu1_mu2_violin_onset$Sex=="F" & (mu1_mu2_violin_onset$Diet=="mu1" | mu1_mu2_violin_onset$Diet=="wt)),]
- t.test(violin_onset_F_mu1$Age~violin_onset_F_mu1$Diet,var.equal=T)$p.value
- violin_onset_M_mu1 <- mu1_mu2_violin_onset[which(mu1_mu2_violin_onset$Sex=="M" & (mu1_mu2_violin_onset$Diet=="mu1" | mu1_mu2_violin_onset$Diet=="wt)),]
- t.test(violin_onset_M_mu1$Age~violin_onset_M_mu1$Diet,var.equal=T)$p.value
- violin_onset_F_mu2 <- mu1_mu2_violin_onset[which(mu1_mu2_violin_onset$Sex=="F" & (mu1_mu2_violin_onset$Diet=="mu2" | mu1_mu2_violin_onset$Diet=="wt)),]
- t.test(violin_onset_F_mu2$Age~violin_onset_F_mu2$Diet,var.equal=T)$p.value
- violin_onset_M_mu2 <- mu1_mu2_violin_onset[which(mu1_mu2_violin_onset$Sex=="M" & (mu1_mu2_violin_onset$Diet=="mu2" | mu1_mu2_violin_onset$Diet=="wt)),]
- t.test(violin_onset_M_mu2$Age~violin_onset_M_mu2$Diet,var.equal=T)$p.value
- onset_mu1_mu2+annotate("text", x=0.90,y=12.4, label="P=5.0x10-9",size=6)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement