Skip to content
Snippets Groups Projects
Commit c7c785e8 authored by Xianzong Meng's avatar Xianzong Meng
Browse files

ANOVA-based statistic

parent c22d4b19
Branches main
No related tags found
No related merge requests found
ANOVA.Rmd 0 → 100644
---
title: "ANOVA"
author: "Alex_Meng"
date: "9/3/2022"
output: html_document
---
#loading packages
```{r}
library(dplyr)
library(tidyverse)
library(reshape2)
```
#import data
```{r}
#EPM
EPM_open_arms_duration <- read_csv('assets/tables/EPM_open_arms_duration.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
EPM_entry_closed_arms <- read_csv('assets/tables/EPM_entry_closed_arms.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
EPM_first_time_entry <- read_csv('assets/tables/EPM_first_time_entry.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
EPM_total_distance <- read_csv('assets/tables/EPM_total_distance.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
#Social contact
SI_total_no_contact <- read_csv('assets/tables/SI_total_no_contact.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
SI_total_mounting <- read_csv('assets/tables/SI_total_mounting.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
SI_total_aggressiveness <- read_csv('assets/tables/SI_total_aggressiveness.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
SI_fighting <- read_csv('assets/tables/SI_fighting.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
#PCR
PCR_CA <- read_csv('assets/tables/PCR_CA.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_dCA1 <- read_csv('assets/tables/PCR_dCA1.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_dCA3 <- read_csv('assets/tables/PCR_dCA3.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_dGRDG <- read_csv('assets/tables/PCR_dGRDG.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_DR <- read_csv('assets/tables/PCR_DR.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_IF <- read_csv('assets/tables/PCR_IF.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_PRL <- read_csv('assets/tables/PCR_PRL.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_PVN <- read_csv('assets/tables/PCR_PVN.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_vCA1 <- read_csv('assets/tables/PCR_vCA1.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_vCA3 <- read_csv('assets/tables/PCR_vCA3.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
PCR_vGRDG <- read_csv('assets/tables/PCR_vGRDG.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
#Oxytocin protein
OXY_CA <- read_csv('assets/tables/OXY_CA.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
OXY_DR <- read_csv('assets/tables/OXY_DR.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
OXY_HIP <- read_csv('assets/tables/OXY_HIP.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
OXY_MFC <- read_csv('assets/tables/OXY_MFC.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
OXY_PVN <- read_csv('assets/tables/OXY_PVN.csv', col_types = cols()) %>% melt() %>% dplyr::rename(group = variable) %>% drop_na()
```
#ANOVA_EPM
```{r}
ano_EPM_open_arms_duration <- aov(value ~ group, data = EPM_open_arms_duration)
pos <-TukeyHSD(ano_EPM_open_arms_duration)
pos
ano_EPM_open_arms_duration <- summary.aov(ano_EPM_open_arms_duration)
ano_EPM_open_arms_duration
print("-------------------------------------------------------------------------")
ano_EPM_entry_closed_arms <- aov(value ~ group, data = EPM_entry_closed_arms)
pos <-TukeyHSD(ano_EPM_entry_closed_arms)
pos
ano_EPM_entry_closed_arms <- summary.aov(ano_EPM_entry_closed_arms)
ano_EPM_entry_closed_arms
print("-------------------------------------------------------------------------")
ano_EPM_first_time_entry <- aov(value ~ group, data = EPM_first_time_entry)
pos <-TukeyHSD(ano_EPM_first_time_entry)
pos
ano_EPM_first_time_entry <- summary.aov(ano_EPM_first_time_entry)
ano_EPM_first_time_entry
print("-------------------------------------------------------------------------")
ano_EPM_total_distance <- aov(value ~ group, data = EPM_total_distance)
pos <-TukeyHSD(ano_EPM_total_distance)
pos
ano_EPM_total_distance <- summary.aov(ano_EPM_total_distance)
ano_EPM_total_distance
```
#ANOVA Social contact
```{r}
ano_SI_total_no_contact <- aov(value ~ group, data = SI_total_no_contact)
pos <-TukeyHSD(ano_SI_total_no_contact)
pos
ano_SI_total_no_contact <- summary.aov(ano_SI_total_no_contact)
ano_SI_total_no_contact
print("-------------------------------------------------------------------------")
ano_SI_total_mounting <- aov(value ~ group, data = SI_total_mounting)
pos <-TukeyHSD(ano_SI_total_mounting)
pos
ano_SI_total_mounting <- summary.aov(ano_SI_total_mounting)
ano_SI_total_mounting
print("-------------------------------------------------------------------------")
ano_SI_total_aggressiveness <- aov(value ~ group, data = SI_total_aggressiveness)
pos <-TukeyHSD(ano_SI_total_aggressiveness)
pos
ano_SI_total_aggressiveness <- summary.aov(ano_SI_total_aggressiveness)
ano_SI_total_aggressiveness
print("-------------------------------------------------------------------------")
ano_SI_fighgting <- aov(value ~ group, data = SI_fighting)
pos <-TukeyHSD(ano_SI_fighgting)
pos
ano_SI_fighgting <- summary.aov(ano_SI_fighgting)
ano_SI_fighgting
```
#ANOVA PCR
```{r}
ano_PCR_CA <- aov(value ~ group, data = PCR_CA)
pos <-TukeyHSD(ano_PCR_CA)
pos
ano_PCR_CA <- summary.aov(ano_PCR_CA)
ano_PCR_CA
print("-------------------------------------------------------------------------")
ano_PCR_dCA1 <- aov(value ~ group, data = PCR_dCA1)
pos <-TukeyHSD(ano_PCR_dCA1)
pos
ano_PCR_dCA1 <- summary.aov(ano_PCR_dCA1)
ano_PCR_dCA1
print("-------------------------------------------------------------------------")
ano_PCR_dCA3 <- aov(value ~ group, data = PCR_dCA3)
pos <-TukeyHSD(ano_PCR_dCA3)
pos
ano_PCR_dCA3 <- summary.aov(ano_PCR_dCA3)
ano_PCR_dCA3
print("-------------------------------------------------------------------------")
ano_PCR_dGRDG <- aov(value ~ group, data = PCR_dGRDG)
pos <-TukeyHSD(ano_PCR_dGRDG)
pos
ano_PCR_dGRDG <- summary.aov(ano_PCR_dGRDG)
ano_PCR_dGRDG
print("-------------------------------------------------------------------------")
ano_PCR_DR <- aov(value ~ group, data = PCR_DR)
pos <-TukeyHSD(ano_PCR_DR)
pos
ano_PCR_DR <- summary.aov(ano_PCR_DR)
ano_PCR_DR
print("-------------------------------------------------------------------------")
ano_PCR_IF <- aov(value ~ group, data = PCR_IF)
pos <-TukeyHSD(ano_PCR_IF)
pos
ano_PCR_IF <- summary.aov(ano_PCR_IF)
ano_PCR_IF
print("-------------------------------------------------------------------------")
ano_PCR_PRL <- aov(value ~ group, data = PCR_PRL)
pos <-TukeyHSD(ano_PCR_PRL)
pos
ano_PCR_PRL <- summary.aov(ano_PCR_PRL)
ano_PCR_PRL
print("-------------------------------------------------------------------------")
ano_PCR_PVN <- aov(value ~ group, data = PCR_PVN)
pos <-TukeyHSD(ano_PCR_PVN)
pos
ano_PCR_PVN <- summary.aov(ano_PCR_PVN)
ano_PCR_PVN
print("-------------------------------------------------------------------------")
ano_PCR_vCA1 <- aov(value ~ group, data = PCR_vCA1)
pos <-TukeyHSD(ano_PCR_vCA1)
pos
ano_PCR_vCA1 <- summary.aov(ano_PCR_vCA1)
ano_PCR_vCA1
print("-------------------------------------------------------------------------")
ano_PCR_vCA3 <- aov(value ~ group, data = PCR_vCA3)
pos <-TukeyHSD(ano_PCR_vCA3)
pos
ano_PCR_vCA3 <- summary.aov(ano_PCR_vCA3)
ano_PCR_vCA3
print("-------------------------------------------------------------------------")
ano_PCR_vGRDG <- aov(value ~ group, data = PCR_vGRDG)
pos <-TukeyHSD(ano_PCR_vGRDG)
pos
ano_PCR_vGRDG <- summary.aov(ano_PCR_vGRDG)
ano_PCR_vGRDG
```
#Oxytocin protein
```{r}
ano_OXY_CA <- aov(value ~ group, data = OXY_CA)
pos <-TukeyHSD(ano_OXY_CA)
pos
ano_OXY_CA <- summary.aov(ano_OXY_CA)
ano_OXY_CA
print("-------------------------------------------------------------------------")
ano_OXY_DR <- aov(value ~ group, data = OXY_DR)
pos <-TukeyHSD(ano_OXY_DR)
pos
ano_OXY_DR <- summary.aov(ano_OXY_DR)
ano_OXY_DR
print("-------------------------------------------------------------------------")
ano_OXY_HIP <- aov(value ~ group, data = OXY_HIP)
pos <-TukeyHSD(ano_OXY_HIP)
pos
ano_OXY_HIP <- summary.aov(ano_OXY_HIP)
ano_OXY_HIP
print("-------------------------------------------------------------------------")
ano_OXY_MFC <- aov(value ~ group, data = OXY_MFC)
pos <-TukeyHSD(ano_OXY_MFC)
pos
ano_OXY_MFC <- summary.aov(ano_OXY_MFC)
ano_OXY_MFC
print("-------------------------------------------------------------------------")
ano_OXY_PVN <- aov(value ~ group, data = OXY_PVN)
pos <-TukeyHSD(ano_OXY_PVN)
pos
ano_OXY_PVN <- summary.aov(ano_OXY_PVN)
ano_OXY_PVN
```
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment