a.menicacci@fr.fournierpharma.com
2005-Nov-25 10:32 UTC
[R] covariance analysis by using R
Hi,
Is anyone has solved MR Xin Meng problem (see below) ?
We have the same analysis configuration : 10 groups (including control one)
with 2 mesures for each (ref at t0 and response at t1).
We expect to compare each group response with control response (group 1)
using a multiple comparison procedure (Dunnett test).
In order to perform this test, we have to normalize our data (as you) to
correct response values by the base line normalized.
Covariance analysis seems to represent the best way to do this. But how to
perform this by using R ?
So, if someone is able to deal with this problem, could you please share
with us your precious knowledge ?
Thanks in advance,
Best Regards.
Alexandre MENICACCI
Bioinformatics - FOURNIER PHARMA
50, rue de Dijon - 21121 Daix - FRANCE
a.menicacci at fr.fournierpharma.com
t??l : 03.80.44.76.17
Original message :
Hello sir:
Here's a question on covariance analysis which needs your help. There're
3
experiments,and x refers to control while y refers to experimental result.
The purpose is to compare the "y" values across the 3 experiments.
experiment_1:
x:0.1 0.2 0.3 0.4 0.5
y:0.5 0.6 0.6 0.7 0.9
experiment_2:
x:1 2 3 4 5
y:3 4 6.5 7.5 11
experiment_3:
x:10 20 30 40 50
y:18 35 75 90 98
Apparently,the control("x") isn't at the similar level so that we
can't
compare the "y" directly through ANOVA. We must normalize
"y" via "x" in
order to eliminate the influence of different level of "x". The method
of
normalize I can get is "covariance analysis",since "x" is
the covariant of
y.
My question is:
How to perform "covariance analysis" by using R? After this
normalization,we can get the according "normalized y" of every
"original
y".
All in all,the "normalized y" of every "original y" is what
I want indeed.
Thanks a lot!
My best regards!
An informal assessment may be useful: PLOT THE DATA.
x y experiment
0.1 0.5 A
0.2 0.6 A
0.3 0.6 A
0.4 0.7 A
0.5 0.9 A
1 3 B
2 4 B
3 6.5 B
4 7.5 B
5 11 B
10 18 C
20 35 C
30 75 C
40 90 C
50 98 C
Save the data as a csv file and read it into an R session:
data.df <- read.csv(file.choose())
data.df
Plot it:
plot(data.df$x, data.df$y)
Not surprisingly, the Cartesian axes obscure the behavior of the Experiment
A so we try logs:
plot(data.df$x, data.df$y, log="xy")
This plot suggests that while Experiments B and C might have a similar
relationship between x and y, Experiment A differs.
Since I know nothing of the physical meaning of these observations I am
unqualified to comment further.
Best wishes.
Charles Annis, P.E.
Charles.Annis at StatisticalEngineering.com
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-----Original Message-----
From: r-help-bounces at stat.math.ethz.ch
[mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of
a.menicacci at fr.fournierpharma.com
Sent: Friday, November 25, 2005 5:32 AM
To: r-help at stat.math.ethz.ch
Subject: Re: [R] covariance analysis by using R
Hi,
Is anyone has solved MR Xin Meng problem (see below) ?
We have the same analysis configuration : 10 groups (including control one)
with 2 mesures for each (ref at t0 and response at t1).
We expect to compare each group response with control response (group 1)
using a multiple comparison procedure (Dunnett test).
In order to perform this test, we have to normalize our data (as you) to
correct response values by the base line normalized.
Covariance analysis seems to represent the best way to do this. But how to
perform this by using R ?
So, if someone is able to deal with this problem, could you please share
with us your precious knowledge ?
Thanks in advance,
Best Regards.
Alexandre MENICACCI
Bioinformatics - FOURNIER PHARMA
50, rue de Dijon - 21121 Daix - FRANCE
a.menicacci at fr.fournierpharma.com
t??l : 03.80.44.76.17
Original message :
Hello sir:
Here's a question on covariance analysis which needs your help. There're
3
experiments,and x refers to control while y refers to experimental result.
The purpose is to compare the "y" values across the 3 experiments.
experiment_1:
x:0.1 0.2 0.3 0.4 0.5
y:0.5 0.6 0.6 0.7 0.9
experiment_2:
x:1 2 3 4 5
y:3 4 6.5 7.5 11
experiment_3:
x:10 20 30 40 50
y:18 35 75 90 98
Apparently,the control("x") isn't at the similar level so that we
can't
compare the "y" directly through ANOVA. We must normalize
"y" via "x" in
order to eliminate the influence of different level of "x". The method
of
normalize I can get is "covariance analysis",since "x" is
the covariant of
y.
My question is:
How to perform "covariance analysis" by using R? After this
normalization,we can get the according "normalized y" of every
"original
y".
All in all,the "normalized y" of every "original y" is what
I want indeed.
Thanks a lot!
My best regards!
______________________________________________
R-help at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html