similar to: Bandwidth selection for kernel regression.

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

Dear R Users, Before running ksmooth( ), a suitable bandwidth selection is needed. I use some functions for this task and receive these results for my data: width.SJ(y,nb=100,method="ste") : 40.25 bcv(y,nb=100) : 40.53 ucv(y) : 41.26 bandwidth.nrd(y) : 45.43 After implementing the function ksmooth(x,y, bandwidth= each of abovementioned bandwidths), I have some NAs

Hello! I do not understand what is meant by: "aic" and "bic" follow a top-down strategy based on the Akaike's and Schwarz's information criteria in the datails to the ADF.test function. What does a "top-down strategy" mean? Probably the respective criterion is minimized and the mode vector contains the lag orders at which the criterion attains it

Apologies in advance if my comments don't help, in which case, no need to respond, but I noted in ?ksmooth: "bandwidth the bandwidth. The kernels are scaled so that their quartiles (viewed as probability densities) are at ? 0.25*bandwidth." So, could this be a source of the discrepancies you cited? Given that ?ksmooth explicitly says: "Note: This function was implemented for

Dear Sir, Madam, or to whom this may concern, my name is Jan Failenschmid and I am a Ph.D. student at Tilburg University. For my project I have been looking into different types of kernel regression estimators and corresponding R functions. While comparing different functions I noticed that stats::ksmooth returned different estimates for the same bandwidth as other kernel regression estimators

Hello all, i hope you don't mind my off topic question. i want to use the Akaike criterion for variable selection in a regression model. Does anyone know some basic literature about that topic? Especially I'm interested in answers to the following questions: 1. Has (and if so how has) the criterion to be modified, if i estimate the transformations of the variables too? 2. How is the

Full_Name: Ommo H?ppop Version: 2.0.1 OS: XP Submission from: (NULL) (84.143.196.187) Hi, Presumably, I've found an error in http://finzi.psych.upenn.edu/R/doc/manual/R-intro.html In Chapter 11.3 Generic functions for extracting model information it says "Select a suitable model by adding or dropping terms and preserving hierarchies. The model with the LARGEST value of AIC

Dear All 1) I'm fitting spatial CAR models using R2Winbugs and although everything seems to go reasonably well (or I think so) the next message appears from WINBUGS 1.4 window: gen.inits() Command #Bugs: gen.inits cannot be executed (is greyed out) The question is if this message means that something is wrong and the results are consequently wrong, or Can I assume it as a simple warning

Hello all, I would like to use the Nadaraya-Watson estimator assuming a Gaussian kernel: So far I sued the library(sm) library(sm) x<-runif(5000) y<-rnorm(5000) plot(x,y,col='black') h1<-h.select(x,y,method='aicc') lines(ksmooth(x,y,bandwidth=h1)) which works fine. What if my data were clustered requiring a bandwidth that varies with x? How can I do that? Thanks in

Hello all, I have been using R's time series capabilities to perform analysis for quite some time now and I am having some questions regarding its reliability. In several cases I have had substantial disagreement between R and other packages (such as gretl and the commercial EViews package). I have just encountered another problem and thought I'd post it to the list. In this case,

Dear R users (r-help@r-project.org), The definition of AIC (Akaike information criterion) for normal error models has just been changed. Please refer to the paper below on this matter. Eq.(22) is the new definition. The essential part is RSS(n+q+1)/(n-q-3); it is close to GCV. The paper is temporarily available at the "Papers In Press" place. Kunio Takezawa(2012): A Revision of

I am hoping to get some advise on the following: I am looking for an automatic variable selection procedure to reduce the number of potential predictor variables (~ 50) in a multiple regression model. I would be interested to use the forward stepwise regression using the partial F test. I have looked into possible R-functions but could not find this particular approach. There is a function

Dear Listers, I got this message when compiling a package: * creating pgirmess-manual.tex ... OK * checking pgirmess-manual.text ... ERROR LaTex errors when creating DVI version. This typically indicates Rd problems. The message is quite explicit but I struggled a lot before understanding that the trouble comes from a single file "selMod.rd" among 44 topics. Even though I have

Hi, I have a problem with a regression I try to run. I did an estimation of the market model with daily data. You can see to output below: /> summary(regression_resn) Time series regression with "ts" data: Start = -150, End = -26 Call: dynlm(formula = ror_resn ~ ror_spi_resn) Residuals: Min 1Q Median 3Q Max -0.0255690 -0.0030378 0.0002787

I'd like to superimpose a ksmooth() onto a barplot(). My data is: > d 2009-06-20 2009-06-21 2009-06-22 2009-06-23 2009-06-24 2009-06-25 2009-06-26 2009-06-27 2009-06-28 2009-06-29 2009-06-30 2009-07-01 2009-07-02 Same Breed (B) 12.64 21.08 13.52 12.51 13.71 9.91 14.24 7.18 11.81 5.92 12.04 17.96

I am working with a model that I have to estimate a nonparametric function. The model is partial linear i.e. Y=X$\beta$ + f(z) + $\epsilon$ I am using the ' double residual methods' Robinson (1988) Speckman (1988) where I estimate a nonparametric function for each of the parametric variables in terms of the nonparametric one i.e. X[,i]=g(Z)+ u this is done because I need the $E(

Hello. I have some further problems with modelling an optimization problem in R: How can I model some optimization problem in R with a linear objective function with subject to some nonlinear constraints? I would like to use "optim" or "constrOptim", maybe with respect to methods like "Simulated Annealing" or "Sequential Quadric Programming" or something

Hi There I have a lm object with 4 parameters and I want to test wether 2 parameters are equal using a Wald test (basically b1=b2 or b1-b2 =0). In the help file from R it says that under ANOVA the optional arguments " Terms" or "L" test whether a linear combination is equal to 0. I tried; >anova(m1, Terms = Beta1-Beta2=0) but I get the error: Object " Beta1"

When I do ARMA(2,2) using one lag of LCPIH data This is eview result > > *Dependent Variable: DLCPIH > **Method: Least Squares > **Date: 08/12/11 Time: 12:44 > **Sample (adjusted): 1970Q2 2010Q2 > **Included observations: 161 after adjustments > **Convergence achieved after 14 iterations > **MA Backcast: 1969Q4 1970Q1 > ** > **Variable Coefficient Std.

All: I am trying to use gam in a scatterplot smoothing problem. The data being smoothed have greater 1000 observation and have multiple "humps". I can smooth the data fine using a function something like: out <- ksmooth(x,y,"normal",bandwidth=0.25) plot(x,out$y,type="l") The problem is when I try to fit the same data using gam out <-