Displaying 20 results from an estimated 1100 matches similar to: "ANOVA and confidence intervals plot"
2007 Mar 27
0
Error when calling residual.plots() on an ANOVA object
Dear *,
I have to perform an ANOVA analysis to study the effects on the
performance of some algorithmic components in the design of my
algorithm.
I have done a full factorial design for the following 5 factors (my
response variable is "best"):
$ model :'data.frame': 3360 obs. of 6 variables:
..$ best : num [1:3360] 0.108 0.573 0.625 1.057 0.451 ...
2007 Feb 09
1
Help in using multcomp.
Hi All,
I am trying use 'multcomp' for multiple comparisons
after my ANOVA analysis. I have used the following
code to do ANOVA:
dat <- matrix(rnorm(45), nrow=5, ncol=9)
f <- gl(3,3,9, label=c("C", "Tl", "T2"))
aof <- function(x) {
m <- data.frame(f, x);
aov(x ~ f, m)
}
amod <- apply(dat,1,aof)
Now, how can I use
2006 Aug 10
1
How to speed up nested for loop computations
Dear all,
here is the result of R.Version():
> R.Version()
$platform
[1] "powerpc-apple-darwin8.6.0"
$arch
[1] "powerpc"
$os
[1] "darwin8.6.0"
$system
[1] "powerpc, darwin8.6.0"
$status
[1] ""
$major
[1] "2"
$minor
[1] "3.1"
$year
[1] "2006"
$month
[1] "06"
$day
[1] "01"
$`svn rev`
2007 Mar 18
1
multcomp
I used the multcomp package sometime back for doing multiple
comparisons. I see that it has been updated and the methods like simint
are no longer supported. When I run the program it prompts to me to use
glht. How do I get the lower and upper conf int and the pValues using
glht? Does anyone have an example?
Thanks ../Murli
[[alternative HTML version deleted]]
2009 Nov 05
1
Newbie question Multcomp
Hello, I'm a totally newbie to R and I'm taking a class using S+.
In the class we use the multcomp command which takes a aov object and
calculates confidence intervals for all pairwise differences by the Fisher
least significant differences method.
How can I do this in R.
Thank you for taking the time with such a basic question. I've been looking
on the net for a few days and I
2013 Oct 12
1
export glht to LaTeX
Hi,
I want to export the result of glht in R into a LaTeX table, such as that result:
Linear Hypotheses:
Estimate Std. Error z value Pr(>|z|)
Group1 - Group2 == 0 -0.14007 0.01589 -8.813 <0.001 "***"
Group1 - Group3 == 0 -0.09396 0.01575 -5.965 <0.001 ***
---
Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05
2007 Nov 21
1
multiple comparison (glht) problem
I am not sure whether there is a bug. When I tested the example given for "glht"
in the help, I entered the following error:
Running commands:
amod <- aov(minutes ~ blanket, data = recovery)
rht <- glht(amod, linfct = mcp(blanket = "Dunnett"),
alternative = "less")
Errors are:
Error in try(coef.(model)) : could not find function
2007 Nov 21
1
question about multiple comparison in ANOVA
I am not sure whether there is a bug. When I tested the example given
for "glht" in the help, I entered the following error:
Running commands:
amod <- aov(minutes ~ blanket, data = recovery)
rht <- glht(amod, linfct = mcp(blanket = "Dunnett"),
alternative = "less")
Errors are:
Error in try(coef.(model)) : could not find function
2012 Jan 02
1
Is using glht with "Tukey" for lme post-hoc comparisons an appropriate substitute to TukeyHSD?
Hello,
I am trying to determine the most appropriate way to run post-hoc
comparisons on my lme model. I had originally planned to use Tukey
HSD method as I am interested in all possible comparisons between my
treatment levels. TukeyHSD, however, does not work with lme. The
only other code that I was able to find, and which also seems to be
widely used, is glht specified with Tukey:
2008 May 28
2
Tukey HSD (or other post hoc tests) following repeated measures ANOVA
Hi everyone,
I am fairly new to R, and I am aware that others have had this
problem before, but I have failed to solve the problem from previous
replies I found in the archives.
As this is such a standard procedure in psychological science, there
must be an elegant solution to this...I think.
I would much appreciate a solution that even I could understand... ;-)
Now, I want to calculate a
2008 Apr 10
1
Tukey in R, extracting values
hey,
how can i extract the values from the CI's when i use following code for a
tukey test?
the output shows three CI's and i know it should work with 'names but i
don't really know how...
thanks
library(multcomp)
data1$soort<-as.factor(data1$soort)
amod<-aov(waarde~soort,data=data1)
g<-glht(amod, linfct=mcp(soort = "Tukey"))
confint(g)
--
View this message in
2006 Oct 24
1
Posthoc tests for 3-way ANOVA analysis
Hi All,
I have performed a 3-way ANOVA analysis for my
experimental data using aov function. My simple R
funtion for this is:
3aof <- function(x){
m <- data.frame(R,S,T, x);
anova(aov(x ~ R+S+T+R*S+R*T+S*T+R*S*T, m) )
}
Now, I am getting P values for all the main and
interactions effects. If I want to perform postdoc
test on one of my main effects, say T, what method I
should use (i have
2005 Nov 25
1
Plotting the diff. between two curves
Dear Rs
I have two vectors A and B
where
A is
V1 d
1 0.000100000 1.123278
2 0.002186431 1.120448
3 0.004351214 1.106661
4 0.006515998 1.107713
5 0.008680781 1.107667
6 0.013010348 1.106353
7 0.019504698 1.104077
8 0.034658181 1.103202
9 0.051976447 1.103200
10 0.073624280 1.094825
11 0.093085682 1.085123
12 0.095250465 1.087325
13 0.132051782 1.086158
14 0.168853098
2011 Mar 01
1
glht() used with coxph()
Hi, I am experimenting with using glht() from multcomp package together with
coxph(), and glad to find that glht() can work on coph object, for example:
> (fit<-coxph(Surv(stop, status>0)~treatment,bladder1))
coxph(formula = Surv(stop, status > 0) ~ treatment, data = bladder1)
coef exp(coef) se(coef) z p
treatmentpyridoxine -0.063 0.939 0.161
2011 Mar 04
2
glht: Problem with symbolic contrast for factors with number-levels
Using a factor with 'number' levels the straightforward
symbolic formulation of a contrast in 'glht' of
the 'multcomp' package fails.
How can this problem be resolved without having to redefine the factor levels?
Example:
#A is a factor with 'number' levels
#B similar factor with 'letter' levels
dat<-data.frame(y=1:4,A=factor(c(1,1,2,2)),
2010 May 17
1
Query on linear mixed model
Hi R Forum
I am a newbie to R and I have been amazed by what
I can get my team to accomplish just by
implementing Scripting routines of R in all my
team's areas of interest..
Recently i have been trying to adopt R scripting
routine for some analysis with longitudanal data..
I am presenting my R script below that I have
tried to make to automate data analysis for
longitudanal data by employing
2013 Jul 25
1
lme (weights) and glht
Dear R members,
I tried to fit an lme model and to use the glht function of multcomp.
However, the glht function gives me some errors when using
weights=varPower().
The glht error makes sense as glht needs factor levels and the model
works fine without weights=. Does anyone know a solution so I do not
have to change the lme model?
Thanks
Sibylle
--> works fine
2011 Jul 26
3
a question about glht function
Hi all:
There's a question about glht function.
My data:data_ori,which inclue CD4, GROUP,time.
f_GROUP<-factor(data_ori$GROUP)
f_GROUP is a factor of 3 levels(0,1,2,3)
result <- lme(sqrt(CD4) ~ f_GROUP*time ,random = ~time|ID,data=data_ori)
glht(result, linfct = mcp(f_GROUP="Tukey") )
Error in `[.data.frame`(mf, nhypo[checknm]) : undefined columns selected
I can't
2009 Mar 22
1
Multiple Comparisons for (multicomp - glht) for glm negative binomial (glm.nb)
Hi
I have some experimental data where I have counts of the number of
insects collected to different trap types rotated through 5 different
location (variable -location), 4 different chemical attractants [A, B,
C, D] were applied to the traps (variable - semio) and all were
trialled at two different CO2 release rates [1, 2] (variable CO2) I also
have a selection of continuous variables
2011 Aug 06
1
multcomp::glht() doesn't work for an incomplete factorial using aov()?
Hi R users,
I sent a message yesterday about NA in model estimates (
http://r.789695.n4.nabble.com/How-set-lm-to-don-t-return-NA-in-summary-td3722587.html).
If I use aov() instead of lm() I get no NA in model estimates and I use
gmodels::estimable() without problems. Ok!
Now I'm performing a lot of contrasts and I need correcting for
multiplicity. So, I can use multcomp::glht() for this.