Différents cas indésirables ont été détectés dans simulate_data() :
dim_x = 0 : Aucun warning. La sortie x est de dimension nulle, ce qui n'empêche pas les autres sorties de prendre des valeurs cohérentes. Je ne sais pas si cela modélise l'absence de variable confondante, ou bien si cela constitue une erreur.
dim_m > 1 ET type_m = "binary" : Une erreur sur les dimensions est renvoyée. Le cas du médiateur binaire multidimensionnel n'a pas encore été implémenté.
type_m = "continuous" ET (sigma_m = 0 OU sigma_m grand) : Un warning est renvoyé. Un NaN est produit sur la dernière sortie, à savoir P(T=1|X,M) = NaN.
n = 0 : Plusieurs warnings sont renvoyés. Les sorties sont des NaN ou des tableaux vides. On pourrait interdire les mauvaises utilisations n <= 0, dim_x <= 0 et dim_m <= 0.
type_m != "binary" : Le branchement sur le cas binaire se fait avec type_m=="binary". On pourrait forcer l'argument "continuous" pour le cas continu, et renvoyer un warning en cas d'argument incorrect. On pourrait utiliser les premières lettres "b" et "c" en guise de raccourci. Ou bien remplacer type_m par un argument booléen m_is_binary.
Les arguments beta_factor représentent les variables oméga de l'article. On pourrait harmoniser les notations.
Différents cas indésirables ont été détectés dans simulate_data() :
dim_x = 0 : Aucun warning. La sortie x est de dimension nulle, ce qui n'empêche pas les autres sorties de prendre des valeurs cohérentes. Je ne sais pas si cela modélise l'absence de variable confondante, ou bien si cela constitue une erreur.
dim_m > 1 ET type_m = "binary" : Une erreur sur les dimensions est renvoyée. Le cas du médiateur binaire multidimensionnel n'a pas encore été implémenté.
type_m = "continuous" ET (sigma_m = 0 OU sigma_m grand) : Un warning est renvoyé. Un NaN est produit sur la dernière sortie, à savoir P(T=1|X,M) = NaN.
n = 0 : Plusieurs warnings sont renvoyés. Les sorties sont des NaN ou des tableaux vides. On pourrait interdire les mauvaises utilisations n <= 0, dim_x <= 0 et dim_m <= 0.
type_m != "binary" : Le branchement sur le cas binaire se fait avec type_m=="binary". On pourrait forcer l'argument "continuous" pour le cas continu, et renvoyer un warning en cas d'argument incorrect. On pourrait utiliser les premières lettres "b" et "c" en guise de raccourci. Ou bien remplacer type_m par un argument booléen m_is_binary.
Les arguments beta_factor représentent les variables oméga de l'article. On pourrait harmoniser les notations.