search for: nbest

I have a patch for libspeex, which optimises some of the loops in vq_nbest and vq_nbest_sign that speeds up encoding - my results: test file: 10s wav file at 16000 Hz, mono encoding with wideband --quality 3, --comp 3 machine: PIII-900Mhz, 256MB RAM before: 2.78s after: 2.38s I'm still trying to grasp the code (I'm just a coder, no background in sound processin...

...actory results. library(leaps) data <- read.table('C:/test_plot_sum2.txt', header = TRUE) attach(data) nox <- data[,1] cnt <- data[,2] vars <- data[,3:48] leaps.setup(x = vars, y = nox, wt = cnt, nvmax = 1) leaps(x = vars, y = nox, wt = cnt, method = "Cp", nbest = 2, names = names(vars)) Error in leaps(x = vars, y = nox, wt = cnt, method = "Cp", nbest = 2, : leaps does not allow more than 31 variables; use regsubsets() regsubsets(x = vars, y = nox, weights = cnt, method = "seqrep", nbest = 1, names = names(vars), nvma...

...ls of my data and or the statistics I am trying to do, or may be arising due to the way leaps works internally. Unfortunately, I am not yet savvy enough to tell why. I can say that this statement works (or at least works to the degree I expect): b <- regsubsets(FUND~.,data=all, intercept=TRUE, nbest=1, nvmax=8, really.big=T, method="forward") It does produce a warning: Warning message: In leaps.setup(x, y, wt = wt, nbest = nbest, nvmax = nvmax, force.in = force.in, : 2 linear dependencies found But I can appreciate that and can reduce my dataset to deal with it. However, if I...

...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 = weights, nbest = nbest, nvmax = nvmax, I think the problem was due to inadequate observations, is there any way to perform subset selection in R if the number of observations are less than the predictors? Thanks! Woo...

...ike to use as a training set (to select a model). I would like to repeat this process until all 25 objects are used as a test set. 2. For each of the training sets I would like to run the following code: library(leaps) forward <- regsubsets(Y ~.,data = training, method = "forward", nbest=1) backward <- regsubsets(Y ~.,data = training, method = "backward", nbest=1) stepwise <- regsubsets(Y ~., data = training, method = "seqrep", nbest=1) exhaustive <- regsubsets(Y ~.,data = training, method = "forward", nbest=1) summary(forward) summary(backwa...

...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 + cov.d + sd.P + sd.v.run + RCS + T + mass.fin, data = DATASET, weights = count, nbest = 10, nvmax = 35, method = "exhaustive") I do however encounter the following warning message which I do not understand: > Reordering variables and trying again: Warning messages: 1: 14 linear dependencies found in: leaps.setup(x, y, wt = wt, nbest = nbest, nvmax = nvmax, force.in...

Hi, ? Is there a fix for this issue??? ---> http://permalink.gmane.org/gmane.comp.audio.compression.speex.devel/2959 ? I am seeing the same thing when I compile speex in visualdsp ? These are the errors I get from using the assembly version of vq_nbest: ? ..\..\..\..\algorithms\voice\speex\src\vq.c [Error ea5004] "C:\Users\coder\AppData\Local\Temp\acc22e8547f000\acc22e8547f001.s":482 Syntax Error in : LOOP vq_loopR2= LC0 = P4; syntax error is at or near text '='. Attempting error recovery by ignoring text until the ';' [...

...some ecological datasets using the regsubsets function in the leaps package. The dataset contains 43 predictor variables plus the response (logcount) all in a dataframe called environment. I am implementing it as follows: library(leaps) subsets<-regsubsets(logcount~.,data=environment,nvmax=10,nbest=2,really.big=FALSE,method="exhaustive") ###the subset regression runs fine when i run it as above and i can get all the usual summaries ###The problem comes when i try and get it to output the variance convariance matric by adding the option vcov=TRUE ##When I do that i get the followin...

Hullo, I am trying to use the leaps package, & keep getting the following error: Error in leaps.setup(x, y, wt = wt, nbest = nbest, nvmax = NCOL(x) + int, : y and x different lengths My data set is attached below. I am rather new to R, & would appreciate any help that could be given. Thanks. http://n4.nabble.com/file/n1755941/lr04.txt lr04.txt -- View this message in context: http://n4.nabble.com/Error-in...

R version: 2.10.0 platform: i486-pc-linux-gnu I'm trying to perform model selection from a data.frame object (creatively named "data") using the leaps function, and I run across the following error: > leaps(data[,3:7], data[,1], nbest = 10) Error in leaps.setup(x, y, wt = wt, nbest = nbest, nvmax = NCOL(x) + int, : character variables must be duplicated in .C/.Fortran Here is a sample of what 'data' looks like: ethanol flask batch delabso delgluc delglyc ph 1 0.00 1 0 1.41 0.0 0.7 1 2 0...

...ves Cp values. To mimic the output of Minitab's MTB > BReg 'Fertility' 'Agriculture'-'Infant.Mortality' ; SUBC> NVars 1 5; SUBC> Best 5; SUBC> Constant. I need to go data(swiss) leapa <-leaps(as.matrix(swiss[,-1]),swiss[,1],method="Cp", nbest=5) leapb <-leaps(as.matrix(swiss[,-1]),swiss[,1],method="adjr2", nbest=5) leapc <-leaps(as.matrix(swiss[,-1]),swiss[,1],method="r2", nbest=5) best <- cbind(leapc$r2,leapb$adjr2,leapa$Cp,leapa$which) colnames(best) <- c( "R.Sq", "Adj.R.Sq","Cp...

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" and have a look at it, but the value of the adjusted R squared is not the same as in the regsubsets output? This is also the case for the other models,...

...orithm returns a best model of each size, the results do not depend on a penalty model for model size: it doesn't make any difference whether you want to use AIC, BIC, CIC, DIC, ... USAGE leaps(x=, y=, wt=rep(1, NROW(x)), int=TRUE, method=c("Cp", "adjr2", "r2"), nbest=10, names=NULL, df=NROW(x), strictly.compatible=TRUE) ARGUMENTS x A matrix of predictors y A response vector wt Optional weight vector int Add an intercept to the model method Calculate Cp, adjusted R-squared or R-squared nbest Number of subsets of each si...

...he computer is hung. I ran regsubsets on a smaller portion of the data, also with linear dependency problems, and got results. However, the hourglass continues its endless spiraling with the full dataset. I am running the following on Windows 7 library(leaps) m_250<-regsubsets(Y~., data=model2, nbest=1, really.big=TRUE) (NOTE: The ~ is a tilda, not a dash, in the regression statement above: Y~.) Does anyone have any opinions on: 1) is R likely to still be running, even after a week, or should i just shut it down? 2) am i doing something wrong with regsubsets? 3) is there a better option th...

...; shape_cb_size entries and subvect_size samples per entry. Ok, that helps. Thanks. >> PS: for me personally, I'd like to know how to change the cb id found >> while keeping the error introduced by this low. > > That's for steganography? I'd say you just look at the nbest list and > chose the 2nd, 3rd, ... best entry instead of the first. Yes. I already introduced into the LSP prediction such an nbest list to do that straightforward. My goal is to make sure if nothing is embedded, the altered Speex version still needs to write exactly the same coefficients as th...

Hi, Jean-Marc Valin wrote: > Well, you could change the order in the encoder as long as you reverse > it in the decoder as well. Ok, I see that in the split_cb_shape_sign_unquant function, that each coefficient is tied to it's position in the nb_subvect exc coefficients. Honestly, I have problems understanding what exactly the codebook search works like. If you have the time to to

Hi, List , i have 6 predictor variables and i want to make possible combinations of these 6 predictors ,all the data is in matrix form , if i am having 6 predictors than possible combination of sets are 64 2 power 6, or 63 ,whatever it may be i want to store the result in another variable to each combination and that i want to put in some model , i want to put every combination in some model

...es and subvect_size samples per entry. > > Ok, that helps. Thanks. > >>> PS: for me personally, I'd like to know how to change the cb id found >>> while keeping the error introduced by this low. >> That's for steganography? I'd say you just look at the nbest list and >> chose the 2nd, 3rd, ... best entry instead of the first. > > Yes. I already introduced into the LSP prediction such an nbest list to > do that straightforward. My goal is to make sure if nothing is embedded, > the altered Speex version still needs to write exactly the...

Hi, I'm using regsubsets from the leaps package to select subsets of variables. I'm calling the function as lp <- regsubsets(x,y,nbest=5,nvmax=9) Then I call plot to see which variables turned up in the models. I use the R^2 scale and see my best model had a R^2 of 0.62. However when I make a linear model using lm() with the same x my R^2 is 0.45. Should'nt I be seeing the same value of R^2? I must be making a mistake somew...

...ive search). The package leaps can help me to generate all possible variable subsets, but I’ll appreciate your guidance as of how to generate one model for each of those possible subsets. I’m new in R! Thanks in advance for your help. library(leaps) subsets <- regsubsets( y ~ . , data=freeny, nbest=100, nvmax=5, method="exhaustive", really.big=T) subset.models <- summary(subsets)$which now fit one lm for each possible combination of variables. [[alternative HTML version deleted]]