similar to: newby trying to solve a system

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

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

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,

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

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

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 R users/developers I requested help to solve the problem of formulating Multivariate Sample selection model by using Full Information Maximum Likelihood (FIML)estimation method. I could not get any response. I formulated the following code of FIML to analyse univariate sample selection problem. Would you please advise me where is my problem library (sem) library(nrmlepln) Selection

Hi there everyone, I have the following model (this is naturally a simplified version just for showing my problem, in case you're wondering this is a translog cost function with the associated cost share equations): C ~ á + â1 log X + â2 log Y + ã1 log Z + ã2 log XX C1 ~ â1 + â2 log YY + ã1 log ZZ Then I have some restrictions on the coefficients, namely that the sum of â equal 1 and the