Displaying 20 results from an estimated 400 matches similar to: "multcomp and glm"
2006 Jul 11
2
Multiple tests on 2 way-ANOVA
Dear r-helpers,
I have a question about multiple testing.
Here an example that puzzles me:
All matrixes and contrast vectors are presented in treatment contrasts.
1. example:
library(multcomp)
n<-60; sigma<-20
# n = sample size per group
# sigma standard deviation of the residuals
cov1 <- matrix(c(3/4,-1/2,-1/2,-1/2,1,0,-1/2,0,1), nrow = 3, ncol=3, byrow=TRUE,
dimnames =
2005 Mar 09
1
multiple comparisons for lme using multcomp
Dear R-help list,
I would like to perform multiple comparisons for lme. Can you report to me
if my way to is correct or not? Please, note that I am not nor a
statistician nor a mathematician, so, some understandings are sometimes
quite hard for me. According to the previous helps on the topic in R-help
list May 2003 (please, see Torsten Hothorn advices) and books such as
Venables &
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
2003 May 05
1
multcomp and lme
I suppose that multcomp in R and multicomp in S-Plus are related and it
appears that it is possible to use multicomp with lme in S-Plus given the
following correspondence on s-news
sally.rodriguez at philips.com 12:57 p.m. 24/04/03 -0400 7 [S] LME summary
and multicomp.default()
Is it possible to use multicomp with lme in R and if so what is the syntax
from a simple readily available
2006 Jul 25
1
Multiple tests on repeated measurements
Dear R-helpers:
My question is how do I efficient and valid correct for multiple tests in a repeated measurement design:
Suppose we measure at two distinct visits with repeated subjects a treatment difference on the same variable.
The treatment differences are assessed with a mixed model and adjusted by two methods for multiple tests:
# 1. Method: Adjustment with library(multcomp)
2006 Feb 07
1
post-hoc comparisons following glmm
Dear R community,
I performed a generalized linear mixed model using glmmPQL (MASS
library) to analyse my data i.e : y is the response with a poisson
distribution, t and Trait are the independent variables which are
continuous and categorical (3 categories C, M and F) respectively, ind
is the random variable.
mydata<-glmmPQL(y~t+Trait,random=~1|ind,family=poisson,data=tab)
Do you think it
2003 May 08
0
multcomp and lme (followup)
I just realized that in the call to `csimint' the argument `asympt=TRUE'
is missing since we need to compute the confidence intervals for a glm
based on the normal approximation.
Torsten
---------------------------------------------------------------------
library(multcomp)
set.seed(290875)
# a factor at three levels
group <- factor(c(rep(1,10), rep(2, 10), rep(3,10)))
# Williams
2004 Feb 20
1
nlme and multiple comparisons
This is only partly a question about R, as I am not quite sure about the
underlying statistical theory either.
I have fitted a non-linear mixed-effects model with nlme. In the fixed
part of the model I have a factor with three levels as explanatory
variable. I would like to use Tukey HSD or a similar test to test for
differences between these three levels.
I have two grouping factors:
2009 Dec 04
2
csimtest function in multcomp package
Hello all,
Quick question: I want to do posthoc contrasts for a linear mixed
effects model. However, when trying to use the csimtest function in
multcomp package I receive an error message saying it cannot find the
function, even after installing and loading package multcomp.
Any pointers would be greatly appreciated
Daniel
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
| : :' :
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
2010 Jun 23
1
Shapefile
Hopefully the attachment will make it this time...
Hi:
I am practicing with the attached shapefile and was wondering
if I can get some help. Haven't used 'rgdal' and 'maptools' much
but it appears to be a great way bring map data into R.
Please take a look at the comments and let me know if I need to
explain better what I am trying to accomplish.
library(rgdal)
2010 Jun 23
1
Plotting Data on a Map
Hi:
I am practicing with the attached shapefile and was wondering
if I can get some help. Haven't used 'rgdal' and 'maptools' much
but it appears to be a great way bring map data into R.
Please take a look at the comments and let me know if I need to
explain better what I am trying to accomplish.
library(rgdal)
library(maptools)
library(ggplot2)
dsn="C:/Documents and
2002 Jun 20
1
new package `multcomp'
New package `multcomp' for general multiple comparisons
written by Frank Bretz, Torsten Hothorn and Peter Westfall
We've uploaded the package `multcomp' to CRAN.
The R package allows for multiple comparisons of
k groups in general linear models. We use the unifying
representations of multiple contrast tests, which include all
common multiple comparison procedures, such as the
2002 Jun 20
1
new package `multcomp'
New package `multcomp' for general multiple comparisons
written by Frank Bretz, Torsten Hothorn and Peter Westfall
We've uploaded the package `multcomp' to CRAN.
The R package allows for multiple comparisons of
k groups in general linear models. We use the unifying
representations of multiple contrast tests, which include all
common multiple comparison procedures, such as the
2002 Oct 14
1
Post hoc Multiple comparison
Dear R-listers
I'm a new R-user who needs some help with a test that I want to do. I
have done a field experiment: four treatments (cont, x, y and xy) at
three sites (A, B and C), the response is count data (0 - 15). I've done
a Poisson regression:
>glm(response~as.factor(treatment)*as.factor(site), family=quasipoisson,
offset(max.response), data=dat)
The "offset" is the
2010 Jun 24
0
rgdal-maptools
The shapefile data can be downloaded from the link below:
download all the six files and save them on your working directory
and make sure the dsn path is set to where the files are saved.
My shapefiles are saved on C:/Data.
https://secure.filesanywhere.com/fs/v.aspx?v=897263875a6472a99baa
Hi:
I am practicing with the attached shapefile and was wondering if I can?
get some help. Haven't
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
2011 Feb 20
2
Problem using F77_CALL(dgemm) in a package
Dear R-devel,
I've written a numerical solver for SOCPs (second order cone programs)
in R, and now I want to move most of the solver code into C for speed.
I've written combined R/C packages before, but in this case I need to
do matrix operations in my C code. As I have never done that before,
I'm trying to write some simple examples to make sure I understand the
basics. I am stuck
2012 Apr 16
0
warning message: coxme with package multcomp
Hi
I'm encountering an error/warning when doing multiple comparisons with
the package multcomp on a coxme model.
My data:
I'm looking at the removal of brood from the nest according to three
treatments I applied on the brood.
The brood and the workers caring about the brood in the nest, belonged
to different colonies.
Factor: treatment (3 levels: tx,uv,meta)
Random effect 1: origin of