Re: [R] (performance) time in Windows vs Linux

2009-06-30 Thread Zeljko Vrba
On Tue, Jun 30, 2009 at 10:20:17AM +0900, Raymond Wan wrote: I. Soumpasis wrote: 2009/6/29 C騷ar Freitas cafanselm...@yahoo.com.br This is true. So I tried the same computer with windows XP and ubuntu 8.10 64bit dual core @3Gz and 4MB RAM Windows 32bit results: user system elapsed

Re: [R] (performance) time in Windows vs Linux

2009-06-29 Thread Zeljko Vrba
On Mon, Jun 29, 2009 at 06:56:55PM +1200, Patrick Connolly wrote: On Mon, 29-Jun-2009 at 02:13AM -0400, milton ruser wrote: | Really? | | In fact I have a quadcore. But how can I know if Linux are really | using only one core, and how can I setup it to use the 4cores? I use GKrellM

Re: [R] (performance) time in Windows vs Linux

2009-06-26 Thread Zeljko Vrba
On Fri, Jun 26, 2009 at 12:23:35PM -0700, Cézar Freitas wrote: I supposed R on Linux should be faster (32 and 64 bit) than windows version. Is this difference because 64 bit R version is slower than 32 bits one? I started the machine in both sittuations and checked free memory. I suspect

Re: [R] confusion on levels() function, and how to assign a wanted order to factor levels, intentionally?

2009-06-16 Thread Zeljko Vrba
On Tue, Jun 16, 2009 at 03:14:01PM +0800, Mao Jianfeng wrote: levels(d$population)-c(YXPy01, KMPy01, YLPy01, GSPy02, BCPy01, LJPy01, GYPt01, YLPd01, CYPd01, CYPd02, CYPd03, BXPd01, NSPt01) I'm not at home with factors myself, but maybe this will do the trick for you: d$population -

Re: [R] GeoXp package

2009-06-10 Thread Zeljko Vrba
On Wed, Jun 10, 2009 at 08:21:06AM +0200, Poizot Emmanuel wrote: Error in fun(...) : GDAL Error 1: libgrass_I.so: Ne peut ouvrir le fichier d'objet partagé: Aucun fichier ou dossier de ce type (sorry for the french :) ) It would have been far more useful had you translated the error

Re: [R] ridiculous behaviour printing to eps: labels all messed up!

2009-06-08 Thread Zeljko Vrba
On Sun, Jun 07, 2009 at 05:00:18PM -0700, maiya wrote: I am completely lost here, can someone help me figure out what is going on here? With an unrelated problem, someone suggested to me never to use dev.copy2eps, but to open a new device with postscript() and draw the figure there. Maybe

Re: [R] Installation of R in Redhat Linux (64 bit)

2009-06-05 Thread Zeljko Vrba
On Fri, Jun 05, 2009 at 08:27:45PM -0700, Santosh wrote: When I run the following: nm /usr/lib64/R/library/stats/libs/stats.so | grep _gfortran_copy_string You probably need to install libraries related to gfortran. __ R-help@r-project.org mailing

[R] ggplot2: How to export several plots with same width?

2009-06-01 Thread Zeljko Vrba
I have three plots and I want the *plot area* to be of the same width on each plot. Since the three plots have different legends, the text width of the legend affects the width of the plot area (longer legend text = narrower plot area). Exporting the three figures to postscript device of same

Re: [R] ggplot2: How to export several plots with same width?

2009-06-01 Thread Zeljko Vrba
On Mon, Jun 01, 2009 at 10:01:40AM +0200, baptiste auguie wrote: I'm not sure it's currently possible with ggplot2 (lattice and latticeExtra offer some workarounds for this). Thanks for the links. I ended up reshaping my dataset and applying faceting.

Re: [R] R's documentation

2009-05-30 Thread Zeljko Vrba
On Fri, May 29, 2009 at 05:20:24PM +0100, Patrick Burns wrote: If you find some documentation that is confusing, then you can write a message about it that states: I think that some kind of a glossary would be helpful. Then I would know whether certain words or phrases are R-specific or

Re: [R] Unable to load R

2009-05-28 Thread Zeljko Vrba
On Thu, May 28, 2009 at 11:42:41AM +0530, anupam sinha wrote: I have checked for the presence of the above mention library and found that the library is present. I have run out of ideas. Can anyone help me out??? I will be greatly indebted. First: how did you check that the library is

Re: [R] Unable to load R

