similar to: Is this a bug (or a feature) in hist(x)$density ??

Hi all, I'm looking to add a "density" smoother on top of a hist when Freq=T. In order to do this I can use the relation between count and density, but I would like to know if there is a way for me to predict it upfront. Here is an example: set.seed(242) z = rnorm(30) hist_z <- hist(z) hist_z$counts / hist_z$density # the relation is 15 # why is this 15 ?? # So I can now do:

I have the following codes but can anyone make it shorter i.e making these FOR loop into one loop ... thanks... par(mfrow=c(2,4)) for(i in 16:23){ hist(data[,i],main=paste(colnames(data)[i],sep=""),ylab="Frequency",xlim=c(1,5),xlab="Score",ylim=c(0,100)) } png("histogram.png") dev.off() par(mfrow=c(2,4)) for(i in

Hi, people. Executing the following command: hist(rpois(100,5), labels=TRUE) yields a graphic in which some labels are truncated (on an X11 device). The truncated labels are those over the highest bars. The hist function should ideally manage enough room for the labels, automatically. (Specifying ylim solves my problem, but yet, hist could be frienlier.) --please do not edit the

Hello dear R help group, I would like to download the tRNA data on: http://gtrnadb.ucsc.edu/download.html And then import it into R. Can anyone direct me as to how to do so? Thanks, Tal ---------------------------------------------- My contact information: Tal Galili E-mail: Tal.Galili@gmail.com Phone number: 972-52-7275845 FaceBook: Tal Galili My Blogs: http://www.talgalili.com (Web and

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]

Hello all, # I am trying to read the text in this URL: u <- http://google.com/complete/search?output=toolbar&q=%d7%a9%d7%9c%d7%95%d7%9d # By using this command: readLines(u) And no matter what variation I tried, I keep getting this output: [1] "<?xml version=\"1.0\"?><toplevel><CompleteSuggestion><suggestion

After reading more, I understand I didn't formulate my last question correctly, so please allow me to rephrase: What I am looking for is a way to save the R console session output. That is, a command that would combine the results of using: ?sink # And ?savehistory My motivation for this is that doing it will allow someone who is a blind user of R to be able to easily export his results to

The following code is giving me problems. I want to export densities of a distribution to a csv file. At the bottom of the code I use the hist function to generate the densities. But hist is returning values greater than 1. I don't understand, why. Any help you can supply is greatly appreciated. # Set word path dir<-"~/Research/MR Distribution Analysis/" setwd(dir)

Hello dear R-help group. I wish to plot a scatter plot using icons (or images) instead of points. Is it possible? and how so? Thanks, Tal ---------------------------------------------- My contact information: Tal Galili E-mail: Tal.Galili@gmail.com Phone number: 972-52-7275845 FaceBook: Tal Galili My Blogs: http://www.talgalili.com (Web and general, Hebrew) http://www.biostatistics.co.il

My question if in the Subject, but if to extend: I am specifically curious about WordPress blogs. But any solution will give me a lead. Thanks, Tal ---------------------------------------------- Contact me: Tal.Galili@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com/ (English) [[alternative HTML version deleted]]

Hello all, I would like some R function that can translate a string to a "URL encoding" (see here: http://www.w3schools.com/tags/ref_urlencode.asp) Is it implemented? (I wasn't able to find any reference to it) Thanks, Tal ----------------Contact Details:------------------------------------------------------- Contact me: Tal.Galili@gmail.com | 972-52-7275845 Read me:

Hello dear R-help mailing list. I wish to be able to have htmlParse work well with Hebrew, but it keeps to scramble the Hebrew text in pages I feed into it. For example: # why can't I parse the Hebrew correctly? library(RCurl) library(XML) u = "http://humus101.com/?p=2737" a = getURL(u) a # Here - the hebrew is fine. a2 <- htmlParse(a) a2 # Here it is a mess... None of

Hello dear R help group, I am trying to read a .txt file, with Hebrew column names, while keeping the column names looking well in R - but without success. I uploaded an example file to: http://www.talgalili.com/files/aa.txt And am trying the command: read.table("http://www.talgalili.com/files/aa.txt", header = T, sep = "\t") This returns me with: X.....ª X...ª......

Hello all, I am using R 2.13.0 with windows 7, after giving my user full privileges to the R folder (as described here<http://www.r-statistics.com/2011/04/how-to-upgrade-r-on-windows-7/> ). This allows me to install new packages just fine. However, when using update.packages(), to update existing packages, I keep getting the following error (for example, when updating the MASS package):

The R poisson random generator rpois appears to have a bug for theta 10 or larger. The sample mean of the pseudo variates is too small: sample mean approx theta - 0.5. I use Version 1.1.1 (August 15, 2000) Of R on a Dell OptiPlex computer with the Windows XP Professional operating system. Has this bug been fixed in later versions? (I found another reported rpois bug, but it appears to be

Hello all, I recently discovered the "comment" command. I see it can only hold a vector of characters. Is there a way (or an alternative), to make it possible to have it keep a list? (for example, to keep different pieces of information like date of creation, information of each variable and so on) The closest solution I can think of is using 'names' on the vector, like this:

Hi all, I'm trying to superimpose a rchisq density line over a histogram with something like: hist(alnlength) lines(density(rchisq(length(alnlength), 4)),col="red") But the rchisq line won't appear anywhere, Anyone knows what I am missing here? Thanks in advance, Albert.

Hi, I'm completely new to R, I am all at sea with the interface and the confusing help files, so would appreciate some help to do a simple task. Need to present the mean and variance of 100 different samples of poisson distributions (N=1000, with fixed lambda) in a file in two columnns, and then produce histograms. So far I have figured out: > N <- 1000 > x <- rpois(N, 3.1) ,

> x <- rnorm(1000) > h <- hist(x,plot=FALSE) > sum(h$density) [1] 2 ----------------------------- shouldn't it be 1?! > h <- hist(x,plot=FALSE, breaks=(-4:4)) > sum(h$density) [1] 1 ----------------------------- now it's 1. why?! -- Sam Steingold (http://sds.podval.org/) on Ubuntu 11.10 (oneiric) X 11.0.11004000 http://www.childpsy.net/ http://www.memritv.org

Hi all, I just failed in loading a saved wordspace (13MB of size), and received these errors: Warning: stack imbalance in 'missing', 52 then 51 Warning: stack imbalance in 'if', 50 then 53 Warning: stack imbalance in 'as.environment', 57 then 59 Warning: stack imbalance in 'ls', 54 then 53 Warning: stack imbalance in '.Internal', 54 then 53 Warning: stack