Displaying 4 results from an estimated 4 matches for "esticon".
2009 Oct 19
1
Defining S3-methods for S4-objects: cannot coerce type 'S4' to vector of type 'integer'
In the 'doBy' package there is an esticon() function for calculating linear contrasts for various model types. I have defined an S3-method 'esticon.mer()' for 'mer' objects from the lme4 package. Building the package and invoking the method gives:
> esticon(fm1, c(1,1))
Confidence interval ( WALD ) level = 0.95
Error in...
2008 Apr 16
2
Post hoc tests with lme
...Med Q3 Max
-1.80200345 -0.64316591 0.05783115 0.70099706 1.63142054
Number of Observations: 36
Number of Groups: 9
========
Now suppose I want to test the following hypothesis
H0: TypeT1 =0 and TypeT2 = 0
I've tried estimable() and glh.test() in package gmodels, esticon() in
package boBy, and linear.hypothesis() in package car, but it seems
none of them would work with objects from lme:
========
> library(gmodels)
> estimable(fm, rbind(c("TypeT1"=1), c("TypeT2"=1)))
Error in FUN(newX[, i], ...) :
`param' has no names and does not m...
2007 Jul 25
2
using contrasts on matrix regressions (using gmodels, perhaps)
Hi,
I want to test for a contrast from a regression where I am regressing the columns of a matrix. In short, the following.
X <- matrix(rnorm(50),10,5)
Y <- matrix(rnorm(50),10,5)
lm(Y~X)
Call:
lm(formula = Y ~ X)
Coefficients:
[,1] [,2] [,3] [,4] [,5]
(Intercept) 0.3350 -0.1989 -0.1932 0.7528 0.0727
X1 0.2007 -0.8505 0.0520
2008 May 12
1
Standard error of combination of parameter estimates
Hi,
Is there a simple command for computing the standard error of a
combination of parameter estimates in a GLM?
For example:
riskdata$age1 <- riskdata$age
riskdata$age2 <- ifelse(riskdata$age<67,0,riskdata$age-67)
model <- glm(death~age1+age2+ldl,
data=riskdata,family=binomial(link=logit))
And we want to find the standard error of the partial linear predictor
for