similar to: queries

Hi folks, I'm trying to call an R function (fisher.test) in my program for like a billion times! Though my program is in Python and I feel that using rpy2 to interface R to python doesn't give me satisfactory performance. So I looked into R code and found out that fisher.test is actually a wrapper around another function called fexact which is implemented in C. Using Cython I managed to

Earlier, I posted the following question: I want to superimpose histograms from three populations onto the same graph, changing the shading of the bars for each population. After consulting the help files and the archives I cannot find out how to do this (seemly) simple graph. To be clear, I want - a single x axis (from -3 to 18) - three groups of bars forming the histograms of each population

I have data like: a <- rnorm(20000) b <- rep(FALSE,20000) b[sample(1:20000,15000)] <- TRUE Using Lattice graphics, I can produce two side-by-side histograms quite easily by: histogram(a | b) However, I would like to produce a "single" histogram with two bars within each bin, one for each group, as the groups are in reality very slightly different from each other. The

Hey everyone, So i cant figure this out. when using histogram() from lattice instead of hist() i get what i want as far as output. But using histogram i can seem to be able to figure out how to get multiple plots on the same panel. So par(mfrow=c(3,2)) for (i in 1:20) hist(rnorm(100),main="",cex.axis=.8) gets me about what i want but i want to use histogram() cause it gives

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

Hi, I want to do something like this: par(mfrow=c(7,7)) for (i in 1:49) hist(RATDACOM[SUBJNO==i],breaks=0.5+(0:6), main="",xlab="",ylab="",xaxt="n",yaxt="n") (Don't think about what RATDACOM and SUBJNO are.) I get an error Error in plot.new() : Figure margins too large. 36 histograms with mfrow=c(6,6) work. But the 36 histograms

Hello, After spending sometime trying to generate histograms of all the numeric variables in a data frame using a for loop, I realized I could use apply to create the histograms. This was one of those few moments when I `see' the alternative to looping. So, I can generate all the histograms with this command: apply(toblo.df[,-c(40,52)],2,hist) # 40 and 52 are ID's When I do this, the

Dear Dr. Sarkar, When I try to run the codes, I found the following problem: > h<- sample(1:14, 319, rep=T) > c<- sample(1:14, 608, rep=T) > n<- sample(1:14, 1140, rep=T) > vt<-c(h, c, n) > ta<-rep(c("h", "c", "n"), c(319, 608, 1140)) > > to<-data.frame(vt,ta) > library(lattice) Attaching package: 'lattice'

Dear R-help listers, I am new. I just want to get helps on how to plot two histograms overlapped in the same plane coordinate. What I did is very ugly. Could you please help me to improve it? I want to got a plot with semi- transparent overlapping region. And, I want to know how to specify the filled colors of the different histograms. I also prefer other solutions other than ggplot2. Many

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. --0-201145637-953833073=:7710 Content-Type: TEXT/PLAIN; charset=US-ASCII Dear all, I've been hacking on histograms the past couple of days. I wanted to do two things, add

Hi, I am having some trouble trying to plot a histogram in more than one colour. What I want to do is, plot two vectors in the same histogram, but with different colours, for instance: > x <- rnorm(1000,20,4); > y <- rnorm(1000,10,2); Then I'd like to have x and y ploted on the same hist (I can do that already doing w <- c(x,y) then hist(w)) but the bars

I would like to make a scatterplot, with a histogram of the x and y variables above and to the right . I can use layout to set up the areas, and hist(x,y) works fine for the upper histogram. However, I need a rotated histogram on the right, and I don't know how to do this. I've seen a solution with a bar plot on the right, but I'd like to use a histogram. Do you know how I can do

Hello, all; I get an unexpected result when trying to plot a probability histogram with R1.9.1 on windows xp: #with the following code: > x <- runif(100,0,1) > hist(x) > hist(x, freq=F) > h <- hist(x, freq=F) > summary(h) # Length Class Mode #breaks 11 -none- numeric #counts 10 -none- numeric #intensities 10 -none- numeric #density 10

Hi I am trying to print out means, STDs and histograms under two sets of factors. I can manage it for one set - see below but not for two sets. That is, I want ot print out the mean STD and Histogram for each ITEM code within each DELIVERABLE code. In addition I can only get to view the histogram for the last item. How do you get R to stop overriding the histogram for eg level 1 for factor 1

Dear all, I am trying to create a lattice plot which consists of 1 xyplot and 2 histograms (each for x and y). My first try was like this: x<-rnorm(1000) y<-rnorm(1000) xy <- xyplot(y~x) hist.x <- histogram(x) hist.y <- histogram(y) print(xy, position=c(0, 0.2, 1, 1), more=TRUE) print(hist.x, position=c(0, 0, 1, 0.33),more=T) print(hist.y, position=c(0.8, 0, 1, 1)) Ok, this is

HI everyone, I'm plotting a histogram in R and within that histogram i need to demonstrate the percentage of another variable (Percentage of MutStatus) within the bins plotted inthe histogram....I don't know how to do that! Data:Validation_Status Mutation_Status TvarRatio Wildtype None 0.08 Wildtype None 0.08 Wildtype None 0.08 Wildtype None 0.08 Wildtype None 0.080139373 Wildtype None

Hi, I'm trying to create a histogram with R. The problem is that the frequency is high for a couple of x-axis categories (e.g. 1500) and low for most of the x-axis categories (e.g. 50) http://r.789695.n4.nabble.com/file/n2237476/LK3_hist.jpg . When I create the histogram, it is not very informative, because only the high frequencies can be seen clearly. Is there any way I could cut the

Dear all Problem: hist()-function, scale = ?percent? I want to generate histograms for changing underlying data. In order to make them comparable, I want to fix the y-axis (vertical-axis) to, e.g., 0%, 10%, 20%, 30% as well as to fix the spaces, too. So the y-axis in each histogram should be identical. Currently, I have 100 histograms and the y-axis scales changes in each. Here is my code:

Hi, I've been passing the hist() function to tapply() to quickly generate histograms based on the list of factors supplied to tapply(). However, I have not figured out how to generate titles for each of the histograms, which paste in the unique values of the list factors as part of the histogram title. I'm hoping that someone can tell me how to do this. Thanks for your time and

Hi, I wonder if anyone has written any code to implement the suggestions of Smart et al (2004) in the Bulletin of the Ecological Society of America for a new way of graphically presenting the results of logistic regression (see www.esapubs.org/bulletin/backissues/085-3/bulletinjuly2004_2column.htm#t ools1 for the full text)? I couldn't find anything relating to this sort of graphical