library(data.table)

bin_sims <- fread("/home/ryan/Desktop/julia_stuff/raksha_seeds/binomial_sims.csv")

bin_sims[,abs_error := abs(value - prop_germ)]

library(tidyverse)

bin_sims[abs_error>0.05, ]


mae <-  bin_sims[,list(MAE=median(abs_error)), by = .(prop_germ, n_seeds)]

mae[n_seeds ==5,] |> 
  ggplot(aes(prop_germ, n_seeds, color = MAE)) + geom_point()

bin_sum
bin_sims[n_seeds ==5,] |> 
  ggplot(aes(prop_germ, value)) + geom_jitter(alpha = 0.02) + 
  xlab("true germination") +
  ylab("predicted proportion germinated") + ggtitle("germinating 5 seeds")


bin_sims[n_seeds ==50,] |> 
  ggplot(aes(prop_germ, value)) + geom_jitter(alpha = 0.02) + 
  xlab("true germination") +
  ylab("predicted proportion germinated") + ggtitle("germinating 50 seeds")

bin_sims[n_seeds ==100,] |> 
  ggplot(aes(prop_germ, value)) + geom_jitter(alpha = 0.02) + 
  xlab("true germination") +
  ylab("predicted proportion germinated") + ggtitle("germinating 100 seeds")


bin_sims[n_seeds ==10,] |> 
  ggplot(aes(prop_germ, value)) + geom_jitter(alpha = 0.02)


bin_sims[n_seeds ==25,] |> 
  ggplot(aes(prop_germ, value)) + geom_jitter(alpha = 0.02)


bin_sims[n_seeds ==50,] |> 
  ggplot(aes(prop_germ, value)) + geom_jitter(alpha = 0.02)
