fix the likelihood computation. fix some other few things

[valse.git] / pkg / data / script_data.R

diff --git a/pkg/data/script_data.R b/pkg/data/script_data.R
index d425a38..7479674 100644 (file)
--- a/pkg/data/script_data.R
+++ b/pkg/data/script_data.R
@@ -1,5 +1,5 @@
-m=11
-p=10
+m=6
+p=6
 
 covY = array(0,dim = c(m,m,2))
 covY[,,1] = diag(m)
@@ -9,7 +9,7 @@ Beta = array(0, dim = c(p, m, 2))
 Beta[1:4,1:4,1] = 3*diag(4)
 Beta[1:4,1:4,2] = -2*diag(4)
 
-#Data = generateXY(100, c(0.5,0.5), rep(0,p), Beta, diag(p), covY)
-
-#Res = valse(Data$X,Data$Y, fast=FALSE, plot=FALSE, verbose = TRUE, kmax=2, size_coll_mod = 100,
-#            selecMod = "BIC")
+#Data = generateXY(200, c(0.5,0.5), rep(0,p), Beta, diag(p), covY)
+#  
+#Res = valse(Data$X,Data$Y, fast=FALSE, plot=FALSE, verbose = TRUE, kmax=3, size_coll_mod = 50, selecMod = "DDSE", mini = 50, maxi=100)
+#plot(Res$tableau[,3], -Res$tableau[,4])