Re: [R] symbols duplicated in plot output

2009-06-29 Thread Dieter Menne
baptiste auguie-5 wrote: > > > pdf() > plot(1, xlab=expression(mu)) > dev.off() > > If I open this pdf in Illustrator CS4, there are two "mu" on top of each > other, giving it a somewhat bold aspect. Other characters not from the > symbol font are just output in one version. I'm guessing this

Re: [R] list subsetting

2009-06-29 Thread Tal Galili
how about ? for(i in 1:3) { b[[i]][[3]] <- 3 } or: for(i in 1:3) { b[[i]][[3]] <- a[[i]] } On Tue, Jun 30, 2009 at 6:54 AM, Ivo Shterev wrote: > > Dear All, > > I am interested to know how assign a list to a list of lists. For example > how to assign the list a = list(3,3,3) to the list

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

2009-06-29 Thread Raymond Wan
Hi Zeljko, Zeljko Vrba wrote: Windows 32bit results: user system elapsed 21.660.02 21.69 Linux 64bit Results user system elapsed 27.242 0.004 27.275 Using wall-clock time metric is not "two different ways" of timing. He could have just as well measured the time using stop-wat

[R] R version-2.9.1 for Linux

2009-06-29 Thread utkarshsinghal
Hi All, I am currently using R version 2.8.1 on linux cent os 4.4 (i386) and want to upgrade to version 2.9.1. It seems to me that version-2.9.1 is it not for my OS. Am I right? Regards Utkarsh [[alternative HTML version deleted]] __ R-hel

[R] lmer (or lme) with heteroscedasticity

2009-06-29 Thread Daryl Morris
Hello, I'm trying to fit a mixed-effects model with a single binary predictor (case/control status in my case), a random intercept (e.g. dependent on radiologist) and also a random slope (a per-radiologist difference between cases and controls). I know how to do that, but what I don't know h

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

2009-06-29 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 > >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 > > 21.660.02

[R] list subsetting

