similar to: all subsets for glm

Hi, I tried several settings by using the "family=gaussian" in "gl1ce", but none of them works. For the case "glm" can work. Here is the error message I got: > glm(Petal.Width~Sepal.Length+Sepal.Width+Petal.Length ,data=iris,family=gaussian()) > gl1ce(Petal.Width~Sepal.Length+Sepal.Width+Petal.Length ,data=iris,family=gaussian()) Error in eval(expr, envir,

Hi, I'm using gl1ce with family=binomial like so: >yy succ fail [1,] 76 23 [2,] 32 67 [3,] 56 43 ... [24,] 81 18 >xx c1219 c643 X1 0.04545455 0.64274145 X2 0.17723669 0.90392792 ... X24 0.80629054 0.12239320 >test.gl1ce <- gl1ce(yy ~ xx, family = binomial(link=logit), bound = 0.5 ) or >omit <- c(2,3) >test.gl1ce

Hi all, I want to know some details about the c routine “lasso” in the functions of “gl1ce()” . However, I have following troubles. First, I can not find the routine in the local directories of this function (or package). Second, if I found the routine, could I call it just like this way, say, fit <- .C("lasso", …,PACKAGE = "lasso2") in my own functions. My system is

Dear all, I was recently contacted by a user about an alledged problem/bug in the latest version of lasso2. After some investigation, we found out that it was a user error which boils down to the following: > x <- matrix(rnorm(200), ncol=2) > var <- "fred" > apply(x, 2, var) Error in get(x, envir, mode, inherits) : variable "fred" of mode "function"

X is a matrix and F is a vector. F2 <- data.frame(cbind(X,F)) F2 V1 V2 V3 F 1 -0.250536332 -1.4755883 1.9580974 -2.136487 2 -0.009856084 0.4953269 0.5486092 -2.744482 3 -0.406962682 0.7729631 0.1861905 -2.891821 4 1.938780097 0.7469251 1.2537781 -1.212992 5 -0.332370358 1.1943637 0.7114278 -1.830441 modF<-formula(F ~ V1 + V2 + V3) #no error message

Hi, Is there an R implementation of least angle regression for binary response modeling? I know that this question has been asked before, and I am also aware of the "lasso2" package, but that only implements an L1 penalty, i.e. the Lasso approach. Madigan and Ridgeway in their discussion of Efron et al (2004) describe a LARS-type algorithm for generalized linear models. Has

Full_Name: Grant Izmirlian Version: 2.2.0 OS: SuSe Linux version 9.2 Submission from: (NULL) (156.40.34.177) Sorry about the last submission, my bug-fix had an error in it because ifelse doesn't vectorize. I'll repost with the correct bug-fix. ------------------------------------------------------------------------------- The option exists to include all parameters, including the

Dear R community: My platform: R 1.7.0 + windows2000. I am trying to install the package "lasso2" which I saw in the following web address: http://cran.us.r-project.org/src/contrib/PACKAGES.html#emplik. However, I failed to install it from R menu "Packages| Install package(s) from CRAN" since I could not find this item in the list. Thanks in advance! Rui [[alternative

Dear R community: My platform: R 1.7.0 + windows2000. I am trying to install the package "lasso2" which I saw in the following web address: http://cran.us.r-project.org/src/contrib/PACKAGES.html#emplik. However, I failed to install it from R menu "Packages| Install package(s) from CRAN" since I could not find this item in the list. Thanks in advance! Rui [[alternative

There are some googlable ways to feed a list of filenames to vim, but I stumble on weird results. With my filelist, I try to do cat list | xargs vim ...to edit the files listed in the file "list". Here's what happens: [root at lasso2 tempdir]# ls -l total 8 -rw-r--r-- 1 root root 0 May 17 18:28 a -rw-r--r-- 1 root root 0 May 17 18:28 b -rw-r--r-- 1 root root 3 May 17

Hi , I want to assign names to the rows in my matrix so that when I use the 'agnes' function from R cluster package , the dendogram that is produced represents the rows of the matrix. This way I would know what elements in the matrix are clustered together. But when I do the following, rownames(distF) <- paste("A",dim(distF)[1]) it gives the following error; Error in

Hi, I'm trying to calculate euclidean distance for my matrix rx1. Its a big matrix with 947 elements. However when i do the following, for(k in 1:947) for(m in (k+1):947){ A =rbind(A,c(k,m,ED(rx1[,k], rx1[,m]))) } I get the following error: ED(rx1[,k], rx1[,m]) : subscript out of bounds I checked for k in 1:100 and it runs fine; its a huge matrix and takes atleast an hour to run before I

Hello, I am trying to shrink the coefficients of a logistic regression for a sparse dataset, I am using the lasso (lasso2) and I am trying to determine the shrinkinage factor by cross-validation. I would like please some of the experts here to tell me whether i'm doing it correctly or not. Below is my dataset and the functions I use w= a b c d e P A 0 0 0 0 0 1 879 1 0 0 0 0 1 3 0 1 0 0 0 7 7

i have a matrix with the following dimensions 136 3 and it looks something like [,1] [,2] [,3] [1,] 402 675 1.802758 [2,] 402 696 1.938902 [3,] 402 699 1.994253 [4,] 402 945 1.898619 [5,] 424 470 1.812857 [6,] 424 905 1.816345 [7,] 470 905 1.871252 [8,] 504 780 1.958191 [9,] 504 848 1.997111...............

I have had problems like this before. Probably there is something important that I don't know about routing. Let me introduce to you "Lasso2", a CentOS 4 www server that has been working perfectly well for years. Now I added a second nw card (eth1), automatically using kudzu. I cannot get this dual nw setup to work. The first nw card (eth0) stopped at once working properly,

Dear all, I am quite puzzled about the bound and absolute.t arguments to the l1ce function in the lasso2 package. (The l1ce function estimates the regression parameter b in a regression model y=Xb+e subject to the constraint that |b|<t for some value t). The doc says: bound numeric, either a single number or a vector: the constraint(s) that is/are put onto the L1 norm of the parameters.

In the package lasso2, there is a Prostate Data. To find coefficients in the prostate cancer example we could impose L1 constraint on the parameters. code is: data(Prostate) p.mean <- apply(Prostate, 5,mean) pros <- sweep(Prostate, 5, p.mean, "-") p.std <- apply(pros, 5, var) pros <- sweep(pros, 5, sqrt(p.std),"/") pros[, "lpsa"] <-

Dear R-help, I am having trouble with the predict function in lasso2. For example: > data(Iowa) > l1c.I <- l1ce(Yield ~ ., Iowa, bound = 10, absolute.t=TRUE) > predict (l1c.I) # this works is fine > predict (l1c.I,Iowa) Error in eval(exper,envir, enclos) : couldn't find function "Yield" And I have similar trouble whenever I use the newdata argument in

Hello list, I've been searching around trying to find whether somebody has written such a package of least angle regression on generalized linear models, like what Lasso2 package does. The extension to generalized linear models is briefly discussed in the comment by D. Madigan and G. Ridgeway. Is such a package available? Thanks, Mike [[alternative HTML version deleted]]

All, Could anyone recommend a package that allows the user to constrain the coefficients from a multiple regression equation? I tried using the gl1ce function in lasso2, but couldn't get it to work. I created a contrived example to illustrate my starting point. data(cars) fmla <- formula(dist ~ speed) gl1c.E <- gl1ce(fmla, data = cars) gl1c.E gl1c.E <- gl1ce(fmla, data =