similar to: axis and tick widths decoupled (especially in rugs!)

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]]

Hi I could only find some discussion on this wrt lattice graphics (which I'm not using). Apologies if I'm missing something obvious. I'd like to produce 3 rug plots under a kernel density plot for a population. The population is subdivided into 3 subpopulations, which I'd like the rug plots to highlight. Naturally, when I do 3 rug plots, they all plot over each other.

Hi all, Is there an equivalent to predict(...,type="linear") of a Proportional hazard model for a Cox model instead? For example, the Figure 13.12 in MASS (p384) is produced by: (aids.ps <- survreg(Surv(survtime + 0.9, status) ~ state + T.categ + pspline(age, df=6), data = Aidsp)) zz <- predict(aids.ps, data.frame(state = factor(rep("NSW", 83), levels =

Hi all, I have a (hopefully quick) question. I've got a fascinating set of fitted surfaces in three dimensions corresponding to local linear multiple regressions. I'd like to add rugs to the X and Y axes (corresponding to my independent variables) in order to get a sense for how many data points I'm working with to graph various portions of the surfaces. The trouble is, I can't

Hi, I have downloaded openssh(http://www.openssh.com/) on my Windows host. Just wanted to know if openssh allows executing a batch script on Windows host. Batch script contains mapping of shares which is required for my scripts. My command run from hpux: $ ssh <windows server> myscript The myscript requires shares to be mapped on login via

There must be a very basic thing I am not getting... I'm working with some entropy functions and the convention is to use log(0)=0. So I wrote a function: llog<-function(x){ if (x ==0) 0 else log(x) } which seems to work fine for individual numbers e.g. >llog(0/2) [1] 0 but if I try whole vectors or tables: p<-c(4,3,1,0) q<-c(2,2,2,2) llog(p/q) I get this: [1] 0.6931472

OK, this is really weird! here's an example code: t1<-c(1,2,3,4) t2<-c(4,2,4,2) plot(t1~t2, xlab="exp1", ylab="exp2") dev.copy2eps(file="test.eps") that all seems fine... until you look at the eps file created, where for some weird reason, if you scroll down to the end, the code reads: /Font1 findfont 12 s 0 setgray 214.02 18.72 (e) 0 ta -0.360 (xp1) tb

Hi! This is going to be a real newbie question, but I can't figure it out. I'm trying to plot densities of various functions of chi-square. A simple chi-square plot I can do with dchisq(). But e.g. chi.sq/degrees of freedom I only know how to do using density(rchisq()/df). For example: plot(1, type="n", xlab="", ylab="", xlim=c(0,2), ylim=c(0,7)) for (i

Hello, I have a simple question on rug(). Currently there is only one color possible for the rug. Is it possible to plot a the rug with different colors, for each rug item ? Thx. Bjoern

Hi All, I'm trying to add a 'rug' representation of my data to a plot created with densityplot(). While I can do this in the simple case, I can't do it properly when I include the "groups" argument. I have an example below. I am running a reasonably new version of R. print(sessionInfo()) R version 2.12.0 Patched (2010-11-07 r53537) Platform: i686-pc-linux-gnu

Dear R list members, I'd like to produce rug plots at the bottom of panels in a trellis display (using the lattice package), but par("usr") doesn't return user coordinates for panels, and consequently rug fails, as the following example (suggested to me by Georges Monette) illustrates: > x <- rnorm(50) > y <- rnorm(50) > f <-

Hi there! I'm not sure I can create a minimal example of my problem, so I'm linking to a minimal .RData file that has only two objects: obs and exp, each is a 6x9 matrix. http://dl.dropbox.com/u/10364753/test.RData link to dropbox file (I hope this is acceptable mailing list etiquette!) Here's what happens: > obs[1, 1] [1] 118 > exp[1, 1] [1] 118 > obs[1, 1]-exp[1, 1] [1]

rug() may add integer axis labels when called with a POSIXlt object as argument. dtimes <- c("09/29/2007 12:54", "09/30/2007 00:14", "10/01/2007 00:14", "10/02/2007 00:14", "10/03/2007 00:14", "10/04/2007 00:14", "10/05/2007 00:14", "10/06/2007 00:14", "10/07/2007

Hello, I'm trying to plot a historgram with a density plot superimposed: hist(x, seq(-1, 1, by = 0.1), prob = T, col = "blue") lines(density(x, bw = 0.1)) rug(x) I don't want to add rug(x) but the histogram will not be plotted unless rug(x) is there. It does not work if the "hist" line alone is present. Can you help? Thanks, Kim [[alternative HTML version

To wit: > plot(1:2,1:2,log="x") > rug(1:2) # should not warn Warning message: some values will be clipped in: rug(1:2) > plot(1:2,1:2,log="x") > rug(c(0,.1,.2)) # should warn, but does not > I believe this fixes the problem: % diff -c /tmp/orig.rug.R /tmp/cberry.rug.R *** /tmp/orig.rug.R Thu Nov 29 12:59:00 2001 --- /tmp/cberry.rug.R Thu Nov 29

I am fairly new to log-linear modelling, so as opposed to trying to fit modells, I am still trying to figure out how it actually works - hence I am looking at the interpretation of parameters. Now it seems most people skip this part and go directly to measuring model fit, so I am finding very few references to actual parameters, and am of course clear on the fact that their choice is irelevant for

Hi! I'm writing my own project in LLVM and I came across a document describing Decoupled software pipelining algorithm written in LLVM by The Liberty Research Group for their VELOCITY C++ compiler: http://liberty.princeton.edu/Research/DSWP/ http://liberty.princeton.edu/Publications/micro38_dswp.pdf and some of spectacular results: http://liberty.princeton.edu/videos/gimp.php

I'm trying to import a table into R the file is about 700MB. Here's my first try: > DD<-read.table("01uklicsam-20070301.dat",header=TRUE) Error: cannot allocate vector of size 15.6 Mb In addition: Warning messages: 1: In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : Reached total allocation of 1535Mb: see help(memory.size) 2: In scan(file, what,

Hello, Can anyone help me out with this problem? After performing logistic regressions and testing the significance with likelihood ratios, I have plotted the results using "termplot". I am wondering, how to get the names of my variables to appear on the x-axis rather than ascending numbers? I have used: termplot(fm, se=T, axes=FALSE, col.se="black") axis(1, 1:4,

Hello dear R community, Recently, the wonderful people<http://www.warwick.ac.uk/statsdept/useR-2011/index.html#About>behind useR!2011 made sure to nudge the speakers to send in their slides from the conference, and updated the conference's website with links to the slides (and videos)<http://www.warwick.ac.uk/statsdept/user-2011/schedule/index.html>of the talks. Since the links