similar to: plot method for rasters and layout

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

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

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

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

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

Hi, I am trying to produce a levelplot using the following: library(lattice) df <- read.table("data", sep="\t", header=TRUE) gr = levelplot(df$z ~ df$x * df$y, xlim=0:1000,ylim=0:1000,aspect="iso",col.regions=heat.colors) print(gr) The produced result is: http://dl.dropbox.com/u/25473/bigplot.png Now this is correct, because each of those points have a very

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

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 everyone! I would like to draw a line surrounding my raster, with i want it in black, and i don't know if its possible to change the colour. A simple example is: r <- raster(nrow=30, ncol=30, xmn=0) r[]<-NA r[393:409]<-99 r[423:439]<-99 r[453:455]<-99 r[456:460]<-30 r[461:469]<-99 r[483:499]<-99 r[513:529]<-99 plot(r,col=terrain.colors(100)) plot( edge(r,

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

I am using this code to calculate the moving average mean.It worked fine but when I wanted to also calculate based on sd(stander deviation) I got the error shown below. I read this documentation of R movingFun and found that sd was mentioned at http://hosho.ees.hokudai.ac.jp/~kubo/Rdoc/library/raster/html/movingFun.html I wonder if needs any other things to conseder or sd is not at all supplied

Just wanted to post this in case someone else was faced with the same problem. It''s database-specific code, worse, it''s specific to a particular * version* of a particular *extension* to a particular database, wires in SQL, and uses eval on a string of data from the database! So, you know, I''m getting coal in my stocking from DHH this year. I can deal. (After the

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

I have several rasters that I want to do some calculations ,basically calculating the moving average. dir2 <- list.files("D:\\2010+2011", "*.bin", full.names = TRUE) saf=stack(dir2) movi <- overlay(stack(saf),fun=function(x) movingFun(x, fun=mean, n=3, na.rm=TRUE)) Error in .overlayList(x, fun = fun, filename = filename,

Dear Everyone, I wrote a Taylor diagram program in R, but it was not very general, so I was happy to find the Taylor.diagram in the plotrix package. On the other hand I can find many problems with the pos.cor=TRUE case, in other words, when we use only the first quarter of the space (positive correlations). (When we choose pos.cor=FALSE, the program seems to me perfect) 1.There is only one line