Displaying 20 results from an estimated 1300 matches similar to: "lme for repeated measurements over time"
2007 Dec 13
2
Function for AUC?
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
2007 Oct 04
2
plot from source file
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')
2008 Apr 05
2
Conditional ploting with logical vector
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
2008 Jan 22
1
Lattice on FreeBSD
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
2008 May 04
1
Validating a mixed-effects model
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
2009 Jul 30
2
lattice shingle plot axis annotation
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
2007 Apr 26
2
Extract p-value from survdiff function
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.
2004 Feb 13
5
predict function
I am using R to do a loess normalisation procedure.
In 1.5.1 I used the following commands to normalise the variable "logratio",
over a 2d surface (defined by coordinates x and y):
> array <- read.table("121203B_QCnew.txt", header=T, sep="\t")
> array$logs555<-log(array$s555)/log(2)
> array$logs647<-log(array$s647)/log(2)
>
2012 Sep 07
2
metafor package: study level variation
Hello. A quick question about incorporating variation due to study in the metafor package. I'm working with a particular data set for meta-analysis where some studies have multiple measurements. Others do not. So, let's say the effect I'm looking at is response to two different kinds of drug treatment - let's call their effect sizes T1 and T2. Some studies have multiple
2008 Jan 02
1
WG: AW: Another problem with encoding
Hello, Peter,
I tried it out: iconv(names(attributes(spss[,'Y6'])[[1]][14]), "UTF-8", "LATIN1", sub='byte') yielded
[1] "<c4>rzte Chirurgie"
and c4 corresponds in most encodings to ?. What can I do next? I wonder whether there is a more comfortable way then to change the
occurences of <..> by the adequate character.
Regards,
Matthias
2008 Jan 09
2
labels to values
I couldn't find out how the following to solve:
There is a column in a data.frame:
> d[, 'Y6']
[1] 6 3 8 11 8 9 6 8 3 5 10 15 NA 9 8 3 8 16 6 6 NA 10 5 2 7 7 6 16 7 15 7 10 12 8 7 12 12 16 7 6 8 8
15 6 [45] NA 8 99 7 12 8 9 16 7 16 8 7 7 1 15 12 8 7 10 7 8 7 8 9 8 6 6 8 6 16 11 5 11 11 1 11 3 7 7 10
10 10 6 11 [89] 16 NA
2002 Sep 09
1
Re: Biplot function of PCA
[was sent a wrong R-help address; manually resent by MM]
Hello
I'am using the 'biplot' and 'biplot.pincomp' functions of the 'mva'
package for my studies.
The biplot represents both the observations and the variables of a
matrix of multivariate data on the same plot.
The observations are represented by their numbers (the line of the data
matrix), but I would need to
2007 Sep 13
2
Multivariate, multilevel regression?
Dear WizaRds,
This is mostly a statistics question, but I'm figuring that R is the right solution (even before I start!)
I have some bio data of heart rate versus time (rats taken from resting to maximal heart rate). I want to regress heart rate on time. The data have been normalized such that resting heart rate is zero at time=0, so that all curves intersect at the origin (and at the origin
2012 Nov 06
3
Survplot, Y-axis in percent
Hi
I am a new fan of R after getting mad with the graphical functional in SPSS.
I have been able to create a nice looking Kaplan Meyer graph using Survplot
function.
However I have difficulties in turning the y axis to percent instead of the
default 0-1 scale.
Further I have tried the function yaxt="n" without any results. Any help in
this matter will be appreciated.
The code is
2011 Nov 09
2
Error in drawing
I have got following error in drawing wavelet fitting. can some one help?
> library(faraway)
> data(lidar)
> newlidar<-lidar[c(1:128),]
> library(wavethresh)
> wds <- wd(newlidar$logratio)
> draw(wds)
Error in plot.default(x = x, y = zwr, main = main, sub = sub, xlab = xlab,
:
formal argument "type" matched by multiple actual arguments
[[alternative HTML
2007 Oct 09
1
Visualize cox proportional hazards
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
2010 Jun 01
2
MacOS X binary for lme4 not available on CRAN
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
2011 Jul 13
3
Colors in R
HI everyone,
I''m trying to assign colors to multiple lines in a graph. Problem is I don''t want to type in as many colors as there are lines....is there a way around this? In brief, I''m plotting the logratio for up to 60 samples and want a different color for each sample. Here is the code I''m using now..
Any help is greatly appreciated..
Best
LT
data <-
2008 Nov 11
0
read.table not clearing colClasses
I am attempting to created colClasses for several tables, then read only specific columns. There are two different table layouts that I am working with. If I use exclusively one layout, the script works perfectly, but if I mix the layouts, it fails to grab the correct columns form layout that is read in second. It appears that colClasses fails to adopt the new structure after the first iteration.
2011 Oct 06
0
TCGA expression data: plotting ....
Hi,
I am new to R.
I am trying to figure out how to graph expression data from the TCGA
database.
If I understand correctly the expression data I have downloaded is from a
microarray using the AgilentG4502A.
I've had trouble reading into R in the level I, level II, and the gene
expression analysis data using
>dat<-read.table("C:\\file.txt", header=T, row.names=1)
for