I'm running SAS on demand for academics. My question how to make the code file= statement work? I'm getting the following error in the LOG: ERROR: HPSPLIT was unable to open the code file for output.   ods graphics on; proc hpsplit data=Wine CVCC seed=123; class Cultivar; model Cultivar = Alcohol Malic Ash Alkan Mg TotPhen Flav NFPhen Cyanins Color Hue ODRatio Proline; code file ="winecode.sas"; grow entropy; prune costcomplexity; run; ... View more

I would like to specify a "spike and slab" prior in PROC MCMC. Let beta be the parameter. I need beta ~ (1 - p)*D + p*Normal(0,var) where D is a point mass at zero and p is a beta(a,b) ... or perhaps more simply, p = 0.5 ... View more

PROC LOGISTIC: The "Error Rate" from FITSTAT does not agree with the "Correct" percent from the CTABLE (with PPROB=.5). In the example below it is easy to see from the PROC FREQ that the error-rate is 0.333, as is given by FITSTAT. (An OBS is correctly classified if P_1 >= 0.5 and Y=1  OR  P_1 < 0.5 and Y=0.) From the classification table the "Correct" is 0.417 and so this gives an error-rate of 0.583 What are the rules for the Classification Table that leads to "Correct" of 0.417?   DATA TEST1; DO ID = 1 to 12; Y = MOD(ID,2); X1 = (ranuni(1) < .6); X2 = (ranuni(1) < .4); X3 = (ranuni(1) < .5); OUTPUT; END; run; PROC LOGISTIC DATA = TEST1 desc; MODEL Y = X1-X3 / CTABLE PPROB=0.5; SCORE DATA = TEST1 OUT=SCORED FITSTAT; run; PROC FREQ DATA = scored; TABLES P_1*Y /norow nocol nopercent; run; ... View more

This documentation link gives the formula for ATT weights (Let ps=propensity score, t is treatment) https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/statug/statug_causaltrt_details06.htm The formulas for the weights are: Control: Let CTOTAL = Sum[ (1-t)*ps/(1-ps) ]. Then a control weight is ( (1-t)*ps/(1-ps) ) / CTOTAL Treated: Let TTOTAL = Sum[ t ]. Then a treated weight is t / TTOTAL I don’t see how the weights from the PROC CAUSALTRT example are in agreement with the formulas. See column ATTwgt2 in the example below. In particular, I think the weights for the treated observations should all equal 0.2 … which is not in agreement with column ATTwgt2. Any help is appreciated! data WORK; input Y T X; datalines; 1 0 1 2 1 1 3 0 1 1 1 1 1 0 3 1 1 3 4 0 4 2 0 4 3 1 4 2 0 4 3 0 4 2 1 4 ; PROC CAUSALTRT DATA= WORK ATT; CLASS X; PSMODEL T(descending)=X; MODEL Y;   OUTPUT OUT=OUTwork2 PS=PS2 IPW=ATTwgt2; proc print data = OUTwork2; run;   Obs    Y      T      X      PS2    ATTwgt2 1      1      0      1      0.50000 2.00000 2      2      1      1      0.50000 2.00000 3      3      0      1      0.50000 2.00000 4      1      1      1      0.50000 2.00000 5      1      0      3      0.50000 2.00000 6      1      1      3      0.50000 2.00000 7      4      0      4      0.33333 1.50000 8      2      0      4      0.33333 1.50000 9      3      1      4      0.33333 3.00000 10     2      0      4      0.33333 1.50000 11     3      0      4      0.33333 1.50000 12     2      1      4      0.33333 3.00000 ... View more