2009-06-29 Thread Ivo Shterev
Dear All, I am interested to know how assign a list to a list of lists. For example how to assign the list a = list(3,3,3) to the list b: > b [[1]] [[1]][[1]] [1] 0 [[1]][[2]] [1] 0 [[1]][[3]] [1] 0 [[2]] [[2]][[1]] [1] 0 [[2]][[2]] [1] 0 [[2]][[3]] [1] 0 [[3]] [[3]][[1]] [1] 0 [[3]][[

Re: [R] Automatically placing a legend in an area with the most white space...

2009-06-29 Thread Frank E Harrell Jr
Also check out the Hmisc package's largest.empty and putKeyEmpty and related functions. Frank Jason Rupert wrote: That fixed it. I thought I was getting the hang of R, but I guess I still have a long way to go. One last question regarding the dynamic placement of the legend in the location

Re: [R] Automatically placing a legend in an area with the most white space...

2009-06-29 Thread Jason Rupert
That fixed it. I thought I was getting the hang of R, but I guess I still have a long way to go. One last question regarding the dynamic placement of the legend in the location on a plot with the largest empty space, is there a way to assure that the legend doesn't go off the plot or is par

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

2009-06-29 Thread Raymond Wan
Hi, I. Soumpasis wrote: 2009/6/29 C騷ar Freitas 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 21.660.02 21.69 Linux 64bit Results user system elapsed 27.242 0.004 27.275

Re: [R] Add ID numbers on a plot

2009-06-29 Thread Jorge Ivan Velez
Hi Steve, Take a look at ?matplot; it might be useful in this particular case. HTH, Jorge On Mon, Jun 29, 2009 at 2:16 PM, Steve Hong wrote: > Dear List, > > I have (for example) 50 observations collected from 50 experimental sites > and want to look at changes of 50 observations as function

Re: [R] Add ID numbers on a plot

2009-06-29 Thread Uwe Ligges
See ?text Uwe Ligges Steve Hong wrote: Dear List, I have (for example) 50 observations collected from 50 experimental sites and want to look at changes of 50 observations as function of time in a graph. I found that I could do that using R-code below: time2 <- 1:25 y1=rnorm(25, mean=0, sd=

Re: [R] Large Stata file Import in R

2009-06-29 Thread Carlos J. Gil Bellosta
Hello, You are dealing with two different problems at the same time: importing Stata data and importing a relatively big file. Can you try to export your data to txt file first and try to import from it directly? Secondly, problems concerning reading big files with R occur quite often and there

Re: [R] Linux Installation Version

2009-06-29 Thread Mark Knecht
On Mon, Jun 29, 2009 at 12:54 PM, Patrick Gedeon wrote: > Hello, > > I would like to download the correct R program RMP for my linux > distribution, but I am not sure which one I need. > > The results of a "cat /proc/version" is: > > Linux version 2.6.23-9.ydl6.1 (r...@build2.terraplex.com) (gcc ve

Re: [R] How to select partially (not completely) unique rows?

2009-06-29 Thread Henrique Dallazuanna
Try this: DF[!duplicated(DF[,c("lat","lon")]),] On Mon, Jun 29, 2009 at 6:55 PM, Mark Na wrote: > Dear R-helpers, > > I know how to use unique to select unique rows, e.g. > > unique.rows<-unique(dataframe) > > but I would like to select those rows that are unique only only TWO of my > datafram

[R] How to select partially (not completely) unique rows?

2009-06-29 Thread Mark Na
Dear R-helpers, I know how to use unique to select unique rows, e.g. unique.rows<-unique(dataframe) but I would like to select those rows that are unique only only TWO of my dataframe's columns (so, two rows with the same value on these two columns would not be kept, even if they had different v

[R] Linux Installation Version

2009-06-29 Thread Patrick Gedeon
Hello, I would like to download the correct R program RMP for my linux distribution, but I am not sure which one I need. The results of a "cat /proc/version" is: Linux version 2.6.23-9.ydl6.1 (r...@build2.terraplex.com) (gcc version 4.1.1 20070105 (Red Hat 4.1.1-52.ydl.1)) #1 SMP Sat Jan 26 20:1

Re: [R] infer haplotypes phasing trios tdthap

2009-06-29 Thread Jing Hua Zhao
Dear Tiago and David, Just for a bit of personal interest, Haplin has been part of CRAN which is potentially helpful. http://cran.r-project.org/web/packages/Haplin/index.html from http://www.uib.no/smis/gjessing/genetics/software/haplin/ Best regards, Jing Hua -Original Message- Fro

[R] Memory issues with 'aov' function

2009-06-29 Thread loxx0076
I am trying to run an ANOVA with a within/between subjects design. It is a 2 (within) x 3 x 3 x 5 design with 100 observations per cell. When I am trying to run this, I get an error message saying that there is insufficient memory, even after allocating max memory for my machine (4 gig). We

[R] Rattle: Trouble opening .csv files

2009-06-29 Thread Thompson, James
I am a novice R user and recently installed Rattle primarily for the decision tree functionality. I immediately run into problems trying to open a .csv file from the "Data" tab. The error message I get in my "R" log window is: Error in if (crv$load.tooltips) loadTooltips() : argume

Re: [R] How to use "subset" in lm function

2009-06-29 Thread Gabor Grothendieck
The dyn package (and also the dynlm package) can be used for this. library(dyn) y <- ts(1:10) dyn$lm(y ~ lag(y, -1)) On Mon, Jun 29, 2009 at 2:50 PM, Hongwei Dong wrote: > Hi, I'm using R to do a time series analysis. In the model, I use the lags > of some variables. such the lags of the variable

Re: [R] ANOVA

2009-06-29 Thread Tal Galili
Hi Sarah I can't see the file. The general code can be seen in the example on ?aov (assuming the design is balanced) On Mon, Jun 29, 2009 at 3:07 PM, Georgina Sarah Humphreys < g.humphrey...@research.gla.ac.uk> wrote: > I have the attached data set (csv) and I want to run an analysis of > v

[R] Help with bootstrapping regression

2009-06-29 Thread Adriana Bejarano
Dear R gurus, I have the following x (conc)-y (rk) data. rk<-c(0.016,0.032,0.048,0.095,0.111,0.143,0.190,0.206,0.222,0.270,0.286,0.302,0.317,0.381,0.397,0.444,0.460, 0.476,0.492,0.508,0.524,0.540,0.556,0.651,0.698,0.714,0.810,0.825,0.841,0.921,0.937,0.952,0.968,0.984,1.000) cc<-c(0.4,0.53,1,1.5

Re: [R] How to use "subset" in lm function

2009-06-29 Thread milton ruser
Like this? x<-1:100 y<-runif(100) df<-data.frame(cbind(x,y)) mycoef.list<-NULL for (i in 10:65) { mycoef<-coef(lm(y~x, data=subset(df, x %in% 1:i)))[2] mycoef.list <-c(mycoef.list, mycoef) } hist(mycoef.list) regards milton brazil=toronto On Mon, Jun 29, 2009 at 2:50 PM, Hongwei Dong wrot

Re: [R] (no subject)

2009-06-29 Thread John Kane
Perhaps? http://finzi.psych.upenn.edu/R/library/quantreg/html/akj.html --- On Mon, 6/29/09, maram salem wrote: > From: maram salem > Subject: [R] (no subject) > To: r-help@r-project.org > Received: Monday, June 29, 2009, 9:05 AM > Hi group, > I found a module for adaptive kernel density est

[R] A graphic art piece I created that only statisticians would understand and the use of R in art

2009-06-29 Thread Thomas Levine
This isn't the typical r-help post, but I think people here will appreciate it. The current Shirt.Woot (http://shirt.woot.com) derby prompt is fake band names, so I chose the name "The Probable Error of a Mean." And I made my submission with R. Here it is: http://www.thomaslevine.org/the_probable_

[R] How to use "subset" in lm function

2009-06-29 Thread Hongwei Dong
Hi, I'm using R to do a time series analysis. In the model, I use the lags of some variables. such the lags of the variables have different length, I just can't use them directly in the lm function. Intuitively, I feel that "subset" might be useful, but I do not know how to use it. Anyone can give

[R] symbols duplicated in plot output

2009-06-29 Thread baptiste auguie
Dear list, I'm manually editing a large collection of pdf files that I produced with R in Adobe Illustrator. In doing so, I've had the surprise to see symbols duplicated on top each other. The following code illustrates this, pdf() plot(1, xlab=expression(mu)) dev.off() If I open this pdf in Ill

[R] Add ID numbers on a plot

2009-06-29 Thread Steve Hong
Dear List, I have (for example) 50 observations collected from 50 experimental sites and want to look at changes of 50 observations as function of time in a graph. I found that I could do that using R-code below: time2 <- 1:25 y1=rnorm(25, mean=0, sd=1) y2=rnorm(25, mean=0, sd=1) ... y50=rnorm

Re: [R] multiple regression w/ no intercept; strange results

2009-06-29 Thread Thomas Lumley
On Mon, 29 Jun 2009, John Hunter wrote: But my question was more numerical: in particular, the R^2 of the model should be equal to the square of the correlation between the fit values and the actual values. No. It is with the intercept and is not w/o it, as my code example shows. Am I corr

Re: [R] a plot of stacked boxes

2009-06-29 Thread Osman Radi
Dear Jim, This almost exactly what I need. However, because I am not trying to plot price change as the color parameter, I need the scale to be positive only (0 to 1) and would like to be able to control the color scheme of the scale. I looked at the help file for map.market() and it's cod

Re: [R] questions about meta-analysis

2009-06-29 Thread Thomas Lumley
On Sat, 27 Jun 2009, Emmanuel Charpentier wrote: Le samedi 27 juin 2009 à 13:02 +0800, sdzhangping a écrit : Dear R users: In the example of meta-analysis (cochrane, package rmeta), I can not found the p-value of Test for overall effect, and some other indices (Z, I, weight and et al). How ca

Re: [R] ggplot2 x axis question

2009-06-29 Thread Christopher Desjardins
Thanks Hadley that worked. Chris On 6/29/09 11:05 AM, hadley wickham wrote: > In that case, try: > > qplot(reorder(factor(model),delta),delta,data=growthm.bic) > > Deepayan: do you think there should also be a numeric method for reorder? > > Hadley > > On Mon, Jun 29, 2009 at 10:39 AM, Christopher

[R] nlsList {nlme} - control arguments problem

2009-06-29 Thread Rick DeShon
Hi All. I'd like to send some control arguments to the nls function when performing a nlsList analysis. I'm fitting a power model to some grouped data and would like to impose lower bounds on the estimates using the "port" algorithm. Obtaining the lower bound constraint works fine with a direct c

Re: [R] Hosmer Lemshaw

2009-06-29 Thread Girish A.R.
Hi, This doesn't quite answer your question, but the following discussion (posted some time back in this forum) would lead you in that direction: http://www.nabble.com/Re%3A-More-Logistic-Regression-Tools--td3793516.html#a3793516 best, -Girish alon.benari wrote: > > Hello R experts, > > Is

Re: [R] ggplot2 x axis question

2009-06-29 Thread hadley wickham
In that case, try: qplot(reorder(factor(model),delta),delta,data=growthm.bic) Deepayan: do you think there should also be a numeric method for reorder? Hadley On Mon, Jun 29, 2009 at 10:39 AM, Christopher Desjardins wrote: > Hi Hadley, > Thanks for the reply and the great graphing package. That

Re: [R] ggplot2 x axis question

2009-06-29 Thread Christopher Desjardins
Hi Hadley, Thanks for the reply and the great graphing package. That code is giving me the following error: > qplot(reorder(model,delta),delta,data=growthm.bic) Error in UseMethod("reorder") : no applicable method for "reorder" Cheers, Chris On 6/28/09 8:21 PM, hadley wickham wrote: > Hi Chris

[R] Hosmer Lemshaw

2009-06-29 Thread Alon Ben-Ari
Hello R experts, Is the a function to preform a Hosmer Lemshaw goodness of fit? I am running R 2.8. Thanks, AlonBen-Ari Department of Anesthesiology and CCM Univ. Pittsburgh [[alternative HTML version deleted]] __ R-help@r-project.org mai

[R] ANOVA

2009-06-29 Thread Georgina Sarah Humphreys
I have the attached data set (csv) and I want to run an analysis of variance on the wingsize data (comparing infected vs non-infected) within and between experiments. Can anyone help me with the command I should use? Many thanks Georgina PhD Student Division of Infection and Immunity B5-29,

[R] Help

2009-06-29 Thread maram salem
Hi group, I found a module for adaptive kernel density estimation for Stata users, but unfortunetly I don't have access to Stata, can I find a similar approach using R? Thank u so much 4 ur time. [[alternative HTML version deleted]] _

Re: [R] Meta-Analysis: Metaplot Labels, SE, and Summary

2009-06-29 Thread Michael Dewey
At 02:37 29/06/2009, Ishwar Bridgelal wrote: Hello, I'm ashamed to admit that I'm an R and meta-analysis neophyte just when my colleagues and I are attempting to run R for that latter purpose. We would like to construct a graph of the 95% confidence intervals of odds ratios (inverse odds ratios,

[R] Stata file Import and Analysis in R

2009-06-29 Thread saurav pathak
Hi I have a stata data set (.dta file) of size 600 MB, I need to import it in R and do a 2SLS multilevel analysis on the data set. I would eb grateful if help is provided for the first part of how to import this big file from Stata to R and then how to open the imported Stata file in R? Kindly hel

[R] Large Stata file Import in R

2009-06-29 Thread saurav pathak
Hi I am using Stata 10 and I need to import a data set in stata 10 to R, I have saved the dataset in lower versions of Stata as well by using saveold command in Stata. My RAM is 4gb and the stata file is 600MB, I am getting an error message which says : "Error: cannot allocate vector of size 3.4

Re: [R] Compute correlation matrix for panel data with specific ordering

2009-06-29 Thread John Kane
Have a look at str(df). Those values are being interpreted as factors not numbers. I don't think this is what you want. --- On Mon, 6/29/09, Serguei Kaniovski wrote: > From: Serguei Kaniovski > Subject: Re: [R] Compute correlation matrix for panel data with specific > ordering > To: r-help

Re: [R] how to call Perl from R

2009-06-29 Thread Gabor Grothendieck
The read.xls function in the gdata package is another example (the call to system is actually in xls2csv which it calls). On Mon, Jun 29, 2009 at 10:02 AM, Marc Schwartz wrote: > On Jun 29, 2009, at 8:58 AM, mau...@alice.it wrote: > >> I read some archived posts about calling R from Perl scripts.

Re: [R] how to call Perl from R

2009-06-29 Thread Marc Schwartz
On Jun 29, 2009, at 8:58 AM, mau...@alice.it wrote: I read some archived posts about calling R from Perl scripts. There seems to be an R package creating the necessaary interface. I'd like to do the opposite. That is to call Perl from an R script. I wonder whether this is possible at all ???

[R] how to call Perl from R

2009-06-29 Thread mauede
I read some archived posts about calling R from Perl scripts. There seems to be an R package creating the necessaary interface. I'd like to do the opposite. That is to call Perl from an R script. I wonder whether this is possible at all ??? What about Bioperl which is a Perl variation built to d

Re: [R] Fold function with several time varying covariates

2009-06-29 Thread Allan Engelhardt
Looks like a problem with the cov.names argument, not cov. Can you make a stand-alone example showing the error? Allan. On 29/06/09 14:31, Jenkins, Helen wrote: > The error message when I use list is: > " s.matrix(!is.na(new.data[, cov.names])) : >invalid subscript type 'list'" > [...] >

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

2009-06-29 Thread I. Soumpasis
2009/6/29 Cézar Freitas > Hi, Ilias. I think that is not ok to compare performance in different > plataforms of different machines. To compare times, is necessary you execute > the code at the two plataforms (Linux and Windows) in the same machine. > But the problem here is other. The advanced us

Re: [R] multiple regression w/ no intercept; strange results

2009-06-29 Thread John Hunter
On Sun, Jun 28, 2009 at 3:38 AM, Dieter Menne wrote: > It seems odd to me that dropping the intercept > would cause the R^2 and F stats to rise so dramatically, and the p > value to consequently drop so much. In my implementation, I get the > same beta1 and beta2, and the R2 I compute using the

Re: [R] Fold function with several time varying covariates

2009-06-29 Thread Jenkins, Helen
The error message when I use list is: " s.matrix(!is.na(new.data[, cov.names])) : invalid subscript type 'list'" **PLEASE NOTE MY EMAIL ADDRESS HAS CHANGED TO h.jenkin...@imperial.ac.uk** Helen Jenkins Dept of Infectious Disease Epidemiology Imperial College, London St Mary's Campus Norfolk Pla

[R] Boxplot with dates

2009-06-29 Thread Mafalda Viana
Dear R users, I am having problems making a boxplot across time. I wish to plot the "Values" (y-axis) against "Date" (x-axis) but would like the Date to be sequential from 2007-11-15 until 2008-01-31 (therefore there will be lots of NA's e.g. from 2007-11-16 to 2007-11-26 nothing will show in the

Re: [R] Exporting large datasets

2009-06-29 Thread Allan Engelhardt
Depends on what you mean by "export". I like save(). Check out the options in the help. It seems to be very efficient and I had no problems working with the 50,000 x 15,000 data sets from the 2009 KDD Cup challenge. Maybe your disk is just slow (network file storage)? Allan. On 29/06/09 07

Re: [R] Stata file Import and Analysis in R

2009-06-29 Thread Liviu Andronic
Hello, On Mon, Jun 29, 2009 at 2:58 PM, saurav pathak wrote: > I have a  stata data set (.dta file) of size 600 MB, I need to import it in > If you're new to R, you can use library(Rcmdr) menus to import the stata data set. Be sure to have at least as much RAM as the size of the data set. Liviu _

Re: [R] Fold function with several time varying covariates

2009-06-29 Thread Allan Engelhardt
On 28/06/09 18:02, Jenkins, Helen wrote: But this just gives me an error message. What error message? __ 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/post

Re: [R] Stata file Import and Analysis in R

2009-06-29 Thread Fredrik Karlsson
Hi, You should install the foreign package and try the read.dta function. I don't know whether the size of your fiel will be a problem, but you would try it at least. /Fredrik On Mon, Jun 29, 2009 at 1:58 PM, saurav pathak wrote: > Hi > I have a  stata data set (.dta file) of size 600 MB, I nee

[R] Select values at random by id value

2009-06-29 Thread James Martin
All, I have data that looks like below. For each id there may be more than one value per day. I want to select a random value for that day for that id. The end result would hopefully be a matrix with the id as rows, date as columns and populated by the random hab value. Thanks to someone on her

[R] (no subject)

2009-06-29 Thread maram salem
Hi group, I found a module for adaptive kernel density estimation for Stata users, but unfortunetly I don't have access to Stata, can I find a similar approach using R? [[alternative HTML version deleted]] __ R-help@r-project.org maili

[R] Stata file Import and Analysis in R

2009-06-29 Thread saurav pathak
Hi I have a stata data set (.dta file) of size 600 MB, I need to import it in R and do a 2SLS multilevel analysis on the data set. I would eb grateful if help is provided for the first part of how to import this big file from Stata to R and then how to open the imported Stata file in R? Kindly hel

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

2009-06-29 Thread Cézar Freitas
--- Em seg, 29/6/09, I. Soumpasis escreveu: De: I. Soumpasis Assunto: Re: [R] (performance) time in Windows vs Linux Para: "Cézar Freitas" Data: Segunda-feira, 29 de Junho de 2009, 5:53 Hi Cezar, I tried your code in a core duo laptop (@2.5Gz) with ubuntu x86_64 with 4GB of RAM. Both R a

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

2009-06-29 Thread Cézar Freitas
Hi, Ilias. I think that is not ok to compare performance in different plataforms of different machines. To compare times, is necessary you execute the code at the two plataforms (Linux and Windows) in the same machine. But the problem here is other. The advanced users explained it well.   So I a

[R] How can I get estimates including standard errors for parameters in constraint for retricted models

2009-06-29 Thread bonnychen
Hello, I am running linear regressions on a restricted model, in which one parameter, say gamma1, equals alpha/(1-alpha), and another parameter, say gamma2, equals beta/(1-alpha), i.e.: y = c0 + gamma1*x1 + gamma2*x2, where gamma1= alpha/(1-alpha), gamma2=beta/(1-alpha). I can get estimates

Re: [R] Exporting large datasets

2009-06-29 Thread jim holtman
What do you want to export it to and in what format? You will probably have a problem trying to write that large an object out at one time. Try instead opening a connecting and the writing much smaller pieces to the file: x <- matrix(runif(1), 100) x[,1] <- 1:100 # put some labels for testi

[R] Odp: re moving value from a vector

2009-06-29 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 29.06.2009 10:58:50: > > Hallo, > > I have a vector of several iterations and I have to remove /filter out all > the values < than 10, > how can I do this? vector[vector>=10] does what you want but word "vector of iterations" frightens me a little

Re: [R] Lattice and high-resolution tiffs

2009-06-29 Thread Dieter Menne
r-help.20.trevva wrote: > > I am trying to produce high-resolution (600dpi+) TIFF figures for use > in a publication. Everything seems to work well when using the > "normal" R-graphics and the relative text size obtained is independent > of the output resolution. However, when I try and make la

Re: [R] re moving value from a vector

2009-06-29 Thread Henrique Dallazuanna
Try this: x[!x < 10] On Mon, Jun 29, 2009 at 5:58 AM, aledanda wrote: > > Hallo, > > I have a vector of several iterations and I have to remove /filter out all > the values < than 10, > how can I do this? > > Thanks a lot > > Alessandra > -- > View this message in context: > http://www.nabble.c

Re: [R] Compute correlation matrix for panel data with specific ordering

2009-06-29 Thread Serguei Kaniovski
I apologize for not being specific enough in my previous posting. Assume you have panel data in the form: df <- data.frame( cbind( rep( c( "AUT" , "BEL" , "DEN" , "GER" ) , 4) , cbind( rep( c( 1999 , 2000 , 2001 , 2002 ) , 4 ) ), sample( 10 , 16 , replace=T) ) ) names(df) <- c( "country" , "y

Re: [R] Saving RGL Graphs

2009-06-29 Thread Duncan Murdoch
Patrick Gedeon wrote: Dear Users: I have created several 3d scatter plots using the rgl package. I wish to save each of them so that I can later open them in R and view/rotate them. What is the command to do this? There is no way to do that. It would be a nice feature, but it requires some

Re: [R] Export rgl plot3d scatter

2009-06-29 Thread Yihui Xie
Will it be so difficult to switch to R during your demonstration? I think it's much more difficult to embed the rgl demo in your document and keep the rotating feature. You may take a look at the example for 'rgl.snapshot', which shows you how to create an animation with ImageMagick, but I still pr

[R] Summary of why R has the 2^31-1 limit?

2009-06-29 Thread Allan Engelhardt
I know it has been discussed before, but is there anywhere a good summary of (1) why R has the 2^31-1 vector length limit on all platforms (specifically 64-bit, of course) and (b) what would be the effort/implications of changing it? I think I have seen it but I couldn't find it and it does no

Re: [R] Saving RGL Graphs

2009-06-29 Thread Liviu Andronic
Hello, On Mon, Jun 29, 2009 at 6:49 AM, Patrick Gedeon wrote: > I have created several 3d scatter plots using the rgl package. I wish to > save each of them so that I can later open them in R and view/rotate them. > What is the command to do this? > Wouldn't it suffice to keep in a text file the c

[R] Lattice and high-resolution tiffs

2009-06-29 Thread r-help . 20 . trevva
Hi, I am trying to produce high-resolution (600dpi+) TIFF figures for use in a publication. Everything seems to work well when using the "normal" R-graphics and the relative text size obtained is independent of the output resolution. However, when I try and make lattice plots, the relative size of

Re: [R] comma-separated thousands in numbers on plot axes?

2009-06-29 Thread Jim Lemon
Alexy Khrabrov wrote: How can I make R separate thousands, millions, etc., on the plot axes, with commas? Hi Alexy, You might try formatting the labels with the "comma" function in the ggplot2 package: axis(1,at=seq(100,500,by=100), labels=comma(seq(100,500,by=100)

[R] re moving value from a vector

2009-06-29 Thread aledanda
Hallo, I have a vector of several iterations and I have to remove /filter out all the values < than 10, how can I do this? Thanks a lot Alessandra -- View this message in context: http://www.nabble.com/removing-value-from-a-vector-tp24250608p24250608.html Sent from the R help mailing list ar

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

2009-06-29 Thread Raymond Wan
Prof Brian Ripley wrote: I meant to write "not so for 'top'" in the final para. Ah, I'm not certain enough to know that "htop" works for threads as well...so I was quick to jump to agreeing with you. :-) I only know it works for multi-cores... Ray ___

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

2009-06-29 Thread Prof Brian Ripley
I meant to write "not so for 'top'" in the final para. On Mon, 29 Jun 2009, Prof Brian Ripley wrote: On Mon, 29 Jun 2009, Raymond Wan wrote: milton ruser wrote: 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?

[R] Exporting large datasets

2009-06-29 Thread DrDimo
Hi all, I am currently working on a particularly large (data wise) project, the problem I am having is with exporting the data (for use later on). My current sample dataset contains about 400million entries (approximately 2 columns x 2 rows). I tried using the write.table command, but 8

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

2009-06-29 Thread Gavin Simpson
On Mon, 2009-06-29 at 19:21 +1200, Patrick Connolly wrote: > On Mon, 29-Jun-2009 at 09:05AM +0200, Zeljko Vrba wrote: > > |> 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? > |> > |> > |>

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

2009-06-29 Thread Raymond Wan
Hi Brian, Thank you for the clarification -- the first part does set the record straight about what I thought about R. I would expect a program to run on a single core by default unless something specifically (and somewhat non-trivial) was done to it. Prof Brian Ripley wrote: If you want

Re: [R] Compute correlation matrix for panel data with, specific ordering

2009-06-29 Thread Dieter Menne
Serguei Kaniovski wifo.ac.at> writes: > Ok I see how to sort the factors, but how do I compute the correlation > matrix in a repeated observations dataset (see the first part of my > question) "Repeated" observations? You must be more specific then. Which function are you using (lme? Error()

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

2009-06-29 Thread Prof Brian Ripley
On Mon, 29 Jun 2009, Raymond Wan wrote: milton ruser wrote: 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 don't know the answer in the context of R -- I didn't know that R can use multiple cores by default

Re: [R] R CMD build for windows package from linux

2009-06-29 Thread Prof Brian Ripley
On Sun, 28 Jun 2009, Jorge Cornejo wrote: Hi, exist any way to create a windows installable package from a Linux R installation ?? Not in general, but see 'Writing R Extensions'. If no compiled code is involved -- you can then zip up the installed Linux package. More generally, you can subm

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

2009-06-29 Thread Patrick Connolly
On Mon, 29-Jun-2009 at 09:05AM +0200, Zeljko Vrba wrote: |> 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 |> >

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

2009-06-29 Thread Raymond Wan
Hi, milton ruser wrote: 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 don't know the answer in the context of R -- I didn't know that R can use multiple cores by default? But in general, I use "htop", who

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 G

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

2009-06-29 Thread Patrick Connolly
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 (install with "aptitude install gkrellm" if you don't have it already). It