I'll look at your program more carefully, but my first impression is that it doesn't look like you are using the log-likelihood function.  It looks like you are summing the (raw) likelihood.   ... View more

For some examples of aligning the Y and Y2 axis, see "How to align the Y and Y2 axes in PROC SGPLOT."   The article discusses situations in which you want the tick marks on the Y axis to be aligned with the tick marks on the Y2 axis, even when Y and Y2 have different ranges and scales. ... View more

For more information about PROC MBC and model-based clustering, see Kessler, D. (2019) "Introducing the MBC Procedure for Model-Based Clustering", https://support.sas.com/resources/papers/proceedings19/3016-2019.pdf  ... View more

Use CONTOUR (you have a typo) ... View more

Dup of Solved: How do i get a block diagonal covariance matrix, V, in SAS Proc Mixed ... - SAS Support Communities ... View more

Thanks. I have alerted the manager for that documentation. For the record, the GEE and GENMOD procedures support the EFFECTPLOT statement, but not the EFFECT statement. ... View more

> You can enrich your dataset with new covariates (spline effects) using PROC TRANSREG and subsequently use the procedure of your choice.   An easier choice is PROC GLMSELECT, which supports an OUTDESIGN= option. See the example at https://blogs.sas.com/content/iml/2019/10/16/visualize-regression-splines.html ... View more

I don't know where that sentence came from, but it is not correct. GENMOD does not support the EFFECT statement so the OP must use GLMMIX.   DOC for EFFECT stmt: SAS Help Center: EFFECT Statement says "This section applies to the following procedures: GLIMMIX, GLMSELECT, HPMIXED, LOGISTIC, ORTHOREG, PHREG, PLS, QUANTLIFE, QUANTREG, QUANTSELECT, ROBUSTREG, SURVEYLOGISTIC, and SURVEYREG." ... View more

Look at the R table of parameter estimates. The R model appears to treat all variables as continuous, whereas in SAS you are treating many variables as factors. Resolve that issue first or you'll never get the same estimates.   Also, I suspect you want to add Sex to the CLASS statement so that the Sex variable is treated as a categorical variable.      ... View more

> do I have to have it?  I included that information so that you can easily check the names of the variables in the input matrix. But you can ignore that suggestion if you prefer an alternative.  ... View more

The changes that we've discussed would make your program look like this: proc iml; use pseudo_data; read all var _NUM_ into rc_migr_dist[colname=varNames]; /* Load matrix with movements shares p_ij */ close pseudo_data; nr = nrow(rc_migr_dist); if upcase(varNames)[1] ^= "MK1" then do; rc_migr_dist = j(nr, 1, 0) || rc_migr_dist; varNames = "MK1" || varNames; end; if nr^=ncol(rc_migr_dist) then print "ERROR: Matrix is not square"; rc_migr_dist = rc_migr_dist / 100; /* ??? Is this necessary? */ print rc_migr_dist[r=varNames c=varNames]; use rclassorgcount; read all var {&SOMETHING_THAT.} into rc_orig_count_vec; /* Load vector with number of humans per movement i, Ni */ close rclassorgcount; /* Calc Matrix_normal_u */ M_norm_u = 0; do i = 1 to &K. - 1; maxval = max(abs(i - &K.), abs(i + (&K_Max. - &K.) - 1)); N_i = rc_orig_count_vec[i]; sum_p = sum( rc_migr_dist[i, i+1:nr] ); /*print sum_p;*/ temp = maxval * N_i * sum_p; M_norm_u = M_norm_u + temp; end; create M_norm_u var {"Matrix_normal_u"}; append from M_norm_u; close M_norm_u; ... View more

I suggest you read in the names of the variables dynamically (instead of using a macro) and then add the first column if necessary: proc iml; use pseudo_data; read all var _NUM_ into rc_migr_dist[colname=varNames]; /* Load matrix with movements shares p_ij */ close pseudo_data; nr = nrow(rc_migr_dist); print varNames; if upcase(varNames)[1] ^= "MK1" then do; rc_migr_dist = j(nr, 1, 0) || rc_migr_dist; varNames = "MK1" || varNames; end; print rc_migr_dist[r=varNames c=varNames]; Be aware that this empirical estimate of the transition matrix will be singular. You are specifying a system that will never allow a transition to the first state. If this is not what you want, then use a small positive probability (such as 1E-3) for the first column and adjust the remaining probabilities so that the rows sum to 1. ... View more

If there are 7 states, then each cell indicates the proportion of accounts that transitioned from State=i into State=j during the time period. If no account transitioned to State=7, then the last column of your 7x7 matrix should contain all zeros.   If "add a nullvector" means "add a column of 0s," then I agree. Otherwise, explain a "nullvector" in this context. In matrix theory, a null vector is a vector v such that A*v=0 for the matrix A. But I don't think you are using that definition(?). ... View more