Re: [R] Problem with retreaving numeric data

2010-02-11 Thread Dieter Menne
dvkirankumar wrote: I am trying to get a value from a dataset in data set its stored as 001 but while i am trying to retreave it its giving as 1 can any one help me to retreave it as 001 I have to guess if you used a boomerang or a sling to retrieve your data. My blind guess is you

Re: [R] Running rscript in windows

2010-02-11 Thread Peter Ehlers
How did you change directory? Does dir() show that your directory is the desktop? Presumably the desktop folder is something like C:/Users/someone/Desktop. Sourcing test.R from there should work. But the simplest way is source(file.choose()) and then navigate to where test.R is located.

Re: [R] Proxy settings

2010-02-11 Thread Van Wyk, Jaap
Thanks, especially to Greg for the help with the proxy settings required for Windows. I have managed to solve the problem, with Greg's latest suggestion. Jacob Jacob L van Wyk Department of Statistics University of Johannesburg (APK) PO Box 524, Auckland Park, 2006 Office ph: 011 559 3080 Fax:

Re: [R] question about boxplot,does it accept ad hoc data

2010-02-11 Thread Dieter Menne
Guanhua Chen wrote: I have a question about boxplot(), does it accept ad hoc data , for example, I only provide quartile information to the function,but get the boxplot as usual, is it possible , or if there is any other function can do this. Thank you very much! boxplot calls bxp

Re: [R] Problem with retreaving numeric data

2010-02-11 Thread Daniel Malter
My guess is that you cannot keep it numeric, but that you have to convert it to a string (the gurus may know better). x=001 x x=001 x Daniel - cuncta stricte discussurus - -Original Message- From: r-help-boun...@r-project.org

Re: [R] Question about rank() function

2010-02-11 Thread Daniel Malter
To follow up on Moshe's post, do rank(fmodel) rank(round(fmodel,4)) rank(fmodel)==rank(round(fmodel,4)) If the two are not identical, you have the explanation Moshe suggested, just that this approach is somewhat more comprehensive as it works on the entire fmodel vector. See the example below:

[R] SAS and RODBC

2010-02-11 Thread Daniel Nordlund
I am using R-2.10.1 binary from CRAN on a WinXP Pro system. I also use SAS v9.2 on the same box. I just started using the SAS ODBC driver that comes with version 9 of SAS. I have been able to set up an ODBC source for SAS datasets using the driver, and then with RODBC I am able to read a

[R] ff package: How to save and open ff(df) files.

2010-02-11 Thread daniel castro
Hello to everyone, I'm a newbie with ff package and I´m starting to use it. I´ve been reading the ff.pdf guide and another documents and questions , but I´m really confused about some procedures I can´t see how to do. I´d want to know if it´s possible (and how) to save a ffdf file(s) and open it

Re: [R] Running rscript in windows

2010-02-11 Thread vikrant
Thanks a lot, really. dir() option showd my desktop files and test.R was shown as test.R.txt when I run source(test.R.txt) it worked. Just for your information, I changed the directory using Menu options File- change dir. I verified the same using getwd() command in R. ' -- View this message

Re: [R] SAS and RODBC

2010-02-11 Thread Prof Brian Ripley
On Thu, 11 Feb 2010, Daniel Nordlund wrote: I am using R-2.10.1 binary from CRAN on a WinXP Pro system. I also use SAS v9.2 on the same box. I just started using the SAS ODBC driver that comes with version 9 of SAS. I have been able to set up an ODBC source for SAS datasets using the

Re: [R] Regarding a error while plotting R chart using qcc package.

2010-02-11 Thread Luca Scrucca
On 11 Feb 2010, at 06:53, Vikrant Shimpi wrote: Dear Luka , I am using qcc package in R to plot SPC charts. BUt while plotting R chart I had a error. My question is it necessary while plotting R Chart the group sample size must be 25 ?. Because when I took group sample size as 1000 it

Re: [R] SAS and RODBC

2010-02-11 Thread Gabor Grothendieck
There is a SASxport package on CRAN that may be of interest. On Thu, Feb 11, 2010 at 3:38 AM, Daniel Nordlund djnordl...@verizon.net wrote: I am using R-2.10.1 binary from CRAN on a WinXP Pro system.  I also use SAS v9.2 on the same box.  I just started using the SAS ODBC driver that comes

