[Genabel-commits] r1051 - pkg/OmicABEL/src

[noreply at r-forge.r-project.org]

Author: dfabregat
Date: 2012-12-13 12:13:33 +0100 (Thu, 13 Dec 2012)
New Revision: 1051

Modified:
   pkg/OmicABEL/src/REML.c
Log:
[BUG FIX] 3 was used instead of wXL.


Modified: pkg/OmicABEL/src/REML.c
===================================================================
--- pkg/OmicABEL/src/REML.c	2012-12-13 11:08:08 UTC (rev 1050)
+++ pkg/OmicABEL/src/REML.c	2012-12-13 11:13:33 UTC (rev 1051)
@@ -81,7 +81,7 @@
 		cf->sigma2[j] = data.sigma;
 		cf->res_sigma2[j] = data.res_sigma;
 		for ( k = 0; k < cf->wXL; k++ )
-			cf->beta_ests[3*j+k] = data.beta[k];
+			cf->beta_ests[cf->wXL*j+k] = data.beta[k];
 	}
 #if DEBUG
 	// Store data - For testing purposes
@@ -311,7 +311,7 @@
 		cf->sigma2[j] = data.sigma;
 		cf->res_sigma2[j] = data.res_sigma;
 		for ( k = 0; k < cf->wXL; k++ )
-			cf->beta_ests[3*j+k] = data.beta[k];
+			cf->beta_ests[cf->wXL*j+k] = data.beta[k];
 	}
 
 #if DEBUG