Re: [R] OdfWeave and Locale

2010-09-01 Thread Prof Brian Ripley
There have been encoding-related problems in odfWeave for years, and notice that it fails its checks on Windows (see the CRAN package check page). However, the problem here appears to be directional quotes, and you may be able to resolve these in the same ways as in Sweave (one of which is

Re: [R] vglm

2010-09-01 Thread Gavin Simpson
On Wed, 2010-09-01 at 11:05 +0800, Dejian Zhao wrote: try fit=vglm(mydata[,Loss]~1,pareto1(location=alpha),trace=TRUE,crit=c) No, please don't. That is not a good example of formula use in R. Several responders have already pointed out that the 'R' way of doing this would be to use a data

Re: [R] wavelet parameters with sowas

2010-09-01 Thread Jim Lemon
On 09/01/2010 01:16 AM, angela wrote: Hello everybody. My name is Angela. I'm doing wavelet using the sowas library. The problem I have is that I don't know how to choose the paremeters to describes differnt wavelet analysis. How I select the noctave, nvoice, w0, s0,...? Hi Angela, I located

Re: [R] wavelet parameters with sowas

2010-09-01 Thread angela
Thank you very much for your suggestion Jim, I'll try that. Regards, Angela -- View this message in context: http://r.789695.n4.nabble.com/wavelet-parameters-with-sowas-tp2401776p2402846.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Detecting Growth Trends

2010-09-01 Thread Lorenzo Isella
Date: Tue, 31 Aug 2010 14:37:16 -0700 (PDT) From: rtist patwarner2...@yahoo.com To: r-help@r-project.org Subject: Re: [R] Detecting Growth Trends Message-ID: 1283290636694-2402347.p...@n4.nabble.com Content-Type: text/plain; charset=us-ascii If the test rejects the null, then it has determined

[R] What solve() does?

2010-09-01 Thread Petar Milin
Hello! Can anyone explain me what solve() function does: Gaussian elimination or iterative, numeric solve? In addition, I would need both the Gaussian elimination and iterative solution for the course. Are the two built in R? Thanks! PM __

[R] Emacs/R graphics / Win32

2010-09-01 Thread SKÁLA Zdeněk
Dear all, apologies for a question that is probably more to the Windows that R/Emacs. Still, please, does anybody else encounter this problem? I use R-2.11.1 with GNU Emacs/ESS (the same behaviour with Xemacs 21.4.22/ESS) under WinXP (32bit). All runs ok but when invoking graphics window (any

Re: [R] vglm

2010-09-01 Thread Dejian Zhao
My previous expression is ok. But I agree using argument 'data' will be a better choice especially when there are many variables from the object specified by 'data'. On 2010-9-1 16:23, Gavin Simpson wrote: On Wed, 2010-09-01 at 11:05 +0800, Dejian Zhao wrote: try

Re: [R] transaction object - how to coerce this data

2010-09-01 Thread jim holtman
try this: x.t - split(x$item, x$basketID) x.t $`1` [1] bread cheese milk $`2` [1] bread cheese eggs $`3` [1] bread cheese beer dput(x.t) structure(list(`1` = c(bread, cheese, milk), `2` = c(bread, cheese, eggs), `3` = c(bread, cheese, beer)), .Names = c(1, 2, 3)) On Wed, Sep 1, 2010

Re: [R] lost in vegan package

2010-09-01 Thread Jari Oksanen
eva epitta at upatras.gr writes: Thank you very very much. I thought I wouldn't never find the right way of calculating species pairs occurrences in vegan package. My data are already transposed. Do you know a way of saving the output? I can only copy and paste the results in an Excel

Re: [R] step

2010-09-01 Thread Jari Oksanen
Gavin Simpson gavin.simpson at ucl.ac.uk writes: On Fri, 2010-08-27 at 11:03 -0300, Silvano wrote: Hi, how can I change the significance level in test F to select variable in step command? I used step(model0, ~x1+x2+x3+x4, direction=c(forward), test='F', alpha=.05)

[R] biclust package

2010-09-01 Thread James Nead
Hi, I wanted to compare the quality of biclusters obtained from the various biclustering algorithms. Is there a function/metric in biclust (or some other package) that will enable me to do this? thanks! [[alternative HTML version deleted]]

[R] [R-pkgs] ez version 2.0

2010-09-01 Thread Mike Lawrence
The ez package was developed to aid those that are new to statistical programming. Over the course of several years of helping colleagues and students learn R, I observed that folks are often initially turned off R because they have difficulty obtaining SPSS-like results quickly (SPSS is the

[R] kruskalmc

2010-09-01 Thread Iasonas Lamprianou
Hi, I would like to find out if the kruskalmc command in the pgirmess package controls for the family-wise type-I error rate. Thank you Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006

Re: [R] Putting legend *outside* plotting area

2010-09-01 Thread jim holtman
Try using 'layout' to create the area for the legend; that way you don't have to worry about how much space to leave. Here is an example where I was putting the legend on the bottom of the chart: pdf('plot.pdf', width=62, height=40) layout(rbind(1, 2), height=c(7, 1)) # put legend on bottom

Re: [R] How to Remove Autocorrelation from Simple Moving Average time series

2010-09-01 Thread jim holtman
Have Excel call your R script. On Mon, Aug 30, 2010 at 10:56 AM, Ramesh Kallol rames...@ambaresearch.com wrote: Hi R experts, I am trying to remove autocorrelation from Simple Moving Average time series. I know that this can be done by using seasonal ARIMA like,            library(TTR)

[R] Save data as .pdf or .JPG

2010-09-01 Thread khush ........
Hi all , I have following script to plot some data. plot( c(1,1100), c(0,15), type='n', xlab='', ylab='', ylim=c(0.1,25) , las=2) axis (1, at = seq(0,1100,50), las =2) axis (2, at = seq(0,25,1), las =2) lines( c(1,304),c(0,0),col='grey', lwd=3 ) lines( c(1,238),c(1,1),col='grey', lwd=3 )

Re: [R] lost in vegan package

2010-09-01 Thread eva
Eva, The printed output can be dumped to a file using command sink(). See its help page for info. The oecosimu result is the output of the basic statistic function amended with item oecosimu with structure: str(mod$oecosimu) List of 6 $ z : num -66.5 $ pval : num 0.01 $

Re: [R] Emacs/R graphics / Win32

2010-09-01 Thread RICHARD M. HEIBERGER
Zdeněk, The way to isolate the problem is to start emacs with emacs -q --no-site-files then manually load ess-site.el and start R. We hope that runs normally. Then add additional lines from your site-start.el and your .emacs until you find the problem line. Followup should be on the ESS

Re: [R] Save data as .pdf or .JPG

2010-09-01 Thread Duncan Murdoch
On 01/09/2010 8:56 AM, khush wrote: Hi all , I have following script to plot some data. plot( c(1,1100), c(0,15), type='n', xlab='', ylab='', ylim=c(0.1,25) , las=2) axis (1, at = seq(0,1100,50), las =2) axis (2, at = seq(0,25,1), las =2) lines( c(1,304),c(0,0),col='grey', lwd=3 )

Re: [R] Reading multiple text files where some files are empty

2010-09-01 Thread Dieter Menne
Paul wrote: Thanks for the responses. 'Try' seems to be the function that I needed. I had to adapt the code suggested below because vapply wasn't recognised, but this seemed to work: x‑lapply(a,function(x){try(read.table(x,colClasses='character'))}) There is the risk with this

Re: [R] how to replace NA with a specific score that is dependant on another indicator variable

2010-09-01 Thread Chris Howden
Hi everyone, I’m looking for a clever bit of code to replace NA’s with a specific score depending on an indicator variable. I can see how to do it using lots of if statements but I’m sure there most be a neater, better way of doing it. Any ideas at all will be much appreciated, I’m

Re: [R] 3D tomography data

2010-09-01 Thread Steve
Like I was saying I want to be able to calculate the nearest neighbour and its vector. I think this can be done using pairdist or K3est in the spatstat package. But I have no idea as to how I prepare my data in a form that the software will recognise. How do I turn my tomography data into pp3

Re: [R] pasting together 2 character arrays

2010-09-01 Thread David Reiner
Always good to have that fourth way, just in case (Array3 - as.vector(outer(Array2[-1], Array1, FUN=function(x, y) paste(y, Array2[1], x, sep=. [1] height.trt0.trt1 height.trt0.trt2 weight.trt0.trt1 weight.trt0.trt2 age.trt0.trt1age.trt0.trt2 [7] sex.trt0.trt1sex.trt0.trt2

Re: [R] Remove Loading required package message

2010-09-01 Thread Dieter Menne
Sébastien Moretti wrote: I use --vanilla --quiet --slave options on command line but I always get Loading required package: ... R 2.5.1 .. and make sure you update your R installation. We are at 2.11.1 currently. Dieter -- View this message in context:

Re: [R] how to replace NA with a specific score that is dependant on another indicator variable

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 9:20 AM, Chris Howden wrote: Hi everyone, I’m looking for a clever bit of code to replace NA’s with a specific score depending on an indicator variable. I can see how to do it using lots of if statements but I’m sure there most be a neater, better way of doing it.

[R] Rd-file error: non-ASCII input and no declared encoding

2010-09-01 Thread Christian Hennig
Dear list, I came across the following error for three of my newly written Rd-files: non-ASCII input and no declared encoding I can't make sense of this. Below I copied in one of the three files. Can anybody please tell me what's wrong with it? Thank you, Christian \name{tetragonula}

Re: [R] how to replace NA with a specific score that is dependant on another indicator variable

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 9:55 AM, David Winsemius wrote: On Sep 1, 2010, at 9:20 AM, Chris Howden wrote: Hi everyone, I’m looking for a clever bit of code to replace NA’s with a specific score depending on an indicator variable. I can see how to do it using lots of if statements but I’m

Re: [R] Rd-file error: non-ASCII input and no declared encoding

2010-09-01 Thread Gavin Simpson
On Wed, 2010-09-01 at 15:09 +0100, Christian Hennig wrote: Dear list, I came across the following error for three of my newly written Rd-files: non-ASCII input and no declared encoding I'm not sure, but is it the m- or n-dash character for the page range in the \source{} section? 2317–2331

[R] opls-da

2010-09-01 Thread Francesca Chignola
Dear all, I would like to apply Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) to a metabolomic dataset, in order to discriminate two groups of samples. I have looked for an available R package and I have found K-OPLS and oscorespls.fit (Orthogonal scores PLSR)

Re: [R] Rd-file error: non-ASCII input and no declared encoding

2010-09-01 Thread Duncan Murdoch
On 01/09/2010 10:09 AM, Christian Hennig wrote: Dear list, I came across the following error for three of my newly written Rd-files: non-ASCII input and no declared encoding I don't think the problem survived emailing, at least I didn't see it. You can display all the non-ASCII

Re: [R] how to replace NA with a specific score that is dependant on another indicator variable

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 10:19 AM, David Winsemius wrote: On Sep 1, 2010, at 9:55 AM, David Winsemius wrote: On Sep 1, 2010, at 9:20 AM, Chris Howden wrote: Hi everyone, I’m looking for a clever bit of code to replace NA’s with a specific score depending on an indicator variable. I can

[R] Autocorrelation with merMCMC object

2010-09-01 Thread berthouly cecile
Hello, *I have some issues with furnction that I've used before and that are not working anymore. I want to check the autocorrelation of an object merMCMC but the function autocorr did not acccept merMCMC object.Is there any ither funtion that I could use?*

Re: [R] Rd-file error: non-ASCII input and no declared encoding

2010-09-01 Thread Christian Hennig
Dear Gavin, very well spotted = solved! Thanks also Duncan! Christian On Wed, 1 Sep 2010, Gavin Simpson wrote: On Wed, 2010-09-01 at 15:09 +0100, Christian Hennig wrote: Dear list, I came across the following error for three of my newly written Rd-files: non-ASCII input and no declared

[R] standardize columns selectively within a dataframe

2010-09-01 Thread Olga Lyashevska
Dear all, I have a dataframe: df-dataframe(a=c(1,2,3),b=c(4,5,6),c=c(7,8,9),d=c(10,11,12)) I want to obtain a new dataframe with columns a and b being standardized ((x-mean(x))/sd(x)); the other two columns (c,d) I want to leave unchanged. What is the best way to achieve this? I have been

Re: [R] standardize columns selectively within a dataframe

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 10:35 AM, Olga Lyashevska wrote: Dear all, I have a dataframe: df-dataframe(a=c(1,2,3),b=c(4,5,6),c=c(7,8,9),d=c(10,11,12)) I want to obtain a new dataframe with columns a and b being standardized ((x-mean(x))/sd(x)); the other two columns (c,d) I want to leave

Re: [R] standardize columns selectively within a dataframe

2010-09-01 Thread Adaikalavan Ramasamy
If you want to scale within columns, you could try cbind( scale(df[,1:2]), df[ ,-c(1:2)] ) a b c d 1 -1 -1 7 10 2 0 0 8 11 3 1 1 9 12 and it is data.frame() btw. On 01/09/2010 15:35, Olga Lyashevska wrote: Dear all, I have a dataframe:

Re: [R] standardize columns selectively within a dataframe

2010-09-01 Thread Olga Lyashevska
Thanks! It is exactly what I was looking for! Cheers Olga __ 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 provide commented, minimal,

[R] invert order

2010-09-01 Thread Sebastian Kruk
Dear R-help list users, I have a huge vector of numbers, how I can invert orden? For example x - 1:1000 I would like to obtain x_r - 1000:1 Thanks, Sebastian. __ R-help@r-project.org mailing list

Re: [R] invert order

2010-09-01 Thread Seeliger . Curt
Sebastian writes: I have a huge vector of numbers, how I can invert orden? For example x - 1:1000 I would like to obtain x_r - 1000:1 Well, 1000:1 works for me. If you need to rev()erse a pre-existing vector, try rev(). cur -- Curt Seeliger, Data Ranger Raytheon Information

[R] Using packages built in linux in windows

2010-09-01 Thread rajesh j
Hi, I built a package in linux and would like to use it in windows. I cant build the package again in windows because it would miss the necessary libraries etc(I'm using Rcpp). So I need to generate a binary version and try to install it in windows. I tried this and I get the error

Re: [R] invert order

2010-09-01 Thread Ivan Calandra
Hi, Maybe sort(x, decreasing=TRUE) would do it HTH, Ivan Le 9/1/2010 16:56, Sebastian Kruk a écrit : Dear R-help list users, I have a huge vector of numbers, how I can invert orden? For example x- 1:1000 I would like to obtain x_r- 1000:1 Thanks, Sebastian.

[R] getting column names of row-by-row sorted matrix

2010-09-01 Thread Murali.Menon
Hi folks, I want to sort a matrix row-by-row and create a new matrix that contains the corresponding colnames of the original matrix. E.g. set.seed(123) a - matrix(rnorm(20), ncol=4); colnames(a) - c(A,B,C,D) a A B C D [1,] -0.56047565 1.7150650

[R] Why is vector assignment in R recreates the entire vector ?

2010-09-01 Thread Tal Galili
Hello all, A friend recently brought to my attention that vector assignment actually recreates the entire vector on which the assignment is performed. So for example, the code: x[10]- NA # The original call (short version) Is really doing this: x- replace(x, list=10, values=NA) # The original

Re: [R] Using packages built in linux in windows

2010-09-01 Thread romain
Hi, The proper venue for Rcpp related questions is the Rcpp-devel mailing list: http://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel Of course you can build a package that uses Rcpp on windows, once all the proper tools are installed, it is just the same as building it on

[R] Where can the benchmark2 function be downloaded from?

2010-09-01 Thread Tal Galili
Hello all, I am looking for the function benchmark2 ( a nice function for comparing the performance of two or more functions). I found online it exists in the butler package http://crantastic.org/packages/butler But for some reason, that package was removed from CRAN:

Re: [R] Using packages built in linux in windows

2010-09-01 Thread rajesh j
This problem doesnt seem to have anything to do with Rcpp. I can't build it on windows because I'm using a library thats only available in linux. Isnt this a resolvable error? On Wed, Sep 1, 2010 at 9:07 PM, rom...@r-enthusiasts.com wrote: Hi, The proper venue for Rcpp related questions is

Re: [R] getting column names of row-by-row sorted matrix

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 11:18 AM, murali.me...@avivainvestors.com wrote: Hi folks, I want to sort a matrix row-by-row and create a new matrix that contains the corresponding colnames of the original matrix. E.g. set.seed(123) a - matrix(rnorm(20), ncol=4); colnames(a) - c(A,B,C,D) a

Re: [R] Why is vector assignment in R recreates the entire vector ?

2010-09-01 Thread Bert Gunter
On Wed, Sep 1, 2010 at 8:09 AM, Tal Galili tal.gal...@gmail.com wrote: Hello all, A friend recently brought to my attention that vector assignment actually recreates the entire vector on which the assignment is performed. So for example, the code: x[10]- NA # The original call (short

Re: [R] Where can the benchmark2 function be downloaded from?

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 11:13 AM, Tal Galili wrote: Hello all, I am looking for the function benchmark2 ( a nice function for comparing the performance of two or more functions). I found online it exists in the butler package http://crantastic.org/packages/butler Have you looked at package

Re: [R] Where can the benchmark2 function be downloaded from?

2010-09-01 Thread Tal Galili
Thank you David! That's exactly what I was looking for. Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com

Re: [R] Using packages built in linux in windows

2010-09-01 Thread Duncan Murdoch
On 01/09/2010 11:34 AM, rajesh j wrote: This problem doesnt seem to have anything to do with Rcpp. I can't build it on windows because I'm using a library thats only available in linux. Isnt this a resolvable error? The resolution looks pretty obvious: port the library to Windows. Duncan

Re: [R] Using packages built in linux in windows

2010-09-01 Thread rajesh j
I thought building a binary took care of that issue. How do I port a linux library to windows? On Wed, Sep 1, 2010 at 9:12 PM, Duncan Murdoch murdoch.dun...@gmail.comwrote: On 01/09/2010 11:34 AM, rajesh j wrote: This problem doesnt seem to have anything to do with Rcpp. I can't build it on

Re: [R] getting column names of row-by-row sorted matrix

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 11:34 AM, David Winsemius wrote: snipped t( apply(a, 1, function(x) colnames(a)[order(x)]) ) [,1] [,2] [,3] [,4] [1,] A C B D [2,] A C B D [3,] D B C A [4,] B A C D [5,] C D B A (apply returns a transposed version.) That is over-broad. apply() only

Re: [R] how to replace NA with a specific score that is dependant on another indicator variable

2010-09-01 Thread Hadley Wickham
first ddply result did I see that some sort of misregistration had occurred; Better with: res -ddply(egraw2, .(category), .fun=function(df) {               sapply(df,                    function(x) {mnx - mean(x, na.rm=TRUE);                                 sapply(x, function(z) if

[R] xyplot (Lattice): Group order in display

2010-09-01 Thread jrflanders
Dear R users, I would like one of my groups in xyplot to appear to lie 'behind' the other groups. I have searched for help and find many, many topics about panel order (e.g., as.table), but that is not what I need. What is odd is that the group that I want to appear 'behind' the others

Re: [R] Using packages built in linux in windows

2010-09-01 Thread Duncan Murdoch
On 01/09/2010 11:51 AM, rajesh j wrote: I thought building a binary took care of that issue. How do I port a linux library to windows? Building a binary package links the package to the library, but if the library doesn't exist, there's nothing to link to. Porting a Linux library to

Re: [R] Why is vector assignment in R recreates the entire vector ?

2010-09-01 Thread Tal Galili
Thank you for the explanation Duncan - very interesting indeed! I wonder if someone in the list might know to answer your question regarding the double duplication. Best, Tal Contact Details:--- Contact me: tal.gal...@gmail.com

Re: [R] Using packages built in linux in windows

2010-09-01 Thread Jeff Newmiller
Normally, a binary package built on one platform will not work on another. While it is possible to cross-compile, in most cases it is best to assemble the necessary development environment oNoten the target platform (including the appropriate libraries compiled for that platform) and build the

Re: [R] Using packages built in linux in windows

2010-09-01 Thread rajesh j
Unfortunately, the library is built with visual studio in windows. It doesnt work with MinGW. On Wed, Sep 1, 2010 at 9:35 PM, Jeff Newmiller jdnew...@dcn.davis.ca.uswrote: Normally, a binary package built on one platform will not work on another. While it is possible to cross-compile, in most

Re: [R] Why is vector assignment in R recreates the entire vector ?

2010-09-01 Thread Matt Shotwell
Tal, For your first example, x is not duplicated in memory. If you compile R with --enable-memory-profiling, you have access to the tracemem() function, which will report whether x is duplicate()d: x - rep(1,100) tracemem(x) [1] 0x8f71c38 x[10] - NA This does not result in duplication of x,

Re: [R] Using packages built in linux in windows

2010-09-01 Thread Ben Bolker
rajesh j akshay.rajesh at gmail.com writes: I thought building a binary took care of that issue. How do I port a linux library to windows? I think we need a little clarification here. When you say a linux library, do you mean a chunk of compiled (non-R) code that the R package (sic) needs

[R] documentation to upgrade R-package from 32 to 64bit

2010-09-01 Thread Hayes, Daniel
Dear all, I am working with the an R-package named GAMLSS (www.gamlss.comhttp://www.gamlss.com) it is currently only functional under the 32-bit version of R (for windows) The author of the package has agreed to help me create 64-bit compatible version. I've been looking through the available

[R] [R-pkgs] testthat: version 0.3

2010-09-01 Thread Hadley Wickham
# testthat Testing your code is normally painful and boring. `testthat` tries to make testing as fun as possible, so that you get a visceral satisfaction from writing tests. Testing should be fun, not a drag, so you do it all the time. To make that happen, `testthat`: * Provides functions that

[R] [R-pkgs] New package: lubridate 0.1

2010-09-01 Thread Garrett Grolemund
Please find below the package announcement for the R package lubridate available from cran. Thank you, Garrett Grolemund Rice University ##lubridate Date-time data can be frustrating to work with in R. R commands for date-times are generally unintuitive and change depending on the type of

Re: [R] How to specify only minimum ylim parameter for boxplot

2010-09-01 Thread mamunbabu2001
Hi, I recently have a similar problem. since your post is about 4 months old I assume you have already got the answer. Try boxplot (data, ylim = c(min_lim, max_lim)) -- View this message in context:

Re: [R] documentation to upgrade R-package from 32 to 64bit

2010-09-01 Thread Matt Shotwell
Try this: http://www.stats.ox.ac.uk/~ripley/Win64/W64porting.html -Matt On Wed, 2010-09-01 at 07:40 -0400, Hayes, Daniel wrote: Dear all, I am working with the an R-package named GAMLSS (www.gamlss.comhttp://www.gamlss.com) it is currently only functional under the 32-bit version of R

[R] how to replace NA with a specific score that is dependant on another indicator variable

2010-09-01 Thread Chris Howden
Hi everyone, I’m looking for a clever bit of code to replace NA’s with a specific score depending on an indicator variable. I can see how to do it using lots of if statements but I’m sure there most be a neater, better way of doing it. Any ideas at all will be much appreciated, I’m

Re: [R] non-linear plot parameters

2010-09-01 Thread Marlin Keith Cox
This is one data set. I have used linear models, but as these data reflect fish consuming energy during a starvation period, there is a lower asymptote. I also realize that data is rather scattered and not the cleanest. I am interested in looking at how fast the animals reach the

Re: [R] documentation to upgrade R-package from 32 to 64bit

2010-09-01 Thread Duncan Murdoch
On 01/09/2010 7:40 AM, Hayes, Daniel wrote: Dear all, I am working with the an R-package named GAMLSS (www.gamlss.comhttp://www.gamlss.com) it is currently only functional under the 32-bit version of R (for windows) The author of the package has agreed to help me create 64-bit compatible

[R] ggplot2 multiple group barchart

2010-09-01 Thread Waller Gregor (wall)
hi there.. i got a problem with ggplot2. here my example: library (ggplot2) v1 - c(1,2,3,3,4) v2 - c(4,3,1,1,9) v3 - c(3,5,7,2,9) gender - c(m,f,m,f,f) d.data - data.frame (v1, v2, v3, gender) d.data x - names (d.data[1:3]) y - mean (d.data[1:3]) pl - ggplot (data=d.data, aes

[R] xyplot (Lattice): Group order in display

2010-09-01 Thread jrflanders
Dear R users, I would like one of my groups in xyplot to appear to lie 'behind' the other groups. I have searched for help and find many, many topics about panel order (e.g., as.table), but that is not what I need. What is odd is that the group that I want to appear 'behind' the others

Re: [R] Why is vector assignment in R recreates the entire vector ?

2010-09-01 Thread Duncan Murdoch
On 01/09/2010 11:09 AM, Tal Galili wrote: Hello all, A friend recently brought to my attention that vector assignment actually recreates the entire vector on which the assignment is performed. So for example, the code: x[10]- NA # The original call (short version) Is really doing this: x-

[R] ggplot2 case study competiton and wiki

2010-09-01 Thread Garrett Grolemund
Hi all, We've moved the ggplot2 wiki to a new location to make it more useful, http://github.com/hadley/ggplot2/wiki. To publicize the new location, we’re organizing a second ggplot2 case study competition. Some great prizes are at stake: * grand prize will be a 32GB ipod touch generously

[R] [Q] Goodness-of-fit test of a logistic regression model using rms package

2010-09-01 Thread GMail (KU)
Hello, I was looking for a way to evaluate the goodness-of-fit of a logistic regression model. After googling, I found that I could use resid(fit, 'gof') method implemented in the rms package. However, since I am not used to the le Cessie-van Houwelingen normal test statistic, I do not know

Re: [R] standardize columns selectively within a dataframe

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 10:42 AM, David Winsemius wrote: On Sep 1, 2010, at 10:35 AM, Olga Lyashevska wrote: Dear all, I have a dataframe: df-dataframe(a=c(1,2,3),b=c(4,5,6),c=c(7,8,9),d=c(10,11,12)) I want to obtain a new dataframe with columns a and b being standardized

Re: [R] standardize columns selectively within a dataframe

2010-09-01 Thread Olga Lyashevska
On Wed, 2010-09-01 at 12:42 -0400, David Winsemius wrote: I suspect you might have tried (df-mean(df))/sd(x) and gotten unsatisfactory results; I know I did. yes, indeed! a few times, but why is that? If you had really wanted to persist and do it from first principles, so to speak, or

[R] problems with unique

2010-09-01 Thread n.via...@libero.it
Hi everyone, I have a question concerning the unique function. I'm using it with two different data frame,which are similar in structure, but as in one case unique works properly in the other doesn't work and I don't understand why. Someone knows why sometimes unique doesn't do things properly??

[R] how to represent error bar in a plot legend?

2010-09-01 Thread Josef . Kardos
I have a simple barplot of 4 mean values, each mean value has an associated 95% confidence interval drawn on the plot as an error bar. I want to make a legend on the plot that uses the error bar symbol, and explains 95% C.I. How do I show the error bar symbol in the legend? I could not find

Re: [R] problems with unique

2010-09-01 Thread Ben Bolker
n.vialma at libero.it n.vialma at libero.it writes: Hi everyone, I have a question concerning the unique function. I'm using it with two different data frame,which are similar in structure, but as in one case unique works properly in the other doesn't work and I don't understand why.

Re: [R] [Q] Goodness-of-fit test of a logistic regression model using rms package

2010-09-01 Thread Frank Harrell
On Wed, 1 Sep 2010, GMail (KU) wrote: Hello, I was looking for a way to evaluate the goodness-of-fit of a logistic regression model. After googling, I found that I could use resid(fit, 'gof') method implemented in the rms package. However, since I am not used to the le Cessie-van

Re: [R] how to represent error bar in a plot legend?

2010-09-01 Thread John Kane
Can you give us a working example of what you are doing? I don't understand why you would need a pch symbol for an error bar. Perhaps if you are using an unboxed legend you could just use arrow() and text() to add the error bar to the legend? --- On Wed, 9/1/10, josef.kar...@phila.gov

Re: [R] [Q] Goodness-of-fit test of a logistic regression model using rms package

2010-09-01 Thread Bert Gunter
... and, furthermore, in most real world situations there are several -- or even lots -- of quite different, incomparable models that give essentially equivalent fits. Distinguishing among the alternatives typically requires focused studies designed for the task. Indeed, as Brian Joiner remarked

Re: [R] xyplot (Lattice): Group order in display

2010-09-01 Thread Deepayan Sarkar
On Wed, Sep 1, 2010 at 9:00 AM, jrflanders jrfland...@gmail.com wrote: Dear R users, I would like one of my groups in xyplot to appear to lie 'behind' the other groups. I have searched for help and find many, many topics about panel order (e.g., as.table), but that is not what I need. What

Re: [R] xyplot (Lattice): Group order in display

2010-09-01 Thread jrflanders
Thank you! It is always so simple... May I throw another question at you? In useOuterStrips, in the latticeExtra package, I am having problems calling custom strip names, e.g. useOuterStrips(xyplot(y-x | LOCATION + EVENT, data = TC6, groups = TYPE), strip=strip.custom(factor. levels = c(Left,

[R] CRAN Mirror definition error.

2010-09-01 Thread Lenart, Joseph
Hello, I was wondering if anyone has run across this problem in the past and might be able to suggest a solution. I did a review of the FAQ's link trying to find a topic that closely related to the experience I was having, however I didn't see anything that might help in diagnosing the

[R] Legend with fill = gray ?

2010-09-01 Thread Jens Oldeland
Hi, I am facing a problem with the legend, I don´t know how to use the fill option in the legend in order to achieve the same standard gray levels that are plotted. Sorry for this easy one, but I really did not find anything so far. It works fine with color: ## C O L O R mat -

Re: [R] Legend with fill = gray ?

2010-09-01 Thread Greg Snow
The phrase does not work can mean many things, it would help us to help you if you spelled out what you expected and how what you saw differs from what was expected. I suspect part of your problem may be in using 0.1:0.3, did you try typing just that at the command line to see the result? The

Re: [R] Legend with fill = gray ?

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 4:17 PM, Jens Oldeland wrote: Hi, I am facing a problem with the legend, I don´t know how to use the fill option in the legend in order to achieve the same standard gray levels that are plotted. Sorry for this easy one, but I really did not find anything so far. It

Re: [R] lattice: plot alignment fails

2010-09-01 Thread Marius Hofert
Okay, I solved it. After using the same arguments for scales (in particular: relation=free), the plots are perfectly aligned. Cheers, Marius On 2010-08-24, at 20:39 , Marius Hofert wrote: Dear expeRts, I would like to have four plots appearing in one figure. The minimal example shows

Re: [R] Legend with fill = gray ?

2010-09-01 Thread Jens Oldeland
Hi, sorry for beeing unprecise. I am looking for a legend that has the same colors as the barplot. For a coloured version I would simply use the same fill=rainbow() option as used in creating the barplot. However, this does not work when there are gray values, because there are not defined in the

[R] [semi-OT] Using fortune() in an email signature

2010-09-01 Thread Stuart Luppescu
Hello, As you can see from my signature in this message, I use the R fortune function to generate a fortune, which is then fed to the signature program, which constructs a named pipe containing the fortune-bearing sig, which is then included in mail messages. The problem is that it's got

Re: [R] [semi-OT] Using fortune() in an email signature

2010-09-01 Thread David Winsemius
On Sep 1, 2010, at 4:49 PM, Stuart Luppescu wrote: Hello, As you can see from my signature in this message, I use the R fortune function to generate a fortune, which is then fed to the signature program, which constructs a named pipe containing the fortune-bearing sig, which is then included

Re: [R] Legend with fill = gray ?

2010-09-01 Thread Greg Snow
Did you try gray.colors(3)? -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: Jens Oldeland [mailto:oldel...@gmx.de] Sent: Wednesday, September 01, 2010 2:43 PM To: Greg Snow Cc:

Re: [R] Legend with fill = gray ?

2010-09-01 Thread Jens Oldeland
No...but it works! I do not know why I did not find it... Many thanks! Hi David, thanks to you as well. Greg Snow schrieb: Did you try gray.colors(3)? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] general question on binomial test / sign test

2010-09-01 Thread Kay Cichini
hello, i did several binomial tests and noticed for one sparse dataset that binom.test(1,1,0.5) gives a p-value of 1 for the null, what i can't quite grasp. that would say that the a prob of 1/2 has p-value of 0 ?? - i must be wrong but can't figure out the right interpretation.. best, kay

Re: [R] general question on binomial test / sign test

2010-09-01 Thread Greg Snow
Try thinking this one through from first principles, you are essentially saying that your null hypothesis is that you are flipping a fair coin and you want to do a 2-tailed test. You then flip the coin exactly once, what do you expect to happen? The p-value of 1 just means that what you saw

Re: [R] [semi-OT] Using fortune() in an email signature

2010-09-01 Thread Dirk Eddelbuettel
On 1 September 2010 at 15:49, Stuart Luppescu wrote: | Hello, As you can see from my signature in this message, I use the R | fortune function to generate a fortune, which is then fed to the | signature program, which constructs a named pipe containing the | fortune-bearing sig, which is then

Re: [R] [semi-OT] Using fortune() in an email signature

2010-09-01 Thread Christian Raschke
Simply using option --slave instead seems to do the trick. $ /usr/bin/R --slave test-fortune.R Similarly to Brian, I'm much more reluctant to help people who don't exist -- who knows, maybe you're a computer program who has just passed the Turing test :-) -- Martin Maechler (about

Re: [R] [semi-OT] Using fortune() in an email signature

2010-09-01 Thread Matt Shotwell
Or using R GNU tools: m...@max:~$ R -e fortunes::fortune() | gawk '/^[^]/ {print}' It's not a question of trying variations, rather of following instructions. -- Brian D. Ripley (about using 'Writing R Extensions') R-help (January 2006) -Matt On Wed, 2010-09-01 at 16:49 -0400, Stuart

  1   2   >