similar to: How do I define moving window fequency

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

I applied bootpred for multinomial logistic reg. (with nnet package). I used same as theta.fit and theta.predict of R for my data. but give me error. Can I do this with response vriable;7 levels predictor variables:5 (1 classifier, 4 continuous)?   Thanks alot Azam   [[alternative HTML version deleted]]

Hi   I want to do leave-one-out cross-validation for multinomial logistic regression in R. I did multinomial logistic reg. by package nnet in R. How I do validation? by which function? response variable has 7 levels   please help me   Thanks alot Azam [[alternative HTML version deleted]]

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

Hi Listers, Is it possible to produce an ordination plot in 2d, where bubbles represent the location of sites (this part is easy enough) and the size of the bubbles is proportional to the sites location in 3d space (I am stuck on this option). So sites that are very near the 2d plane of the xy axes would be larger while sites that are actually further away in 3 d space would be proportionally

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   Pleas help me to increase the memory in R.   In order increase memory, I read the FAQ and follow the instruction as below     Close R, then right-click on your R program icon (the icon on your desktop or in your programs directory). Select ``Properties'', and then select the ``Shortcut'' tab. Look for the ``Target'' field and after the closing quotes around the

Hi   I did validation for prediction by logistic regression according to following:   validationsize <- 23 set.seed(1) random<-runif(123) order(random) nrprofilesinsample<-sort(order(random)[1:100]) profilesample <- data[nrprofilesinsample,] profilevalidation <- data[-nrprofilesinsample,] salich<-profilesample$SALIC.H.1 salic.lr<-glm(salich~wetnessindex, profilesample,

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

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.