similar to: Hmisc break points error

Dear Author and Maintainer of Hmisc contributed package, Dear R developers, i noticed a problem with histbackback routine when i installed the last version of R (2.4.1) (With my previous version, 2.3.1, there was no such a problem) If, e.g., you simply edit the command > histbackback(rnorm(20), rnorm(30)) a warning message is printed on the console. Warning messages: 1: argument

Hi, I'm trying to recreate a sensitivity-specificity graph using the histbackback function. The only problem is that these graphs are typically drawn with vertical rather than horizontal bar plots (and the histbackback function only seems to work with horiz=TRUE argument, using "horiz=FALSE" doesn't work). Does anyone know if: 1) there's a different graphing function that

When `length( skewed.probs ) > 200' uniform samples are generated in R-devel. R-3.5.1 behaves as expected. `epsilon` can be a lot bigger than illustrated and still the uniform distribution is produced. Chuck > set.seed(123) > > epsilon <- 1e-10 > > ## uniform to 200 then small > p200 <- prop.table( rep( c(1, epsilon), c(200, 999-200))) > ## uniform to 201

Hi all, I want to put several figures in a one figure for easy comparison, so i need to use the same methods to plot these figures. The following is an example. I also list my method, but it does not work. #Example data x<- 1:10; y<- 1:10; z<- outer( x,y,"+");z2<- outer( x,y,"-") #Quick view them image.plot(x,y,z) #relatively larger value image.plot(x,y,z2)

How do you label Xtick and Ytick marks with Currency symbols: $2000 instead of 2000? I would like to add dollar symbols to tickmarks on boxplots, histograms and back-to-back histograms. My Examples (requiring the lattice and Hmisc packages): data(case0102, package="Sleuth2") str(case0102) boxplot(Salary~Sex, case0102) histogram(~ Salary | Sex, data=case0102) require(Hmisc) #

Hi, I would like to do coloring of map regions based on the region values "weight". The approach I am taking is first to break regions into equal intervals, classIntervals(spdf$weight,4)$brks #4 intervals in this case and coloring all regions within the interval with the same color col = brewer.pal(4,"RdYlGn")) The max "weight" is as well the boundary of the

Greetings R users, My goal is to generate quartile groups of each variable in my data set. I would like each experiment to have its designated group added as a subsequent column. I can accomplish this individually with the following code: brks <- with(data_variables, cut2(var2, g=4)) #I don't want the actual numbers, I need a numbered group data$test1=factor(brks,

Question: after fitting a gamma function to some data, how do I get predicted values? I'm a SAS programmer, I new R, and am having problems getting my brain to function with the concept of "object as class ...". The following is specifics of what I am doing: I'm trying to determine the pdf from data I have created in a simulation. I have generated frequency counts

Dear all, I would like to produce a map with information about the patenting activity in German districts, by coloring districts with different degrees of patenting activity in different colors. I work with the packages maptools, maps and spdep. The map data is read from an external .shp file (+ the corresponding .shx and .dbf files). Plotting a map with the IDs or the patenting indicator itself

Basic question. For the below data, i would like to but each of the values in a bin that represents their value. So the below would hopefully put .1 in the 0-.1 bin, .2 in the .11-.2 bin and so forth. The outlying values would then be put into and outer category representing everything >1. Im using the breaks to inform some code for making a clorepleth map that represents probabilities,

I am not as expert as John, but I thought it worth pointing out that the variable substitution technique gives up one set of constraints for another (b=0 in this case). I also find that plots help me see what is going on, so here is my reproducible example (note inclusion of library calls for completeness). Note that NONE of the optimizers mentioned so far appear to be finding the true best

I've seen a number of problems like this over the years. The fact that the singular values of the Jacobian have a ration larger than the usual convergence tolerances can mean the codes stop well before the best fit. That is the "numerical analyst" view. David and Jeff have given geometric and statistical arguments. All views are useful, but it takes some time to sort them all out and

Hi all, I can not seem to get what I want using the Lattice package to generate an array of histograms of spatialgrid dataframe data. I can use the sp package and spplot to generate an array of maps that display an array of spatialgrid dataframe data -- that's good. I have:

I have two .txt files with lists of generated data I want to plot them on the same histogram. I'm aware of the histbackback function, but that is not really what I am looking for, I just want them on the same histogram but so they are still grouped separately (different colours, different norm curves etc). Does anyone know how to do this? Thank you. [[alternative HTML version deleted]]

Hello, recently I read about the SweaveListingUtils package and now I want to try it out. However, I can not make it work... Below a minimal example. The problem seems to be the following line (generated by SweaveListingPreparations()?): \ifthenelse{\boolean{Sweave@gin}}{\setkeys{Gin}{width=0.6\textwidth}}{}% If I comment out this line, it works. What can I do about this? I am using the

Hello all, After more digging I was able to find out how to do this. The answer came from an example here: https://stackoverflow.com/questions/3541713/how-to-plot-two-histograms-together-in-r yr_1997<-data.frame(bias=ann_bias$bias1997) yr_1998<-data.frame(bias=ann_bias$bias1998) yr_1999<-data.frame(bias=ann_bias$bias1999) yr_2000<-data.frame(bias=ann_bias$bias2000)

Hmisc 3.1-2 has been uploaded to the CRAN incoming directory. Change log 3.2-1 1/25/2007: Hmisc function 'ecdf' has been renamed 'Ecdf' to deconflict it with the existing 'ecdf' function in base. Fixed Bug in format.df that would create numbers with many trailing zeros. Added arguments 'math.row.names' and

Hmisc 3.1-2 has been uploaded to the CRAN incoming directory. Change log 3.2-1 1/25/2007: Hmisc function 'ecdf' has been renamed 'Ecdf' to deconflict it with the existing 'ecdf' function in base. Fixed Bug in format.df that would create numbers with many trailing zeros. Added arguments 'math.row.names' and

Glad you found an answer, though it looks more self-educational than efficient (see suggestions below). In the future, follow the recommendations of the Posting Guide: use plain text, and provide a reproducible example. Some elaborations on what "reproducible" means are [1][2][3]. One issue here was that you did not include sample data to work with (I have assumed below that ann_bias has

Hola, El problema con la propuesta de Olivier es que los intervalos son diferentes para cada variable. La forma sencilla y rápida es: spplot(zm["part88"], col.regions=plotclr, at=class$brks) Pero para que quede más elegante hay que dar algunos pasos más: ## Intervalos en forma character op <- options(digits=4) tab <- print(class) options(op) intChar <- names(tab) ## Indice