Thank you. I tried revising with Arrays to try the matrix muliplication and found that I am still hitting a wall. See the code below for the change. proc model data=R_&y._&q._set2 OUTPARMS=M_rho_Initial_&y._&q. ;
parms M ;
exogenous M;
label M='M_factor_t' ;
array current [7] End_1-End_7;
array long [7] Endl_1-Endl_7;
array Scaled [7] scale_1-scale_7;
Bounds -1 <= M <=1;
eq.Esq = COUNT*(
( ( ( probnorm((long[1] - sqrt(&rho.)*M)/sqrt(1-&rho.)) - current[1] )**2 )*Scaled[1] )+
( ( ( probnorm((long[2] - sqrt(&rho.)*M)/sqrt(1-&rho.)) - current[2] )**2 )*Scaled[2] )+
( ( ( probnorm((long[3] - sqrt(&rho.)*M)/sqrt(1-&rho.)) - current[3] )**2 )*Scaled[3] )+
( ( ( probnorm((long[4] - sqrt(&rho.)*M)/sqrt(1-&rho.)) - current[4] )**2 )*Scaled[4] )+
( ( ( probnorm((long[5] - sqrt(&rho.)*M)/sqrt(1-&rho.)) - current[5] )**2 )*Scaled[5] )+
( ( ( probnorm((long[6] - sqrt(&rho.)*M)/sqrt(1-&rho.)) - current[6] )**2 )*Scaled[6] )+
( ( ( probnorm((long[7] - sqrt(&rho.)*M)/sqrt(1-&rho.)) - current[7] )**2 )*Scaled[7] )
);
solve M ;
run;
quit; I am multiplying through the Matrices and getting the following error: ERROR: The Newton method Jacobian matrix of partial derivatives of the
equations with respect to the variables to be solved is singular
at observation 1, for iteration 2. The system of equations cannot
be solved. I have set this up such that the sum of squared errors is fit on the multiplied matrices, but the matrix should multiply to a single summed value instead of a column; I do not think this is where the issue is. I could use further guidance on the Proc Model use and steps, as well as input on the approach, if any help can be provided.
... View more
