featuring : Ajout du script pour faire les boxplots

figure : Ajout des boxplots des estimations des paramètres

src/figure.R

+distrib_alpha <-  ggplot(data = data.frame(alpha) ) + 
+  geom_boxplot(aes(y = log(alpha))) + 
+  geom_hline(yintercept = log(true_v[1])) +
+  labs(title="Distribution de log(alpha)" , subtitle = "sigma = 0.542") +
+  theme(plot.title = element_text(hjust = 0.5) , plot.subtitle = element_text(hjust = 0.5)) +
+  ylim(-15, 0)
+
+
+distrib_mu <-  ggplot(data = data.frame(mu) ) + 
+  geom_boxplot(aes(y = log(mu))) +
+  geom_hline(yintercept = log(true_v[2])) +
+  labs(title="Distribution de log(mu)" , subtitle = "sigma = 0.542") +
+  theme(plot.title = element_text(hjust = 0.5) , plot.subtitle = element_text(hjust = 0.5)) +
+  ylim(-32, -20)
+
+distrib_omega_alpha <-  ggplot(data = data.frame(omega_alpha) ) + 
+  geom_boxplot(aes(y = omega_alpha)) +
+  geom_hline(yintercept = true_v[4]**2) +
+  labs(title="Distribution de omega_alpha" , subtitle = "sigma = 0.542") +
+  theme(plot.title = element_text(hjust = 0.5) , plot.subtitle = element_text(hjust = 0.5)) +
+  ylim(0,28)
+
+
+distrib_omega_mu <-  ggplot(data = data.frame(omega_mu)) + 
+  geom_boxplot(aes(y = omega_mu)) + 
+  geom_hline(yintercept = true_v[5]**2) +
+  labs(title="Distribution de omega_mu" , subtitle = "sigma = 0.542") +
+  theme(plot.title = element_text(hjust = 0.5) , plot.subtitle = element_text(hjust = 0.5)) +
+  ylim(0,30)
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+

src/identifiabilite.R

@@ -9,7 +9,7 @@ control <- saemix::saemixControl(map = FALSE,fix.seed = T, fim = FALSE, ll.is =
 
 alpha_pop = exp(-6.34)
 mu_pop = exp(-26.8)
-sigma = 0.542
+sigma = 0.05
 omega_alpha = 3.37
 omega_mu = 3.78