similar to: R-alpha: Minor plotting problems

In R0.61, In hist(), should the line rel.freqs <- counts/(sum(x) * diff(breaks)) computing the relative frequencies or height of the rectangles in a histogram not be rel.freqs <- counts/(sum(counts) * diff(breaks)) instead, or do I misunderstand something? Thanks, Philippe -- -------------------------------------------------------- Philippe Lambert Tel:

Hi I have recently been analyzing birthweight data from a clinic. The data has obvious defects in that there is digit preference on certain weights making them overrepresented. This shows as spikes in the histogram on certain well rounded weights like 2, 2.5, 3, etc. I would like to show this to government officials but can't figure out how I should present the finding in an easy to

Hello everybody, This is a bug involving functions in core R package: graphics::hist.default, grDevices::nclass.FD, and base::pretty.default. It is not yet on Bugzilla. I cannot submit it myself, as I do not have an account. Could somebody else add it for me, perhaps? That would be much appreciated. Kind regards, Sietse Sietse Brouwer Summary ------- Floating point errors can cause a data

> x = rnorm(100) > b = seq(min(x) - 1, max(x) + 1, length = 11) > b [1] -3.4038769 -2.7451072 -2.0863375 -1.4275678 -0.7687980 -0.1100283 [7] 0.5487414 1.2075111 1.8662808 2.5250506 3.1838203 > > invisible(hist(x, breaks = b, include.lowest = TRUE, plot = FALSE)) Warning message: argument 'include.lowest' is not made use of in: hist.default(x, breaks = b,

Hi Folks, I'd like to know how hist() decides how many cells to use when it ignores my "suggestion" to use say 'hist(...,breaks=50)'. More specifically, I have the results of 10000 simulations, each returning an 8-vector, therefore 8 variables each with 10000 values. Some of these 8 have somewhat skew distributions. Say one of these 8 variables is X. I ask for H <-

Hey everybody, I try to make a graph with two different plots. First I make a boxplot of my data. It is a collection off correlation values of different pictures. For example: 0.23445 pica 0.34456 pica 0.45663 pica 0.98822 picb 0.12223 picc 0.34443 picc etc. Ok, I make this boxplot and I get for every picture the boxes. After this I want to know, how many correlations per picture exist. So I

Dear all I have a dataset that looks like: fr.wt site 1 4400 glen 2 235 glen 3 225 glen ' ' ' ' ' ' ' ' ' 82 550 glen 83 550 kom 84 550 kom ' ' ' ' ' ' ' ' ' 191 820 kom 192 2000 soet ' ' ' ' ' ' I need to do a series of histograms for each of the codes, levels or factors in

Hi Group, I've a vector of 1000 numeric values for which I want to draw a histogram. I've read this vector into R with no variable name.I mean only the 1000 values, which makes V1 the name of the variable by default?? Then I tried > hist(V1, breaks = "Sturges", +      freq = NULL, probability = !freq, +      include.lowest = TRUE, right = TRUE, +      density = NULL, angle =

Hi, Is this a bug in truehist()? > library(MASS) > x <- rep(1, 10) > truehist(x) Error in pretty(data, nbins) : invalid 'n' value Thanks, Gad > R.version platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2 minor 4.1 year 2006 month 12 day 18 svn

Hi, Is this a bug in truehist()? > library(MASS) > x <- rep(1, 10) > truehist(x) Error in pretty(data, nbins) : invalid 'n' value Thanks, Gad > R.version platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2 minor 4.1 year 2006 month 12 day 18 svn

Dear developers, the current implementation of hist.default() calculates 'density' (and 'intensities') as dens <- counts/(n*h) where h has been calculated before as h <- diff(fuzzybreaks) which results in 'fuzzy' values for the density, see e.g. > tmp <- hist(1:10,breaks=c(-2.5,2.5,7.5,12.5),plot=FALSE) > print(tmp$density,digits=15) [1]

Full_Name: Reinhold Koch Version: 1.6.1 OS: redhat 8.0 Submission from: (NULL) (131.152.84.111) I came across rather weird behavior of the breaks in hist: hist(1:3) gives the expected result, besides an unnecessary gap between 2nd and 3rd column hist(1:4) always merges up the first two columns, also if I resort to hist.default(1:4,breaks=1:4). hist.default(1:4, include.lowest=F) gives an

Hy! I would like to know how run a montecarlo simulation with R. Thank you!!!! Francesca Matalucci __________________________________________________________________ Accesso Internet Gratis per utenti Excite! Attivalo subito! http://www.excite.it/hitech/accesso Il Mio Excite. Personalizza la tua Home page Excite come vuoi tu! http://www.excite.it AAA/Relazioni. Sfoglia gli annunci e trova la

Dear all, I have done further modifications on the histogram functions that I reported earlier this year, and I hope this can be of use and perhaps included in the distribution. I have been using this stuff a couple of months myself, and while it is nothing sophisticated, it has it's applications. :-) I did a few small modifications today to make it a bit more compact. I have modified the

I have a number of files containing anywhere from a few dozen to a few thousand integers, one per record. The statement "refdata18 = read.csv("K:\\MerchantData\\RiskModel\\Capture.Week.18.csv", header = TRUE,na.strings="")" works fine, and if I type refdata18, I get the integers displayed, one value per record (along with a record number). However, when I try "

Hello, Is there any test  for bimodality in R that x <- c(rnorm(1000,0,1),rnorm(1000,3,1)) hist(x,nclass=100) Thank you in advance for any help. Regards, Samor [[alternative HTML version deleted]]

Buenas. Sé que en R hay multitud de datasets y me haría falta alguno que trataran de variables relacionadas con salud, sobre todo para aprender más acerca de cómo realizar una regresión logística binomial o multinomial. Gracias..

Hello all, I hope someone of you can help me out, I have searched other posts as well but I can't find any solution to the problem I'm dealing with. I want to make a histogram from the data Telephone Lines MDGdataset <-read.csv("MDG_dataset_2010.csv", header=T) MDGdatasetAdapted <- subset(MDGdataset, select = c(Country_Code, Country_Name, Year,

Hi there, I want to install locfit in R for Windows NT. I have recently downloaded the file locfit.zip from http://www.ci.tuwien.ac.at/R/bin/windows/windows-NT/contrib/ but it doesn't include any explanation on how to do the installation. Have anybody done this? If so I'd really appreciate your help. Cheers, Oscar -- ------------------------------------------------------- Oscar

Here's a number of items that came up in connection with my course for medical ph.d. students this week: boxplot(): - deals ungracefully with empty groups and factor levels not present in grouping variable. - no indication of what variable is being plotted data.entry() (mostly W95 rseptbeta problems) - Entry of first variable name doesn't take before 2nd