similar to: AIC bug?

I would like to figure out the equations for calculating "AIC" in both "step() function" and "AIC () function". They are different. Then I just type "step" in the R console, and found the "AIC" used in "step() function" is "extractAIC". I went to the R help, and found: "The criterion used is AIC = - 2*log L + k *

R 2.3.0, all packages up to date Linux, SuSE 10.0 Hi I want to calculate AIC or BIC from several results from mle calculation. I found the AIC function, but it does not seem to work with objects of class mle - If I execute the following: ml1 <- mle(...) AIC(ml1) I get the following error messale: Error in logLik(object) : no applicable method for "logLik" Therefore I am using the

Dear R users, According Brockwell & Davis (1991, Section 9.3, p.304), the penalty term for computing the AIC criteria is "p+q+1" in the context of a zero-mean ARMA(p,q) time series model. They arrived at this criterion (with this particular penalty term) estimating the Kullback-Leibler discrepancy index. In practice, the user usually chooses the model whose estimated index is

Perhaps, I should complement my previous e-mail question: I am trying to select variables through a quadratic discriminant analysis. It seems to me that a robust criterion is the value of the discriminatory power (i.e. the ratio of the separation of the class means to within-class variance) How can I obtain this ratio from V & R's qda() output ? Thanks (I hope this question is relevant) --

Dear all, I have written a small set of functions for drawing graphical representations of R objects and evaluations of R expressions. It is mainly intended for pedagogical purposes. Here are some samples: http://panini.u-paris10.fr/~sloiseau/sampleGraphics.pdf The idea is to use some graphical conventions in order to make clear the abstract syntax and to explain function, extraction,

Do you know any backend that implement instructions as a flag modifier in instruction ? Thank, Vincent Lejeune ----- Mail original ----- > De : Jakob Stoklund Olesen <stoklund at 2pi.dk> > À : Vincent Lejeune <vljn at ovi.com> > Cc : "llvmdev at cs.uiuc.edu (LLVMdev at cs.uiuc.edu)" <llvmdev at cs.uiuc.edu> > Envoyé le : Lundi 6 août 2012 20h06 >

Hello, I'm a Korean researcher who have been started to learn the "R" package. I want to make gam model and AIC value of the model to compare several models. I did the GAM model, but there were error for AIC. SO, how can I do? pleas help me!!! I did like below; > a.fit <- gam(pi~ s(t1r), family = gaussian(link="log")) >

Hi, R600 hardware (Radeon gfx card) does neither have a NEG nor an ABS instruction ; however any sources operand can be negated/abs'd by setting a bit for every source operand in the final bytecode (but not DST). A good way of modeling this behavior in LLVM is by using TargetFlag on operand. Currently the R600 LLVM backend in Mesa lower NEG and ABS DAG instruction to a MOV + TargetFlag using

Hello, Does logistic regression really provide better results than lda or qda ? (my purpose is not classification but highlighting of discriminant variables). If this is true, where could I get an R script performing stepwise logistic regression ? Thanks -- Daniel AMORESE Lab. M2C "Morphodynamique Continentale et C?ti?re" UMR CNRS 6143 Caen/Rouen Centre de G?omorphologie UCBN

Hello, I have random errors that I can not analyze nor fix: when several users connect to the same Samba share, they get sometimes error messages "the local drive name is already in use. This connection has not been restored", although the connection is still usable afterwards. This sometimes breaks a file operation. Clients are running Windows 2000, and the server is a Solaris 8

Hello All; I need to run a multiple regression analysis and use Akaike's Information Criterion for model selection. I understand that this command will give the AIC value for specified models: AIC(object, ..., k = 2) with "..." meaning any other optional models for which I would like AIC values. But, how can I specify (in the place of "...") that I want R to

On Aug 6, 2012, at 11:00 AM, Vincent Lejeune <vljn at ovi.com> wrote: > Ok. > > I tried to do it using a pass after register allocation, lowering NEG/ABS instructions. > However I met a problem : apparently getNextOperandForReg() can returns a MachineOperand before the one I'm processing. > > The following code snippet : > > > void

This is an informal poll. If you are primarily a Ruby programmer, - What was your primary language before you started coding in Ruby? Else, - What''s your current programming language of choice? I''ll start: Perl -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi, We have some problems with Xen 4.0.1 and resinzing of lvm volumes. When I resize a lvm volume with "lvresize -L+50G -n /dev/data/windows-data" the Windows Dom U don''t scan the new size of the lvm volume. The only way to see the extend is to halt and create the DomU We have try with sign and unsigned drivers gplpv, with Windows 2008 and Windows 2003 server :

Here it is :  Starting program: /home/vlj/llvmbin/bin/llc -march=r600 -mcpu=cayman /home/vlj/shader [Thread debugging using libthread_db enabled] Using host libthread_db library "/lib64/libthread_db.so.1". llc: /home/vlj/llvm/include/llvm/ADT/SmallVector.h:143: const T& llvm::SmallVectorTemplateCommon<T, <template-parameter-1-2> >::operator[](unsigned int) const [with T

I'm trying to send data to a server with openssh 7.9p1, but it's hanging somewhere. the client stop at the line : Jul 7 11:52:16 TOTO sshd[19553]: debug3: channel 0: will not send data after close and after 5 minutes the client closes the connection, why ? This is the trace of the server openssh : ( DEBUG3 level) Jul 7 11:52:15 TOTO sshd[31175]: debug3: fd 6 is not O_NONBLOCK Jul

Dear all, I have tried to modify the code of MASS:::stepAIC.default(), dropterm() and addterm() to use AICc instead of AIC for model selection. The code is appended below. Somehow the calculations are still not correct and I would be grateful if anyone could have a look at what might be wrong with this code... Here is a working example: ## require(nlme) model1=lme(distance ~ age + Sex, data =

Hello, I encounter a problem using the UPS MGE Evolution 850 VA USB with Debian Etch and the nut package available for this version. I got this error: ==syslog== Can't connect to UPS [MGE850VA] (usbhid-ups-auto): No such file or directory upsd[2767]: Startup successful upsmon[2769]: Startup successful upsd[2767]: Connection from 127.0.0.1 upsd[2767]: Client MGE at 127.0.0.1 logged

Greetings, I have recently come into some confusion over weather or not AIC results for comparing among models requires that they be nested. Reading Burnham & Anderson (2002) they are explicit that nested models are not required, but other respected statisticians have suggested that nesting is a pre-requisite for comparison. Could anyone who feels strongly regarding either position

I would like test AICc as a criteria for model selection for a glm using stepAIC() from MASS package. Based on various information available in WEB, stepAIC() use extractAIC() to get the criteria used for model selection. I have created a new extractAIC() function (and extractAIC.glm() and extractAIC.lm() ones) that use a new parameter criteria that can be AIC, BIC or AICc. It works as