Displaying 20 results from an estimated 3000 matches similar to: "BIC() in "stats" {was [R-sig-ME] how to extract the BIC value}"
2006 Jul 04
1
lmer print outs without T
Hi,
I have been having a tedious issue with lmer models with lots of
factors and lots of levels. In order to get the basic information at
the beginning of the print out I also have to generate these enormous
tables as well. Is there a method command to leave off all of the
effects and correlations? Or, do I have to go to string commands?
2006 Oct 18
1
lmer- why do AIC, BIC, loglik change?
Hi all,
I am having issues comparing models with lmer. As an example, when
I run the code below the model summaries (AIC, BIC, loglik) differ between
the summary() and anova() commands. Can anyone clear up what's wrong?
Thank you!
Darren Ward
library(lme4)
data(sleepstudy)
fm1<-lmer(Reaction ~ Days + (1|Subject), sleepstudy)
summary(fm1)
fm2<-lmer(Reaction ~ Days +
2006 Aug 11
1
help:coerce lmer.coef to matrix
Hi,
Thanks for your response, it nearly worked! But it only wrote one coloumn
of data and not the three columns I need.
Using fixef(m1) doesnt give the same results as coef(m1) when you are
using more than one random effect. I need the coefficients for each
individual so I use coef(m1) to get this which results in an object of
class lmer.coef, 3 columns by 700 rows.
as.data.frame() wont work on
2007 Jun 25
1
conflict between lme4 and RMySQL packages (PR#9753)
Full_Name: Dale Barr
Version: 2.5.1 (patched)
OS: Ubuntu linux x86_64
Submission from: (NULL) (138.23.70.108)
When RMySQL is loaded in before lme4, the summary() function for lmer objects in
the lme4 packages produces the following error:
Error in printMer(object) : no slot of name "status" for this object of class
"table"
When RMySQL is loaded AFTER lme4, however, no such
2007 Jul 31
1
Extracting random parameters from summary lme and lmer
LS,
I'm estimating multilevel regression models, using the lme-function
from the nlme-package. Let's say that I estimated a model and stored
it inside the object named 'model'. The summary of that model is
shown below:
Using summary(model)$tTable , I receive the following output:
> summary(model)$tTable
Value Std.Error DF t-value
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
2007 Aug 16
3
function to find coodinates in an array
Dear list,
I am looking for a function/way to get the array coordinates of given
elements in an array. What I mean is the following:
- Let X be a 3D array
- I find the ordering of the elements of X by ord <- order(X) (this
returns me a vector)
- I now want to find the x,y,z coordinates of each element of ord
Can anyone help me?
Thanks!
Ana
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
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
2009 May 19
0
Automatic Differentiation for R
Martin (see below) gives a good explanation of the difference between AD and symbolic
differentiations. I'm of the opinion we can use both. However, the real issue as far
as I'm concerned (from an optimizer's point of view, which may also be that of ODE and
PDE folk) is that right now none of the offerings that we have are easy to use. Indeed,
usability is one of the key issues in my
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
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 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
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]]
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
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())));
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
2011 Dec 20
2
Extract BIC for coxph
Dear all,
is there a function similar to extractAIC based on which I can extract the
BIC (Bayesian Information Criterion) of a coxph model?
I found some functions that provide BIC in other packages, but none of them
seems to work with coxph.
Thanks,
Michael
[[alternative HTML version deleted]]
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.