Displaying 20 results from an estimated 898 matches for "bics".
Did you mean:
bias
2005 Apr 18
2
Why no BIC.default function?
I'm using R 2.0.1.
I looked in the email archives but didn't see anything on this topic.
I've noticed a surprising (to me) difference between AIC and BIC:
> methods("AIC")
[1] AIC.default* AIC.logLik*
> methods("BIC")
[1] BIC.gls* BIC.lm* BIC.lme* BIC.lmList* BIC.logLik* BIC.nls*
The BIC.gls BIC.lm BIC.lme BIC.lmList and BIC.nls functions appear
2010 May 18
1
BIC() in "stats" {was [R-sig-ME] how to extract the BIC value}
>>>>> "MM" == Martin Maechler <maechler at stat.math.ethz.ch>
>>>>> on Tue, 18 May 2010 12:37:21 +0200 writes:
>>>>> "GaGr" == Gabor Grothendieck <ggrothendieck at gmail.com>
>>>>> on Mon, 17 May 2010 09:45:00 -0400 writes:
GaGr> BIC seems like something that would logically go into stats
2006 Mar 06
2
[Q] BIC as a goodness-of-fit stat
Dear R-List
I have a question about how to interpret BIC as a goodness-of-fit statistic.
I was trying to use "EMclust" and other "mclust" library and found that BIC
was used as a goodness-of-fit statistic.
Although I know that smaller BIC indicates a better fit, it is not clear to
me how good a fit is by reading a BIC number. Is there a standard way of
interpreting a BIC
2009 Apr 05
1
Which model to keep (negative BIC)
Hi,
My questions concern the function 'mclustBIC' which compute BIC for a range of clusters of several models on the given data and the other function 'mclustModel' which choose the best model and the best number of cluster accordind to the results of the previous cited function.
1) When trying the following example (see ?mclustModel), I get negative BIC computed by
2008 May 13
1
R help: problems with step function
Dear List Members,
I have encountered two problems when using the step function to
select models. To better illustrate the problems, attached is an
R image which includes the objects needed to run the code attached.
lm.data.frame have factor variables with 3 levels.
The following run shows the first problem. AICs (* and **) are different.
I noticed that the Df for rs13482096:rs13483699 is 4,
2009 Jun 28
1
ggplot2 x axis question
Hi,
I have 45 models that I have named: 1, 2, 3, ... , 45 and I am trying to
plot them in order of ascending BIC values. I am however unclear as to
how I can get the models to line up on the x-axis by BIC and not by
numeric order. For example, if model 5 has a lower BIC than 1, I want it
to be the first point on the left hand side of the curve. This seems to
work in plot:
plot(1:45,
2009 Jun 28
1
ggplot2 x axis question
Hi,
I have 45 models that I have named: 1, 2, 3, ... , 45 and I am trying to
plot them in order of ascending BIC values. I am however unclear as to
how I can get the models to line up on the x-axis by BIC and not by
numeric order. For example, if model 5 has a lower BIC than 1, I want it
to be the first point on the left hand side of the curve. This seems to
work in plot:
plot(1:45,
2008 Sep 24
2
Error message when calculating BIC
Hi All,
Could someone help me decode what this error means ?
> BIC(nb.80)
Error in log(attr(object, "nobs")) :
Non-numeric argument to mathematical function
>
BTW, nb.80 is a negative binomial glm model created using the MASS
library with the call at the bottom of the message
In the hopes of trying to figure this out I tried the following
workaround but it did not work
2008 Dec 19
0
What BIC is calculated by 'regsubsets'?
The function 'regsubsets' appears to calculate a BIC value that is
different from that returned by the function 'BIC'. The latter is
explained in the documentation, but I can't find an expression for the
statistic returned by 'regsubsets'.
Incidentally, both of these differ from the BIC that is given in Ramsey
and Schafer's, The Statistical Sleuth. I assume
2009 Sep 05
2
About BIC
Hello,
I am working on getting optimal lags by using BIC, But I don't know
how to calculate BIC. Is there any code or useful function for it?
Thanks and regards,
Dan Zhao
2010 Dec 26
1
Calculation of BIC done by leaps-package
Hi Folks,
I've got a question concerning the calculation of the Schwarz-Criterion
(BIC) done by summary.regsubsets() of the leaps-package:
Using regsubsets() to perform subset-selection I receive an regsubsets
object that can be summarized by summary.regsubsets(). After this
operation the resulting summary contains a vector of BIC-values
representing models of size i=1,...,K.
My problem
2003 Nov 21
1
: BIC for gls models
Hi all,
I would like to know how the BIC criterion is calculated for models estimated using gls( ) function. I read in Pinheiro & Bates (2000) p84 that
BIC = -2logL + npar*log(N) (for the ML method), or
BIC = -2logLR + npar*log(N-p) (for the REML method)
but when I use any of these formulae I don't obtain the result given by R.
Thanks in advance for any help.
Eve CORDA
Office national
2011 Sep 01
4
Question about BIC of two different regression models? how should we compare two regression models?
Hi All,
In order to compare two different logistic regressions, I think I need to compare them based on their BIC values, but I am not sure if the smaller BIC would mean a better model or the reverse is true?
Thanks a lot,Andra
[[alternative HTML version deleted]]
2008 Jul 15
2
sem & testing multiple hypotheses with BIC
I'm coming from the AMOS world and am wondering if there is a simple
way to do multiple hypothesis testing in the manner of BIC analyses in
AMOS using the sem package in R. I've read the documentation, but
don't see anything in there except for basic BIC scores. Perhaps
someone has devised a simple way to compare the relative likelihood of
all possible path-fittings within a
2005 Oct 16
1
BIC doesn't work for glm(family=binomial()) (PR#8208)
Full_Name: Ju-Sung Lee
Version: 2.2.0
OS: Windows XP
Submission from: (NULL) (66.93.61.221)
BIC() requires the attribute $nobs from the logLik object but the logLik of a
glm(formula,family=binomial()) object does not include $nobs. Adding
attr(obj,'nobs') = value, seems to allow BIC() to work.
Reproducing the problem:
library(nmle);
BIC(logLik(glm(1~1,family=binomial())));
2004 Jul 01
3
BIC vz SBIC vz SIC
DeaRs,
I have a doubt about:
BIC (Bayesian Information Criterion)
SBIC (Schwartz Bayesian Informarion Criterion)
SIC (Schwartz Information Criterion)
In many references these are know as the same (eg. stepAIC() function) but I
just found a SAS8.2 output that show either the BIC and SIC values for a
logistic regression.. simillary values but different.
1) question: What are the differences?
2006 Oct 05
1
lmer BIC changes between output and anova
list,
i am using lmer to fit multilevel models and trying to use anova to compare the models. however, whenever i run the anova, the AIC, BIC and loglik are different from the original model output- as below. can someone help me out with why this is happening? (i'm hoping the output assocaited with the anova is right!).
thank you,
darren
> unconditional<-lmer(log50 ~ 1 + (1 |
2008 May 14
0
Problems with step function
Dear List Members,
I have encountered two problems when using the step function to
select models. To better illustrate the problems, an R image
(step.add1.test.RData)
which includes the objects needed to run the code (step.add1.test.R) can be
found at
www.biostat.wisc.edu/~pwang/r-help/<http://www.biostat.wisc.edu/%7Epwang/r-help/>
lm.data.frame have factor variables with 3 levels.
The
2010 Aug 17
2
how to selection model by BIC
Hi All:
the package "MuMIn" can be used to select the model based on AIC or AICc.
The code is as follows:
data(Cement)
lm1 <- lm(y ~ ., data = Cement)
dd <- dredge(lm1,rank="AIC")
print(dd)
If I want to select the model by BIC, what code do I need to use? And when
to select the best model based on AIC, what the differences between the
function "dredge" in
2007 Sep 17
1
Stepwise logistic model selection using Cp and BIC criteria
Hi,
Is there any package for logistic model selection using BIC and Mallow's Cp
statistic? If not, then kindly suggest me some ways to deal with these
problems.
Thanks.
--
View this message in context: http://www.nabble.com/Stepwise-logistic-model-selection-using-Cp-and-BIC-criteria-tf4464430.html#a12729613
Sent from the R help mailing list archive at Nabble.com.