2009-05-28 Thread Zeljko Vrba
On Thu, May 28, 2009 at 02:02:52PM +0530, anupam sinha wrote: Actually I did a search using locate but could not find the file. But when Locate reports useful results only if its database is up-to-date. Try running one of ldd /usr/lib/R/bin/exec/R ldd /usr/lib64/R/bin/exec/R and see what it

[R] How to exclude a column by name?

2009-05-27 Thread Zeljko Vrba
Given an arbitrary data frame, it is easy to exclude a column given its index: df[,-2]. How to do the same thing given the column name? A naive attempt df[,-name] did not work :) __ R-help@r-project.org mailing list

Re: [R] How to exclude a column by name?

2009-05-27 Thread Zeljko Vrba
On Wed, May 27, 2009 at 12:52:41PM +0200, Paul Hiemstra wrote: This piece of code does the trick. Most important is the which() command: df = data.frame(a = runif(10), b = runif(10)) df[,-which(names(df) == a)] Thanks to you and Linlin. It did not occur to me to use which(); I thought

[R] Problem with adding labels in ggplot2

2009-05-27 Thread Zeljko Vrba
I apologize for not pasting a complete example, but the data-set is too large, so I hope someone can help me just by description of symptoms. I define a generic plot object name (note the missing y=.. in aes()) to plot different y-values against the same set of x-values. p.b4.generic.wg -

Re: [R] Problem with adding labels in ggplot2

2009-05-27 Thread Zeljko Vrba
On Wed, May 27, 2009 at 01:40:52PM +0200, Zeljko Vrba wrote: I apologize for not pasting a complete example, but the data-set is too large, so I hope someone can help me just by description of symptoms. -snip- I have solved the problem by introducing an artificial variable in the original

[R] Selection and aggregation in one operation?

2009-05-26 Thread Zeljko Vrba
I have a large data-frame with measurements such as: id i v1 v2 v3 1 1 1.1 1.2 1.3 1 2 1.4 1.5 1.6 1 3 1.5 1.7 1.8 2 1 2.1 2.2 2.3 2 2 2.7 2.5 2.6 2 3 2.4 2.8 2.9 For each unique value of 'id' (which in the real data-set is a combination of three variables) I want to compute the median

[R] postscript problems (landscape orientation)

2009-05-21 Thread Zeljko Vrba
I use the following function to export some figures to .eps: p.eps - function(p, fname, title = NULL, width, height) { postscript(file=fname, onefile=FALSE, paper=special, width=width, height=height, horizontal=FALSE) print(p + opts(title = title)) dev.off() } Whenever I have

Re: [R] Inverse

2009-05-21 Thread Zeljko Vrba
On Thu, May 21, 2009 at 08:36:27PM +1000, Kon Knafelman wrote: I have coded the following polynomial function y= function(x) x^3-2*x^2+1 I need to find the inverse of this, but the code i am using now isnt returning what i want it to. This function is not injective, so the inverse

Re: [R] postscript problems (landscape orientation)

2009-05-21 Thread Zeljko Vrba
On Thu, May 21, 2009 at 12:32:28PM +0100, Prof Brian Ripley wrote: Using Acrobat Reader to view PostScript! It is a PDF viewer. Ah, sorry, I explicitly convert the PS with ghostscript's ps2pdf. suspect you need to track down where conversion to PDF is happening and disable

Re: [R] postscript problems (landscape orientation)

2009-05-21 Thread Zeljko Vrba
On Thu, May 21, 2009 at 02:14:01PM +0200, Peter Dalgaard wrote: I think the trick is jade:~/ env | grep GS_ GS_OPTIONS=-dAutoRotatePages=/None Thanks, I found that myself. However, when using ps2pdf from Miktex 2.7, I get the following error: Unrecoverable error: typecheck in

Re: [R] What's the best way to tell a function about relevant fields in data frames

2009-05-12 Thread Zeljko Vrba
On Tue, May 12, 2009 at 12:18:59PM +0200, Titus von der Malsburg wrote: Is there a canonical way to tell a function which fields in a data frame are relevant? What other alternatives are possible? What are the pros and cons of the alternatives? Why not simply rearrange your data frames to

Re: [R] From two colors to 01 sequences

2009-05-12 Thread Zeljko Vrba
On Tue, May 12, 2009 at 12:20:56PM +0100, Paul Smith wrote: I have got several pdf files with rows of colored rectangles: red rectangles should be read as 0; green rectangles as 1. No other color exists. Is there some way to have R reading the colored rectangles to a matrix or data frame

[R] ggplot2: recommended workaround for broken legend.position=top

