similar to: regsubsets selection criterion

Hello, I am analysing several samples whose sizes are from 9 to 110. I would like to test their distribution with R, whether they are normal or not. I wonder which test for normality from R should I use . Thank you for help. Samuel Samuel BERTRAND Doctorant Laboratoire de Biomecanique LBM - ENSAM - CNRS UMR 8005 151, bd de l'Hopital 75013 PARIS Tel. +33 (0) 1 44 24 64 53 Fax +33 (0) 1

Hi, I want to make model selection with regsubsets. My code is: a<-regsubsets(Gesamt ~ CommunistSocialist + CountrySize + GNI + Lifeexp + Schoolyears + ExpMilitary + Mortality + PopPoverty + PopTotal + ExpEdu + ExpHealth, data=olympiadaten, nbest=2) summary(a) plot(a,scale="adjr2") (output attached) The problem is now, that I want to fit the best model again "manually"

Dear R-help, I have used the regsubsets function from the leaps package to do subset selection of a logistic regression model with 6 independent variables and all possible ^2 interactions. As I want to get information about the statistics behind the selection output, I?ve intensively searched the mailing list to find answers to following questions: 1. What should I do to get the statistics

Hi, I have 12 independent variables and one dependent variable. Now I want to select the best adj. R squared model by using the regsubsets command, so I code: > plot(regsubsets(Gesamt ~ CommunistSocialist + CountrySize + GNI + Lifeexp + Schoolyears + ExpMilitary + Mortality + + PopPoverty + PopTotal + ExpEdu + ExpHealth, data=olympiadaten, nbest=1, nvmax=12), scale='adjr2') Then I

Hi, I have found a minor problem with leaps(). In 1.3.1 under Windows 2000 I seem to only be able to obtain values for one statistic at a time. That is choosing method=c("Cp","adjr2","r2") just gives Cp values. To mimic the output of Minitab's MTB > BReg 'Fertility' 'Agriculture'-'Infant.Mortality' ; SUBC> NVars 1 5; SUBC>

Hi i need to create a model from 250 + variables with high collinearity, and only 17 data points (p = 250, n = 750). I would prefer to use Cp, AIC, and/or BIC to narrow down the number of variables, and then use VIF to choose a model without collinearity (if possible). I realize that having a huge p and small n is going to give me extreme linear dependency problems, but I *think* these model

Colleagues, please help me with the simple question. How can I find R2 and p while doing best subsets regression? Also how can I see B and p for coefficients? Maybe you can advice me detailed manual, because in manual for the package this is not described. Thank you

Hello, I have a question about doing ALL possible subsets regression with a general linear model. My goal is to produce cumulative Akaike weights for each of 7 predictor variables-to obtain this I need R to: 1. Show me ALL possible subsets, not just the best possible subsets 2. Give me an AIC value for each model (instead of a BIC value). I have tried to do this in library(RcmdrPlugin.HH),

Hello I am trying to use all subsets regression on a test dataset consisting of 11 trails and 46 potential predictor variables. I would like to use Mallow's Cp as a selection criterion. The leaps function would provide the required output but does not work with this many variables (see below). The alternative function regsubsets should be used, but I am not able to define the function in

Hi, I'm new to R and this mailing list, so I will attempt to state my question as appropriately as possible. I am running R version 2.7 with Windows XP and have recently been exploring the use of the function regsubsets in the leaps package in order to perform all-subsets regression. So, I'm calling the function as:

Hello dear R-help members, I recently became interested in using biglm with leaps, and found myself somewhat confused as to how to use the two together, in different settings. I couldn't find any example codes for the leaps() package (except for in the help file, and the examples there are not as rich as they could be). That is why I turn to you in case you could share some good tips and

I'm interested in writing some model selection functions (for linear regression models, as a start), which incorporate the PRESS criterion since it, to my knowledge, is not currently implemented in any available model selection procedure. I thought it would be simplest to build on already existing functions like regsubsets in package leaps. It's easy enough to calculate the PRESS

Dear List members I am using the regsubsets function to select a few predictor variables using Mallow's Cp: > sel.proc.regsub.full <- regsubsets(CO2 ~ v + log(v) + v.max + sd.v + tad + no.stops.km + av.stop.T + a + sd.a + a.max + d + sd.d + d.max + RPA + P + perc.stop.T + perc.a.T + perc.d.T + RPS + RPSS + sd.P.acc + P.dec + da.acc.1 + RMSACC + RDI + RPSI + P.acc + cov.v + cov.a +

I am generating a bunch of examples for my students and when I type R BATCH demo10.R (file demo10.R reproduced below) I get an error: ------------- Error message ------ > plot(mods$size,mods$Cp, + main="Cp versus talla modelos", + xlab=expression(p), + ylab=expression(C[p])) Error in plot.window(xlim, ylim, log, asp, ...) : need finite ylim values In

Hi, I'd like to run regsubsets for model selection by exhaustive search. I have a list with 20 potential explanatory variables, which represent the real and the imaginary parts of 10 "kinds" of complex numbers: x <- list(r1=r1, r2=r2, r3=r3, ..., r10=r10, i1=i1, i2=i2, i3=i3, ..., i10=i10) Is there an easy way to constrain the model search so that "r"s and

Hi, I would like to perform a forward selection procedure on a data set with 6 observations and 10 predictors. I tried to run it with regsubsets (I set nvmax=number of observations) but I keep getting these warning messages: Warning messages: 1: 5 linear dependencies found in: leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, 2: nvmax reduced to 5 in: leaps.setup(x, y, wt =

Hi, does anyone know if there is a way to easily extract p-values from the regsubsets() function? Thanks, James Stegen -- James C. Stegen NSF Postdoctoral Fellow in Bioinformatics University of North Carolina Chapel Hill, NC 919-962-8795 stegen at email.unc.edu http://www.unc.edu/~stegen/index.html

Hi, does anyone know if there is a way to easily extract p-values from the regsubsets() function? Thanks, James Stegen p.s. this is a reposting due to me not putting in a useful subject heading... -- James C. Stegen NSF Postdoctoral Fellow in Bioinformatics University of North Carolina Chapel Hill, NC 919-962-8795 stegen at email.unc.edu http://www.unc.edu/~stegen/index.html

How would I rewrite my code so that I can implement the use of multicore on an Rstudio server to run regsubsets using the "exhaustive" method? The data has 1200 variables and 9000 obs so the code has been shortened here: model<-regsubsets(price~x + y + z + a + b + ...., data=sample, nvmax=500, method=c("exhaustive")) Our server is a quad core 7.5 gb ram, is that