Displaying 20 results from an estimated 6000 matches similar to: "Extracting level-1 variance from lmer()"
2008 Sep 11
1
plot of all.effects object
All,
I'm trying to plot an all.effects() object, as shown in the help for
all.effects and also Crawley's R book (p.178, 2007). The data has a repeated
measures structure, but I'm using all.effects for the simple lm() fit here.
Below is a reproducible example that yields the error message.
fm.ex = lm(dv ~ time.num*drug*X, data = dat.new)
fm.effects = all.effects(fm.ex, xlevels =
2008 Sep 19
1
Type I SS and Type III SS problem
Dear all:
I m a newer on R.? I have some problem when I use?anova function.? I use anova function to get Type I SS results, but I also need to get Type III SS results.? However, in my code, there is some different between the result of Type I SS and Type III SS.? I don?t know why the ?seqe? factor disappeared in the result of Type III SS.? How can I do??
Here is my example and result.
2003 Dec 17
1
repeated measures aov problem
Hi all,
I have a strange problem and rigth now I can't figure out a
solution.
Trying to calculate an ANOVA with one between subject factor (group)
and one within (hemisphere). My dependent variable is source
localization (data). My N = 25.
My data.frame looks like this:
> ML.dist.stack
subj group hemisphere data
1 1 tin left 0.7460840
2 2 tin left
2000 Jul 05
1
Tukey.aov with split-plot designs
I am using R 1.1 with Redhat 6.2 and RW 1.001 with Win98 (the upkey doesn't
work on my IBM either as has been previously reported by others).
The function aov doesn't return either the residuals or the residual
degrees of freedom for split-plot designs.
If you use the following code from Baron and Li's "Notes on the use of R
for psycology experiments and questionnaires"
2002 Oct 09
3
proc mixed vs. lme
Dear All,
Comparing linear mixed effect models in SAS and R, I found the following
discrepancy:
SAS R
random statement random subj(program); random = ~ 1 |
Subj
-2*loglik 1420.8 1439.363
random effects
variance(Intercept) 9.6033 9.604662
2002 Apr 18
1
Help with lme basics
In Baron and Li's "Notes on the use of R for psychology experiments and questionnaires" http://cran.r-project.org/doc/contrib/rpsych.htm they describe a balanced data set for a drug experiment:
"... a test of drug treatment effect by one between-subject factor: group (two groups of 8 subjects each) and two within-subject factors: drug (2 levels) and dose (3 levels). "
2017 Dec 26
1
identifying convergence or non-convergence of mixed-effects regression model in lme4 from model output
Hi R community!
I've fitted three mixed-effects regression models to a thousand
bootstrap samples (case-resampling regression) using the lme4 package in
a custom-built for-loop. The only output I saved were the inferential
statistics for my fixed and random effects. I did not save any output
related to the performance to the machine learning algorithm used to fit
the models (REML=FALSE).
2010 Sep 16
1
ANOVA - more sophisticated contrasts
dear list,
i am using a multifactorial design with two treatments (factor A: drugs,
three levels; factor B: theraphy, two levels) and a time factor (three
levels, different timepoint). hypothetically, i measured the same subjects
for all treatements and timepoints, so its a repeated measurement design.
now i ran an anova in R and also some Tukey post-hoc tests using glht. but
what i am actually
2001 Feb 27
2
Remove columns by name data[-c("subj","drug")]
Is there an easy way to remove data frame columns
by name instead of by index? The following gives
the idea
remove<-c("subj","drug")
data[-remove]
I found a solution with a few evals and substitutes,
similar to that used in reshapeLong, but there must
be an easier way out.
Dieter
---------------------------------------
Dr. Dieter Menne
Biomed Software
72074 T?bingen
Tel
2006 Oct 08
1
Simulate p-value in lme4
Dear r-helpers,
Spencer Graves and Manual Morales proposed the following methods to
simulate p-values in lme4:
************preliminary************
require(lme4)
require(MASS)
summary(glm(y ~ lbase*trt + lage + V4, family = poisson, data =
epil), cor = FALSE)
epil2 <- epil[epil$period == 1, ]
epil2["period"] <- rep(0, 59); epil2["y"] <- epil2["base"]
2006 Feb 20
1
Extracting variance components from lmer
Hi All.
I need a bit of help extracting the residual error variance from the VarCorr
structure from lmer.
#Here's a 2-way random effects model
lmer.1 <- lmer(rating ~ (1|person)+(1|rater), data = dat)
#Get the structure
vc.fit <- VarCorr(lmer.1)
#results in.....
$person
1 x 1 Matrix of class "dpoMatrix"
(Intercept)
(Intercept) 0.7755392
$rater
1 x 1 Matrix
2009 Oct 19
1
Reposting various problems with two-way anova, lme, etc.
Hi,
I posted the message below last week, but no answers, so I'm giving it
another attempt in case somebody who would be able to help might have missed
it and it has now dropped off the end of the list of mails.
I am fairly new to R and still trying to figure out how it all works, and I
have run into a few issues. I apologize in advance if my questions are a bit
basic, I'm also no
2006 Aug 03
3
between-within anova: aov and lme
I have 2 questions on ANOVA with 1 between subjects factor and 2 within factors.
1. I am confused on how to do the analysis with aov because I have seen two examples
on the web with different solutions.
a) Jon Baron (http://www.psych.upenn.edu/~baron/rpsych/rpsych.html) does
6.8.5 Example 5: Stevens pp. 468 - 474 (one between, two within)
between: gp
within: drug, dose
aov(effect ~ gp * drug *
2005 Jun 28
1
How to extract the within group correlation structure matrix in "lme"
Dear R users,
I fitted a repeated measure model without random effects by using lme. I will use the estimates from that model as an initial estimates to do multiple imputation for missing values of the response variable in the model. I am trying to extract the within group correlation matrix or covariance matrix.
here is my code:
f = lme(y ~x0+x1+trt+tim+x1:tim +tim:trt,random=~-1|subj,
2006 Sep 04
1
Problem with Variance Components (and general glmm confusion)
Dear list,
I am having some problems with extracting Variance Components from a random-effects model:
I am running a simple random-effects model using lme:
model<-lme(y~1,random=~1|groupA/groupB)
which returns the output for the StdDev of the Random effects, and model AIC etc as expected.
Until yesterday I was using R v. 2.0, and had no problem in calling the variance components of the
2011 Jun 08
1
using stimulate(model) for parametric bootstrapping in lmer repeatabilities
Hi all,
I am currently doing a consistency analysis using an lmer model and
trying to use parametric bootstrapping for the confidence intervals.
My model is like this:
model<-lmer(y~A+B+(1|C/D)+(1|E),binomial)
where E is the individual level for consistency analysis, A-D are
other fixed and random effects that I have to control for.
Following Nakagawa and Scheilzeth I can work out the
2007 Jul 09
2
ANOVA: Does a Between-Subjects Factor belong in the Error Term?
I am executing a Repeated Measures Analysis of Variance with 1 DV (LOCOMOTOR
RESPONSE), 2 Within-Subjects Factors (AGE, ACOUSTIC CONDITION), and 1
Between-Subjects Factor (SEX).
Does anyone know whether the between-subjects factor (SEX) belongs in the
Error Term of the aov or not? And if it does belong, where in the Error Term
does it go? The 3 possible scenarios are listed below:
e.g.,
1.
2005 Oct 10
1
lmer / variance-covariance matrix random effects
Hello,
has someone written by chance a function to extract the
variance-covariance matrix from a lmer-object? I've noticed the VarCorr
function, but it gives unhandy output.
Regards,
Roel de Jong
2007 Jan 02
1
How to extract the variance componets from lme
Here is a piece of code fitting a model to a (part) of a dataset, just
for
illustration. I can extract the random interaction and the residual
variance
in group meth==1 using VarCorr, but how do I get the other residual
variance?
Is there any way to get the other variances in numerical form directly -
it
seems a litte contraintuitive to use "as.numeric" when extracting
estimates,
2012 May 01
1
VarCorr procedure from lme4
Folks
In trying to use lmer for a hierarchical model, I encountered the
following message:
Error in UseMethod("VarCorr") :
no applicable method for 'VarCorr' applied to an object of class "mer"
foo.mer <- lmer(y ~ TP + (TP|M),data=joe.q)
> head(joe.q[,1:5])
TP M AB Trt y
1 1 Jan A NN 19.20002
2 1 Jan A NN 19.06378
3 1 Jan A NN