From 311459b34e4649b5d873c9c846ccbbc3846b25ba Mon Sep 17 00:00:00 2001
From: ottiP <ottavia.prunas@edu.unito.it>
Date: Wed, 19 Jan 2022 14:55:09 -0500
Subject: [PATCH] plot to make AIC histograms

---
 create_figures_ms/AIC_computation.R | 30 +++++++++++++++++++++++++++++
 1 file changed, 30 insertions(+)
 create mode 100644 create_figures_ms/AIC_computation.R

diff --git a/create_figures_ms/AIC_computation.R b/create_figures_ms/AIC_computation.R
new file mode 100644
index 0000000..ff37f49
--- /dev/null
+++ b/create_figures_ms/AIC_computation.R
@@ -0,0 +1,30 @@
+compute_AIC <- function(ll,p){
+  
+  AIC <- 2*p+2*ll
+  return(AIC)
+  
+}
+
+result1<- readRDS("../Results/result100_Jul_nointerA.rds")
+min_ll1 <- sapply(result1,'[[','minimum')
+parms1 <- sapply(result1,'[[','estimate')
+p1 <- nrow(parms1)
+AIC1 <- compute_AIC(min_ll1,p1)
+
+result2<- readRDS("../Results/result100_Jul_allinter_TestA.rds")
+min_ll2 <- sapply(result2,'[[','minimum')
+parms2 <- sapply(result2,'[[','estimate')
+p2 <- nrow(parms2)
+AIC2 <- compute_AIC(min_ll2,p2)
+
+
+
+
+pdf("hist_12.pdf")
+hist(AIC1-AIC2, col="gray", xlab="AIC difference",ylab="",main=paste("AIC difference between  model 1 and model 2"))
+dev.off()
+
+pdf("hist_13.pdf")
+hist(AIC2-AIC3, col="gray", xlab="AIC difference",ylab="",main=paste("AIC difference between  model 2 and model 3"))
+dev.off()
+