[R] package for factor analysis

2005-12-14 Thread Birgit Kessler
I have Windows XP Professional Version 2002 and the R-Version 2.1.1. I want to do factor analysis with R. In Google Search I find the reference to the RScaLAPACK-Package, but my R-Version 2.1.1. told me that RScaLAPACK is invalid package, bevor 2.0.0. installed . Which package can I take for

[R] tool for cluster analysis

2005-12-14 Thread Birgit Kessler
I have Windows XP Professional Version 2002 and the R-Version 2.1.1. I did cluster analysis with the cluster package and the agnes (method = “ward”). The results are satisfactory. But the dendrogram of agnes is confused to work with the results. Is there a tool, I can get a clear arrangement of

Re: [R] package for factor analysis

2005-12-14 Thread Roger Bivand
On Wed, 14 Dec 2005, Birgit Kessler wrote: ?factanal (in the stats package, which is part of every base distribution of R) I have Windows XP Professional Version 2002 and the R-Version 2.1.1. I want to do factor analysis with R. In Google Search I find the reference to the

[R] package for factor analysis

2005-12-14 Thread Birgit Kessler
I have Windows XP Professional Version 2002 and the R-Version 2.1.1. I want to do factor analysis with R. In Google Search I find the reference to the RScaLAPACK-Package, but my R-Version 2.1.1. told me that “RScaLAPACK is invalid package, bevor 2.0.0. installed” . Which package can I take for

[R] glmmADMB: Generalized Linear Mixed Models using AD Model Builder

2005-12-14 Thread Hans Julius Skaug
Dear R-users, Half a year ago we put out the R package glmmADMB for fitting overdispersed count data. http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html Several people who used this package have requested additional features. We now have a new version ready. The major new feature is

Re: [R] tool for cluster analysis

2005-12-14 Thread TEMPL Matthias
I have Windows XP Professional Version 2002 and the R-Version 2.1.1. I did cluster analysis with the cluster package and the agnes (method = ward). The results are satisfactory. But the dendrogram of agnes is confused to work with the results. Is there a tool, I can get a clear

[R] [R-pkgs] Updated version of gdata, gtools, gplots and gmodels

2005-12-14 Thread Jain, Nitin
Hello, We have submitted the updated version of gdata, gmodels, gplots and gtools to CRAN. Summary of the changes is attached at the end. Best, Nitin __ Nitin Jain, PhD [EMAIL PROTECTED] Non Clinical Statistics Pfizer, Inc. (Groton, CT) Bldg: 260, # 1451 Ph: (860) 686-2526

[R] Problem with RSPerl

2005-12-14 Thread madhurima bhattacharjee
Hello, I am using RSPerl of omegahat.org to call R modules from Perl . I have tried running the test,pl script provided by RSPerl package. It gives the following error: perl: symbol lookup error: /usr/lib/perl5/site_perl/5.8.5/i386-linux-thread-multi/auto/R/R.so: undefined symbol: tryEval I

Re: [R] package for factor analysis

2005-12-14 Thread Jacques VESLOT
you can try : plot(hclust(dist(ach, manh), ward), hang=-1) you can choose distance method and clustering method. you can find a lot of useful functions for multivariate analysis in the ade4 package ; but, for clustering, the cluster package is quite nice. Birgit Kessler a écrit : I have

Re: [R] correct C function usage

2005-12-14 Thread Ido M. Tamir
On Tuesday 13 December 2005 22:35, you wrote: On Tue, 13 Dec 2005, Ido M. Tamir wrote: Hello, I am not sure if I am interfacing with C correctly and _safely_ or if there is a better way esp. with regards to terminating the returned array. You need to pass the length to the C routine and

Re: [R] About help on 'mahalanobis'

2005-12-14 Thread Martin Maechler
HenrikB == Henrik Bengtsson [EMAIL PROTECTED] on Wed, 14 Dec 2005 17:05:32 +1100 writes: HenrikB Hi, HenrikB help on 'mahalanobis' (in the stats package in Rv2.2.0) now says: HenrikB Description: HenrikB Returns the Mahalanobis distance of all rows in 'x' and the vector

Re: [R] Age of an object?

