Creates a l\*m\*n x (m^l)\*(n^m) dimensional matrix if y_zx_dependent = FAlSE, creates l\*m\*n x (m^l)\*((n^m)+1) dimension matrix otherwise.
Create_P(l, m, n, Rxy, y_zx_depenendent = FALSE)|Z|
|X|
|Y|
Matrix of response function indicators
If set to true P will be calculated assuming Y is dependent on z and x, otherwise will assume y is only dependent on x.