[R] help with multivariate analysis
spencer.graves at pdf.com
Sat Dec 17 00:48:01 CET 2005
Dear Dr. Vidali:
1. I would first plot the data a number of different ways. Maybe
you've already done plenty of this, but I have had a tendancy to rush to
the model I thought was the "best" without adequate data integrity
checks. I typically want to start with normal probability (using
"qqnorm"). Outliers, the need for transformations, mixtures, etc., will
jump off the page when I do this. Sometimes, I find that I didn't get
the data read properly. Also, I may consider log(heart.frequency) and
1/heart.frequency as well as heart.frequency. In most cases, I will
base the subsequent analysis on whichever looks most nearly normally
2. I will also want to make scatterplots of the two variables as
well as plotting each vs. time (with one line for each patient and
separate plots for the two groups, but on the same scale). This can
tell me how much correlation to expect as well as what kinds of effects
I want to model.
3. Have you consulted Pinheiro and Bates (2000) Mixed-Effects Models
in S and S-Plus (Springer)? The first half of that book describes how
to analyze data sets like you describe when effects are linear and
additive. The latter half describes what to do with nonlinear effects.
4. If you still would like further assistance from this group,
please submit another post, but PLEASE do read the posting guide!
"www.R-project.org/posting-guide.html". Anecdotal evidence suggests
that posts more consistent with that guide tend to get more useful
hope this helps.
Matteo Vidali wrote:
> dear R users,
> I need some help for multivariate analysis.
> I have 2 anaesthetic treatment groups (20 patients/group) where I
> register heart frequency and pressure for 60 min (repeated measures
> every 5 minutes). I would like to perform a test to check if treatments
> are different in controlling freq and pressures during the anaesthesia,
> but i would like to have also an overall measure and not only multiple p
> for different time intervals. I also think I should choose a test in
> which time is meaningful since the measures are not simple repeated
> measurements but measurements taken at specific time points.
> 1 million dollar question.... how to do in R?
> thanks in advance
> Dr Matteo Vidali
> Dep. of Medical Sciences
> University of East Piedmont "A. Avogadro"
> R-help at stat.math.ethz.ch mailing list
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Spencer Graves, PhD
Senior Development Engineer
PDF Solutions, Inc.
333 West San Carlos Street Suite 700
San Jose, CA 95110, USA
spencer.graves at pdf.com
More information about the R-help