diff --git a/melpca.cc b/melpca.cc
index 3c703ee7aada9e918f1a32e891beff32b2895dc3..6bf883463c4b31eda2ddb2e21ebef2dfd73f2518 100644
--- a/melpca.cc
+++ b/melpca.cc
@@ -35,7 +35,7 @@ namespace Melodic{
 	if(opts.paradigmfname.value().length()>0)
 	{
 		basicGLM tmpglm;
-		tmpglm.olsfit(in,melodat.get_param(),Identity(melodat.get_param().Ncols()));
+		tmpglm.olsfit(in,melodat.get_param(),IdentityMatrix(melodat.get_param().Ncols()));
 		std_pca(tmpglm.get_residu(),weights,Corr,tmpE,tmpD);
 //		std_pca(in,weights,Corr,tmpE,tmpD,melodat.get_param());
 	}