similar to: Re : Re: setting a number of values to NA over a data.frame.

Hi Hadley, thank you for providing this "scagnostics" primer.... I was trying to do some basic testing, and I see that I probably missed some points : first it's not clear for me if the argument of "scagnostics" should be raw data or "processed" data (results of calling "splom" or whatever...). If the first, I thought (from Wilkinson & al.) that if

Hi list, My attention was drawn to the foreach package by recent posts...I decided to have a look... I'm using R.2.9.1 on Windows, I have downloaded the foreach package today (v 1.2.1), together with iterators (v. 1.0.1) and codetools (v.0.2-2). Full of hope I try the most simple thing of all out of the package vignette : > x <- foreach(i = 1:3) %do% sqrt(i) and get : > Erreur

Hello eveRybody (and probably special regards to Gabor G....), I recently upgraded from R-2.9.2-win32 to R-10.1.0-win32, and I'm using Windows XP Professional 2002 with service pack 3 on a PC with IntelCore Duo CPU@ 3.00 GHz. Last time I used it (some weeks ago, before upgrading), one of my scripts (including the two lines cited below) used to run OK. When I ran it yesterday evening (after

Hello, list why not add the smart proposal by Greg Snow as a built-in function in {stats}, just changing the "x234" and "newc" lines to allow for more distributions to be generated ? Or do I miss an already existing function to do that ? Regards. Olivier # slight modification of the original code by Greg Snow [mailto:Greg.Snow at intermountainmail.org] # on April 04, 2007

Sorry if this answer was already given, or if I miss the point, but did you have a look to "lowerTriangle" and "upperTriangle" functions in the "gdata" package ? # example # A<-matrix(rnorm(9),3,3) # B<-B<-matrix(NA,dim(A)[1],dim(A)[2]) # lowerTriangle(B)<-lowerTriangle(A) # upperTriangle(B)<-lowerTriangle(A) # diag(B)<-diag(A) Hope this helps,

Dear list, (sorry if I post to the wrong place...), Though having spent some time on it, I cannot find an answer by myself to the following behaviour of the read.image function (adimpro package) : I'm running R.2.6.2 on Windows XP. The home directory is C:\Program files\R\R-2.6.2 Version 0.4.4 of adimpro is loaded. ImageMagick 6.3.0 is installed, its directory is C:\Program

Hi John, Unless I miss a point, why dont you try something like : # some fake data > fake<-as.data.frame(cbind(seq(1,10,by=1),c(rep(1,4),rep(0,4),rep(2,2)))) V1 V2 1 1 1 2 2 1 3 3 1 4 4 1 5 5 0 6 6 0 7 7 0 8 8 0 9 9 2 10 10 2 # change 0 by NA > fake[fake==0]<-NA # or fake$V2[fake$V2==0]<-NA if you

Dear All, Are there any existing scripts to plot Taylor Diagram (definition see http://www.ipsl.jussieu.fr/~jmesce/Taylor_diagram/taylor_diagram_definition.html ) ? Thanks a lot! Linda [[alternative HTML version deleted]]

Hi, I'm looking for software that can perform kriging on systems with dimensionality higher than 3, say d=5. Are anyone aware of packages in R that can do this? Thanks, Eivind Sm??rgrav ------------------------------------------------------------------- The information contained in this message may be CONFIDENTIAL and is intended for the addressee only. Any unauthorised use,

Hi, I am not sure to get the issue, but assuming your data are arranged as in your example with dates in column x$V1 and signals as x$V2 , I think that you could use "rle" in the following way : test<-rle(x$V2) testmat<-matrix(NA,length(test$values),2) testmat[,1]<-x[c(1,cumsum(test$length)[1:(length(test$values)-1)]+1),1] testmat[,2]<-test$values And you'll get your

Hi R-listers, I run R version 2.13.0, on a i386-pc-mingw32/i386 platform under Windows XP. I perform daily update of my installed packages. I've got the most recent Java for this platform installed, as I am told by the JavaUpdater when checking. (and I've read the posting guide, so that I'm not going to get fired by Uwe as poor wwreith [joke])... But : since one of the last daily

Dear List.... sorry to bother you R-gurus with such an "unstatistical" question... but I face a problem using Crimson Editor with R 2.01.1 that I never had using R 2.00.1. I already posted on the Crimson Editor forum but it seems to be VERY few R-users there.... I successfully used R v2.00.1until now (under Windows XP professionnal, version 2002, Service Pack 2, P4 processor CPU 1.8

Hello there, any ideas on how to save all the objects on my workspace inside a list object? For example, say my workspace is as follows ls() [1] "x" "y" "z" and suppose I want to put these objects inside a list object, say object.list <- list() without having to explicitly write down their names as in object.list$x = x object.list$y = y object.list$z = z Is

Hi list, I thought I was used to layouts, but today I am facing a problem I cannot overcome : On my R installation (Windows 7 Pro, SP1, R version 2.13.0, daily update of packages), I am not able to put raster plots in user defined layouts : > layout.matrix<-matrix(c(1,2,3,4,5,5),2,3) > layout(mat=layout.matrix) > layout.show(5) works fine, I get the correct frames in

Hi listeRs, Maybe I have missed something, but I am facing a problem I don't understand even after reading archives and manuals (…maybe not carefully enough !) : On my PC (64-bit OS, Intel Xeon CPU , two L5506@2.13 GHz processors, 12 Go RAM) under Windows 7 SP1, I run both 32-bit (i386-pc-mingw32/i386) and 64-bit (x86_64-pc-mingw32/x64) version 2.14.2 (2012-02-29) of R, packages

Hello to everybody, ... and sorry if the question has already been answered. I am dealing with complex numbers and need a way of plotting them on angular plots (2D and 3D). Is there a package that already manage this, or can somebody in this R-world who could help me with his skill ? Thanks to all. Have a nice week-end ! Olivier Eterradossi Ph.D., "PsychoSensory Properties of

Dear R-gurus, being very new to R, (as well as lazy and not too smart !) I have some problems (and get lost in the docs) trying to write something to find the 9 values (A1,B1,C1,A2,B2,C2....C3) which are solutions of a 12 equations system of the form : > x1-(A1/(A1+B1+C1)) = 0 > y1-(B1/(A1+B1+C1))= 0 > z1-(C1/(A1+B1+C1)) = 0 > 3 same equations with subscript 2 > 3 same

Hi again, this is a follow-up of my previous post : I've noticed that > attr(first.list[[1]]$size.ecdf,"call") gives ecdf(test.moms[, "m.pxs"]) and that there is no "test.moms" field in my records (a local variable named test.moms was used inside a function, but it was stored under another name in the list, outside the function) May this be the answer (and

Hi there I have some timeseries data which I plot in a OHLC Plot. In the same plot I'd like to have the EMA of this timeseries. I tried to add the EMA point to OHLC with lines(), but this doesn't work. Has anyone an idea how to handle it? Regards, Michael Zak

Hi, I'm having trouble with quantmod's addTA plotting functions. They seem to work fine when run from the command line. But when run inside a function, only the last one run is visible. Here's an example. test.addTA <- function(from = "2010-06-01") { getSymbols("^GSPC", from = from) GSPC.close <- GSPC[,"GSPC.Close"] GSPC.EMA.3