Hi all,   I am trying to understand why in the example for GEE (page 3) the Type 1 option is being used rather than the Type 3 option.    https://support.sas.com/documentation/cdl/en/statug/63962/HTML/default/viewer.htm#statug_genmod_sect064.htm   Also it seems to contradict with what is described here: https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_genmod_sect034.htm   Which says "This type of analysis is not available for GEE models, since the deviance is not computed for this type of model."   Thanks, Doug   ... View more

Hi All,   So I have data as follows:  Time = X  Result = Y   Time (Months)   Result 0 7.17 3 5.98 6 7.31 12 73.87 18 80.92 30 39.11 36 37.93 42 33.3 48 35.2 78 33.26   I need to fit a cubic spline to the data, and then find Y-values for 8 equal time points between 19.2 and 44.4 months. And then do a regression on those 8 data points to find a slope.   I have been playing around with proc glmselect, and I get 5 knots and their parameters...but I can't seem to apply the equation correctly after the first knot, I am not sure how to get SAS to show me the other parts of the predictive equation (i.e. the second intercept)   This is how far I got: ods select ANOVA ParameterEstimates SplineKnots TPFSplineDeatils; ods output ParameterEstimates=param SplineKnots=spline ; proc glmselect data=patient1; effect spl = spline( time/ details naturalcubic basis=tpf(noint) knotmethod=percentiles(5) ); model result = spl / selection=none ; /* fit model by using spline effects */ output out=SplineOut predicted=fit ; /* output predicted values for graphing */ quit; ... View more

So currently I am assessing agreement between two evaluators. They are assessing disease response and the table would like below:                      Evaluater 2  CR     CRc      CRi        PR        SD       PD Evaluater 1 CR CRc CRi PR SD PD   To calculate a weighted Kappa is assigned an Ordinal value to each response: 1, 2, 3, 4, 5, 6. Which sas that uses the following weights as default: 1.00, 0.8, 0.6, 0.4, 0.2.   However, I would really like for CRc and CRi to have the same weights for distance from CR. Does anyone know how I should reorder the ordinal values if even possible, my issue is that they can't have the same ordinal value --- is the only way to combine CRc & CRi? I don't know if that would be clinically meaningful. ... View more

Looking for general information about 1-sided CMH tests. Also if anyone has any code with regards to this test that would be helpful.   can use variables strata*treatment*outcome. ... View more