turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- I get the same BIC and -2 Res Log!?

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

06-16-2015 05:55 PM

Hi guys!

I need to know why I get the same BIC and -2 Res log likelihood values. Shouldnt that be impossible when the equation for BIC is -2ReslogL + 2p, p = number of parameters.

Im using this code for a dataset of macro economic data

PROC IMPORT OUT=oil DATAFILE= "W:\oil.xlsx"

DBMS=xlsx REPLACE;

GETNAMES=YES;

RUN;

proc mixed data=oil covtest noitprint;

class country year quarter;

model hcip=oilprice interest exchange unemploy output_g/solution residual;

random intercept/subject=country;

random intercept/subject=year(country);

repeated /type=ar(1) subject=year;

run;

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

06-16-2015 08:51 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

06-17-2015 08:55 AM

What about the other information criteria? Are they messed up as well?

It could be that the -2 log likelihood value is so large that adding in 2p gets swamped in the printed output. Is that a possibility? To check, use ODS output to get the IC into a dataset, and look at the values in a long format.

Steve Denham

Message was edited by: Steve Denham Removed calculated value that depended on number of fixed effects

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

06-17-2015 09:26 AM

By the way, the formula for BIC is -2LL + p.ln(n). You gave the formula for AIC. Also note that p is the number of variance-covariance parameters with REML estimation (what you are doing). The fixed effect parameters have nothing to do with it (you may already know this, but other readers may not). Are all your variance-covariance parameters 0? Steve also has a good point. Would be good to see all your output.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

06-17-2015 03:04 PM

I'll add that procedures that minimize some sort of information criterion (LOESS, GLMSELECT, SEVERITY, any model selection procedure) sometimes add or substract constants. Not relevant for this problem, I don't think, but some sources might give F(n, p; x) + C, whereas SAS reports F(n,p; x).

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

06-17-2015 04:31 PM

Good point about statistical software. This can lead to a lot of confusion in comparing output from different programs with the same data and same model. In MIXED and GLIMMIX, however, one does get the complete likelihood (or restricted likelihood) in the information-criterion statistics, including the term(s) that does not depend on the the parameters. For REML estimation, only variance-covariance parameters are used in the calculation because the fixed effect are eliminated at each ML step. I hope to see a compete output from the OP to help figure out this issue.