X-Git-Url: https://git.auder.net/css/rpsls.css?a=blobdiff_plain;f=test%2Fscript_data.R;fp=test%2Fscript_data.R;h=102961c7baf7fdf3a33ceb29aa4316506356bb09;hb=17b9fa5f6004bb55e915d8916f1c9a1f128a65ce;hp=0000000000000000000000000000000000000000;hpb=923ed737d0fa335b858204b813c964432488abbe;p=valse.git diff --git a/test/script_data.R b/test/script_data.R new file mode 100644 index 0000000..102961c --- /dev/null +++ b/test/script_data.R @@ -0,0 +1,15 @@ +m=6 +p=6 + +covY = array(0,dim = c(m,m,2)) +covY[,,1] = diag(m) +covY[,,2] = diag(m) + +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(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])