diff --git a/data/nasaTLX.csv b/data/nasaTLX.csv
index b674f14..c5b5362 100644
--- a/data/nasaTLX.csv
+++ b/data/nasaTLX.csv
@@ -25,4 +25,10 @@ qwerty, 5, 2, 15, 6, 3
 circle, 13, 4, 7, 12, 13
 dvorak, 9, 11, 14, 14, 14
 qwerty, 0, 0, 5, 7, 2
-circle, 15, 15, 11, 15, 11
\ No newline at end of file
+circle, 15, 15, 11, 15, 11
+dvorak, 16, 14, 6, 16, 15
+qwerty, 2, 3, 8, 5, 4
+circle, 14, 12, 8, 10, 12
+dvorak, 14, 13, 9, 16, 19
+qwerty, 3, 4, 5, 7, 3
+circle, 17, 14, 7, 12, 13
\ No newline at end of file
diff --git a/data/results.csv b/data/results.csv
index 001e4ba..ea9e753 100644
--- a/data/results.csv
+++ b/data/results.csv
@@ -8,3 +8,5 @@ ab, 21, circle-dvorak-qwerty, 12.1, 0.034, 5.4, 0.029, 5.8, 0.105
 mz, 24, qwerty-dvorak-circle, 15.8, 0.063, 9.8, 0.043, 10.3, 0.065
 oa, 21, qwerty-circle-dvorak, 17.5, 0.011, 9.5, 0.011, 11.6, 0.010
 dc, 22, dvorak-qwerty-circle, 15.2, 0.033, 7.0, 0.097, 5.1, 0.009
+pt, 23, dvorak-circle-qwerty, 16.3, 0.027, 6.5, 0.012, 6.2, 0.005
+rn, 24, circle-qwerty-dvorak, 14.8, 0.052, 8.8, 0.034, 9.5, 0.022
diff --git a/doc/report.Rnw b/doc/report.Rnw
index cfc5b16..7a986aa 100644
--- a/doc/report.Rnw
+++ b/doc/report.Rnw
@@ -163,6 +163,80 @@ dev.off()
 \end{figure}
 
 \subsection{Inferential Statistics}\label{inferential-statistics}
+Independent var: 
+- QWERTY
+- DVORAK
+- Circle
+
+Dependent var:
+- WPM
+- TER
+- Nasa-TLX
+
+%Anova RM for WPM
+<<echo=FALSE>>=
+library(tidyr)
+library(dplyr)
+
+# Add participant ID
+results$id <- 1:nrow(results)
+
+# --- WPM Long Format ---
+wpm_long <- results %>%
+  select(id, qwerty_wpm, dvorak_wpm, circle_wpm) %>%
+  pivot_longer(
+    cols = -id,
+    names_to = "layout",
+    values_to = "wpm"
+  )
+
+wpm_long$id <- factor(wpm_long$id)
+
+wpm_long$layout <- factor(wpm_long$layout,
+                         levels=c("qwerty_wpm","dvorak_wpm","circle_wpm"),
+                         labels=c("QWERTY","DVORAK","CIRCLE"))
+
+# --- RM ANOVA for WPM ---
+anova_wpm <- aov(wpm ~ layout + Error(id/layout), data=wpm_long)
+
+# Print ANOVA table
+summary(anova_wpm)
+@
+
+%Anova RM for TER
+<<echo=FALSE>>=
+
+# --- TER Long Format ---
+ter_long <- results %>%
+  select(id, qwerty_ter, dvorak_ter, circle_ter) %>%
+  pivot_longer(
+    cols = -id,
+    names_to = "layout",
+    values_to = "ter"
+  )
+
+ter_long$id <- factor(ter_long$id)
+
+ter_long$layout <- factor(ter_long$layout,
+                         levels=c("qwerty_ter","dvorak_ter","circle_ter"),
+                         labels=c("QWERTY","DVORAK","CIRCLE"))
+
+# --- RM ANOVA for TER ---
+anova_ter <- aov(ter ~ layout + Error(id/layout), data=ter_long)
+
+summary(anova_ter)
+@
+
+% Post-Hoc analysis with bonferroni correction for WPM
+<<echo=FALSE>>=
+suppressMessages(library(emmeans))
+
+suppressMessages(emm_wpm <- emmeans(anova_wpm, ~ layout))
+
+posthoc <- pairs(emm_wpm, adjust = "bonferroni")
+
+print(posthoc)
+@
 
 \subsubsection{Objective Measures}\label{objective-measures-1}
 
diff --git a/figures/nasa_boxplots.pdf b/figures/nasa_boxplots.pdf
index 87c1cc5..1e02a5b 100644
Binary files a/figures/nasa_boxplots.pdf and b/figures/nasa_boxplots.pdf differ
diff --git a/figures/ter_plot.pdf b/figures/ter_plot.pdf
index 45f4099..a6e0804 100644
Binary files a/figures/ter_plot.pdf and b/figures/ter_plot.pdf differ
diff --git a/figures/wpm_plot.pdf b/figures/wpm_plot.pdf
index 98daf20..e15ec23 100644
Binary files a/figures/wpm_plot.pdf and b/figures/wpm_plot.pdf differ