Mike: On Sat, May 25, 2013 at 7:24 PM, C W <tmrs...@gmail.com> wrote: > Thomas, thanks for the cool trick. I always thought browser() was the > only thing existed, apparently not.
Which you would have known had you read the docs! See section 9 on Debugging in the R Language Definition Manual -- Bert > > Mike > > On Sat, May 25, 2013 at 9:09 PM, Thomas Stewart > <tgs.public.m...@gmail.com> wrote: >> Mike- >> >> You can use the traceback function to see where the error is: >> >>> bob <- matrix(rnorm(100*180), nrow=180) >>> yyy <- rnorm(180) >>> fit1 <- cv.glmnet(bob, yyy, family="mgaussian") >> Error in rep("(Intercept)", nrow(a0)) : invalid 'times' argument >>> traceback() >> 6: predict.multnet(object, newx, s, type, exact, offset, ...) >> 5: predict.mrelnet(glmnet.object, type = "nonzero") >> 4: predict(glmnet.object, type = "nonzero") >> 3: lapply(X = X, FUN = FUN, ...) >> 2: sapply(predict(glmnet.object, type = "nonzero"), length) >> 1: cv.glmnet(bob, yyy, family = "mgaussian") >> >> So, thee error is in the predict.multnet function. If you peak at that >> function, you see where the function falls apart. It seems that the >> function wants a0 to be a matrix but in this example it is a vector. I'm >> not familiar enough with the package to offer advice on how to fix this. >> >> -tgs >> >> >> >> On Sat, May 25, 2013 at 4:14 PM, C W <tmrs...@gmail.com> wrote: >>> >>> Dear list, >>> I am using glmnet. I have no idea what this error is telling me. >>> Here's my code, >>> >>> > bob <- matrix(rnorm(100*180), nrow=180) >>> > yyy <- rnorm(180) >>> > fit1 <- cv.glmnet(bob, yyy, family="mgaussian") >>> Error in rep("(Intercept)", nrow(a0)) : invalid 'times' argument >>> >>> In fact, I peeked inside cv.glmnet() using, >>> > glmnet:cv.glmnet >>> >>> Can't even find the error message in the code. I am clueless at the >>> moment. >>> Thanks in advance, >>> Mike >>> >>> ______________________________________________ >>> R-help@r-project.org mailing list >>> https://stat.ethz.ch/mailman/listinfo/r-help >>> PLEASE do read the posting guide >>> http://www.R-project.org/posting-guide.html >>> and provide commented, minimal, self-contained, reproducible code. >>> >>> >> > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.