Re: [R] R2HTML

2010-02-11 Thread Regina Schinner
I am using R 2.10.1 for Mac OS. But Bart's Tip to use 'HTML(summary(iris))' completely solved the problem. Thanks to everybody for helping. Regina Am 10.02.2010 um 20:01 schrieb Greg Snow: When I copy/paste/run your code below I get a file with the summary output in a nice table. Telling

[R] ZOO object colnames refering to Dates

2010-02-11 Thread Research
, version of virus signature database 4857 (20100211) __ The message was checked by ESET Smart Security. http://www.eset.com __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

[R] Sweve/cacheSweave

2010-02-11 Thread Unger, Kristian
Hi there I have a problem with using Sweave in combination with the option driver = cacheSweave. Whichever code I try to run - when it comes to converting the tex file into pdf it comes up with the same errors (\csname \endcsname errors). Does anybody have an idea what it going wrong?

[R] Biclustering package

2010-02-11 Thread Alex Roy
Hello, I am looking for R package which can perform biclustering a part from biclust package. thanks Alex [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] Sweave, lty = 3 line incorect in pdf output

2010-02-11 Thread Achim Zeileis
Ken, I can't contribute much, just that I tried to look into it and can't replicate it. I tried it both with .Stex and .Rnw yielding identical results, just like plotting it interactively (on my default X11) and printing it via dev.copy2pdf(). All of the approaches seemed to yield the same

Re: [R] Biclustering package

2010-02-11 Thread Gábor Csárdi
Alex, the isa2 package implements the biclustering algorithm discussed in Bergmann S, Ihmels J, and Barkai N. Iterative signature algorithm for the analysis of large-scale gene expression data. Phys Rev E Stat Nonlin Soft Matter Phys 2003 Mar; 67(3 Pt 1) 031902 Best, Gabor On Thu, Feb 11, 2010

Re: [R] Unexpected output in first iteration of for-loop

