similar to: how to overlap plots

Hi all, I want to draw a histgram for each row of a matrix and compare them. However the plot I got does not have the same y range and x range, which makes it difficult to make the comparison. Is there a easy way to fix the x range and y range in a xy plot for several plots, instead of specifying them for each plot. The following is my code for generalizing the matrix and draw the histogram.

Hello, R experts: I got data like this: group duplicate treatment A Y 5 A Y 3 A N 6 B Y 2 B N 4 B Y 1 How to sort the data and calculate the average treatment value for each group in two level of duplicate. Results like this: group duplicate treatment A Y 4 A N

Hello, We all know linux distros will use most of the free memory available for caching stuff, I/O buffers etc (hence sites like linuxatemyram.com). In many virt scenarios we have some sort of memory deduplication mechanisms, such as KSM on KVM (don't know how it's called in Vmware). My situation is one where I need to cram up as many virtual machines as I can on a single hypervisor

Hello, We all know linux distros will use most of the free memory available for caching stuff, I/O buffers etc (hence sites like linuxatemyram.com). In many virt scenarios we have some sort of memory deduplication mechanisms, such as KSM on KVM (don't know how it's called in Vmware). My situation is one where I need to cram up as many virtual machines as I can on a single hypervisor

Hello, R experts, I tried to retrieve all biological process GO terms at level 3 starting "biological process" as level 1 using the code as bellows: 1 library(GO) 2 library(GOstats) 3 level2<-getGOChildren("GO:0008150")$"GO:0008150"$Children 4 for ( i in 1:length(level2)) { 5 level3 <- getGOChildren(level2[i])$level2[i]$Children 6 for ( j in

Dear R experts: I tried two ways to install Package:XML on windows xp for R 1.9.1, all failed. Messages are given as bellows: 1> download from CRAN > install.packages("XML", CRAN = getOption("CRAN"), + contriburl = contrib.url("http://cran.r-project.org"), + available = NULL, destdir = NULL, + installWithVers = FALSE) trying URL

Dear R-users: I am runing R 1.6.2 with Windows XP. I try to calculate Pearson correlation and Spearman correlation of any pairwise columns of 8000 x 80 data matrix with missing values and randomize the matix 1000 times and calculate this two correlations again. The code bellow for Pearson is working fine but for Spearman got the error bellow for randomized data matrix and R console is stop

Hello, I have written a function that demonstrates the CLT by generating samples following the exponential distribution, calculating the means, plotting the histogram, and drawing the limiting normal curve as an overlay. I have the title of each histogram state the sample size and rate (1/theta) for the exponential (the output is actually 4 histograms), but I can't get the greek letter theta

I am new to R. I want to graph group data using a "Traditional Bar Chart with Standard Error Bar", like the kind shown here: http://samiam.colorado.edu/~mcclella/ftep/twoGroups/twoGroupGraphs.html Is there a simple way to do this? So far, I have only figured out how to plot the bars using barplot. testdata <- scan(, list(group=0,xbar=0,se=0)) 400 0.36038 0.02154 200 0.35927

Hi, there I'm trying to plot what is returned from a call to tapply, and can't figure out how to do it. My guess is that it has something to do with the inclusion of row names when you ask for the values you're interested in, but if anyone has any ideas on how to get it to work, that would be stellar. Here's some example code: y1<-rnorm(40, 2) x1<-rep(1:2, each=20)

Has anyone written utilities to do rescaled range analysis in R? Jeff -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch

Hi Folks, I have data that looks like this: freq gender xBar 1000 m 2.32 1000 f 3.22 2000 m 4.32 2000 f 4.53 3000 m 3.21 3000 f 3.44 4000 m 4.11 4000 f 3.99 I want to plot two lines (with symbols) for the two groups "m" and "f". I have tried the following: plot(xBar[gender=="m"]~freq[gender=="f"]) followed by

Dear all, I am a newbie in R and could not find help on this problem. I am trying to plot an histogram with probabilities in the y axis. This is the code I am using: #TLC uniform n=30 mi=1; mx=6 nrep=1000 xbar=rep(0,nrep) for (i in 1:nrep) {xbar[i]=mean(runif(n,min=mi,max=mx))} hist(xbar,prob=TRUE,breaks="Sturges",xlim=c(1,6),main=paste("n =",n), xlab="Média",

Hello R experts, I'd like to studentize a matrix (tmp1) by column using apply() function and skip some specific values such as zeros in the example below to tmp2 but not tmp3. I used the script below and only can get a matrix tmp3. Could you please help me to studentize the matrix (tmp1) without changing the zeros and generate a new matrix tmp2? Thanks, Joshua tmp1 [,1] [,2] [,3] [,4]

Greets, My data looks like: > p3.18 s xbar subgroup 1 0.84 12.2 1 2 1.64 11.2 2 3 2.07 10.6 3 4 2.49 12.2 4 5 0.84 11.2 5 ... Using the command > qcc(p3.18$xbar,type="xbar",sizes=5,center=mean(p3.18$xbar),std.dev=mean(p3.18$s)/0.94,title="X-bar Chart for Paper Sheet Length Data") I get the x-bar chart I expect. However,

Hi Gents, I need to get the control limits from qcc function. As follows: qcc(MDI, type = "xbar.one") Call: qcc(data = MDI, type = "xbar.one") xbar.one chart for MDI Summary of group statistics: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.3266 0.4249 0.4371 0.4333 0.4451 0.4858 Group sample size: 1 Number of groups: 383 Center

Dear Rusers, I would like to color different sections of a histogram different colors. I have an example that was done by "brute force" given below. Has anyone implemented something like this in general? If not, any suggestions/pointers on how to write a general function to do so would be most appreciated. Alan-

I am having issues with the following: (muhat = 1/n^2(sum of all the xi's) ) essentially if xbar = the sample mean, muhat = sample mean but square the n. Question: Use R to run a Monte Carlo simulation which compares the finite-sample performance of xbar and muhat. Specifically generate 1000 samples n=30 from a standard normal distribution. For each sample calculate xbar and muhat. I have

Hi, I've had a quick look through the package list, and unless I've missed something, I can't seem to find anything that will do I-MR / Xbar-R / Xbar-S control charts ? Assuming there is something out there, can anyone point me in the right direction ? Thanks ! TIm

Hello R users, I have gene expression data of two groups of genes (large and small). Gene expression intensities of those genes are classified into 1 to 10 levels. What I want is to make a random set of genes that have the same levels as the small group from large group using sample(). I used smallvec to hold the number of genes in each levels (1 to 10) for small group, largevec for large group.