see ?findIntervals: vec refers to the second argument, i.e. bins in your call. That one has to be sorted. Why do you think that mean-sd is larger than min? Most things in this world are not normally distributed.
On 04/02/2014 06:28 PM, Alessandra Carioli wrote: > Dear Bloggers, > I’ve been working with a for loop to create maps of dependent variables to > speed up a bit the work but I keep getting this error message > > Error in findInterval(varofint2, bins, rightmost.closed = T) : 'vec' must be > sorted non-decreasingly > > I am not sure what to do (I have tried sorting the variables, just to see > what happens but I still get the message). > Any help on the matter is very much appreciated! > > Ale > > R code: > > library(spdep) > library(RColorBrewer) > > var.data<- data.frame(var1,var2,var3,var4,var5,var6,var7,var8,var9) > var.names<-c(“name variable 1”,…,”name variable 9" > > m <- length(var.data) > var.color <- matrix(nrow=n,ncol=m, data=rep(0,n*m)) > for (i in 1:m) { > varofint <- var.data[,i] > mean <- mean(varofint) > sd <- sd(varofint) > min <- min(varofint) > max <- max(varofint) > bins <- c(min,mean-sd, mean, mean+sd, max) > lb <- length(bins) > med <- median(1:lb) > colpal <- brewer.pal(lb, "PRGn")[-med] > var.color[,i] <- colpal[findInterval(varofint, bins, rightmost.closed=T)] > #this is what does not work > plot(shape.shp, col=var.color[,i], bg="grey5",lty=0) #shape.shp name of > shape file > title(var.names[[i]], cex=1) > legend("bottomright",fill=colpal,legend=paste(round(bins[-length(bins)],2), > "-", > round(bins[-1],2)),cex=0.7, bg="white") > } > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics (ifgi), University of Münster Heisenbergstraße 2, 48149 Münster, Germany. Phone: +49 251 83 33081 http://ifgi.uni-muenster.de GPG key ID 0xAC227795
signature.asc
Description: OpenPGP digital signature
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo