*-http://www.pinggu.org/bbs/thread-992234-1-1.html

  clear
  input city ind famly
  110101                    1                       1
  110101                    2                       1
  110101                    3                       1
  110101                    1                       2
  110101                    2                       2
  110101                    1                       3
  110101                    2                       3
  110101                    3                       3
  110101                    4                       3
  110101                    1                       4
  110101                    2                       4
  end

  drop famly

  qui tab city
  local N_city = r(r)
  qui sum ind
  local N_ind = r(N)
  gen family = .
  set trace on
  forvalues c = 1/`N_city'{
     local j = 0
     forvalues i = 1/`N_ind'{
           if ind[`i']==1{
             replace family = `++j' in `i'
           }
           else{
             replace family = `j' in `i'
           }            
         }
  }