2005-12-14 Thread Martin Maechler
Trevor == Trevor Hastie [EMAIL PROTECTED] on Tue, 13 Dec 2005 12:51:34 -0800 writes: Trevor It would be nice to have a date stamp on an object. Trevor In S/Splus this was always available, because objects were files. [are you sure about always available? In any case,

Re: [R] Age of an object?

2005-12-14 Thread Philippe Grosjean
Martin Maechler wrote: Trevor == Trevor Hastie [EMAIL PROTECTED] on Tue, 13 Dec 2005 12:51:34 -0800 writes: Trevor It would be nice to have a date stamp on an object. Trevor In S/Splus this was always available, because objects were files. [are you sure about always

Re: [R] correct C function usage

2005-12-14 Thread Martin Maechler
Ido == Ido M Tamir [EMAIL PROTECTED] on Wed, 14 Dec 2005 10:37:07 +0100 writes: Ido On Tuesday 13 December 2005 22:35, you wrote: On Tue, 13 Dec 2005, Ido M. Tamir wrote: Hello, I am not sure if I am interfacing with C correctly and _safely_ or if there is a

[R] Design library

2005-12-14 Thread Mohammad Reza
Dear friends Hello I'm a PhD student in Edinburgh University and interested in survival analysis. By chance, I found design library in Splus software. I have some questions about it and it's highly appreciated if you can help me. In its help (Design library help) I found the topics(not the

Re: [R] Age of an object?

2005-12-14 Thread Robert Kinley
Trevor == Trevor Hastie [EMAIL PROTECTED] on Tue, 13 Dec 2005 12:51:34 -0800 writes: Trevor It would be nice to have a date stamp on an object. Indeed it would. If I had a better brain and memory, and didn't have to run so many projects in parallel I could probably cope OK.

Re: [R] help with writing function

2005-12-14 Thread Petr Pikal
Hi or use outer x-1:10 x-sample(x) mat-outer(x,x,-) # result with zeroes matrix(mat[!mat==0],9,9) # get rid of them will give you a matrix with columns with intended result. HTH Petr On 13 Dec 2005 at 11:53, Robert Burrows wrote: Date sent: Tue, 13 Dec 2005 11:53:00 -0500

[R] Fitting binomial lmer-model, high deviance and low logLik

2005-12-14 Thread Ivar Herfindal
Hello I have a problem when fitting a mixed generalised linear model with the lmer-function in the Matrix package, version 0.98-7. I have a respons variable (sfox) that is 1 or 0, whether a roe deer fawn is killed or not by red fox. This is expected to be related to e.g. the density of red

Re: [R] help with writing function

2005-12-14 Thread Oarabile Molaodi
Thanks for the help Petr it works, and all who contributed managed with the suggested function below as I needed to remove the zeros and the upper triangular portion of the matrix.: IED - function(risk) { n - length(risk) mrisk - matrix( rep(risk, n), ncol=n, byrow=TRUE ) diff

Re: [R] superimpose density line over hist

2005-12-14 Thread Martin Maechler
Romain == Romain Francois [EMAIL PROTECTED] on Tue, 13 Dec 2005 15:40:59 +0100 writes: Romain A few comments : Romain - your code should be reproductible, otherwise it is useless. (that Romain recommandation is on the posting guide) Romain - that

Re: [R] Age of an object?

2005-12-14 Thread Prof Brian Ripley
On Sat, 14 Jan 2006, Philippe Grosjean wrote: Martin Maechler wrote: Trevor == Trevor Hastie [EMAIL PROTECTED] on Tue, 13 Dec 2005 12:51:34 -0800 writes: Trevor It would be nice to have a date stamp on an object. Trevor In S/Splus this was always available, because objects were

Re: [R] Fitting binomial lmer-model, high deviance and low logLik

2005-12-14 Thread Doran, Harold
If you suspect a local maxima, have you tried different starting to values to see if the likelihood is maximized in the same place? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ivar Herfindal Sent: Wednesday, December 14, 2005 5:34 AM To:

Re: [R] Fitting binomial lmer-model, high deviance and low logLik

2005-12-14 Thread vmuggeo
Hi, I am not able to explain fully your results..However note that the deviance obtained in GLM with binary data (i.e Bernoulli 0/1) is meaningless..you should group your observations to get a valid GoF-type statistic. Point estimates are OK, of course. regards, vito Hello I have a problem

Re: [R] superimpose density line over hist

2005-12-14 Thread Romain Francois
Le 14.12.2005 11:54, Martin Maechler a écrit : Romain == Romain Francois [EMAIL PROTECTED] on Tue, 13 Dec 2005 15:40:59 +0100 writes: Romain A few comments : Romain - your code should be reproductible, otherwise it is useless. (that Romain

[R] Plotting model over data

2005-12-14 Thread Kare Edvardsen
Dear R-list. How would I plot model result over actual data? Here are the lme results: Linear mixed-effects model fit by REML Data: D-set AIC BIClogLik 884.3326 911.428 -432.1663 Random effects: Formula: ~C | ID Structure: General positive-definite, Log-Cholesky

Re: [R] Age of an object?

2005-12-14 Thread Kjetil Brinchmann Halvorsen
Philippe Grosjean wrote: Martin Maechler wrote: Trevor == Trevor Hastie [EMAIL PROTECTED] on Tue, 13 Dec 2005 12:51:34 -0800 writes: Trevor It would be nice to have a date stamp on an object. One way to do this with important objects is to use the comment function (in package base)

Re: [R] Design library

2005-12-14 Thread Ido M. Tamir
I'm a PhD student in Edinburgh University and interested in survival analysis. By chance, I found design library in Splus software. I have if yes, how can I find a book to introduce with these concepts? would you please introduce me some references? best wishes In R there is a pcakge called

Re: [R] Age of an object?

2005-12-14 Thread Kjetil Brinchmann Halvorsen
Philippe Grosjean wrote: Martin Maechler wrote: Trevor == Trevor Hastie [EMAIL PROTECTED] on Tue, 13 Dec 2005 12:51:34 -0800 writes: Trevor It would be nice to have a date stamp on an object. Following up on my post of a few minutes ago, I tried to write an timestamp function

Re: [R] Plotting model over data

2005-12-14 Thread Petr Pikal
Hi ?predict over ?order #(ed) data could be one way. HTH Petr On 14 Dec 2005 at 13:36, Kare Edvardsen wrote: Date sent: Wed, 14 Dec 2005 13:36:42 +0100 From: Kare Edvardsen [EMAIL PROTECTED] To: R-help r-help@stat.math.ethz.ch Subject:

[R] [R-pkgs] New Quantian release with over 800 CRAN/BioC packages

2005-12-14 Thread Dirk Eddelbuettel
[ Reposting this here as this Quantian release contains - all CRAN packages as of December 10, 2005 * except the Windows only ones (BRugs, mimR, rcom, RWinEdt) * and three that would not build (Rlsf, ROracle, rJava) for lack of Java, Oracle or rlsf - all BioC packages

[R] unable to force the vector format

2005-12-14 Thread Charles Plessy
Dear all, I am so ashamed to pollute the list with a trivial question, but it is a long time I have not used R, and I need a result in the next one or two hour... I have a table which I have loaded with read.table, and I want to make the mean of its columns. slides - read.table(slides.txt)

[R] Kalman Filter Forecast using 'SSPIR'

2005-12-14 Thread Sumanta Basak
Dear R Users, I am new to state-space modeling. I am using SSPIR package for Kalman Filter. I have a data set containing one dependent variable and 7 independent variables with 250 data points. I want to use Kalman Filter for forecast the future values of the dependent

Re: [R] Age of an object?

2005-12-14 Thread Roger Bivand
On Wed, 14 Dec 2005, Martin Maechler wrote: Kjetil == Kjetil Brinchmann Halvorsen [EMAIL PROTECTED] on Wed, 14 Dec 2005 08:59:24 -0400 writes: Kjetil Philippe Grosjean wrote: Martin Maechler wrote: Trevor == Trevor Hastie [EMAIL PROTECTED] on Tue, 13 Dec 2005

Re: [R] Translation Dictionary (Rosetta Stone)

2005-12-14 Thread Gabor Grothendieck
There is R for Octave users at: http://cran.r-project.org/doc/contrib/R-and-octave-2.txt On 12/14/05, Doran, Harold [EMAIL PROTECTED] wrote: Dear List: A while back I found a web site called the Rosetta Stone for Computer Algebra Systems showing how to do similar operations in different

[R] R for Windows server

2005-12-14 Thread Marco Geraci
Hi, I need to run a burdensome simulation. I'm using a Dell PowerEdge 2850, 3.8 Ghz CPU, 2 Gb RAM, Windows server 2003, 32 bit. The task manager indicates that Rgui.exe is using 25% of the CPU speed. There is no other process that is slowing down R. Does anyone know if it's a limitation of

Re: [R] Age of an object?

2005-12-14 Thread Duncan Temple Lang
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 In addition to what and how to associate time information with an object, I imagine Trevor was also expecting that it would be done automatically rather than require the user to set the timestamp explicitly. To that end, the addTaskCallback()

Re: [R] Age of an object?

2005-12-14 Thread Duncan Temple Lang
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Roger Bivand wrote: On Wed, 14 Dec 2005, Martin Maechler wrote: Kjetil == Kjetil Brinchmann Halvorsen [EMAIL PROTECTED] on Wed, 14 Dec 2005 08:59:24 -0400 writes: Kjetil Philippe Grosjean wrote: Martin Maechler wrote: Trevor

[R] concatenating expressions and standard text

2005-12-14 Thread manuel.martin
Hi all, is it possible to concatenate expressions and basic text when for instance labeling axis of a plot? I would like to see something like the concatenation of expression(C[0]) and for case 1 on my x axis. Obviously a plot(x, y, xlab=paste(expression(C[0]), in case1)) will not work. Thank

Re: [R] R for Windows server

2005-12-14 Thread Prof Brian Ripley
Is hyperthreading on? In that case the maximum is 50% Are there two processors? In that case the maximum is 25% On Wed, 14 Dec 2005, Marco Geraci wrote: I need to run a burdensome simulation. I'm using a Dell PowerEdge 2850, 3.8 Ghz CPU, 2 Gb RAM, Windows server 2003, 32 bit. The

[R] Memory shortage running Repeated Measures (nlme)

2005-12-14 Thread Evandro do Nascimento Silva
I forgot tomention R version (it is 2.0.1) and OS (Windows XP Professional 2002) in my previous message. Cheers, Evandro Silva __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Age of an object?

2005-12-14 Thread Martin Maechler
Roger == Roger Bivand [EMAIL PROTECTED] on Wed, 14 Dec 2005 15:14:47 +0100 (CET) writes: Roger On Wed, 14 Dec 2005, Martin Maechler wrote: Kjetil == Kjetil Brinchmann Halvorsen [EMAIL PROTECTED] on Wed, 14 Dec 2005 08:59:24 -0400 writes: Kjetil Philippe

Re: [R] R for Windows server

2005-12-14 Thread Philippe Grosjean
Hello, You don't tell us much info (what processor(s) are exactly in your server?). You must know the following: 1) Windows server prioritizes client accesses over local running programs (but you say there is no other process slowing down R... is it deconnected from the network?), 2) Rgui

Re: [R] Design library

2005-12-14 Thread Frank E Harrell Jr
Mohammad Reza wrote: Dear friends Hello I'm a PhD student in Edinburgh University and interested in survival analysis. By chance, I found design library in Splus software. I have some questions about it and it's highly appreciated if you can help me. In its help (Design library help) I

Re: [R] Age of an object?

2005-12-14 Thread Seth Falcon
On 14 Dec 2005, [EMAIL PROTECTED] wrote: This does treat any existing comment rather brutally, could stamp rather be: stamp - paste(Sys.time(), comment(obj), value) probably enhanced with some field separators to let the inspector grab just its chunk? Something like DCF? This does treat

Re: [R] concatenating expressions and standard text

2005-12-14 Thread Kristel Joossens
Can you give an example of the expression you would like to have on you horizontal axis? Simple expressions like ex1 - expression(1+ 0:9) seem to work well. Kristel manuel.martin wrote: Hi all, is it possible to concatenate expressions and basic text when for instance labeling axis of a

[R] Looking for a sort of tapply() to data frames

2005-12-14 Thread January Weiner
Hi, I read about the by() function, but it does not seem to do the job I need. Here is the problem: Say - I have a data frame, with three columns. The first one contains strings that describe the data points, with repeats (for example, days of a week). The other two contain numbers. Something

[R] Problem with dir.create (R2.2.0 Windows XP 2002 SP 2)

2005-12-14 Thread hadley wickham
I've run into a problem with dir.create on R2.2.0 Windows XP 2002 SP 2. setwd(d:/) print(dir.create(d:\\otis-sim\\rdata, recursive=T)) print(dir.create(d:\\otis-sim\\, recursive=T)) Both return false and fail to create the directories. setwd(c:/) print(dir.create(d:\\otis-sim\\rdata,

[R] Memory shortage running Repeated Measures (nlme)

2005-12-14 Thread Evandro do Nascimento Silva
I forgot to mention R version (which is 2.0.1) and OS (Windows XP Professional 2002) in my previous message. Cheers, Evandro Silva __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Age of an object?

2005-12-14 Thread Peter Dalgaard
[Got stuck in the mail filter: Too many recipients. Cancelled and resending.] Duncan Temple Lang [EMAIL PROTECTED] writes: In addition to what and how to associate time information with an object, I imagine Trevor was also expecting that it would be done automatically rather than require

Re: [R] Looking for a sort of tapply() to data frames

2005-12-14 Thread Thomas Lumley
On Wed, 14 Dec 2005, January Weiner wrote: Hi, I read about the by() function, but it does not seem to do the job I need. Here is the problem: by() will work, you just need to use the right function in it. You want by(df[,-1], df$Day, function.that.means.each.column) so all you need to do

[R] Translation Dictionary (Rosetta Stone)

2005-12-14 Thread Doran, Harold
Dear List: A while back I found a web site called the Rosetta Stone for Computer Algebra Systems showing how to do similar operations in different symbolic programs. This has been very helpful for moving between Mathematica, Yacas and Maxima. I wonder if anyone maintains something akin to this

[R] Append tables

2005-12-14 Thread Bill Hunsicker
R Help: I have read a number of tables into R with identical headings and I would now like to make a single table that has all the data appended under this single heading line. for example: t1 - read.csv(f1,header=TRUE) t2 - read.csv(f2,header=TRUE) all - c(t1,t2) #all is now twice as

[R] The fastest way to select and execute a few selected functions inside a function

2005-12-14 Thread Ales Ziberna
Dear useRs? I have the following problem! I have a function that calls one or more functions, depending on the input parameters. I am searching for the fastest way to select and execute the selected functions and return their results in a list. The number of possible functions is 10, however

Re: [R] Append tables

2005-12-14 Thread Ferdinand Alimadhi
?rbind Bill Hunsicker wrote: R Help: I have read a number of tables into R with identical headings and I would now like to make a single table that has all the data appended under this single heading line. for example: t1 - read.csv(f1,header=TRUE) t2 - read.csv(f2,header=TRUE) all -

Re: [R] Append tables

2005-12-14 Thread Don MacQueen
Does rbind() do what you want? At 11:47 AM -0500 12/14/05, Bill Hunsicker wrote: R Help: I have read a number of tables into R with identical headings and I would now like to make a single table that has all the data appended under this single heading line. for example: t1 -

Re: [R] Labeling a range of bars in barplot?

2005-12-14 Thread Dan Bolser
Marc Schwartz (via MN) wrote: On Tue, 2005-12-13 at 10:53 +, Dan Bolser wrote: Hi, I am plotting a distribution of (ordered) values as a barplot. I would like to label groups of bars together to highlight aspects of the distribution. The label for the group should be the range of values

Re: [R] unable to force the vector format

2005-12-14 Thread Charles Plessy
On Wed, Dec 14, 2005 at 02:46:35PM +0100, Petr Pikal wrote : Hi is colMeans(your.df) what you want? BTW searching column mean in CRAN gives me the answer on first hit Thank you all so much for all the answers I got. rowMeans is what I needed. I will bookmark the search page of

[R] ANCOVA Post-hoc test

2005-12-14 Thread Nicolas Poulet
Hello, Despite my search, I didn't find a post-hoc test for an ANCOVA. I used the functions aov() and lm() to run the ANCOVA then I tried TukeyHSD() but it didn't work (because of the covariable is a continuous variable?). Furthermore, I would like to plot the adjusted values (i.e. the values

Re: [R] getting faster results

2005-12-14 Thread Elizabeth Lawson
I am running algorithms for Bayesian hierarchical models with 20 or more variables or writing alogrithms that solve for partially linear models using wavlelet estimation. Soemtimes just creating the predicted vaules from the Bayesian analysis takes over a day to run! I am just wondering if R

Re: [R] Age of an object?

2005-12-14 Thread Dave Roberts
This would be extraordinarily helpful, but I have not thought of a graceful way to do it. Everything in R now has a class attribute, but a timestamp for such simple things as vectors seems like overkill. On the other hand, those of us writing packages could implement this pretty easily for

Re: [R] unable to force the vector format

2005-12-14 Thread Ido M. Tamir
slides - read.table(slides.txt) slides [1:5,] V1V2 V3 V4 V5 V6 V7 V8 1 PLB00090AA02 0.147 0.018 0.046 0.064 -0.018 -0.008 -0.063 2 PLB00090BC08 0.171 0.011 -0.001 0.009 0.052 0.032 -0.065 3 PLB00090CG02 0.029 -0.014 -0.042 0.006 0.024 -0.009

[R] suggestions for nls error: false convergence

2005-12-14 Thread Rajarshi Guha
Hi, I'm trying to fit some data using a logistic function defined as y ~ a * (1+m*exp(-x/tau)) / (1+n*exp(-x/tau) My data is below: x - 1:100 y - c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,1,1,1,2,2,2,2,2,3,4,4,4,5, 5,5,5,6,6,6,6,6,8,8,9,9,10,13,14,16,19,21,

Re: [R] getting faster results

2005-12-14 Thread Berton Gunter
What is almost certainly critical here are the algorithms you use (to do the integrations in the Bayesian models) and perhaps a fast BLAS. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The business of the statistician is to catalyze the scientific learning process.

Re: [R] getting faster results

2005-12-14 Thread Liaw, Andy
It can depend on a whole list of things why the computations were taking so long, and just saying what methods you are using to do the computation is far from sufficient, as there are efficient ways and inefficient ways to do any computations, and we don't know how you did any of them. With the

Re: [R] getting faster results

2005-12-14 Thread Elizabeth Lawson
I am not sure... But my employer would like to get me a new machine and considering that R is the program I use most often, he was wondering what type of machine would be the best to get form me. jim holtman [EMAIL PROTECTED] wrote: If you job is CPU intensive (it should be unless it is

[R] Need help for a statistical problem

2005-12-14 Thread james brown
Hello Dear I need to select among 700 objects a good representative sample. These objects could be residential houses, commercial buildings, trucks, etc. How to get a good sample size and select a set of objects that is very representative. The second part of my question is to find a statistical

Re: [R] Age of an object?

2005-12-14 Thread Brahm, David
Trevor Hastie [EMAIL PROTECTED] wrote: It would be nice to have a date stamp on an object. In S/Splus this was always available, because objects were files. The g.data package on CRAN stores R objects in individual files, like the (old) S-Plus model. The timestamp on a file tells you the last

[R] SAS.xpt/STATA.dta, field descriptions, and dbWriteTable

2005-12-14 Thread Brett Magill
Hello all, I have a large database (~100MB in 13 relational files) that made its way to me in two formats. The files were sent in both SAS xport and STATA, thanks, I am told, to stat-transfer. The two files, ostensibly, contain the same data, just in different formats. My goal: Move these

[R] memory tops out at 1.84gb on OS X 10.4 machine w/ 5GB ram

2005-12-14 Thread Ken Termiso
Hi all, Sorry if this is a dumb question, but I am on 10.4 with R2.2, and when loading a big text file (~500MB) with scan(file, what=character) I am throwing malloc errors that say I am out of memory...I have 5GB on this machine, and Activity Monitor tells me R is only up to ~1.84GB both times

Re: [R] Labeling a range of bars in barplot?

2005-12-14 Thread Gabor Grothendieck
Note that if I follow this correctly then you could remove the loop. In particular note that 1. st is just the cumulative sum of new.break.points but summed from the end: st - rev(cumsum(rev(my.new.breaks))) 2. segments and text both take vector arguments and 3. averaging over the groups

Re: [R] The fastest way to select and execute a few selected func tions inside a function

2005-12-14 Thread Liaw, Andy
See if the following helps you any: flist - list(min=min, mean=mean, median=median, max=max, n=length) x - rnorm(1000) system.time(for (i in 1:1e4) sapply(flist, function(f) f(x))) [1] 5.90 0.01 6.05 NA NA system.time(for (i in 1:1e4) sapply(flist[c(mean, median)], function(f) f(x))) [1]

Re: [R] Problem with dir.create (R2.2.0 Windows XP 2002 SP 2)

2005-12-14 Thread Liaw, Andy
No problem for me (Version 2.2.0 (2005-10-06 r35749), WinXPPro SP2): setwd(q:/) print(dir.create(q:/Andy/what/is/this, recursive=TRUE)) [1] TRUE setwd(c:/) print(dir.create(q:/Andy/what/is/that/thing, recursive=TRUE)) [1] TRUE Andy From: hadley wickham I've run into a problem with

Re: [R] Problem with dir.create (R2.2.0 Windows XP 2002 SP 2)

2005-12-14 Thread Michael H. Prager
I hope this isn't an insulting question, but does drive D: exist on that machine? Mike on 12/14/2005 11:09 AM hadley wickham said the following: I've run into a problem with dir.create on R2.2.0 Windows XP 2002 SP 2. setwd(d:/) print(dir.create(d:\\otis-sim\\rdata, recursive=T))

Re: [R] Problem with dir.create (R2.2.0 Windows XP 2002 SP 2)

2005-12-14 Thread hadley wickham
I hope this isn't an insulting question, but does drive D: exist on that machine? Yes, it definitely does! I'd expect setwd(d:/) to give an error if it didn't Hadley __ R-help@stat.math.ethz.ch mailing list

Re: [R] dendrogram: how to obtain leaf height

2005-12-14 Thread Dr. Herwig Meschke
try height.of.leafs - dendrapply(hcd, function(e) attr(e, height)) unlist(height.of.leafs) Regards Herwig -- Dr. Herwig Meschke Wissenschaftliche Beratung Hagsbucher Weg 27 D-89150 Laichingen phone +49 7333 210 417 / fax +49 7333 210 418 email [EMAIL PROTECTED] On 12 Dec 2005 at 21:15,

Re: [R] Problem with dir.create (R2.2.0 Windows XP 2002 SP 2)

2005-12-14 Thread Michael H. Prager
on 12/14/2005 4:43 PM hadley wickham said the following: I hope this isn't an insulting question, but does drive D: exist on that machine? Yes, it definitely does! I'd expect setwd(d:/) to give an error if it didn't Hadley H Oops. You are right, of course. I checked that on

Re: [R] Problem with dir.create (R2.2.0 Windows XP 2002 SP 2)

2005-12-14 Thread Prof Brian Ripley
One comment is that Windows does not like trailing directory separators as in print(dir.create(d:\\otis-sim\\, recursive=T)) and so this will give a warning as in print(dir.create(d:\\otis-sim\\, recursive=T)) [1] FALSE Warning message: 'd:\otis-sim\' already exists Beyond that, the code is

Re: [R] Problem with dir.create (R2.2.0 Windows XP 2002 SP 2)

2005-12-14 Thread Peter Dalgaard
hadley wickham [EMAIL PROTECTED] writes: I hope this isn't an insulting question, but does drive D: exist on that machine? Yes, it definitely does! I'd expect setwd(d:/) to give an error if it didn't What is the file system? I seem to recall that the top directory of a (V)FAT drive is

[R] 32 vs 64 bit

2005-12-14 Thread Paul Gilbert
I ran R on an AMD Athlon 64 with a 32 bit Linux (by accident, I moved the hard disk from another machine). Now, after install 64 bit Linux and 64 bit R, I am finding some of my R tests are quite a bit slower, I think because memory demands cause a lot of swap. Does anyone know if extra memory

Re: [R] Problem with dir.create (R2.2.0 Windows XP 2002 SP 2)

2005-12-14 Thread Prof Brian Ripley
On Wed, 14 Dec 2005, Peter Dalgaard wrote: hadley wickham [EMAIL PROTECTED] writes: I hope this isn't an insulting question, but does drive D: exist on that machine? Yes, it definitely does! I'd expect setwd(d:/) to give an error if it didn't What is the file system? I seem to recall

Re: [R] 32 vs 64 bit

2005-12-14 Thread Prof Brian Ripley
On Wed, 14 Dec 2005, Paul Gilbert wrote: I ran R on an AMD Athlon 64 with a 32 bit Linux (by accident, I moved the hard disk from another machine). Now, after install 64 bit Linux and 64 bit R, I am finding some of my R tests are quite a bit slower, I think because memory demands cause a lot

Re: [R] unable to force the vector format

2005-12-14 Thread jim holtman
Not too sure what you want the means of, but try ?colMeans ?rowMeans On 12/14/05, Charles Plessy [EMAIL PROTECTED] wrote: Dear all, I am so ashamed to pollute the list with a trivial question, but it is a long time I have not used R, and I need a result in the next one or two hour... I

[R] X11 png jpeg cledit false when running from script

2005-12-14 Thread Jon Dressel
Hello, When launching R from a script, when queried it reports the following capabilities as false: jpeg, png, x11 and cledit. When R is run from a terminal session, all of these capabilities are reported true. R is running on FC4. As the purpose of this script is to output the png file to a

Re: [R] unable to force the vector format

2005-12-14 Thread Petr Pikal
Hi is colMeans(your.df) what you want? BTW searching column mean in CRAN gives me the answer on first hit HTH Petr On 14 Dec 2005 at 22:12, Charles Plessy wrote: Date sent: Wed, 14 Dec 2005 22:12:16 +0900 From: Charles Plessy [EMAIL PROTECTED] To:

Re: [R] Age of an object?

2005-12-14 Thread Henrik Bengtsson
If R would have timestamps telling when any object was last modified, we could extend R with a 'GNU make'-style functionality (or syntax) together with some fancy caching to persistent storage (files, data bases, ...). That would really nice! As B.R. and M.M. writes, timestamping is most

Re: [R] The fastest way to select and execute a few selected functions inside a function

2005-12-14 Thread Ferdinand Alimadhi
Is this what you want? myf-function(FUN, x){ res- vector( mode=list) for(i in FUN) res[[i]]-do.call(i,args=list(x)) return(res) } myf(FUN=mean,x=1:10) myf(FUN=c(mean,max),x=1:10) myf(FUN=c(mean,max, sum, median),x=1:10) Ales

[R] no subject

2005-12-14 Thread Marco Blanchette
Dear all, I am still fairly new to R and try to analyze large tables of data generated from genomic experiment. Currently, I am trying to plot pair of experiments coming from different file, trying to look at the behavior of individual feature in pair of experiment. My problem is that I have

Re: [R] Age of an object?

2005-12-14 Thread Martin Maechler
Kjetil == Kjetil Brinchmann Halvorsen [EMAIL PROTECTED] on Wed, 14 Dec 2005 08:59:24 -0400 writes: Kjetil Philippe Grosjean wrote: Martin Maechler wrote: Trevor == Trevor Hastie [EMAIL PROTECTED] on Tue, 13 Dec 2005 12:51:34 -0800 writes: Trevor It would be nice

[R] Delete missing values

2005-12-14 Thread John Sorkin
I am trying to delete rows containing missing values from a groupeddata object. Several of the columns are character (sexChar, HAPI, rs2304785) the rest are numeric. For some reason I am excluding all rows with missing values. Your suggestions for corrections would be appreciated. This did not

Re: [R] Merging lists on common key (was no subject)

2005-12-14 Thread Marc Schwartz
On Wed, 2005-12-14 at 18:14 -0800, Marco Blanchette wrote: Dear all, I am still fairly new to R and try to analyze large tables of data generated from genomic experiment. Currently, I am trying to plot pair of experiments coming from different file, trying to look at the behavior of

Re: [R] X11 png jpeg cledit false when running from script

2005-12-14 Thread Marc Schwartz
On Wed, 2005-12-14 at 15:53 -0800, Jon Dressel wrote: Hello, When launching R from a script, when queried it reports the following capabilities as false: jpeg, png, x11 and cledit. When R is run from a terminal session, all of these capabilities are reported true. R is running on FC4.

Re: [R] Delete missing values

2005-12-14 Thread Marc Schwartz
On Wed, 2005-12-14 at 21:34 -0500, John Sorkin wrote: I am trying to delete rows containing missing values from a groupeddata object. Several of the columns are character (sexChar, HAPI, rs2304785) the rest are numeric. For some reason I am excluding all rows with missing values. Your

Re: [R] no subject

2005-12-14 Thread Jacques VESLOT
try with merge() : table1 = data.frame(CGID=c(CG_1,CG_3,CG_2, CG_4, CG_5), diff=c(3,5,6,4,3)) table2 = data.frame(CGID=c(CG_2,CG_3,CG_4, CG_1, CG_5), diff=c(4,6,3,9,10)) newtable- merge(table1, table2, by=CGID) matplot(merge(table1,

[R] residual and null deviance of an lme object with correlation structure

2005-12-14 Thread Corey Bradshaw
Hello, I am attempting to calculate the residual and null deviance of an lme object that includes a corAR1 correlation structure. I tried deviance(lme.object) and it only returned NULL. Can anyone help? Thank you. Corey Bradshaw __

Re: [R] X11 png jpeg cledit false when running from script

2005-12-14 Thread Prof Brian Ripley
It could be as simple as not having DISPLAY set when running the script. capabilities() first checks there is a display specified, then that it can be opened. SELinux (which my FC3 box runs) does from time to time get in the way (a month or so ago it stopped scripts accessing Xvfb), but

Re: [R] Need help for a statistical problem

2005-12-14 Thread Uwe Ligges
james brown wrote: Hello Dear I need to select among 700 objects a good representative sample. These objects could be residential houses, commercial buildings, trucks, etc. How to get a good sample size and select a set of objects that is very representative. The second part of my

[R] Job Ad: Centocor RD, USA

2005-12-14 Thread Pikounis, Bill [CNTUS]
(Apologies if this came through in HTML format, as my email client server configuration does not send Plain Text reliably.) LEAD STATISTICIAN Centocor R D, an operating company of Johnson Johnson, seeks a highly motivated statistician/data analyst to work in its Nonclinical Statistics group.