Re: [R] Backward logistic regression

2009-07-03 Thread Tal Galili
If the daya set is not too big, I believe the easiest will be: stepAIC from the MASS package... On Thu, Jul 2, 2009 at 11:21 PM, Angela Radulescu angela.radule...@gmail.com wrote: Hi, I have performed a binary logistic regression using the glm ( family = binomial(logit) )command on 6

Re: [R] How to use current value of variable in function definition?

2009-07-03 Thread Mark Wardle
2009/7/2 Peter Dalgaard p.dalga...@biostat.ku.dk: Mark Wardle wrote: [...] Otherwise, it sounds like a recipe for obfuscation! What are you trying to do really? Not really obfuscation, there are legitimate reasons to want this sort of behaviour. See demo(scoping) for a related example

[R] Query

2009-07-03 Thread MINI GHOSH
Dear R-User, Is it possible to simulate Impulsive differential equation using R? If yes, please let me know how? Thanking you, Best Regards, Mini Ghosh See the Web#39;s breaking stories, chosen by people like you. Check out Yahoo! Buzz. http://in.buzz.yahoo.com/

Re: [R] Package 'operators'

2009-07-03 Thread Romain Francois
My point is your problem is not entirely related to this package in particular but more should I trust additional R packages, ie should I use package cron when R can handle date objects, should I use ggplot when i could do it all in grid, ... operators comes with unit tests that allow me to

Re: [R] Package 'operators'

2009-07-03 Thread Barry Rowlingson
On Thu, Jul 2, 2009 at 11:57 PM, seeliger.c...@epamail.epa.gov wrote: Thanks for replying, Francois. To directly answer your question, the difference between using base R functions and a library comes down to code correctness and stability, as well as future support.  Portions of the

[R] Wiki down?

2009-07-03 Thread Mario Valle
From yesterday I cannot connect anymore to wiki.r-project.org Is the problem known? Thanks! mario -- Ing. Mario Valle Data Analysis and Visualization Group| http://www.cscs.ch/~mvalle Swiss National Supercomputing Centre (CSCS) | Tel: +41 (91) 610.82.60 v.

[R] Lo's modified R/S

2009-07-03 Thread mihaela
Hi all, Does anyone know if Lo's modified R/S statistic is implemented in R? I found rsFit function from fArma package (LrdModelling), but I guess this is to calculate the classical R/S statistic. Thank you, M. -- View this message in context:

Re: [R] Problem with groupedData and lme

2009-07-03 Thread Dieter Menne
MUHC-Research wrote: I'm currently having trouble with the implementation of a groupedData object in the lme() function. Executing the following function applyScalingSimp - function(input.population) { ## GA is a time value varInOrder - c(GA,weight,grouping,sex)

Re: [R] Wiki down?

2009-07-03 Thread Philippe Grosjean
Yes, it is known. It is a major power supply fealure during all night... and UPS were exhausted. The site is now back again. Sorry for the inconvenience. Best, Philippe Grosjean Mario Valle wrote: From yesterday I cannot connect anymore to wiki.r-project.org Is the problem known? Thanks!

Re: [R] sum of list elements

2009-07-03 Thread Patrick Burns
My guess is that the original poster is expecting a matrix not a scalar. If so, then: Reduce('+', a) [,1] [,2] [,3] [1,]300 [2,]030 [3,]003 Patrick Burns patr...@burns-stat.com +44 (0)20 8525 0696 http://www.burns-stat.com (home of The R Inferno and

[R] Presumably simple question about sorting/ordering

2009-07-03 Thread Colin Beale
Hi, I'm stuck with a fairly basic re-ordering problem. I want to extract part of a matrix as a vector, and reorder it to match with an erratic sequence of x,y coordinates. Here's an example that shows what I want and how close I've got (but for some reason my mind just can't get the final

Re: [R] mapping states with colors

2009-07-03 Thread Paul Hiemstra
Hi, You could have a look at the sp-package, it provides R-classes for spatial data. In addition, it also provides very nice plotting facilities. To read your data into these sp classes you can use the rgdal package. There is also an R mailing list specifically for geographic data

Re: [R] A problem on zoo object

2009-07-03 Thread Bogaso
Thanks Gabor for this reply. However can please clarify one more thing? How I want to create a list, wherein each member of that list is the monthly observations. For example, 1st member of list contains daily observation of 1st month, 2nd member contains daily observation of 2nd month etc.

Re: [R] to creates an array

