Displaying 20 results from an estimated 400 matches similar to: "adjusted p-values with TukeyHSD?"
2011 Aug 15
1
Get significant codes from a model output fit with GEE package
Does anyone know how could I get the significant codes from mixed model
output fitted with a GEE package?
The output I got is the following:
GEE: GENERALIZED LINEAR MODELS FOR DEPENDENT DATA
gee S-function, version 4.13 modified 98/01/27 (1998)
Model:
Link: Logit
Variance to Mean Relation: Binomial
Correlation Structure: Exchangeable
Call:
gee(formula = bru
2012 Jul 27
1
Understanding the intercept value in a multiple linear regression with categorical values
Hi!
I'm failing to understand the value of the intercept value in a
multiple linear regression with categorical values. Taking the
"warpbreaks" data set as an example, when I do:
> lm(breaks ~ wool, data=warpbreaks)
Call:
lm(formula = breaks ~ wool, data = warpbreaks)
Coefficients:
(Intercept) woolB
31.037 -5.778
I'm able to understand that the value of
2008 Mar 04
2
Asking, are simple effects different from 0
Hello, R-i-zens. I'm working on an data set with a factorial ANOVA
that has a significant interaction. I'm interested in seeing whether
the simple effects are different from 0, and I'm pondering how to do
this. So, I have
my.anova<-lm(response ~ trtA*trtB)
The output for which gives me a number of coefficients and whether
they are different from 0. However, I want the
2006 Mar 09
1
bugs in simtest (PR#8670)
# R for Windows will not send your bug report automatically.
# Please copy the bug report (after finishing it) to
# your favorite email program and send it to
#
# r-bugs at r-project.org
#
######################################################
This report is joint from Richard Heiberger <rmh at temple.edu>
and Burt Holland <bholland at temple.edu>.
Burt Holland is the coauthor
2004 Feb 03
1
output from multcomp and lm
Dear R-users
I analysed the same data set by two different ways;
analysis of covariance by using lm and anova functions
and multiple comparison by using simtest function in
the multcomp library.
The output from the analysis of covariance is;
> y<-lm(D~Cond+Q1,data=x)
> anova(y)
Analysis of Variance Table
Response: D
Df Sum Sq Mean Sq F value Pr(>F)
Cond 2
2005 Jul 15
1
Adjusted p-values with TukeyHSD (patch)
Dear R-developeRs,
Attached follows a patch against svn 34959 that adds the
printing of p-values to the TukeyHSD.aov function in stats package. I
also updated the corresponding documentation file and added a 'see also'
reference to the simint function of the multcomp package.
As it was already brought up in a previous thread [1] in R-help,
one can obtain the adjusted
2003 Jan 30
1
TukeyHSD and BIBD
Hi,
the function TukeyHSD gives incorrect results for balanced incomplete block
designs, as the example below shows, but I can only half fix it. There are
two problems,
1. It uses model.tables to estimate treatment means,
2. It uses the wrong standard error
The first problem can be fixed using dummy.coef, if the lines
> TukeyHSD.aov
function (x, which = seq(along = tabs), ordered = FALSE,
2012 Jun 13
1
Tukey Kramer with ANOVA (glm)
Hello,
I am performing a BACI analysis with ANOVA using the following glm:
fit1<-glm(log(Cucs_m+1)~(BA*Otter)+BA+Otter+ID+Primary, data=b1)
The summary(aov(fit1)) shows significance in the interaction; however, now I
would like to determine what combinations of BA and Otter are significantly
different (each factor has two levels). ID and PRIMARY substrates are
categorical and included in
2002 Jun 26
2
contrast matrix in package multcomp
Hi,
I've got a problem building a contrast matrix for the Dunnet contrast in
package multcopm. The following works fine:
> summary(simtest(adiff ~ trial))
Simultaneous tests: Dunnett contrasts
Data: adiff by trial
Contrast matrix:
trial1 trial2 trial3 trial4 trial5
trial2-trial1 -1 1 0 0 0
trial3-trial1 -1 0 1 0 0
2004 Jan 18
1
multcomp, simint, simtest and computation duration
Dear R-listers,
I am trying to compute simultaneous confidence intervals with simint from the package multcomp. 230 measures (abundance) have been taken in 23 sites (factor) of a data.frame (donnees: a file can be sent on request, saved with save(donnees,file="donnees")). I would like to get all pairwise comparisons with :
mc<- simint(ren~ID,type="Tukey",data=donnees)
I
2008 Oct 29
1
problem with "simtest"
Hello all
I am working with the package multcomp but I have problems with the function
simtest; the program say that can not find this function, nevertheless I doesn't have any problem with the function glht that it is in the same package.
Someone knows what could be the problem?
Thank you
[[alternative HTML version deleted]]
2003 Nov 05
2
Multiple comparisons with a glm
I've never seen anything written about multiple comparisons,
as in the multcomp package or with TukeyHSD, but using a glm.
Do such procedures exist? Are they sensible?
Are there any packages in R that implement such comparisons?
Thank you.
--
Ken Knoblauch
Inserm U371
Cerveau et Vision
18 avenue du Doyen Lepine
69675 Bron cedex
France
Tel: +33 (0)4 72 91 34 77
Fax: +33 (0)4 72 91 34 61
2007 May 21
2
more simplified output from glht object
Hi,
I use glht to make multcomp, using Tukey, from a glm model.
It is possible to get a more simplified output of result? Somethink like
ordering by letters.
Thanks
Ronaldo
--
Human kind cannot bear very much reality.
-- T. S. Eliot, "Four Quartets: Burnt Norton"
--
> Prof. Ronaldo Reis J?nior
| .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Ecologia
| : :' :
2004 Aug 13
5
simtest for Dunnett's test
Hi!
I use simtest fonction of multcomp package to compile a Dunnett's test.
I have 10 treatments and one control group, so i create a matrix with:
m<-matrix(0,10,11)
m[1,1]<--1
m[1,2]<-1
m[2,1]<--1
m[2,3]<-1
m[3,1]<--1
m[3,4]<-1
m[4,1]<--1
m[4,5]<-1
m[5,1]<--1
m[5,6]<-1
m[6,1]<--1
m[6,7]<-1
m[7,1]<--1
m[7,8]<-1
m[8,1]<--1
m[8,9]<-1
2004 May 20
4
pmvt problem in multcomp
Hi, all:
Two examples are shown below.
I want to use the multiple comparison of Dunnett.
It succeeded in upper case "example 1".
However, the lower case "example 2" went wrong.
In "example 2", the function pmvt return NaN, so I cannot show
this simtest result. Is there any solution?
(I changed the variable "maxpts" to a large number in front of
the
2011 May 01
1
Simulation Questions
I have the following script for generating a dataset. It works like a champ except for a couple of things.
1. I need the variables "itbs" and "map" to be negatively correlated with the binomial variable "lunch" (around -0.21 and -0.24, respectively). The binomial variable "lunch" needs to remain unchanged.
2. While my generated variables do come out
2011 Apr 11
1
Help on calculating a variable using random numbers
I'm new to R, but I'm trying to write a program for a dissertation that generates a dataset as follows...
subject=1:1000
treat=rbinom(1*1000,1,.13)
gender=rbinom(1*1000,1,.5)
eth=runif(1*1000, min=1, max=4)
cogat=rnorm(1*1000, 100, 16)
map=rnorm(1*1000, 200, 9)
simtest=data.frame (subject=subject, treat=treat, gender=gender, eth=round(eth,digits=0),
2005 Dec 07
1
contrasts for lm
I would like estimate a number of contrasts from a one-way ANOVA model. I see
that the lm command has a contrasts option, but I can't figure out how to use
it! Any help that can be offered would be greatly apreciated.
Here is my model statement:
Model<-lm(log2PM~P+T+P*T)
where P has 16 levels, T(treatment) has 12 levels and I am interested in
looking at different treatment
2003 May 14
1
Multiple comparison and lme (again, sorry)
Dear list,
As a reply to my recent mail:
> simint and TukeyHSD work for aov objects.
> Can someone point me to similar functions for lme objects?
Douglas Bates wrote
There aren't multiple comparison methods for lme objects because it is
not clear how to do multiple comparisons for these. I don't think the
theory of multiple comparisons extends easily to lme models. One
could
2005 May 23
0
using lme in csimtest
Hi group,
I'm trying to do a Tukey test to compare the means of a factor
("treatment") with three levels in an lme model that also contains the
factors "site" and "time":
model = response ~ treatment * (site + time)
When I enter this model in csimtest, it takes all but the main factor
"treatment" as covariables, not as factors (see below).
Is it