similar to: Likelihood ration test on glm

Does anyone know of a package that includes the Modified Sims test [Gewerke, 1983, Sims, 1972]? This test is used in econometrics and is a kind of alternative to the Granger test [Granger, 1969], which is in the package lmtest. Thanks in advance, chris Refernces: Gewerke, J., R. Meese, and W. Dent (1983), "Comparing Alternative Tests of Causality in Temporal Systems: Analytic Results and

I mainly program in Common Lisp and use R for statistical analysis. While in R I miss the power and ease of use of Lisp, especially its many primitives such as find, member, cond, and (perhaps a bridge too far) loop. Has anyone created a package that includes R analogs to a subset of Lisp functions? Chris Elsaesser, PhD Principal Scientist, Machine Learning SPADAC Inc. 7921 Jones

New to R; please excuse me if this is a dumb question. I tried to RTFM; didn't help. I want to do a series of regressions over the columns in a data.frame, systematically varying the response variable and the the terms; and not necessarily including all the non-response columns. In my case, the columns are time series. I don't know if that makes a difference; it does mean I have to call

I need a package that can compute missing values of n-dimensional vectors for n > 2. This is a kind of interpolation, complicated in dimensions higher than 2. The idea is that I have a set of fully specified vectors (i.e., with no missing values) and I get a new vector that has one or more missing attributes; I need to fill in the missing values with values that are, based on the 'training

Hi! I have written an R-package (http://tocsy.agnld.uni-potsdam.de/wavelets/index.html) in R 2.4.1 that requires the ts() function. Users using R 2.5.0 now have a problem installing this package. I checked the package using R 2.5.0: _______________________________________________________ * Installing *source* package 'sowas' ... ** libs gcc -std=gnu99

Dear All, I'm trying to use waldtest to test poolability (parameter stability) between two logistic regressions. Because I need to use robust standard errors (using sandwich), I cannot use anova. anova has no problems running the test, but waldtest does, indipendently of specifying vcov or not. waldtest does not appear to see that my models are nested. H0 in my case is the the vector of

Hi all, I have been working through the examples in one of the vignettes associated with the 'mlogit' package, 'Kenneth Train's exercises using the mlogit package for R.' In spite of using the code unchanged, as well as the data used in the examples, I have been unable to run a Wald test to test two models. Specifically, I have run the following command, where mc and mi2 are

Hello. I am still a newbie in R. Excuse me if I am asking something obvious. My efforts to get an answer through browsing the mailing archives failed. I want to perform an augmented Dickey-Fuller test and to obtain AIC and BIC and to be able to impose some linear restrictions on the ADF regression so as to decide the correct order of autoregression. However I could find no obvious way to impose

Here's a more general question following up on the specific question I asked earlier: Can anybody recommend an R command other than mle.aic() (from the wle package) that will give back a ranked list of submodels? It seems like a pretty basic piece of functionality, but the closest I've been able to find is stepAIC(), which as far as I can tell only gives back the best submodel, not a

Hi, I am trying to run a panel regression using 88 observations and 9 variables. In-built Hausman Test did not work, then I found a code for auxiliary regression method for the Hausman test. The panel models are: fe=plm(gd ~ l+g+o+c+g1+h+n+r, model = "within", data = new.frame,index = c("id")) re=plm(gd ~ l+g+o+c+g1+h+n+r, model = "random", data = new.frame,index =

Hello all, I am trying to figure out an optimal linear model by using stepwise regression which requires partial f-test, I did some Googling on the Internet and realised that someone seemed to ask the question before: Jim Milks <jrclmilks at joimail.com> writes: > Dear all: > > I have a regression model that has collinearity problems (between > three regressor variables). I

Dear R users, I am using lme and nlme to account for spatially correlated errors as random effects. My basic question is about being able to correct F, p, R2 and parameters of models that do not take into account the nature of such errors using gls, glm or nlm and replace them for new F, p, R2 and parameters using lme and nlme as random effects. I am studying distribution patterns of 50 tree

Hi, I know that is not possible make a stepwise procedure using REML in R, I can use ML for this. For nested design it may be very dangerous due the difference in variance structure, mainly in a splitplot design. ML make significative variables that REML dont make. I read an article that is made a stepwise procedure using GENSTAT. from article: "Terms were dropped from a model in a

Dear R-helpers, I have run the code below which I expected to make an object called dd1, but that object does not exist. So, in summary, my problem is that my function is meant to make an object (dd1), and it does indeed make that object (I know that the last line of the function prints it out) but then, after the function has run, the object has disappeared. It's late on a Friday so I may

Hey, I'm trying to implement a GARCH model with Johnson-Su innovations in order to simulate returns of financial asset. The model should look like this: r_t = alpha + lambda*sqrt(h_t) + sqrt(h_t)*epsilon_t h_t = alpha0 + alpha1*epsilon_(t-1)^2 + beta1 * h_(t-1). Alpha refers to a risk-free return, lambda to the risk-premium. I've implemented it like this: #specification of the model

Hi all-- We are planning an intervention study for adolescent alcohol use, and I am planning to use simulations based on a hurdle model (using the hurdle() function in package pscl) for sample size estimation. The simulation code and power code are below -- note that at the moment the "power" code is just returning the coefficients, as something isn't working quite right. The

I am calculating cox propotional hazards models with the coxph function from the survival package. My data relates to failure of various types of endovascular interventions. I can successfully obtain the LR, Wald, and Score test p-values from the coxph.object, as well as the hazard ratio as follows: formula.obj = Surv(days, status) ~ type coxph.model = coxph(formula.obj, df) fit =

Hello, I'm moving from R 2.2.1 (Winows XP) to R 2.5.1 and have problems with installing "myfuncs", which worked OK in 2.2.1 R CMD INSTALL myfuns # gives installing to '' ---------- Making package myfuncs ------------ adding build stamp to DESCRIPTION installing R files installing man source files installing indices installing help >>>

Shown below is most of the FORTRAN subroutine named HCLUST. My question concerns the DO loop labeled as '10'. What happened to its CONTINUE statement? I will assume that after FLAG(I)=.TRUE. is executed that control returns to DO 10 I=1,N. Am I correct? Dave ---------------------------- C Initializations C DO 10 I=1,N C We do not initialize MEMBR in order to be able to

Stephen, thanks for your effort. The more appropriate list for this discussion is probably R-devel (as far as I understand it) so I've moved the discussion there. Related topics have already been discussed in the past. Specifically, I remember contributions by Paul Johnson ("rockchalk" package) and John Fox ("effects" and "car" package) as their packages