Estimating Unobserved Complementarities between Entrepreneurs and Venture Capitalists Matlab Code

Jump to navigation Jump to search

Main Project here: Estimating Unobserved Complementarities between Entrepreneurs and Venture Capitalists

Important contractions

  • GMM: Generalized Method of Moments
  • MSM: Method of Simulated Moments
  • MLE: Maximum Likelihood Estimation
  • GA: Genetic Algorithm
  • SE: Standard Error

How to run

First, you need a gurobi license, which can be obtained (free) here.

In master.m, edit the options section to reflect what you want the code to do. Then run it.

Location/structure of the code

The old codebase is located in


The new codebase is located in


The new-new codebase is located in


Some of this code is just library code (prefix mtimesx).


options section

task: Can take the values {'data', 'monte', 'monte_data'}.

estimator: Can take the values {'MSM'}. (other estimators removed in readjusted code)

use_solver: Can take the values {'ga'}. (other solvers removed in readjusted code)

error_type (currently hard coded as 1 and isn't fully written to support 2): 1 for match specific errors, with the error distribution following an exchangeable structure. var(e) = sig^2, and cov(e,e') = 1/4*sig^2. 2 for agent specific errors, with the error structure of match <i, j> as sig*ei*ej.


Does the majority of the work for this problem. Runs the GA, saves the results to 'empirics_match_specific_1st_stage_ga', then runs it again with different globals, and saves it to 'empirics_match_specific_2nd_stage_ga'.


Constraints on GA. For [c,ceq] = nonlinearcons_msm(x), GA constrains x such that c ≤ 0 and ceq = 0. c and ceq are row vectors when there are multiple constraints. ceq is unused for our purposes.


This is the fitness function. Takes a vector and returns a scalar. GA minimizes this function.


Generates the moments needed for the GMM.

Location/Structure of Data

Pro tip: running the command

 whos -file filename

from a matlab session will tell you the contents of any .mat file.

There are two of both psdata and dyad_tech_mkt_data, corresponding to 4 industries or 5 industries. You can find them in

 E:\McNair\Projects\MatchingEntrepsToVC\OriginalCode\FVEIC4 data


 E:\McNair\Projects\MatchingEntrepsToVC\OriginalCode\FVEIC5 data

The one wanted should be copied into the code directory:


psdata.mat contains

  • vc: the number of VC's in each market, in the form of (10, 12, 13, ...): the first market has 10 VC's, the second market has 12 VC's....
  • firm: the number of firms in each market, (13, 14, ...): the first market has 13 downstream, the second 14 downstream... this is a one (VC) to many (firms) market, and the number of firms in each market is at least as many as the number of firms
  • m_id: unused, I'm guessing the market id

dyad_tech_mkt_data.mat variables used include:

  • pvc_exp_n, mean_pvc_exp_n, std_pvc_exp_n
  • lnfpat, mean_lnfpat:
  • mean_m_dist_1000, std_m_dist_1000:
  • mean_exp_sector, std_exp_sector:
  • m_match: contains the matching outcome: for each dyad (firm VC pair), 0 means not matched and 1 means matched. The variable asg in the code contains the matching outcomes of all markets. The matching outcomes are delineated by firms using the position data N1 (vc' from psdata)and N2 (firm' from psdata).

Location/Structure of Output






Simulated moments matrix dimension mismatch. This bug (also in the adjusted code) does not allow the "optimal stage" to complete.

 Assignment has more non-singleton rhs dimensions than non-singleton subscripts
 Error in gmm_2stage_estimation (line 65)
       EV(:, :, m) = EV(:, :, m) + temp2;
 Error in master (line 513)

This is caused because temp2 is 166x166 in monte_data, but nm = M = 1, so EV is 1x1. nm is set to size(M0, 1), but M0 is sometimes a scalar?

This is solved by commenting out mkt_resample.m line 107:

 % M0 = moments(M, 1, S, K, pv', asg, H, HK, G, FF, Hu, Hd, z1, z2, psa, ps1, ps2, N1, N2) 

because it returns a scalar, but should return a 166x1 vector.

Better fix notes: The problem is narrowed down to the section of mkt_resample.m before the call to moments. Somehow the parameters to moments are changed such that it returns a scalar.

Solution found (I think): monte_M cannot be set to 1. For values greater than 1, the program seems to work. It doesn't crash, anyway.


Ran first iteration with task='monte', use_solver='ga', estimator='MSM', gurobi output disabled (took 42 sec to finish 1st iteration)

Unobserved comp profile monte.png

Monte_data and monte run much faster than data, because data is so large.


There is a lot of unused code. For example, monte_std and monte_x seem to only be assigned, not used.