2010-02-11 Thread Chaehan So
Ok, you're right - may I rephrase: How should I modify the assignment of result - latentVariableNames so it produces the output without the first line? I thought result - NULL should do the job, but it didn't because of the following names(result) assignment (which I understand, but not how I can

Re: [R] Sweave, lty = 3 line incorect in pdf output

2010-02-11 Thread Ken Knoblauch
Hi Achim, Thanks for trying. Maybe, it's a Mac-thing. I'll wait to see if I get other replies. It's strange, as the quartz graphic is correct but the pdf produced by Sweave doesn't contain the dotted line, with the caveat, on my Mac for the moment. Ken Quoting Achim Zeileis

Re: [R] ZOO object colnames refering to Dates

2010-02-11 Thread Achim Zeileis
really appreciated. Thanks in advance, Costas __ Information from ESET Smart Security, version of virus signature database 4857 (20100211) __ The message was checked by ESET Smart Security. http://www.eset.com __ R-help@r-project.org

[R] gdata

2010-02-11 Thread Casimir
Hi Using R 2.10.1 on a mac os 10.6.2, I have have a problem with gdata package. When I use the command read.xls, I get this error-message: Erreur dans xls2sep(xls, sheet, verbose = verbose, ..., method = method, : Unable to read translated csv file

Re: [R] ZOO object colnames refering to Dates

2010-02-11 Thread Gabor Grothendieck
4857 (20100211) __ The message was checked by ESET Smart Security. http://www.eset.com __ 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

Re: [R] Sweave, lty = 3 line incorect in pdf output

2010-02-11 Thread Ken Knoblauch
Some follow-up after Achim's message, If I open the Sweave created pdf, lty3Test-001.pdf, attached to the first mail, with the display command from imageMagick, indeed, the dotted line is present, but faintly, as with the hardcopy printout, i.e., not respecting the lwd = 3 argument. But, the

Re: [R] gdata

2010-02-11 Thread Gabor Grothendieck
The most recent version of gdata uses certain perl packages not used by the older versions and perhaps that is causing the problem. You could try installing the latest version of active perl or if that fails you could try an older version of gdata: http://www.activestate.com/activeperl/

Re: [R] AFTREG with ID argument

2010-02-11 Thread Philipp Rappold
Göran, thanks! One more thing that I found: As soon as you have at least one NA in the independent vars, the trick that you mentioned does not work anymore. Example: testdata start stop censor groupvar var1 1 01 01 0.1284928 2 12 01 0.4896125

Re: [R] R wiki link ?

2010-02-11 Thread Martin Maechler
PhGr == Philippe Grosjean phgrosj...@sciviews.org on Mon, 08 Feb 2010 19:55:08 +0100 writes: PhGr Prof. John C Nash wrote: Is this a transient problem, or has the link to the R wiki on the R home page (www.r-project.org) to http://wiki.r-project.org/ been corrupted? I

Re: [R] Sweve/cacheSweave

2010-02-11 Thread Duncan Murdoch
Unger, Kristian wrote: Hi there I have a problem with using Sweave in combination with the option driver = cacheSweave. Whichever code I try to run - when it comes to converting the tex file into pdf it comes up with the same errors (\csname \endcsname errors). Does anybody have an idea

Re: [R] Running rscript in windows

2010-02-11 Thread Duncan Murdoch
vikrant wrote: Thanks a lot, really. dir() option showd my desktop files and test.R was shown as test.R.txt when I run source(test.R.txt) it worked. Just for your information, I changed the directory using Menu options File- change dir. I verified the same using getwd() command in R. '

Re: [R] Sweve/cacheSweave

2010-02-11 Thread Duncan Murdoch
On 11/02/2010 6:06 AM, Duncan Murdoch wrote: Unger, Kristian wrote: Hi there I have a problem with using Sweave in combination with the option driver = cacheSweave. Whichever code I try to run - when it comes to converting the tex file into pdf it comes up with the same errors (\csname

[R] ACF and PACF

2010-02-11 Thread sathiya_mtm
Hi helpers, can you help me in plotting acf and pacf functions in R. I am using the code acf(variable name) but it is not working. Expecting your reply. Thanks -- View this message in context: http://n4.nabble.com/ACF-and-PACF-tp1477149p1477149.html Sent from the R help mailing list

[R] (linux) display messages in R startup screen

2010-02-11 Thread Wolfgang Raffelsberger
Dear guRus, I'm maintaining central R installation (Linux) for our institute. The startup screen of R traditionally echoes some messages (starting with R version ...Copyright (C).. The R Foundation for Statistical Computing...) and I'd like to display for example right after these messages

[R] Coda or R or input error? (Error in lowess(xp, yp[, k]) :...)

2010-02-11 Thread Raisinmeister
I have been using R with the packages locfit and coda to analyse output from a population genetics program (msvar) which I have been running on an external server (http://cbsuapps.tc.cornell.edu/ index.aspx). One of the first steps is to read the msvar output tables into R and then use coda to

[R] Zero-inflated Negat. Binom. model

2010-02-11 Thread Luciano La Sala
Dear R crew: I am sorry this question has been posted before, but I can't seem to solve this problem yet. I have a simple dataset consisting of two variables: cestode intensity and chick size (defined as CAPI). Intensity is a count and clearly overdispersed, with way too many zeroes. I'm

Re: [R] (linux) display messages in R startup screen

2010-02-11 Thread Dirk Eddelbuettel
On 11 February 2010 at 12:08, Wolfgang Raffelsberger wrote: | I'm maintaining central R installation (Linux) for our institute. The | startup screen of R traditionally echoes some messages (starting with R | version ...Copyright (C).. The R Foundation for Statistical | Computing...) and I'd

[R] cluster/distance large matrix

2010-02-11 Thread Bart Thijs
Hi all, I've stumbled upon some memory limitations for the analysis that I want to run. I've a matrix of distances between 38000 objects. These distances were calculated outside of R. I want to cluster these objects. For smaller sets (egn=100) this is how I proceed: A-matrix(scan(file,

[R] [R-pkgs] GWRM package

2010-02-11 Thread Antonio José Sáez Castillo
I'd like to announce the availability of the GWRM package. It provides a set of functions for fitting Generalized Waring Regression Models. It also includes the dataset and the example of Rodriguez-Avi, J; Conde-Sanchez, A; Saez-Castillo, A.J., Olmo-Jimenez, M. J. and Martinez Rodriguez, A.

Re: [R] Sweve/cacheSweave

2010-02-11 Thread Unger, Kristian
Thank you very much Duncan! Adding the line \usepackage{Sweave} fixed the problem. Best wishes Kristian _ Dr Kristian Unger Bioinformaticist Chernobyl Tissue Bank Division of Surgery and Cancer Imperial College London Rm G02, Ground floor G Block Hammersmith Hospital Du Cane Road,

Re: [R] cluster/distance large matrix

2010-02-11 Thread Christian Hennig
Dear Bart, a strange thing in your question is that the term Ward's method usually refers to a method based on the k-means criterion, which, in its standard form, is not based on dissimilarities, but on objects*variables-data. So I wonder how and why you want to use Ward's method on a

Re: [R] ACF and PACF

2010-02-11 Thread Uwe Ligges
On 11.02.2010 12:48, sathiya_mtm wrote: Hi helpers, can you help me in plotting acf and pacf functions in R. I am using the code acf(variable name) but it is not working. For me it is, hence citing the footer of al R-help messages: PLEASE do read the posting guide

[R] WMF conversion

2010-02-11 Thread Terry Therneau
Could you save the original graphics as a pdf, and include the pdf in Powerpoint? -- begin included message -- Unfortunately, after placed in the PowerPoint and the PowerPoint is converted to PDF via MS Office's built in conversion utility, the resulting image have diagonal

Re: [R] Total least squares linear regression

2010-02-11 Thread Terry Therneau
This is also called Deming regression and perhaps many other things. It is a lively topic in the validation of competing assay methods in the laboratory. I have a function 'deming.R' that does a generalized form of this, based on the (very nice) article below, the code is attached. (The

[R] Odp: Simplify the code

2010-02-11 Thread Petr PIKAL
Hi I did not see the answer yet, probably nobody is much interested in such function, so I will try. r-help-boun...@r-project.org napsal dne 10.02.2010 21:57:42: Hi all, I wrote the following code for the function that is attached. I want to simply the code. Why. If it works there is

[R] Question re model development

2010-02-11 Thread Stratos Laskarides
Dear Sir/Madam I'm trying quite a long shot here... I have used R previously, and would just like to know if any software/model has been designed which will allow me to project the future risk profile of a group of insured lives (by factors such as age, gender etc.), based on past and exisiting

[R] for loop the can be skipped

2010-02-11 Thread Fabrizio Pollastri
Hello, there is compact way to skip a for loop when the requested number of iterations is 0 ? For example, something equivalent to the following code: for (i in some_function(from=1,to=iteration_cycles)) { ... } Where the loop is skipped when iteration_cycles==0 Thanks in advance. F:

Re: [R] for loop the can be skipped

2010-02-11 Thread Duncan Murdoch
On 11/02/2010 9:28 AM, Fabrizio Pollastri wrote: Hello, there is compact way to skip a for loop when the requested number of iterations is 0 ? For example, something equivalent to the following code: for (i in some_function(from=1,to=iteration_cycles)) { ... } Where the loop is skipped

Re: [R] Biclustering package

2010-02-11 Thread Alex Roy
Thank you very much Gabor. Alex On Thu, Feb 11, 2010 at 10:55 AM, Gábor Csárdi csa...@rmki.kfki.hu wrote: Alex, the isa2 package implements the biclustering algorithm discussed in Bergmann S, Ihmels J, and Barkai N. Iterative signature algorithm for the analysis of large-scale gene

Re: [R] cluster/distance large matrix (fwd)

2010-02-11 Thread Thomas Lumley
On Thu, 11 Feb 2010, Christian Hennig wrote: It is well know that hierarchical methods are problematic with too large dissimilarity matrices; even if you resolve the memory problem, the number of operations required is enormous. There is at least one exception to this. Single-linkage

Re: [R] Suprising behavior of paste or cat?

2010-02-11 Thread Russell Pierce
If the list would provide additional reports I'd appreciate it. I replicated the effect this morning on another Windows XP machine... R version 2.10.1 (2009-12-14) i386-pc-mingw32 locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3]

[R] R ANOVA gives diferent results than SPSS

2010-02-11 Thread Protzko
I guess my subject says it all. But I loaded a dataset in spss and used the foreign package to read and save it in R. Running an anova (using the aov command) gives a different F and p value in R than it does in SPSS. ANy idea what is going on? -- View this message in context:

[R] Using Rscript in combination with eval() results in error

2010-02-11 Thread Janko Thyson
Dear List, I'm having trouble running a .bat file which calls an R-Script via the command line (using Rscript). I put the following line in a file called test.bat: Rscript --vanilla test.R Then I tried to launch test.bat via Windows' CMD (I plan to make this a scheduled Windows task). The

[R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Johan Jackson
Hi all, First off, it is surprising that there are no examples of how to use read.table() under ?read.table ! I am trying to read in a flat file of type 'raw'. It has 1000 rows and 600K columns. I have the RAM to accomplish this, but can't get the data into R using read.table: x -

Re: [R] R ANOVA gives diferent results than SPSS

2010-02-11 Thread Thomas Lumley
On Thu, 11 Feb 2010, Protzko wrote: I guess my subject says it all. But I loaded a dataset in spss and used the foreign package to read and save it in R. Running an anova (using the aov command) gives a different F and p value in R than it does in SPSS. ANy idea what is going on? Yes.

Re: [R] Unexpected output in first iteration of for-loop

2010-02-11 Thread Don MacQueen
You have a mistake in how you're setting up the object named result before the loop. You set result - latentVariableNames. It is a vector of length 6, when you call the function with LV. The printout from the first iteration shows this. But then you rbind result with a data frame that has

Re: [R] histogam plots

2010-02-11 Thread li li
Thank you very much for your reply. That was very helpful. I also want to add two density curves on top of each histgram. One is the density curve of a standard normal random variable. The other is the denstity curve according to the histgram. I was trying to use the function panel.mathdensity,

Re: [R] Creating a new Access database with R

2010-02-11 Thread Ryusuke Kenji
I am facing the same problem as well, I would like to code with following concept but wondering how to cope it. if (*mdb file exist) { add new row/col } else { add new *mdb file } -- From:

Re: [R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Ivan Calandra
Hi! |colClasses| character. A vector of classes to be assumed for the columns. I'm not an R expert and I don't know what your flat file raw is, but the colClasses argument is to define whether the column will be treated as containing factors, logical, integer etc... For more on

Re: [R] Simplify the code

2010-02-11 Thread li li
Thank you for your reply. I was just trying to make the code less complicated. I got your point. Thanks! 2010/2/11 Petr PIKAL petr.pi...@precheza.cz Hi I did not see the answer yet, probably nobody is much interested in such function, so I will try. r-help-boun...@r-project.org napsal dne

Re: [R] Zero-inflated Negat. Binom. model

2010-02-11 Thread Achim Zeileis
On Thu, 11 Feb 2010, Luciano La Sala wrote: Dear R crew: I am sorry this question has been posted before, but I can't seem to solve this problem yet. Just for the others reader who might not recall that you asked virtually the same question last week. This was my reply:

Re: [R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Ivan Calandra
Well, it's too complicated for me! Here are what I would do (limited since I'm still a newbie) 1) the syntax seems correct, it should work. The problem is somewhere else, coming from your own file. Did you try skipping the colClasses argument? To see how it looks like... If you can import it

Re: [R] AFTREG with ID argument

2010-02-11 Thread David Winsemius
On Feb 11, 2010, at 5:58 AM, Philipp Rappold wrote: Göran, thanks! One more thing that I found: As soon as you have at least one NA in the independent vars, the trick that you mentioned does not work anymore. Example: testdata start stop censor groupvar var1 1 01 0

[R] Access variables by string

2010-02-11 Thread Philipp Rappold
Dear all, I have two probably very easy questions: (1) Is there a way to access certain variables by their string-based name representation? Example: numbers - c(one, two, three) varname - numbers print(varname[2]) (2) I need this functionality for a customized na.exclude() function that I

Re: [R] Suprising behavior of paste or cat?

2010-02-11 Thread David Winsemius
On Feb 11, 2010, at 10:17 AM, Russell Pierce wrote: If the list would provide additional reports I'd appreciate it. I replicated the effect this morning on another Windows XP machine... Worked as expected on a Mac OSX 10.5.8, R 2.10.1, [R.app GUI 1.30 (5534) x86_64-apple-darwin9.8.0] --

Re: [R] Question re model development

2010-02-11 Thread David Winsemius
On Feb 11, 2010, at 9:10 AM, Stratos Laskarides wrote: Dear Sir/Madam I'm trying quite a long shot here... I have used R previously, and would just like to know if any software/model has been designed which will allow me to project the future risk profile of a group of insured lives (by

Re: [R] Integral of function of dnorm

2010-02-11 Thread Trafim Vanishek
This is exactly what I mean. I need to find integrate(dnorm(mean=8,sd=1)*log(dnorm(mean=8,sd=1)), -Inf, Inf) Which doesn't work like that, because it says: Error in dnorm(mean = 8, sd = 1) : element 1 is empty; the part of the args list of '.Internal' being evaluated was: (x, mean, sd,

Re: [R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Don MacQueen
The error message says there is no method for converting from 'character' to 'raw'. Apparently, R is seeing character data in the file, and is trying to convert it to raw, since you specified raw, and it can't. See, for example, as('aa','raw') Error in as(aa, raw) : no method or default

Re: [R] Access variables by string

2010-02-11 Thread Peter Dalgaard
Philipp Rappold wrote: Dear all, I have two probably very easy questions: (1) Is there a way to access certain variables by their string-based name representation? Example: numbers - c(one, two, three) varname - numbers print(varname[2]) (2) I need this functionality for a

Re: [R] Access variables by string

2010-02-11 Thread Don MacQueen
For your first question, use the get() function -Don At 5:18 PM +0100 2/11/10, Philipp Rappold wrote: Dear all, I have two probably very easy questions: (1) Is there a way to access certain variables by their string-based name representation? Example: numbers - c(one, two, three) varname

Re: [R] R ANOVA gives diferent results than SPSS

2010-02-11 Thread Jeremy Miles
I've always found exactly the same results. If you post code that allows us to reproduce this, I suspect someone would be able to shed light on it. And output too. J On 11 February 2010 06:47, Protzko prot...@gmail.com wrote: I guess my subject says it all. But I loaded a dataset in

Re: [R] Suprising behavior of paste or cat?

2010-02-11 Thread William Dunlap
Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Russell Pierce Sent: Wednesday, February 10, 2010 9:21 PM To: r-help@r-project.org Subject: [R] Suprising behavior of

Re: [R] R ANOVA gives diferent results than SPSS

2010-02-11 Thread Peter Dalgaard
Thomas Lumley wrote: On Thu, 11 Feb 2010, Protzko wrote: I guess my subject says it all. But I loaded a dataset in spss and used the foreign package to read and save it in R. Running an anova (using the aov command) gives a different F and p value in R than it does in SPSS. ANy idea

Re: [R] Integral of function of dnorm

2010-02-11 Thread Charles Annis, P.E.
Here's a suggestion: Plot the function: x - seq(3, 13, length=101) plot(x, y=dnorm(x, mean=8,sd=1)*log(dnorm(x,mean=8,sd=1))) x - seq(7.5, 8.5, length=101) plot(x, y=dnorm(x, mean=8,sd=1)*log(dnorm(x,mean=8,sd=1))) x - seq(7.9, 8.1, length=101) plot(x, y=dnorm(x,

Re: [R] Creating a new Access database with R

2010-02-11 Thread Paul-
As a workaround, you can keep an empty mdb file on your filesystem. When you need a new database, you can copy and rename the empty file. -- View this message in context: http://n4.nabble.com/Re-Creating-a-new-Access-database-with-R-tp1477400p1477533.html Sent from the R help mailing list

[R] aggregate function / custom column names?

2010-02-11 Thread Chuck White
This question is about column names returned by the aggregate function. Consider the following example df - data.frame( id = c(rep('11',30),rep('22',30),rep('33',30)), value = c(rnorm(30,2,0.5), rnorm(30,3,0.5), rnorm(30,6,0.5)) ) aggregate(df[,c(value),drop=FALSE], by=list(id=df$id),

Re: [R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Johan Jackson
Hi Don and all, I guess we're getting somewhere. Thanks. The file (first three columns, first five rows) looks like this: X10 X20 X30 00 0001 00 02 02 00 00 00 00 01 01 00 00 00 I guess R is reading 00 as a character? But here's the weird thing: this data (a raw matrix in R)

Re: [R] aggregate function / custom column names?

2010-02-11 Thread Gabor Grothendieck
Try this: aggregate(list(Max = df$value), df['id'], max) id Max 1 11 2.610491 2 22 3.796836 3 33 6.562515 or if using value rather than Max is ok then just: aggregate(df['value'], df['id'], max) idvalue 1 11 2.610491 2 22 3.796836 3 33 6.562515 On Thu, Feb 11, 2010 at 12:18 PM,

Re: [R] Access variables by string

2010-02-11 Thread Benilton Carvalho
On Thu, Feb 11, 2010 at 4:18 PM, Philipp Rappold philipp.rapp...@gmail.com wrote: Dear all, I have two probably very easy questions: (1) Is there a way to access certain variables by their string-based name representation? Example: numbers - c(one, two, three) varname - numbers

[R] Fwd: Recall function: evaluation nested too deeply: infinite recursion / options(expressions=)?

2010-02-11 Thread Allan Clark
hello all i searched the archives and couldn't get a solution to the following question. i have the following function: F=function(z,v) { if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1)) } and while testing whether it works (ie F(z,v) is approx =

Re: [R] Integral of function of dnorm

2010-02-11 Thread Peter Ehlers
Here's another suggestion: evaluate -0.5 * (log(2 * pi) + 1) ## -1.418939 Or you can do the integral: f - function(x, mu = 8, s = 1){ dnorm(x, mu, s) * log(dnorm(x, mu, s))} integrate(f, lower = -2, upper = 18) #-1.418939 with absolute error 2.8e-10 Try other values for lower=

Re: [R] Integral of function of dnorm

2010-02-11 Thread Greg Snow
Try: tmpfun - function(x) dnorm(x,mean=8,sd=1)*log(dnorm(x,mean=8,sd=1)) integrate( tmpfun, -Inf, Inf) Also you may want to look at the log argument to dnorm rather than taking the log of the function. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain

Re: [R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Johan Jackson
Hi Ivan, Thanks for the reply. Damn IT! My original post was screwed up. HERE is what I did: x - read.table(data,header=TRUE,colClasses=rep('raw',60)) #returns error: no method or default for coercing character to raw I've read the ?read.table and the colClasses argument. I'm still

Re: [R] Fwd: Recall function: evaluation nested too deeply: infinite recursion / options(expressions=)?

2010-02-11 Thread jim holtman
You might be impacted by FAQ 7.31. Your statement if(v == -.5) might be suspect. You might want to use 'all.equal' On Thu, Feb 11, 2010 at 11:40 AM, Allan Clark allan.cl...@uct.ac.za wrote: hello all i searched the archives and couldn't get a solution to the following question. i

Re: [R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Greg Snow
The read.table function does not know how to convert the character representation that it reads into raw variables. Try using 'integer' for the colClasses to read the data in as integers, then convert those back to raw (if that is really what you need). -- Gregory (Greg) L. Snow Ph.D.

Re: [R] histogam plots

2010-02-11 Thread Dennis Murphy
Hi: ## Histogram + density plots in lattice and ggplot2 # example data frame x - data.frame(x = rnorm(600), g = factor(rep(1:6, each = 100))) # lattice library(lattice) histogram(~ x | g, data = dd, type = 'density', panel = function(x, subscripts, ...) {

[R] Recall function: evaluation nested too deeply: infinite recursion / options(expressions=)?

2010-02-11 Thread Allan Clark
hello all i searched the archives and couldn't get a solution to the following question. i have the following function: F=function(z,v) { if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1)) } and while testing whether it works (ie F(z,v) is approx =

Re: [R] R ANOVA gives diferent results than SPSS

2010-02-11 Thread Protzko
a one-way ANOVA should be a one-way ANOVA I guess, model is simple enough I thought. The F value seems pretty clear, I'm doing nothing fancy here, just trying to figure out how to do in R what I'm doing in SPSS. -- View this message in context:

Re: [R] Access variables by string

2010-02-11 Thread Dan davison
Philipp Rappold philipp.rappold at gmail.com writes: Dear all, [...] (2) I need this functionality for a customized na.exclude() function that I am building, which should only exclude rows that have NA in certain columns. Maybe there is already a function which does exactly what I

[R] Dual Category X-Axis (Multi-Level)

2010-02-11 Thread Powell, Eric
Hello, I was wondering if anyone knows of a way to produce multilevel x-axis groupings similar to Excel's pivotcharts. (If you're not sure what I'm talking about, here's an example: http://peltiertech.com/WordPress/chart-with-a-dual-category-axis/). I realize that I could use mtext to

Re: [R] Creating a new Access database with R

2010-02-11 Thread Dieter Menne
Paul- wrote: As a workaround, you can keep an empty mdb file on your filesystem. When you need a new database, you can copy and rename the empty file. Creating a new database is not part of (R)ODBC because there are too many differences between implementations. You you use some RDCOM

Re: [R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Greg Snow
The other possibility is that you could create the function to convert from character to raw (possibly wrapping as.raw around as.integer) so that read.table knows what to do. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111

Re: [R] histogam plots

2010-02-11 Thread Dieter Menne
Dennis Murphy wrote: ## Histogram + density plots in lattice and ggplot2 There was a typo in your example (should be dd - data.frame), but anyway. ggplot2 needs 10 times as long. While Bill Venables (?) would argue what are you going to do with 3 seconds in your life, it adds up in

[R] Reading raw intensity data from getGEO's gse object

2010-02-11 Thread Dipen
Hello, I am trying to download a SOFT matrix file for GSE series to extract raw intensity data. gse- getGEO('GSE', AnnotGPL=T) However, I am not able to use the function Table(gse), which returns the following error Error in function (classes, fdef, mtable) : unable to find an

Re: [R] trouble with read.table and colClasses='raw'

2010-02-11 Thread Peter Ehlers
Johan Jackson wrote: Hi Don and all, I guess we're getting somewhere. Thanks. The file (first three columns, first five rows) looks like this: X10 X20 X30 00 0001 00 02 02 00 00 00 00 01 01 00 00 00 I guess R is reading 00 as a character? But here's the weird thing: this

Re: [R] Suprising behavior of paste or cat?

2010-02-11 Thread Russell Pierce
Thank you for your input so far r-help denizens. Neither David nor Peter were able to replicate my result. Has anybody other than me been able to generate the failure I'm describing? So far I've experienced it on 3 machines, Windows XP/P4/2.1.10, Windows XP/Atom/2.1.10/2.1.11(release), Windows

Re: [R] R ANOVA gives diferent results than SPSS

2010-02-11 Thread Greg Snow
A couple of possibilities: The data is not the same, e.g. something in the file was interpreted differently by the 2 programs, one of the programs may have stopped reading at an unrecognized value, while the other skipped it and went on. Or it used to be common to encode missing values as

Re: [R] Find each time a value changes

2010-02-11 Thread Tim Clark
It was brought to my attention that the rle() answer to this question was not posted. The following gives the correct answer once the last value is deleted. x-seq(1:100) y1-rep(1,10) y2-rep(2,10) y-c(y1,y2,y1,y1,y1,y2,y1,y2,y1,y2) xy-cbind(x,y) print(xy) print(str(xy)) # SEE WHAT RLE GIVES

Re: [R] Suprising behavior of paste or cat?

2010-02-11 Thread Duncan Murdoch
I don't think you have said how you are examining the output files. Is it possible that your text editor is assuming that the files are UCS-2 (Unicode), even though R is writing ASCII? Duncan Murdoch On 11/02/2010 1:44 PM, Russell Pierce wrote: Thank you for your input so far r-help

Re: [R] Creating a new Access database with R

2010-02-11 Thread Barry Rowlingson
On Thu, Feb 11, 2010 at 6:23 PM, Dieter Menne dieter.me...@menne-biomed.de wrote: Paul- wrote: As a workaround, you can keep an empty mdb file on your filesystem. When you need a new database, you can copy and rename the empty file. Creating a new database is not part of (R)ODBC because

Re: [R] Dual Category X-Axis (Multi-Level)

2010-02-11 Thread Greg Snow
Using mtext is probably one of the most straight forward ways. You can also use something like: par(mfrow=c(1,3), mar=c(5.1,0,4.1,0), oma=c(0,4.1,0,1.1)) one - 1:3 two - 2:5 three - 6:4 names(one) - LETTERS[1:3] names(two) - letters[2:5] names(three) - state.abb[6:4] tmp.ylim -

Re: [R] Find each time a value changes

2010-02-11 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Tim Clark Sent: Thursday, February 11, 2010 10:53 AM To: r-help@r-project.org Subject: Re: [R] Find each time a value changes It was brought to my attention that the rle()

  1   2   >