02-08-2010 07:59 PM
02-10-2015 01:10 PM
I'm having this same issue. I've been using surveylogistic mostly in an analysis of HCUP NIS data, but when I went to analyze the effect of year on a couple of continuous variables with surveyreg, I'm getting infinite t-values, and in one domain, the reference category for year has values when it should have none.
proc surveyreg data=sga.sga4yrs nomcar;
class yearcat insure;
model los = yearcat female insure conganom / solution;
All the domain variables are dichotomous categorical variables. Weight is supplied by HCUP. LOS is continuous. Yearcat is a categorical (4 categories) variable indicating which year the observation is from (reversed the order so the earlier year would be the reference cat). If I run yearcat as continuous, I still get problems.
Does anything pop out at people? Any more information needed?
02-10-2015 01:42 PM
1) it would be better to start a new question
2) in the appropriate part of the forum such as statistical procedures
3) post some of your output such as the data summary, fit statistics and estimated regression coefficients
02-10-2015 04:09 PM
Maybe I'm just fortunate, but I have never seen a t-value of infinity in SAS output. I have seen missing t-values, is that what you have, or is it truly showing as infinity?
02-11-2015 07:59 AM
Well, to get a t value of infinity, that would imply that the standard error is zero. I have seen that in complete collinearity situations.
02-11-2015 08:56 AM
Hey, I'm the old guy with memory problems. Anyway, I have seen F and t values reported as "Infnty". These are almost always associated with df=0 (mixed models).
02-11-2015 09:59 AM
Steve, when I run the models, I do get 0 df for whatever the first variable in the class statement is. I'm going to try to attach some results. I would paste them in, but for some reason my browser isn't letting me paste into this text box, and I have no other options at the moment.
If I remove the weight statement, I suddenly get results that at least appear normal (no infinite t-values, no values for my reference category).
02-13-2015 10:15 AM
Yeah, a look at the output should clear this up some. I am thinking that it has to do with the default parameterization with the WEIGHT statement as opposed to without the statement.
02-13-2015 03:05 PM
Everything looks fine except the second portion (unweighted and t values of Infnty) and the first where all the yearcat coefficients appear (none set to zero). I notice that the output in the second portion does not include the clusters (assume that it is the second of the possible domains). Could you possible split those 3M obs into a separate dataset and do some exploring there? I know you said you have subsetted and reduced the model before, but it looks like that particular domain is giving you an unweighted analysis, and I don't see any reason for it.
As far as the first, since a generalized inverse was used, and the estimates are identical to within roundoff errors, I suspect. Try printing the X'WX matrix and its inverse (X and I options to the MODEL statement), and maybe something will jump out.