From 8795d8557c76862344880f91f09ccc61c4bf770c Mon Sep 17 00:00:00 2001
From: unknown <isabelle.sanchez@inra.fr>
Date: Wed, 11 May 2022 09:32:47 +0200
Subject: [PATCH] update KBO_EM: ppnew=max(sum(),0.5)
---
R/utils_functions.R | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/R/utils_functions.R b/R/utils_functions.R
index 9f59033..8d7d317 100644
--- a/R/utils_functions.R
+++ b/R/utils_functions.R
@@ -616,9 +616,9 @@ KBO_EM<-function(param,kappaOpt,Y,Tps,N,dix){
C= sum((p%*%(matrice_C*Z))[1:(N-1)])
Yb= sum((p%*%(matrice_Yb*Z))[1:(N-1)])
- newm0= (C*Yb-B*Ya)/(C^2-A*B)
- newmm= (C*Ya-A*Yb)/(C^2-A*B)
- ppnew=sum(p%*%Z)/N
+ newm0=(C*Yb-B*Ya)/(C^2-A*B)
+ newmm=(C*Ya-A*Yb)/(C^2-A*B)
+ ppnew=max(sum(p%*%Z)/N , 0.5)
# Outputs
resultat=list("m0"=newm0,
--
GitLab