similar to: lasso constraint

Hi, First, when I try the example Prostate with bound 0.44 (as in the manual), I got a different result: > l1c.P <- l1ce(lpsa ~ ., Prostate, bound=0.44) > l1c.P .... Coefficients: (Intercept) lcavol lweight age lbph svi 1.0435803 0.4740831 0.1953156 0.0000000 0.0000000 0.3758199 lcp gleason pgg45 0.0000000 0.0000000

All, For my understanding, I wanted to see if I can get glmnet predictions using both the predict function and also by multiplying coefficients by the variable matrix. This is not worked out. Could anyone suggest where I am going wrong? I understand that I may not have the mean/intercept correct, but the scaling is also off, which suggests a bigger mistake. Thanks for your help. Juliet Hannah

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.

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

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,

In an .Rd example for a package, I want to use data from another package, but avoid loading the entire package and avoid errors/warnings if that other package is not available. If I don't care about loading the other package, I can just do: if (require("ElemStatLearn", quietly=TRUE)) { data(prostate) # rest of example } I'd rather just be able to do something like:

Hi Everyone, I have a question regarding the construction of 3D graphs in 'R', BUT these graphs also need to illustrate movement (with time) of the prostate gland (using radiological techniques). I am not sure how to do this in 'R' although I'm sure there is some way of doing it. Below, I have copied and pasted some of the data with which I'm working on. The data

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

Describe fails for me with a message similar to what was an issue in 2008 and got fixed according to posts. R version 2.15.0 (2012-03-30) Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i386-pc-mingw32/i386 (32-bit) # output truncated > options(chmhelp = FALSE, help_type = "text") > .help.ESS <- help >

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

Dear all, I did post this more or less identical mail in a follow up to another question I posted, but under another heading. I try again, but now under the correct header. upon running this code (from the Hmisc library-latex function) I believe the call to summary.formula is allright and produces wonderful tables, but the latex command results in a correct formatted table but where all the

Dear Asterisk Community, Does your company provide inbound 800# origination? If so, please read this message and e-mail us a quote for monthly co-lo hosting of our asterisk server and per-minute inbound 800# origination. The Prostate Cancer Research and Education Foundation (PC-REF) is a non-profit organization dedicated to helping prostate cancer sufferers and their loved ones. We have

I have the following files list: > list.files() [1] "Prostate-Cancer_cvs_Dir" [2] "Prostate_Cancer-miRNAs&Genes.Pathway.xml" [3] "Prostate_Cancer_Pathways-miRNAs-GeneTargets-Dir" [4] "Prostate_Cancer_Pathways-miRNAs-GeneTargets-Dir.zip" [5] "Prostate-miRNAs.OrganTargets.txt"

Full_Name: Kem Phillips Version: 2.7.1 (2008-06-23) OS: Windows Xp professional Submission from: (NULL) (98.221.200.108) The Hmisc function describe fails, giving the error message: Error in formatDateTime(dd, atx, !timeUsed) : could not find function "format.dates" Loading the chron package, where function dates apparently resides, does not fix the problem. Note

Hi, I am working on a project about hospital efficiency. Due to the high multicolinearlity of the data, I want to fit the model using ridge regression. However, I believe that the data from large hospital(indicated by the number of patients they treat a year) is more accurate than from small hosptials, and I want to put more weight on them. How do I do this with lm.ridge? I know I just need

Dear All, I'm reading Frank Harrell's wonderful Regression Modeling Strategies book and ran into a problem following the example in Chapter 8. I'm working on platform: Ubuntu 8.04 (i486-pc-linux-gnu) R version: 2.7.2 (2008-08-25) and my command sequence was: library(chron) library(Hmisc) load("prostate.sav") describe(prostate) The last command returned the error

Hi, I've lost my mind on it... I have to scatterplot two vectors, grouped by a third variable, with two different dimensions according to whether each cell line in the plot is sensitive or resistant to a given drug, and with a different color for each of 9 tissues of origin. Here's what I've done:

Hi everybody, this is a real dummy thing. I sorted a matrix based on a given column, and what I get is right, until it comes to columns of negative and positive values; than, "order" orders everything from max to min in the negative values, and then AGAIN from max to min in the positive values!!! Why isn't everything order from max to min, and that's it? Thank you!!! Attached

Sir, I found your description of the dataset about nodal involvement in prostate cancer. It comes from the book biostatistics casebook. I like to use the dataset for doing logistics regression. Can you tell me where I can find the dataset. Thanks and greetings Wim van Baarle wvbaarle@wxs.nl [[alternative HTML version deleted]]

Dear Rexperts, I am trying to add a '+' identifying the mean in a boxplot using the following sizelist <- split(size, grp) centers <- boxplot(sizelist, style.bxp = "att", medpch = "o", ylab = "Prostate Volume (cm3)") points(centers, unlist(lapply(sizelist, mean)), pch = "+") But, I get error Error in xy.coords(x, y) :