Re: [R] venn diagramm

2011-05-07 Thread Breheny, Patrick
I've never actually used it with 5 subsets, but the 'venn' function in the gplots package claims to be able to make a Venn diagram with up to 5 subsets. --Patrick From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of khush ..

Re: [R] Overlapping x axes using Lattice

2011-05-03 Thread Breheny, Patrick
I'm not clear on what you're looking for here. Your x-axis is numeric, why are you converting it to a factor? If you keep it numeric, the labels don't overlap. Or perhaps you don't want it to be numeric, in which case why not just change the aspect ratio of the plot until they no longer overl

Re: [R] Axis trouble

2011-05-03 Thread Breheny, Patrick
The expression 0:g_range[2] is not meaningful. The : operator is for integers, while your data is continuous. Likely, you want something along the lines of axis(2, las=1, at=pretty(vecAVG)) ___ Patrick Breheny Assistant Professor Department of Biostatistics Department of

Re: [R] Lasso with Categorical Variables

2011-05-02 Thread Breheny, Patrick
Clementonia, It sounds like you are looking for the group lasso (Yuan & Lin, 2006). There are two packages on CRAN that have implemented this idea: grpreg and grplasso. The syntax of each is similar to lars (in particular requiring a numeric design matrix as produced by model.matrix), except

Re: [R] Putting x-axis in opposite order

2011-04-28 Thread Breheny, Patrick
Try this: plot(xx, yy, type="l", xlim=rev(range(xx))) ___ Patrick Breheny Assistant Professor Department of Biostatistics Department of Statistics University of Kentucky -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On B

Re: [R] changing a specific column name

2011-04-28 Thread Breheny, Patrick
> colnames(d)[ind] <- 'new name' > colnames(d) [1] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width" "new name" -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of pdb Sent: Thursday, April 28, 2011 7:39 AM To: r-help@r-pr

Re: [R] density plot of simulated exponential distributed data

2011-04-27 Thread Breheny, Patrick
There is an extensive statistical literature on how to correct for boundary bias in kernel density estimates. See, for example, An Improved Estimator of the Density Function at the Boundary S. Zhang, R. J. Karunamuni and M. C. Jones Journal of the American Statistical Association Vol. 94, No. 4

Re: [R] glmnet package: penalty.factor option

2011-04-27 Thread Breheny, Patrick
Have you tried: pen <- rep(1,ncol(X)) pen[c(1,5,7)] <- 0 fit <- glmnet(X,y,penalty.factor=pen) -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Brian H. Chen Sent: Wednesday, April 27, 2011 12:01 PM To: r-help@r-project.org Subject