2009-05-10 Thread Zeljko Vrba
Searching the mail archives I found that using legend.position as in p.ring.3 + opts(legend.position=top) is a known bug. I tried doing p.ring.3 + opts(legend.position=c(0.8, 0.2)) which works, but the legend background is transparent, i.e. I see the plot background through the legend. Adding

Re: [R] Generating plots with two lines, each line plotted with a different color

2009-05-07 Thread Zeljko Vrba
On Thu, May 07, 2009 at 08:00:59AM -0400, John Sorkin wrote: I would like to have the two lines generated by my plot have different colors, in addition to the two different line types specified by the plot command below. I would appreciate advice on how I can specify line colors.

Re: [R] Do you use R for data manipulation?

2009-05-06 Thread Zeljko Vrba
Sorry for reply to the wrong person, I lost the original email. Farrel Buchinsky wrote: Is R an appropriate tool for data manipulation and data reshaping and data organizing? I think so but someone who recently joined our group thinks not. The new recruit believes that python or another

Re: [R] problem with ggplot2 boxplot, groups and facets

2009-05-05 Thread Zeljko Vrba
On Tue, May 05, 2009 at 09:21:49PM -0400, stephen sefick wrote: If you would provide a reproducible example I could tinker. The code I posted is real. I attach here the relevant part of the data-set as CSV. Best regards, Zeljko. ,sp,wg,n,v.realtime 9968,GP,1,1,28.924 55802,GP,1,2,23.5566

[R] Calculating relative ratios in a data frame..

2009-05-04 Thread Zeljko Vrba
I have a data-set that is structured as follows: spwgn v.realtime v.cputime v.tcputime v.idletime v.nswtch 9 0 11 28.61300 28.6128.6039 0.0e+00 407 1563 0 12 15.20270 30.3828.5981 9.80523e-01 483 3128 0 1

Re: [R] Calculating relative ratios in a data frame..

2009-05-04 Thread Zeljko Vrba
On Mon, May 04, 2009 at 04:44:37PM -0400, jim holtman wrote: Here is a way -- you are dividing by zero in the scaling: Thanks to both of you who answered. I know about division by zero, but I basically ignore rows having n==1 when plotting. __

Re: [R] running R on netbooks/minis?

2009-05-03 Thread Zeljko Vrba
On Sun, May 03, 2009 at 03:27:43AM -0500, Erin Hodgess wrote: Is it possible to run R on a netbook/mini, please? There should be no reason not to be possible, if the notebook uses an OS that R supports. __ R-help@r-project.org mailing list

[R] Reccomendation for graphics package

2009-05-01 Thread Zeljko Vrba
Hello, What would you recommend for producing publication-quality plots with R? Built-in graphics, trellis, ggplot2, or something else? Basic requirements: - I need to draw line-, box-, density-plots, bar-charts and histograms - error bars on bar- and box-plots - easy tiling of multiple plots

Re: [R] Reccomendation for graphics package

2009-05-01 Thread Zeljko Vrba
On Fri, May 01, 2009 at 01:06:34PM -0500, hadley wickham wrote: It should be trivial with ggplot2 too, but it's hard to provide concrete advice without a concrete problem. Elementary problem: qplot(wg, v.realtime, data=df.best.medians$gv1, facets = . ~ n, colour=sp) produces a nice plot.

Re: [R] Reccomendation for graphics package

2009-05-01 Thread Zeljko Vrba
On Fri, May 01, 2009 at 09:38:13PM +0200, Zeljko Vrba wrote: On Fri, May 01, 2009 at 01:06:34PM -0500, hadley wickham wrote: It should be trivial with ggplot2 too, but it's hard to provide concrete advice without a concrete problem. Elementary problem: qplot(wg, v.realtime, data

Re: [R] Reccomendation for graphics package

2009-05-01 Thread Zeljko Vrba
On Fri, May 01, 2009 at 03:23:18PM -0500, hadley wickham wrote: If you have a categorical x axis, you need to specify the group aesthetic which defines what group of points should form a line. It's hard to tell what that should be from your example, maybe sp? Yes, adding group=sp works.

Re: [R] Batch importing data with respective naming

2009-04-28 Thread Zeljko Vrba
On Tue, Apr 28, 2009 at 07:09:04PM -1000, Taylor Hermes wrote: I tried the following: Add this before for(): peak - list() for (i in 1:100) { peak[[i]] - read.table(paste(i,--one--hist.txt, sep=), sep=,, header=TRUE) } __