similar to: Conditional ploting with logical vector

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

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

Hello (R-)Experts I hope someone can help with this problem concerning axis annotation of a lattice shingle plot. I want a plot with three shingles to display some laboratory value over time. In the first panel over the first few days, then in the next panel some months, and in the last panel some years. In the following minimal example the axis annotation will be in days, but I'd like to

Hi I constructed a mixed-effects model from longitudinal repeated measurements of lab values in 22 patients seperated into two groups with the groups as fixed effect using lme. I thought about using the jackknife procedure, i. e., removing any one subject and calculating the fixed effect, to assess the stability of the fixed effect and thereby validate the model. I suppose this has been done in

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

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')

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

Fedora Linux 8 R 2.6.1 Intel CPU I am trying to get started with R under Linux. It has been suggested to me that I run R under EMACS. I have looked at the CRAN website, and several other places to try to find an RPM for EMACS and, if I need it for ESS. I have not succeeded in finding one. Can someone tell me where to go to get the RPMs? If anyone would have a suggestion for an editor other than

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

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 all, Suppose: y<-rnorm(100) x1<-rnorm(100) lm.yx<-lm(y~x1) To predict from a new data source, one can use: # works as expected dum<-data.frame(x1=rnorm(200)) predict(lm.yx, newdata=dum) Suppose lm.yx has been run and we have the lm object. And we have a dataframe that has columns that don't correspond by name to the original regressors. I very! naively assumed that doing

Hi all, I am using R1.5.0 under Unix, I have a couple of questions here. 1. My program is running out of memory. I am writing a program to grow a list of trees using rpart() on a subset of a large dataset(5807x693) with a different response for every tree. I saw that after each tree was constucted, 116 MB of data was being added to the Vcells. I have no idea what this data is. My dataset is

Hi list, I have installed SIPp into my server. But not able to used it properly. how to configure with my server ? how to see logs on webpage ? how to start call testing .... when i start SIPp then found verious hits on myserver. *CLI:- * [Dec 27 17:37:54] NOTICE[28001]: chan_sip.c:20785 handle_request_invite: Call from '' to extension 'service' rejected because extension not

Bugzilla: https://bugzilla.redhat.com/show_bug.cgi?id=2168506 This small set covers the new /dev/mapper/VG-LV "--key" ID format in the libguestfs LUKS-on-LVM inspection test. Thanks, Laszlo Laszlo Ersek (3): update common submodule LUKS-on-LVM inspection test: rename VGs and LVs LUKS-on-LVM inspection test: test /dev/mapper/VG-LV translation common

Bugzilla: https://bugzilla.redhat.com/show_bug.cgi?id=2168506 This small set covers the new /dev/mapper/VG-LV "--key" ID format in the LUKS-on-LVM virt-inspector test. Thanks, Laszlo Laszlo Ersek (3): update common submodule inspector: rename VGs and LVs in LUKS-on-LVM test inspector: test /dev/mapper/VG-LV translation in LUKS-on-LVM test common

Hi, I'm doing some preliminary work to get Clang to compile for PS3 targets. As an intermediate step in that direction, could someone do me the favor of reviewing and submitting the enclosed patch, or giving me feedback on it? Basically, I need to be able to differentiate the triple for a PS3 target. These are the triples currently used in the gcc-based compiler from the PS3 devkit:

Hi, all. I got a lot of help from the this mailing list, it save me a lot of time. Here is my question. x<-rep(c(2,3,4),20) y<-sample(rep(c(0,1),30)) tr<-rep(c(1:5),6) data1<-data.frame(x,y,tr) data1<-data1[order(data1$tr),] > data1 x y tr1 2 0 16 4 1 111 3 0 116 2 0 121 4 0 126 3 1 131 2 0 136 4 1 141 3 0 146 2 0 151 4 0 156 3 0 12 3 0 27 2 0

http://r.789695.n4.nabble.com/file/n2260087/%E6%8D%95%E8%8E%B7.png I am going to plot my data set like this, with means and 25% & 75% quentiles. I've tried "boxplot", but the output is not what I want. Should I use other functions? Thanks -- View this message in context: http://r.789695.n4.nabble.com/ploting-dots-with-quentiles-tp2260087p2260087.html Sent from the R help

I'm trying to do something very simple... I wan to plot a histogram where the y axis represent the percentage of the total sample that each bin represents. I know how to plot a histogram with the counts and density... but can't find anything that gives me perenct of sample on the y axis. Any help is appriciated Below is the script I'm working with par(mfrow=c(1,2))