similar to: lattice shingle plot axis annotation

Hi list Maybe someone can help with the following problem (thanks in advance): In a function I have a plot and want to add symbols/text only when indicated by a logical vector (which was generated by the function before, not manually like in the following example): plot(1:10, 1:10) lv <- c(T,T,T,F,F,F,T,T,T,F) text(1:10, 1:10, '*', pos = 3) In this example the asterisk is plotted at

Hello Is there an easy way, i.e. a function in a package, to calculate the area under the curve (AUC) for drug serum levels? Thanks for any advice -- Armin Goralczyk, M.D. -- Universit?tsmedizin G?ttingen Abteilung Allgemein- und Viszeralchirurgie Rudolf-Koch-Str. 40 39099 G?ttingen -- Dept. of General Surgery University of G?ttingen G?ttingen, Germany -- http://www.chirurgie-goettingen.de

I'm interested in a version of tapply() that operates with shingles instead of factors. For instance: x <- c(1,1,2,2,3,3) y <- c(1,1,1,0,0,0) s <- shingle(x,intervals=cbind(c(0.5,1.5),c(2.5,3.5))) # the following function should exist! tapply.shingle(x,s,mean) # returns the vector c(0.75,0.25) I've written such a function as follows: tapply.shingle <-

I know that I can get a count of histogram bins in base R with plot=FALSE. However, I'd like to do the same thing with lattice. The problem is that I've set up shingles, and I'd like to get the count within each bin within each shingle. plot=FALSE doesn't seem to do it. [[alternative HTML version deleted]]

Dear all, back in 2002 Martin Henry H. Stevens wrote (https://stat.ethz.ch/pipermail/r-help/2002-May/019851.html) > How do I control the text in strips? Specifically, I want to put in the > ranges generated in shingle(x) where x is continuous. with an answer from Deepyan Sarkar (see strip.new towards the end of this message). I assume that the answer worked back then, but I've tried

How do I control the text in strips? Specifically, I want to put in the ranges generated in shingle(x) where x is continuous. Thanks, Hank Martin Henry H. Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.muohio.edu/~botcwis/bot/henry.html

Hello! Suppose I have a set of values: a <- c(1:10, 5:10) Suppose I also have a set of intervals: b <- cbind(c(0,2.5,4.5,6.5), c(5.5,7.5,9.5,11)) I can create a shingle that counts how many values are in each interval: c <- shingle(a, b) I can display the shingle to see the counts: summary(c, showValues=FALSE) The display looks like this: Intervals: min max count 1 0.0 5.5

Hi R-Users, I have a time series of bivariate observations (x,y,t). I've developed a few panel routines to explore certain aspects of the data, and now would like to break the time series up into smaller chunks to explore the temporal trends. I would like to know if anyone has any experience breaking up time series with lattice. Base graphics offers coplot, utilizing the

Dear group, I have noticed that the MacOS X binary for lme4 is not available on CRAN at the moment. I am aware that it may be possible to install from source but I am not very familiar with that procedure and would rather avoid it. But I need the package for a statistics course next week. So can anybody give an update about the situation with the MacOS X binaries? Will it be resolved within the

Hi list I tried to install package lattice on FreeBSD, but install.package() does not seem to find it in the repositories, even trying different mirrors. Could it be that lattice package is not available for FreeBSD? -- Armin Goralczyk, M.D. -- Universit?tsmedizin G?ttingen Abteilung Allgemein- und Viszeralchirurgie Rudolf-Koch-Str. 40 39099 G?ttingen -- Dept. of General Surgery University of

Hi all, By default a call to xyplot from the Lattice package when using 2 factors [eg xyplot( dv~iv | XY * AB ) ] yields the following shingle structure: |_A_|_A_|_B_|_B_| |_X_|_Y_|_X_|_Y_| However, I'm wondering if it is possible to merge the upper shingle within levels of that factor, as in: |___A___|___B___| |_X_|_Y_|_X_|_Y_| Mike -- Mike Lawrence Graduate Student, Department of

Hi: Thanks a lot to Deepayan Sarkar, author of lattice I think, who solved my first query. I am afraid that I have another one. I am plotting several mutipanels boxplots (with one conditioning variable) on one page. The x, y and conditioning variable are all continuous variables. The x and conditioning variables are transformed to shingles before being plotted. The plot looks nice but there

Hi, My problem is simple: since having updated the lattice package, I cannot load lattice anymore. If I type in the command 'library(lattice)' the loading fails with the following message: --- cut here --- Error in library.dynam(lib, package, package.lib) : shared library 'lattice' not found In addition: Warning messages: 1: In loadNamespace(package, c(which.lib.loc,

Hi all, I've got some problems when changing the trellis settings for the lattice plots. The plots look exactly as I want them to when calling show.settings() as well as when plotting them in the graphical window. But when printing to a pdf file, none of the settings are used!? Does anyone know what might have happened? Because the when changing the trellis settings, these should remain in

I am using R 2.4.0 and lattice to produce some xyplots conditioned on a factor and a shingle. The shingle merely chops up the data along the x-axis, so it is easy to identify which part of the shingle a panel is in by looking at the x-axis markings. I only want to have a strip at the top for the factor. Is this possible? I looked into calculateGridLayout() and it seems to me that there

Dear list Maybe someone can help with the following problem: I have a source file containing the following code for plotting: pdf(file = 'data/mario/boxplot.pdf') bwplot(sig100 ~ target | group, data = fish, main = 'Boxplot der Signale in 100 Zellen nach Gruppe', xlab = '', ylab = 'Anzahl der Signale') dev.off() pdf(file = 'data/mario/xyplot.pdf')

There are now several R geographic user groups, and a few have mailing lists on the R mailing list system. Thanks to Martin M, there's also a pointer to a page I'm maintaining to list/describe the groups. The page is at http://macnash.telfer.uottawa.ca/RUG.html Contact me if you have a listing. I'm prepared to wikify it if there is sufficient interest. John Nash

Hi list, I want to use the p-value from the survdiff function (package survival) to reuse within a function in a Kaplan-Meier plot. The p-value is somehow not a component of the value list ?! Thanks in advance -- A. Goralczyk G?ttingen, Ger.

Hi: Yesterday I downloaded R and got it up and running on my PC without a hitch, not so when I tried to do the same on my Mac powerbook, using Mac OSX: version 10.4.11.. I have virtually no experience, so I can't tell what exactly went wrong. At first it appeared that the software loaded correctly, it seemed like I was able to install R , and I even successfully unpacked the

Hello all I would like to visualize the hazard ratios of a cox proportional hazards model. I have seen some good examples in the New England Journal of Medicine: http://content.nejm.org/cgi/content/short/353/26/2747/F3 http://content.nejm.org/cgi/content/short/350/26/2654/F2 (I hope these are open access, but I am not sure) And something similar in the book 'R Graphics' by Paul