From ac4ae0b6d2d5eef56ba476d93a0f9fe37d8ea97a Mon Sep 17 00:00:00 2001
From: Michel Juillard <michel.juillard@mjui.fr>
Date: Sat, 22 Oct 2011 15:28:41 +0200
Subject: [PATCH] make sure that dataset_.missing.aindex always exist, even
 when there are no missing observations

---
 matlab/utilities/dataset/describe_missing_data.m | 2 +-
 matlab/utilities/dataset/initialize_dataset.m    | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/matlab/utilities/dataset/describe_missing_data.m b/matlab/utilities/dataset/describe_missing_data.m
index 7b7fedf4a..595ee3c6d 100644
--- a/matlab/utilities/dataset/describe_missing_data.m
+++ b/matlab/utilities/dataset/describe_missing_data.m
@@ -59,7 +59,7 @@ for obs=1:T
     if rows(tmp(:))
         i(obs) = { tmp(:) };
     else
-        i(obs) = { [] };
+        i(obs) = { variable_index(:) };
     end
 end
 
diff --git a/matlab/utilities/dataset/initialize_dataset.m b/matlab/utilities/dataset/initialize_dataset.m
index 95d81cbb0..4f9995209 100644
--- a/matlab/utilities/dataset/initialize_dataset.m
+++ b/matlab/utilities/dataset/initialize_dataset.m
@@ -78,7 +78,7 @@ if dataset_.missing.state
     dataset_.missing.number_of_observations = n;
     dataset_.missing.no_more_missing_observations = s;
 else
-    dataset_.missing.aindex = [];
+    dataset_.missing.aindex = num2cell(repmat(1:dataset_.info.nvobs,dataset_.info.ntobs,1),2);
     dataset_.missing.vindex = [];
     dataset_.missing.number_of_observations = [];
     dataset_.missing.no_more_missing_observations = [];