From 5c6aa44be3904df994460721589a9f067cd00a46 Mon Sep 17 00:00:00 2001
From: Matias-Lopez-13 <112888781+Matias-Lopez-13@users.noreply.github.com>
Date: Fri, 21 Jun 2024 13:40:08 -0300
Subject: [PATCH] Create plot_p

---
 R/plot_p | 36 ++++++++++++++++++++++++++++++++++++
 1 file changed, 36 insertions(+)
 create mode 100644 R/plot_p

diff --git a/R/plot_p b/R/plot_p
new file mode 100644
index 0000000..c4b83d7
--- /dev/null
+++ b/R/plot_p
@@ -0,0 +1,36 @@
+plot_p <- function(obs_support,total_obs){
+  ## Test to make sure that obs_support is less than or equal to total_obs
+  stopifnot("The number of observations in favor of the working hypothesis must be less than or equal to the total number of observations"=obs_support<=total_obs)
+  obs_oppose <- obs_support+1
+  stopifnot("Observations are already compatible with the null. The number of observations in favor of the working hypothesis must be greater than or equal to half of the total number of observations"=obs_support >= (total_obs/2))
+  ## We assume odds=1 here
+ thep <- dFNCHypergeo(x=obs_support, m1 = obs_support, m2 = obs_oppose,
+                      n = total_obs, odds = 1)
+
+
+
+k_values <- 0:total_obs  # Possible number of observations favoring working theory
+
+# Probability of getting each k value in a sample of observations from null
+probabilities <- dhyper(k_values, m = obs_support, n = obs_oppose, k = total_obs)
+
+
+null <- data.frame(k_values, probabilities)
+
+plot_null<-ggplot(null, aes(x = factor(k_values), y = probabilities, 
+               fill = factor(k_values == obs_support))) +
+  geom_bar(stat = "identity", color = "black") +
+  scale_fill_manual(values = c("skyblue", "red"), guide = FALSE)+
+  labs(title = "Null distribution",
+       x = "Observations favoring working theory",
+       y = "Frequency of occurrence")+ 
+  scale_y_continuous(labels = scales::percent_format(accuracy = 1),
+                     breaks =c(0.05,0.1))+
+  theme_classic()+
+    geom_hline(yintercept=0.05,linetype="dashed")+
+  geom_hline(yintercept=0.1,linetype="dotted")
+
+  return(plot_null)
+
+}
+