Hi Samuel, The help file explains which is which:
A list with components devcvm Average drop in CV deviance for each lambda value ncallcvm=ncallcvm Average number of features with non-zero wts in the CV, for each lambda value se.devcvm Standard error of average drop in CV deviance for each lambda value devcv Drop in CV deviance for each lambda value ncallcv Number of features with non-zero wts in the CV, for each lambda value folds Indices for CV folds call Call to this function So the following code will do: > plot(aa[[2]], aa[[1]], type="b", ylim=c(0,25)) > segments(aa[[2]], aa[[1]]-aa[[3]], aa[[2]], aa[[1]]+aa[[3]]) Can't figure out how he made the horizontal bars at the ends of each error bar vary in length. How do you like uniCox? Have seen my previous post regarding the error message? https://stat.ethz.ch/pipermail/r-help/2010-November/260857.html ...Tao ________________________________ To: jim holtman <jholt...@gmail.com> Cc: r-help@r-project.org Sent: Wed, December 8, 2010 5:59:40 AM Subject: Re: [R] UniCox in R Thank you for your reply. My question was how to create Fig. 2? I am using UniCox; the aa=uniCoxCV gives list of aa$ devcvm aa$ncallcvm aa$se.devcvm, aa$devcv, aa$ ncallcv and aa$ folds, which data I should plot here to create the graph in Fig. 2? Many thanks, Sam --- On Wed, 8/12/10, jim holtman <jholt...@gmail.com> wrote: From: jim holtman <jholt...@gmail.com> Subject: Re: [R] UniCox in R Cc: r-help@r-project.org Date: Wednesday, 8 December, 2010, 5:25 Do a little reading on how to use the graphics commands. I would look at plot lines segments a combination of those will let you easily create the output. It would seem you need some data object that has the dimensions of the bar lengths you want to create, but given that, it is not much of a problem to solve. > > > > > Hello, > > I am interested in Figure 2 in > > http://www-stat.stanford.edu/~tibs/ftp/cus.pdf > > Can anyone tell please how to create this plot? > > Many thanks > Samuel > > > > > > > [[alternative HTML version deleted]] > > ______________________________________________ > 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. > -- Jim Holtman Data Munger Guru What is the problem that you are trying to solve? [[alternative HTML version deleted]] [[alternative HTML version deleted]] ______________________________________________ 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.
