similar to: plot 3d info in 2d

Suppose that I'm working on Hadley's diamond dataset and I want to review the relationship between price, colour and carat. I might run the following: library(ggplot2) #plot scatter and add some hex binning q<-qplot(carat,price,data=diamonds, geom=c("hex"), main="Variability of Diamond Prices by Carat and Colour") #facet to get one scatter for

Dear R-users, I am trying to make a plot with ggplot-geom_tile(), but cannot remove some unwanted (white) lines through my plot. Below a reproducible example: ##### library(ggplot2) tot=as.data.frame(rep(seq(-50,50,5),each=21)) names(tot)="precip" temp=rep(seq(-5,5,0.5),21) tot$temp=temp disc=array(dim=c(21,21)) for(i in 1:21){ for(y in 1:21){ temp<-

Hello everyone, I was asked to repost this again, sorry for any inconvenience. I'm looking replacement for ddply function from plyr package. Function allows to apply function by category stored in any column/columns. Regular loops or lapplys slow down greatly because my unique combination count exceeds 9000. Is there any available solution which allow me to apply function by category?

Hi, I have a problem on my mac when trying in R to produce png images. I am getting this warnings with the ArrayQualityMetrics package: > arrayQualityMetrics(rma_fatBody, outdir="normData", force =T) The report will be written into directory 'normData'. KernSmooth 2.23 loaded Copyright M. P. Wand 1997-2009 (loaded the KernSmooth namespace) libpng warning: Application built

Dear All, when I'm running a PCA with prcomp(USArrests, scale = TRUE) I get the right principal components, but with the wrong sign infront Rotation: PC1 PC2 PC3 PC4 Murder 0.5358995 -0.4181809 0.3412327 0.64922780 Assault 0.5831836 -0.1879856 0.2681484 -0.74340748 UrbanPop 0.2781909 0.8728062 0.3780158 0.13387773 Rape 0.5434321 0.1673186 -0.8177779 0.08902432 instead of PC1 PC2 PC3 PC4

Hi All, I get the message failed with message ‘Dimension out of range’ when using cuhre in package R2Cuba. Does anyone know what this mean? Or would I need to email the package author? The funny thing is it does give a result and comparing it to "adaptIntegrate" in package cubature, the two numbers are very close. Thanks, Sachin [[alternative HTML version deleted]]

Hello I have the following problem when loading the package hydroGOF on Windows 7 running R.12.2 library(hydroGOF) Error : package 'hydroTSM' does not have a name space Error: package/namespace load failed for 'hydroGOF' The same command does not result in error on R.13.1 at my Mac running Lion. Have I done something wrong? Many thanks Ed

Hi Users, I want to fit the maximums of a data series with a cubic spline. How do I go about this in R. I failed to figure out how I can use the mgcv library to do this. Thanks ---------------------------- ZABLONE

Dear All, The following code save my graphs as pdf: pdf("j:/mix.pdf", width = 18, height = 16) grid.newpage() pushViewport(viewport(layout = grid.layout(3,1))) vplayout <- function(x, y) viewport(layout.pos.row = x, layout.pos.col = y) print(Aplot, vp = vplayout(1, 1)) print(Bplot, vp = vplayout(2, 1)) print(Cplot, vp = vplayout(3, 1)) dev.off() How can I save it in

Hi list   I want to define a moving window based on ferquency in a grid file. is there a package for this? I found packages for moving window mean, median or quantile.   thanks  [[alternative HTML version deleted]]

Dear R users, Just to clarify. I am not offering to pay someone to do my Dissertation. These 4-5 commands on Kalman Filter would be only a tiny part of my 10,000 words dissertation. A part that even after trying for a few days, I am still stuck on. I am offering ?50, just to say thanks. Regards -- View this message in context:

Dear all, I would like to use the following function fitdist(data, distr, method=c("mle", "mme", "qme", "mge"), start=NULL, fix.arg=NULL, ...) for many different distr values like distr=c("norm","lnorm","pois") (just a small example) and take back into a list the parameter name which is what is inside distr plus what the

Hi When we use the sav_gol command in R , it shows an error which says: " error in as.matrix". We've downloaded the necessary packages. Kindly help us with this issue. If there is any other function to perform Savitzky-Golay smoothing in R, please let me know. With Regards Reynolds

Hi, Sorry for my ignorance, but how can I see if a function is 'compiled code' or plain R? E.g. the daisy function from the cluster package.

Hi everybody, my problem is the following: I have a function which produces for an input number a dataframe with a fixed number of columns and 0 to 10 rows. Now i want to apply this function to a vector with different inputs and merge all these result-dataframes in one data.frame. Can i give a data.frame-reference to the function, so that the result of the function is not returned but appended

Dear Sir, Please help me to open my rda files. It's urgent to see the content and even if I looked for solutions i didn't find one. When I try to load my .rda file I receive this message: "R> load(file = "SpecPhantomExp2WaterSupp Set2.rda") Error: bad restore file magic number (file may be corrupted) -- no data loaded In addition: Warning message: file

I am trying to convert a matlab code to R. Most of the conversion statements are simple, but problem is with return. Since it's a recursive function it's a bit complicated. Anyway in matlab we can have return statment like [price, notional] = functionname() and also for the same function... notional(1) = functionname() how to do this in R? i have tried using data frame but still

Hi   I have a problem for saving grid file in R. I want to save a grid file as "tif".   I write this: writeGDAL(predict.grid, 'predict.tif')   after running, I don't recieve any error. But when I want to open the 'file.tif', it doesn't open. for example in paint program: I recieve this: Paint cannot read this file. This is not valid bitmap file, or its format

Hi   I want to make a grid map using diversity indices. I have a soil grid map with 7 class in r. how can I make a diversity map based on soil map in r?  I found something for calculation diversity in a table (.csv).   Thanks [[alternative HTML version deleted]]

Hi there, I have a complex math equation which does not have a closed form solution. It is - y <- (p*exp(-a*d)*(1-exp((d-p)*(a-x[1]))))/((p-d)*(1-exp(-p*(a-x[1])))) For this equation, I have all the values except for x[1]. So I need to solve this problem numerically. Can anyone suggest an optimization package that I can use to estimate the value for x[1]? Thanks in advance, Diviya