Changeset 14855c0 in tol

Timestamp:

May 13, 2010, 7:18:59 AM (16 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

9421997

Parents:

97a3b1b

Message:

Comparación de resultados con los datos reales

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@2164 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol_pkg/BysVecLinReg/test/test_0003/test.tol (modified) (2 diffs)

tol_pkg/BysVecLinReg/test/test_0003/test.tol

@@ -43 +43 @@
 
 VMatrix avg = Tra(Constant(1,sampleLength,1/sampleLength)*mcmc);
-VMatrix vec.A = BysVecLinReg::vmatrix2vec(
+VMatrix vec.A_ = BysVecLinReg::vmatrix2vec(
   SubRow(avg,Range(1,MV::_.m * MV::_.d,1)));
-VMatrix vec.V = BysVecLinReg::vmatrix2vec(
+VMatrix vec.V_ = BysVecLinReg::vmatrix2vec(
   SubRow(avg,Range(MV::_.m * MV::_.d+1,VRows(avg),1)));
 
-VMatrix A_ = BysVecLinReg::vec2vmatrix(vec.A,MV::_.A);
-VMatrix V_ = BysVecLinReg::vec2vmatrix(vec.V,MV::_.V);
+VMatrix A_ = BysVecLinReg::vec2vmatrix(vec.A_,MV::_.A);
+VMatrix V_ = BysVecLinReg::vec2vmatrix(vec.V_,MV::_.V);
+//Instant residuals correlation
+VMatrix V.COR_ = Mat2VMat(NormDiag(VMat2Mat(V_)));
+
 VMatrix E_ = Y - A_*X;
 VMatrix E_t = Tra(E_);
 VMatrix E01_ = CholeskiSolve(V.L,E_,"PtLLtP");
+
+VMatrix S2_ = MMtSqr(E_);
+VMatrix S2.COR_ = Mat2VMat(NormDiag(VMat2Mat(S2_)));
+
+VMatrix cmpA = BysVecLinReg::vmatrix2vec(A) |
+               BysVecLinReg::vmatrix2vec(A_); 
+VMatrix cmpV = BysVecLinReg::vmatrix2vec(V) |
+               BysVecLinReg::vmatrix2vec(V_); 
 VMatrix cmpE = BysVecLinReg::vmatrix2vec(E) |
                BysVecLinReg::vmatrix2vec(E_); 
@@ -60 +71 @@
 Real stdErr  = VMatAvr(E01 ^2);
 Real stdErr_ = VMatAvr(E01_^2);
+Real quality.adjust = Min(1,Exp((stdErr-stdErr_)*N*d));
 
-Real quality.adjust = Min(1,Exp((stdErr-stdErr_)*N*d));
 
 VMatrix C.dist = C*vec.A - c;