wolffd@0: function f = family(A,i,t)
wolffd@0: % FAMILY Return the indices of parents and self in sorted order
wolffd@0: % f = family(dag,i,t)
wolffd@0: %
wolffd@0: % t is an optional argument: if present, dag is assumed to be a 2-slice DBN
wolffd@0: 
wolffd@0: if nargin < 3 
wolffd@0:   f = [parents(A,i) i];
wolffd@0: else
wolffd@0:   if t == 1
wolffd@0:     f = [parents(A,i) i];
wolffd@0:   else
wolffd@0:     ss = length(A)/2;
wolffd@0:     j = i+ss;
wolffd@0:     f = [parents(A,j) j] + (t-2)*ss;
wolffd@0:   end
wolffd@0: end