2009-07-03 Thread Jim Lemon
barbara.r...@uniroma1.it wrote: Is there a command as mat.or.vec(nr,nc) to create an array that I must calculate with more cicle? Hi Barbara, The approach that I usually use is something like this: new.matrix-matrix(NA,nrow=3,ncol=3) This creates a 3x3 matrix named new.matrix that is

[R] Correlation Network Diagram?

2009-07-03 Thread rory . winston
Hi all On page 39 of this paper [1] by Andrew Lo there is a very interesting correlation network diagram (sorry I dont have a more convenient link to the type of diagram I'm talking about). does anyone know of any package in R that can generate these types of diagrams? Cheers -- Rory [1]

[R] Correlation Network Diagram - Correction

2009-07-03 Thread rory . winston
Sorry - In my previous mail, I said page 39, when I actually meant page 36. Cheers -- Rory [[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

[R] Function to eliminate blank space within strings?

2009-07-03 Thread Agustin Lobo
Is there an specific function to eliminate blank space within strings? Thanks -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: agustin.l...@ija.csic.es http://www.ija.csic.es/gt/obster

[R] Confidence Limits for a Cross-Product Ratio

2009-07-03 Thread francogrex
Data from Fisher's paper: Confidence Limits for a Cross-Product Ratio. y col1 col2 [1,] 103 [2,]2 15 fisher.test(y) Fisher's Exact Test for Count Data data: y p-value = 0.0005367 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence

Re: [R] Function to eliminate blank space within strings?

2009-07-03 Thread Romain Francois
On 07/03/2009 10:34 AM, Agustin Lobo wrote: Is there an specific function to eliminate blank space within strings? Thanks gsub( [[:space:]]+, , string ) -- Romain Francois Independent R Consultant +33(0) 6 28 91 30 30 http://romainfrancois.blog.free.fr |- http://tr.im/qJ8V : RGG#153: stars

Re: [R] Function to eliminate blank space within strings?

2009-07-03 Thread jim holtman
How much blank space do you want to remove? all to one, all to none, leading, trailing? You can use regular expressions: x -mytest string # leading sub(^ +, , x) [1] mytest string # trailing sub( +$, '', x) [1]mytest string # multiple to one

[R] fix() and edit() not working with Rcmdr and german LANG-variable

2009-07-03 Thread Jannis Vajen
Dear Mailinglist, I just set up an R 2.9.1 environment with Rcmdr 1.4-6 on Ubuntu Jaunty 9.04. As I'm from Germany my $LANG variable is set to de_DE.UTF-8. Now, when I open up Rcmdr and try to edit a new datamatrix there is no edit window appearing: Datenmatrix - edit(as.data.frame(NULL)) ERROR:

Re: [R] Function to eliminate blank space within strings?

2009-07-03 Thread Duncan Murdoch
Agustin Lobo wrote: Is there an specific function to eliminate blank space within strings? No, but gsub could do that. Duncan Murdoch __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] Remove all spaces from a string so it can be used by assign()

2009-07-03 Thread Kurt Smith
Hi I have a string 56 Fe [1] that I would like to use as a variable name by using assign however I think the spaces and brackets might be causing R some trouble. How can I change the string so that it just becomes 56Fe1 and can be used as a variable name. Thank You Kurt

Re: [R] Remove all spaces from a string so it can be used by assign()

2009-07-03 Thread Romain Francois
On 07/03/2009 12:56 PM, Kurt Smith wrote: Hi I have a string 56 Fe [1] that I would like to use as a variable name by using assign however I think the spaces and brackets might be causing R some trouble. Not really; R assign( 56 Fe [1], 10 ) R `56 Fe [1]` [1] 10 How can I change the

[R] fix() and edit() not working with Rcmdr and german LANG-variable

2009-07-03 Thread Jannis Vajen
Dear Mailinglist, I just set up an R 2.9.1 environment with Rcmdr 1.4-6 on Ubuntu Jaunty 9.04. As I'm from Germany my $LANG variable is set to de_DE.UTF-8. Now, when I open up Rcmdr and try to edit a new datamatrix there is no edit window appearing: Datenmatrix - edit(as.data.frame(NULL)) ERROR:

Re: [R] save the result into a word file

2009-07-03 Thread Allan Engelhardt
Some suggestions in increasing complexity: 1. Cut-and-paste from the R console. (Cut-and-paste has got to be the best computer invention ever.) 2. If you want to generate the output automatically from an unattended script, try help(cat) and help(file) 3. Try help(Sweave) for even more

Re: [R] Remove all spaces from a string so it can be used by assign()

2009-07-03 Thread Romain Francois
On 07/03/2009 02:13 PM, Kurt Smith wrote: Are you sure it would work? It works when I physically enter 56 Fe [1] but fails when I try to enter that anything other then directly. instrument.input[1,6] [1] 56 Fe [ 1 ] 43 Levels: 10840277.06 109014288.37 11055630.67 11456522.47 ... CPS

Re: [R] Remove all spaces from a string so it can be used by assign()

2009-07-03 Thread Kurt Smith
Thanks for the help, I really appreciate it. However, I have one more small issue and I don't really want to post a second message in such a short space of time, so I will put it up here if someone is willing to help. Basically I need to transfer two vectors into a data.frame/matrix (with a

[R] R on flash

2009-07-03 Thread Raphael Saldanha
Hi! Is there some R version for flashdisk, pendrive etc ? -- Regards, Raphael Saldanha saldanha.plan...@gmail.com [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Odp: R on flash

2009-07-03 Thread Petr PIKAL
Quicvk and dirty Insall R to PC and copy the appropriate directory to flash disk. Regards Petr r-help-boun...@r-project.org napsal dne 03.07.2009 14:34:59: Hi! Is there some R version for flashdisk, pendrive etc ? -- Regards, Raphael Saldanha saldanha.plan...@gmail.com

Re: [R] Error: cannot allocate vector of size 332.3 Mb

2009-07-03 Thread Steve Ellis
R-helpers, I thank Jonathan Greenberg and David Winsemius for their responses. I will keep R64.app in mind but I found that by deleting some large objects that I didn't need I was able to do my computations using R 2.9.1. (This is consistent with Winsemius's experiment on a 10GB

[R] Computing Time Intervals On A Series

2009-07-03 Thread rory . winston
Hi I have a dataset that looks like this (dput'd below): head(x, 20) time status 1 2009-07-02 10:32:37 1 2 2009-07-02 10:32:43 0 3 2009-07-02 10:32:43 1 4 2009-07-02 10:32:44 0 5 2009-07-02 10:32:44 1 6 2009-07-02 10:32:48 0 7 2009-07-02 10:32:48 1 8 2009-07-02 10:32:54 0 9 2009-07-02 10:33:04

Re: [R] bigglm() results different from glm()

2009-07-03 Thread utkarshsinghal
Hi Sir, Thanks for making package available to us. I am facing few problems if you can give some hints: Problem-1: The model summary and residual deviance matched (in the mail below) but I didn't understand why AIC is still different. AIC(m1) [1] 532965 AIC(m1big_longer) [1] 101442.9

Re: [R] Correlation Network Diagram - Correction

2009-07-03 Thread Allan Engelhardt
Try the network package from the statnet suite of packages. See also http://www.jstatsoft.org/v24/i02 with page 27 having similar graphs. The statnet tutorial in the same issue is recommended. Home page is http://csde.washington.edu/statnet/. Allan On 03/07/09 11:43,

Re: [R] Correlation Network Diagram?

2009-07-03 Thread Ronggui Huang
Maybe not address your question directly. There is a package for correlation visualization in r-forge (http://r-forge.r-project.org/projects/corrplot/). 2009/7/3 rory.wins...@gmail.com: Hi all On page 39 of this paper [1] by Andrew Lo there is a very interesting correlation network diagram

[R] is AIC always 100% in evaluating a model?

2009-07-03 Thread alexander russell
Hello, I'd like to say that it's clear when an independent variable can be ruled out generally speaking; on the other hand in R's AIC with bbmle, if one finds a better AIC value for a model without the given independent variable, versus the same model with, can we say that the independent variable

Re: [R] Computing Time Intervals On A Series

2009-07-03 Thread Gabor Grothendieck
Try this: d - diff(DF$time) d[rep(c(FALSE, TRUE), length = length(d))] Time differences in secs [1] 0 0 0 10 1 6 23 2 1 d[rep(c(TRUE, FALSE), length = length(d))] Time differences in secs [1] 6 1 4 6 0 2 4 8 1 0 On Fri, Jul 3, 2009 at 9:21 AM, rory.wins...@gmail.com wrote: Hi I

Re: [R] save the result into a word file

2009-07-03 Thread Liviu Andronic
Hello, On Fri, Jul 3, 2009 at 2:15 PM, Allan Engelhardtall...@cybaea.com wrote: 1. Cut-and-paste from the R console.  (Cut-and-paste has got to be the best computer invention ever.) First option is probably interesting when using in Word a monospace font for R results. For tables, you can also

Re: [R] Odp: R on flash

2009-07-03 Thread Prof Brian Ripley
We were not told the OS, but if this is Windows, see the rw-FAQ 'Can I run R from a CD or USB drive?' If it is a Unix-alike: the scheme below will not work and you need to edit R_HOME in the 'R' script (including the copy in R_HOME/bin) if you relocate the R file tree. That can be a problem

Re: [R] save the result into a word file

2009-07-03 Thread Ronggui Huang
How about export it to a HTML and then copy to MS word. I have packaged some code to export lm object to HTML in AJS style. You can have a look at ASAtable (http://r-forge.r-project.org/R/?group_id=328). Ronggui 2009/7/3 Suyan Tian st...@mail.rockefeller.edu: Hi, everyone: I forget on how to

Re: [R] Computing Time Intervals On A Series

2009-07-03 Thread David Winsemius
On Jul 3, 2009, at 9:21 AM, rory.wins...@gmail.com wrote: Hi I have a dataset that looks like this (dput'd below): head(x, 20) time status 1 2009-07-02 10:32:37 1 2 2009-07-02 10:32:43 0 3 2009-07-02 10:32:43 1 4 2009-07-02 10:32:44 0 5 2009-07-02 10:32:44 1 6 2009-07-02 10:32:48 0 7

Re: [R] Plot two graphs with different ranges in one

2009-07-03 Thread John Kane
Better yet use two graphs on the same page. It is damned hard to read and understand a two axes display. --- On Thu, 7/2/09, Mark Knecht markkne...@gmail.com wrote: From: Mark Knecht markkne...@gmail.com Subject: Re: [R] Plot two graphs with different ranges in one To: Rolf Turner

Re: [R] Computing Time Intervals On A Series

2009-07-03 Thread David Winsemius
I noticed that Gabor's interpretation of your desires was the opposite of what mine was, and also that you were thinking tapply might be a useful strategy. Here is a tapply implementation that lets you choose: x$wanted - c(0,-diff(x$status)) x$dif.tm - c(diff(x$time),NA) sum(x[x$wanted==1,

Re: [R] save the result into a word file

2009-07-03 Thread jim holtman
One technique that I use a screen capture that you can paste into the document. This works nicely since it will also capture the color/size of the fonts. Works well for the graphics that you want to include. If you have MS OneNote, you can use this. There is a freeware package (MWSnap) that I

[R] lattice, Hmisc: vertical placement of superpose.symbol in Dotplot()

2009-07-03 Thread John G. Bullock
I often use Dotplot() from the Hmisc package to construct dotplots. Typically, I plot multiple symbols on each line of a dotplot. Sometimes these symbols overlap. In these cases, I would like to eliminate the overlap by slightly adjusting the vertical placement of the plotted symbols. Is there

Re: [R] Correlation Network Diagram - Correction

2009-07-03 Thread Sigbert Klinke
Hi, the igraph package will be able to do the job. Yours sincerely Sigbert Klinke __ 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

[R] Color of ecdf plots

2009-07-03 Thread Lars Bergemann
Hi. I have the following two ecdf plots in one graph: plot( ecdf(), do.points=FALSE, verticals=TRUE, main=paste(Ecdf of distances ,DIM,sep=), col=red ); lines( ecdf(), do.points=FALSE, verticals=TRUE ); How do I change the color of the resulting graph?

Re: [R] Function to eliminate blank space within strings?

2009-07-03 Thread Sigbert Klinke
Hi, Is there an specific function to eliminate blank space within strings? Thanks gsub( , , txt) Sigbert __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] (no subject)

2009-07-03 Thread Nikolay Vasilev
Dear All, I did a new installment of R on my office computer. However, I encounter one technical problem with my R Help. When I write down ?, followed by the function I am interested in, the next message appears in Help Menu:This program cannot display the webpage. Please, find attached

Re: [R] Error

2009-07-03 Thread jim holtman
did you load the package with library(Rdonlp2) before running the example code? On Fri, Jul 3, 2009 at 9:53 AM, Yoonus CAyoonus...@ifmr.ac.in wrote: Dear Sir, I installed successfully using the install.packages(Rdonlp2, repos= http://R-Forge.R-project.org;). I tried on the example which is

Re: [R] Plot Fonts in Windows vs Mac OSX

2009-07-03 Thread Winston Chang
On Thu, Jul 2, 2009 at 7:11 PM, Sean Carmody seancarm...@gmail.com wrote: I have been plotting the same charts using png on a Windows machine and on a Mac OSX and the quality of the resulting images, particularly in relation to the fonts, look far superior in the plots produced on the Mac. Is

[R] Error

2009-07-03 Thread Yoonus CA
Dear Sir, I installed successfully using the install.packages(Rdonlp2, repos= http://R-Forge.R-project.org;). I tried on the example which is defined on the Rdonlp help file ,p - c(-10,10)par.l - c(0,0); par.u - c(100,100)lin.u - 1; lin.l - 1A - t(c(1,1))fn - function(x){ x[1]^2+x[2]^2}ret -

[R] Histogram

2009-07-03 Thread aledanda
Hallo, How can I have the control of the width of each bin in the histogram? I would like to plot my data modifying the number and the width of the bins. Is there a good pdf or manual where i can learn to use the plot/hist in R very well? I find the help not always clear Thanks a lot Ale --

[R] normalised curve fitting with error bars

2009-07-03 Thread Pooja Jain
Dear List, My data consist of nine columns and about 50,000 rows. It looks like this. -9.0225 3.46464 2.80926 -0.3847 3.73735 1.1058 -2.989361.38901 -8.1846 -2.4315 -5.1189 1.8225 3.3798 1.7874 4.693 -3.9286 1.4266 5.7849 -3.4894 -4.0305 3.7879 3.5195 2.9186 2.8685

Re: [R] block diagrams

2009-07-03 Thread jim holtman
Check out the R-Graphviz package on Bio-Conductor. On Fri, Jul 3, 2009 at 5:43 AM, rajesh jakshay.raj...@gmail.com wrote: Hi, I need to draw something like the attached Image.Where can I do this in R? -- Rajesh.J __ R-help@r-project.org mailing

Re: [R] Histogram

2009-07-03 Thread Uwe Ligges
aledanda wrote: Hallo, How can I have the control of the width of each bin in the histogram? I would like to plot my data modifying the number and the width of the bins. Is there a good pdf or manual where i can learn to use the plot/hist in R very well? See ?hist and its argument breaks.

Re: [R] help system (no subject)

2009-07-03 Thread Uwe Ligges
Please provide a sensible subject line as you are asked in the posting guide which you obviously have not read before posting. Please also specify your OS, the R version and in this case which help system your are going to use. HTML? Plain text? Compiled html (if under Windows)? PDF? My

Re: [R] Histogram

2009-07-03 Thread Mark Knecht
On 7/3/09, aledanda danda.ga...@gmail.com wrote: Hallo, How can I have the control of the width of each bin in the histogram? I would like to plot my data modifying the number and the width of the bins. Is there a good pdf or manual where i can learn to use the plot/hist in R very

Re: [R] Color of ecdf plots

2009-07-03 Thread Uwe Ligges
Lars Bergemann wrote: Hi. I have the following two ecdf plots in one graph: plot( ecdf(), do.points=FALSE, verticals=TRUE, main=paste(Ecdf of distances ,DIM,sep=), col=red ); lines( ecdf(), do.points=FALSE, verticals=TRUE ); How do I change the color

Re: [R] normalised curve fitting with error bars

2009-07-03 Thread David Winsemius
On Jul 3, 2009, at 11:56 AM, Pooja Jain wrote: Dear List, My data consist of nine columns and about 50,000 rows. It looks like this. -9.0225 3.46464 2.80926 -0.3847 3.73735 1.1058 -2.98936 1.38901 -8.1846 -2.4315 -5.1189 1.8225 3.3798 1.7874 4.693 -3.9286 1.4266 5.7849 -3.4894

Re: [R] normalised curve fitting with error bars

2009-07-03 Thread Pooja Jain
Hi David, Sorry for not being absolutely clear. x-axis will be the values ranging from (-50 to 50) and y-axis will be the density distribution as you would get after plotting every column fitted to the normal distribution. However, what I need here is to reduce the nine column file to a

Re: [R] normalised curve fitting with error bars

2009-07-03 Thread John Kane
Reading in your data == xx - t(read.table(textConnection(x), header=FALSE, as.is=TRUE)); xx mm - apply(xx, 1, mean) mhi - apply(xx,1, max) mlow - apply(xx, 1, min) plot(mm, ylim=c(min(xx)-1, max(xx)+1)) points(mhi, col=red)

[R] Two questions about the cloud function in the lattice package

2009-07-03 Thread Martin Eklund
Hi, I have two questions regarding the cloud function in the lattice package: 1) Is there a way to not print the surrounding frame (i.e. the square surrounding the entire plot)? 2) Is there a way to italicize the text displayed with the key argument? Some sample code: data(iris)

Re: [R] normalised curve fitting with error bars

2009-07-03 Thread David Winsemius
On Jul 3, 2009, at 12:52 PM, Pooja Jain wrote: Hi David, Sorry for not being absolutely clear. x-axis will be the values ranging from (-50 to 50) and y-axis will be the density distribution as you would get after plotting every column fitted to the normal distribution. However, what I

Re: [R] Two questions about the cloud function in the lattice package

2009-07-03 Thread David Winsemius
On Jul 3, 2009, at 1:02 PM, Martin Eklund wrote: data(iris) cloud(Sepal.Length~Petal.Length*Petal.Width,data=iris, groups=Species,screen=list(z=20,x=-70), col=c(red,blue,green), pch=c(20,20,20), #try adding par.box=list(col=white), perspective=FALSE,

[R] What command lists everything in a package?

2009-07-03 Thread Mark Knecht
Hi, Two easy questions I'm sure. 1) As an example if I use the code require(zoo) then once it's loaded is there a command that lists everything that zoo provides so that I can study the package? Certainly help(zoo) gives me some clues about what zoo does but I'd like a list. Maybe

[R] windrose (circular package) odd table and windrose plot

2009-07-03 Thread Alessandro Bigi
R.version _ platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2

[R] A fast version of ccf () accepting missing values ?

2009-07-03 Thread bruno Piguet
Dear R-Helpers, I need to compute cross-correlation on two signals wich may contain missing values. One cannot pass Na.action=na.pass to the ccf() function. So, I wrote a naive function of my own (see below). Unsurprisingly, this function is not very fast. Do you think that it is

Re: [R] What command lists everything in a package?

2009-07-03 Thread Gabor Grothendieck
library(help = zoo) will list all the help files and, at the bottom, all the vignettes. It does not list demos so to get those you have to issue a second command: demo(package = zoo) On Fri, Jul 3, 2009 at 1:21 PM, Mark Knechtmarkkne...@gmail.com wrote: Hi,   Two easy questions I'm sure. 1)

Re: [R] Color of ecdf plots

2009-07-03 Thread David Winsemius
Does not work for me with rnorm(20) as an argument. (And supplying a dummy value for DIM.) The error after fixing the missing DIM issue was: Error in plot.stepfun(x, ..., ylab = ylab, verticals = verticals, pch = pch) : argument 4 matches multiple formal arguments The help page for

Re: [R] Two questions about the cloud function in the lattice package

2009-07-03 Thread Deepayan Sarkar
On Fri, Jul 3, 2009 at 10:02 AM, Martin Eklundmartin.ekl...@farmbio.uu.se wrote: Hi, I have two questions regarding the cloud function in the lattice package: 1) Is there a way to not print the surrounding frame (i.e. the square surrounding the entire plot)? There is an example in the

[R] exporting interactive rggobi plots

2009-07-03 Thread Veerappa Chetty
Hi , I have created a parallel coordinate plot using rrgobi. I like to make a powerpoint slide. I am not able to copy and paste. Is there any way to do this? Thanks for any help. Chetty -- Professor of Family Medicine Boston University Tel: 617-414-6221, Fax:617-414-3345 emails:

Re: [R] What command lists everything in a package?

2009-07-03 Thread Duncan Murdoch
On 7/3/2009 1:21 PM, Mark Knecht wrote: Hi, Two easy questions I'm sure. 1) As an example if I use the code require(zoo) then once it's loaded is there a command that lists everything that zoo provides so that I can study the package? ls(package:zoo) will list all the exported items in

Re: [R] What command lists everything in a package?

2009-07-03 Thread David Winsemius
On Jul 3, 2009, at 1:21 PM, Mark Knecht wrote: Hi, Two easy questions I'm sure. snipped as already answered 2) Related to the above, how do I tell what packages are currently loaded at any given time so that I don't waste time loading things that are already loaded? search() tells me

Re: [R] What command lists everything in a package?

2009-07-03 Thread Henrique Dallazuanna
To see all functions in a package: library(zoo) ls(package:zoo) On Fri, Jul 3, 2009 at 2:21 PM, Mark Knecht markkne...@gmail.com wrote: Hi, Two easy questions I'm sure. 1) As an example if I use the code require(zoo) then once it's loaded is there a command that lists everything that

Re: [R] What command lists everything in a package?

2009-07-03 Thread Mark Knecht
Thanks to everyone for your answers. They were great and I think I have what I was looking for. Cheers, Mark On 7/3/09, Duncan Murdoch murd...@stats.uwo.ca wrote: On 7/3/2009 1:21 PM, Mark Knecht wrote: Hi, Two easy questions I'm sure. 1) As an example if I use the code

Re: [R] exporting interactive rggobi plots

2009-07-03 Thread David Winsemius
On Jul 3, 2009, at 1:35 PM, Veerappa Chetty wrote: Hi , I have created a parallel coordinate plot using rrgobi. I like to make a powerpoint slide. I am not able to copy and paste. Is there any way to do this? Thanks for any help. Save it (specifics may depend on OS) and then in PPT (this

Re: [R] productivity tools in R?

2009-07-03 Thread cls59
Duncan Murdoch-2 wrote: The Windows GUI was written before the Mac GUI, and was written using a fairly strange toolkit, which are two reasons the Mac looks prettier. But the Windows GUI has a few advantages over the Mac: (snip) You can copy from the console, prompts and all, and

Re: [R] exporting interactive rggobi plots

2009-07-03 Thread baptiste auguie
I've never tried it but I would have thought the package DescribeDisplay can help you, http://www.ggobi.org/describe-display/ baptiste 2009/7/3 David Winsemius dwinsem...@comcast.net On Jul 3, 2009, at 1:35 PM, Veerappa Chetty wrote: Hi , I have created a parallel coordinate plot using

Re: [R] productivity tools in R?

2009-07-03 Thread Gabor Grothendieck
On Fri, Jul 3, 2009 at 2:14 PM, cls59ch...@sharpsteen.net wrote: Duncan Murdoch-2 wrote: The Windows GUI was written before the Mac GUI, and was written using a fairly strange toolkit, which are two reasons the Mac looks prettier. But the Windows GUI has a few advantages over the Mac:

Re: [R] productivity tools in R?

2009-07-03 Thread Duncan Murdoch
On 7/3/2009 2:14 PM, cls59 wrote: Duncan Murdoch-2 wrote: The Windows GUI was written before the Mac GUI, and was written using a fairly strange toolkit, which are two reasons the Mac looks prettier. But the Windows GUI has a few advantages over the Mac: (snip) You can copy from the

[R] Variable names in lattice XY-plot

2009-07-03 Thread ukoenig
Hi, how can I get a more descriptive text instead of the variable names in my XY-lattice plot, according to the table below? Variabletext acet = Acetylaspartate Thalamus chol = Choline Thalamus acetp= Acetylaspartate parieoc ino = Inositole Thalamus I could not find a

Re: [R] productivity tools in R?

2009-07-03 Thread Gabor Grothendieck
On Fri, Jul 3, 2009 at 3:07 PM, Duncan Murdochmurd...@stats.uwo.ca wrote: On 7/3/2009 2:14 PM, cls59 wrote: Duncan Murdoch-2 wrote: The Windows GUI was written before the Mac GUI, and was written using a fairly strange toolkit, which are two reasons the Mac looks prettier. But the Windows

Re: [R] Warning when trying to access a variable out of scope?

2009-07-03 Thread Greg Snow
Does this do what you want? b.test - 3 f - function(a,b) { + a+b.test + } f(10,20) [1] 13 environment(f) - baseenv() f(10,20) Error in f(10, 20) : object 'b.test' not found -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org

Re: [R] productivity tools in R?

2009-07-03 Thread cls59
Duncan Murdoch-2 wrote: cls59 wrote: I recently noticed that the Mac GUI contains a way to do block selection- so you can select just the code you want to copy and exclude the prompts. Just hold down the option key and the text cursor will turn into a crosshair which allows you to

Re: [R] productivity tools in R?

2009-07-03 Thread Max Kuhn
As for the editor debate... I'm not touching that one with a ten hundred foot pole :p You don't have to: in emacs, C-x M-p will do it for you. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] Variable names in lattice XY-plot

2009-07-03 Thread baptiste auguie
?strip.custom p - xyplot(acet+chol+ino+acetp ~ zp, group=grp, data=data, type=l, scales=list(relation=free), auto.key=list(title= Neurotransmitters, border=TRUE)) update(p, strip=strip.custom(factor.levels=letters[1:4])) HTH, baptiste 2009/7/3

[R] Multi-line plots - max Y?

2009-07-03 Thread Mark Knecht
Hi all, OK, I'm struggling a bit with getting good basic plots so I put together this little example. It works - sort of - but it seems a bit hokey to me. Can someone suggest how a pro would do it? The goal is to take data from 4 columns of a data.frame (data read in by read.csv) and

Re: [R] Variable names in lattice XY-plot

2009-07-03 Thread Deepayan Sarkar
On Fri, Jul 3, 2009 at 1:02 PM, baptiste auguiebaptiste.aug...@googlemail.com wrote: ?strip.custom p - xyplot(acet+chol+ino+acetp ~ zp, group=grp, data=data, type=l, scales=list(relation=free), auto.key=list(title= Neurotransmitters, border=TRUE)) update(p,

Re: [R] Multi-line plots - max Y?

2009-07-03 Thread Mark Knecht
On Fri, Jul 3, 2009 at 1:38 PM, Dylan Beaudettedylan.beaude...@gmail.com wrote: Hi, how about something like this: a - 1:10 b - cumsum(a) c - cumsum(b) d - cumsum(c) X- data.frame(a,b,c,d) plot(b ~ a, data=X, type=l, col=blue, ylim=c(0,max(X))) lines(c ~ a, data=X, col=green)

Re: [R] Multi-line plots - max Y?

2009-07-03 Thread Dylan Beaudette
Hi, how about something like this: a - 1:10 b - cumsum(a) c - cumsum(b) d - cumsum(c) X- data.frame(a,b,c,d) plot(b ~ a, data=X, type=l, col=blue, ylim=c(0,max(X))) lines(c ~ a, data=X, col=green) lines(d ~ a, data=X, col=red) legend('topleft', legend=c('a', 'b', 'c'), col=c('blue', 'green',

Re: [R] Variable names in lattice XY-plot

2009-07-03 Thread ukoenig
Thanks a lot Baptiste and Deepayan! That was very helpful-things can be so easy... -Udo Quoting Deepayan Sarkar deepayan.sar...@gmail.com: On Fri, Jul 3, 2009 at 1:02 PM, baptiste auguiebaptiste.aug...@googlemail.com wrote: ?strip.custom p - xyplot(acet+chol+ino+acetp ~ zp,

[R] The time series analysis functions/packages don't seem to like my data

2009-07-03 Thread Ted Byers
I have hundreds of megabytes of price data time series, and perl scripts that extract it to tab delimited files (I have C++ programs that must analyse this data too, so I get Perl to extract it rather than have multiple connections to the DB). I can read the data into an R object without any

Re: [R] Plot Fonts in Windows vs Mac OSX

2009-07-03 Thread Sean Carmody
Thanks for the reponse Winston. Your chart (http://stdout.org/~winston/X/r-antialias/pch.html) illustrates my problem perfectly. I did experiment with Cairo, with no success in improving the display. In the blog post you linked to, the following comment appeared: If you don't have anti-aliasing

Re: [R] Correlation Network Diagram?

2009-07-03 Thread Mark Difford
Hi Rory, There are several. Have a look at the gR Task Views. There you will also find a link to the statnet suite, where you will find links to a dedicated set of jstatsoft articles. Regards, Mark. Rory Winston wrote: Hi all On page 39 of this paper [1] by Andrew Lo there is a very

Re: [R] The time series analysis functions/packages don't seem to like my data

2009-07-03 Thread Mark Knecht
On Fri, Jul 3, 2009 at 2:48 PM, Ted Byersr.ted.by...@gmail.com wrote: I have hundreds of megabytes of price data time series, and perl scripts that extract it to tab delimited files (I have C++ programs that must analyse this data too, so I get Perl to extract it rather than have multiple

Re: [R] is AIC always 100% in evaluating a model?

2009-07-03 Thread Ben Bolker
alexander russell-2 wrote: Hello, I'd like to say that it's clear when an independent variable can be ruled out generally speaking; on the other hand in R's AIC with bbmle, if one finds a better AIC value for a model without the given independent variable, versus the same model with,

[R] Error using the Rdonlp2‏ Package

2009-07-03 Thread Lars Bishop
Dear experts, I'm attempting to solve a constrained optimization problem using the Rdonlp2 package. I created a Lagrange function (L=f(x)-lambda(g(x)-c)), where x is a vector of 16 parameters. This is what I'm using as objective function in the code below. In addition, I set bounds on these

Re: [R] Simple one

2009-07-03 Thread Ben Bolker
JoK LoQ wrote: Just a quickly beginner's question. I wanna find the mean only from the values from a column related to specific values from another one. Like, theres a 'region' column, i want the mean of the value on 'profit' column only from south sells from 'region' column

Re: [R] The time series analysis functions/packages don't seem to like my data

2009-07-03 Thread Ted Byers
Hi Mark Thanks for replying. Here is a short snippet that reproduces the problem: library(PerformanceAnalytics) thedata = read.csv(K:\\Work\\SignalTest\\BP.csv, sep = \t, header = FALSE, na.strings=) thedata x = as.timeseries(thedata) x table.Drawdowns(thedata,top = 10)

  1   2   >