diff --git a/figure/identif/sigma_0.05/alpha_0.05.png b/figure/identif/sigma_0.05/alpha_0.05.png new file mode 100644 index 0000000000000000000000000000000000000000..d5dcd990061581e7f59f5bf30f7c1d7bb55bf14f Binary files /dev/null and b/figure/identif/sigma_0.05/alpha_0.05.png differ diff --git a/figure/identif/sigma_0.05/mu_0.05.png b/figure/identif/sigma_0.05/mu_0.05.png new file mode 100644 index 0000000000000000000000000000000000000000..7849d9b7c297700bc5b066db2b78188ad8327c09 Binary files /dev/null and b/figure/identif/sigma_0.05/mu_0.05.png differ diff --git a/figure/identif/sigma_0.05/oalpha_0.05.png b/figure/identif/sigma_0.05/oalpha_0.05.png new file mode 100644 index 0000000000000000000000000000000000000000..c1b394cd34431dd86fdb6c455fc6eb560051a692 Binary files /dev/null and b/figure/identif/sigma_0.05/oalpha_0.05.png differ diff --git a/figure/identif/sigma_0.05/omu_0.05.png b/figure/identif/sigma_0.05/omu_0.05.png new file mode 100644 index 0000000000000000000000000000000000000000..8b3d7a74cae61b3fe5aa2b6b06339cd3076c6103 Binary files /dev/null and b/figure/identif/sigma_0.05/omu_0.05.png differ diff --git a/figure/identif/sigma_0.542/alpha_0.542.png b/figure/identif/sigma_0.542/alpha_0.542.png new file mode 100644 index 0000000000000000000000000000000000000000..0f4f98dc09df6c7edfbf44731a640103ab159f87 Binary files /dev/null and b/figure/identif/sigma_0.542/alpha_0.542.png differ diff --git a/figure/identif/sigma_0.542/mu_0.542.png b/figure/identif/sigma_0.542/mu_0.542.png new file mode 100644 index 0000000000000000000000000000000000000000..e18968c6a723d38507f139b8159149c241b313e7 Binary files /dev/null and b/figure/identif/sigma_0.542/mu_0.542.png differ diff --git a/figure/identif/sigma_0.542/oalpha_0.542.png b/figure/identif/sigma_0.542/oalpha_0.542.png new file mode 100644 index 0000000000000000000000000000000000000000..7c3041d93f7776f2765389d1d48e88dbc0f7aff3 Binary files /dev/null and b/figure/identif/sigma_0.542/oalpha_0.542.png differ diff --git a/figure/identif/sigma_0.542/omu_0.542.png b/figure/identif/sigma_0.542/omu_0.542.png new file mode 100644 index 0000000000000000000000000000000000000000..8a3c7ebbb37ea0cbc241d55575bd4e602ad21d79 Binary files /dev/null and b/figure/identif/sigma_0.542/omu_0.542.png differ diff --git a/src/figure.R b/src/figure.R new file mode 100644 index 0000000000000000000000000000000000000000..a6be3208a6d708a03d54ed15c92f6c8d2947e057 --- /dev/null +++ b/src/figure.R @@ -0,0 +1,47 @@ +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) + + + + + + + + + + + + + + + + + + diff --git a/src/identifiabilite.R b/src/identifiabilite.R index 1b0650004e6dea5d63d8b65a7fa44fa4c868deef..3eb65361a9b29db4ba8505d895402e07200e8ff1 100644 --- a/src/identifiabilite.R +++ b/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