similar to: I might be dumb : a simple question about "foreach"

Hi again, Awfully sorry John, I should have been sleeping and did not see your full post.... here is a way, unless I miss the point again : fake<-as.data.frame(cbind(seq(1,10,by=1),c(rep(1,4),rep(0,4),rep(2,2)))) # from my previous post # one moree column this time ! fake3<-cbind(fake,fake$V2) index<-c(2,3) fake3[,index][fake3[,index]==0]<-NA not nice, but seems to do the job.

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

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

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]]

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

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,

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,

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

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

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

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

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

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

Can somebody help me in finding package/Example in R which could do EMA(Exponential Mean Average). I installed "TTR" package but the 'EMA" function which I was trying to use is giving the following error. Error: Could not find function "EMA" Thanks & Regards, Thomas [[alternative HTML version deleted]]

Hi Does anyone know of any R routines to send emails from R, under Windows? I thought about writing such a facility using the R(D)COM package to drive e.g. MS Outlook, but I don't want to reinvent the wheel. I have found a function Sys.mail in the library syskern, but this only works under Unix by shelling out a mail command. Thanks, David