diff --git a/matlab/partial_information/PCL_Part_info_irf.m b/matlab/partial_information/PCL_Part_info_irf.m
index e28c062cd..4f291098f 100644
--- a/matlab/partial_information/PCL_Part_info_irf.m
+++ b/matlab/partial_information/PCL_Part_info_irf.m
@@ -29,7 +29,10 @@ function  y=PCL_Part_info_irf( H, varobs, ivar, M_, dr, irfpers,ii)
 % The jump variables have dimension NETA
         
         
-    [junk,OBS] = ismember(varobs,M_.endo_names,'rows');
+    OBS = [];
+    for i=1:rows(varobs)
+        OBS = [OBS find(strcmp(deblank(varobs(i,:)), cellstr(M_.endo_names))) ];
+    end
     NOBS = length(OBS);
     
         G1=dr.PI_ghx;
diff --git a/matlab/partial_information/PCL_Part_info_moments.m b/matlab/partial_information/PCL_Part_info_moments.m
index b68b3d17d..d0f7f8c11 100644
--- a/matlab/partial_information/PCL_Part_info_moments.m
+++ b/matlab/partial_information/PCL_Part_info_moments.m
@@ -1,4 +1,4 @@
-function  AutoCOR_YRk=PCL_Part_info_irmoments( H, varobs, dr,ivar)
+function  AutoCOR_YRk=PCL_Part_info_moments( H, varobs, dr,ivar)
 % sets up parameters and calls part-info kalman filter
 % developed by G Perendia, July 2006 for implementation from notes by Prof. Joe Pearlman to 
 % suit partial information RE solution in accordance with, and based on, the 
@@ -31,7 +31,10 @@ function  AutoCOR_YRk=PCL_Part_info_irmoments( H, varobs, dr,ivar)
     warning_old_state = warning;
     warning off
 
-    [junk,OBS] = ismember(varobs,M_.endo_names,'rows');
+    OBS = [];
+    for i=1:rows(varobs)
+        OBS = [OBS find(strcmp(deblank(varobs(i,:)), cellstr(M_.endo_names))) ];
+    end
     NOBS = length(OBS);
     
     G1=dr.PI_ghx;