Displaying 20 results from an estimated 10000 matches similar to: "lme v. aov?"
2003 Apr 02
2
replication of latin squares --- again
Dear all, this is a newbie's question.
I have a 4x4 latin square replicated 3 times. That is:
operators
batches 1 2 3 4
1 A B C D
2 B C D A
3 C D A B
4 D A B C
operators
batches 1 2 3 4
5 A B C D
6 B C D A
7 C D A B
8 D A B C
operators
batches 1 2 3 4
9 A B C D
2003 Nov 16
1
SE of ANOVA (aov) with repeated measures and a bewtween-subject factor
Hallo!
I have data of the following design:
NSubj were measured at Baseline (visit 1) and at 3
following time points (visit 2, visit 3, visit 4).
There is or is not a treatment.
Most interesting is the question if there is a
difference in treatment between the results of visit 4
and baseline. (The other time points are also of
interest.) The level of significance is alpha=0.0179
(because of an
2008 Aug 05
1
Mixed model with multiple response variables?
Hi,
I have a data set collected from 10 measurements (response variables)
on two groups (healthy and patient) of subjects performing 4 different
tasks. In other words there are two fixed factors (group and task),
and 10 response variables. I could analyze the data with aov() or
lme() in package nlme for each response variable separately, but since
most likely there are correlations among the 10
2003 May 23
1
variance components
Dear All,
I need to calculate the variance components in a mixed effect model (one
fixed and one random effect) with REML (maximizing the proportion of
the likelihood that does not depend on the fixed effects). In S+ there is
the varcomp function, but I would like to do it in R. Is there a way to do
that?
Thanks!
Katalin
___
Katalin Csillery
Division of Biological Sciences
University of
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 *
2006 Aug 14
2
lme() F-values disagree with aov()
I have used lme() on data from a between-within subjects experiment. The correct
ANOVA table is known because this is a textbook example (Experimental Design by
Roger Kirk Chapter 12: Split-Plot Factorial Design). The lme() F-values differ from
the known results. Please help me understand why.
d<-read.table("kirkspf2.dat",header=TRUE)
for(j in 1:4) d[,j] <- factor(d[,j]) ### Make
2006 Oct 06
1
Once again: aov vs. lme/lmer
First of all, I apologize for asking a question that has appeared
recurrently in this mailing list. However, I have googled for it, have
looked at the mailing list archives, and also looked at Pinheiro & Bates book
(although not very thoroughly, I must confess), to no avail.
Here is the question: I am trying to obtain with lme or lmer the same exact
numerical results (p-values) that I obtain
2008 Aug 25
1
aov, lme, multcomp
I am doing an analysis and would like to use lme() and the multcomp
package to do multiple comparisons. My design is a within subjects
design with three crossed fixed factors (every participant sees every
combination of three fixed factors A,B,C). Of course, I can use aov() to
analyze this with an error term (leaving out the obvious bits):
y ~ A*B*C+Error(Subject/(A*B*C))
I'd also like
2006 Aug 10
5
Variance Components in R
Hi,
I'm trying to fit a model using variance components in R, but if very
new on it, so I'm asking for your help.
I have imported the SPSS database onto R, but I don't know how to
convert the commands... the SPSS commands I'm trying to convert are:
VARCOMP
RATING BY CHAIN SECTOR RESP ASPECT ITEM
/RANDOM = CHAIN SECTOR RESP ASPECT ITEM
/METHOD = MINQUE (1)
/DESIGN
2005 Oct 27
2
Extracting Variance Components
Dear List,
Is there a way to extract variance components from lmeObjects or
summary.lme objects without using intervals()? For my purposes I don't
need the confidence intervals which I'm obtaining using parametric
bootstrap.
Thanks,
Mike
[[alternative HTML version deleted]]
2003 Jul 16
1
Q re Linear Models
Hello,
I apologize in advance if what I'm about to ask is trivial or has been
answered before. In the latter case I would appreciate a pointer to the
right list/location
I'm trying to model the following experimental design with groupedData
and lme in R:
Subjects were measured on two tasks (continuous outcome variable is
"Prob", tasks are coded as within-subjects factor
2012 Jan 30
1
Linear Mixed Model set-up
Hello,
I have some data covering contaminant concentrations in fish over a time
period of ~35 years. Each year, multiple samples of fish were taken (with
varying sample sizes each year). Ultimately, I want an estimation of the
variance between years, and the variance within years + random effects. I
used a linear mixed model to estimate these variances, but after reading a
number of different
2001 Dec 23
1
aov for mixed model (fixed and random)?
I'm starting to understand fixed and random effects, but I'm
puzzled a bit. Here is an example from Hays's textbook (which is
great at explaining fixed vs. random effects, at least to dummies
like me), from the section on mixed models. You need
library(nlme) in order to run it.
------
task <- gl(3,2,36) # Three tasks, a fixed effect.
subj <- gl(6,6,36) # Six subjects, a random
2003 Oct 27
2
variance component analysis for nested model
Given a set of data:
> names(data)
[1] "city" "house" "visit" "value"
I am looking for a way to compute the variance components of the
nested model (ie, visit 1 at house 2 at city 3 isn't related to visit
1 and house 2 at city 4), but different houses in the same city may be
related, and different visits to the same house are probably
2006 Apr 25
5
Heteroskedasticity in Tobit models
Hello,
I've had no luck finding an R package that has the ability to estimate a
Tobit model allowing for heteroskedasticity (multiplicative, for example).
Am I missing something in survReg? Is there another package that I'm
unaware of? Is there an add-on package that will test for
heteroskedasticity?
Thanks for your help.
Cheers,
Alan Spearot
--
Alan Spearot
Department of Economics
2008 Mar 15
1
Anova
Hi all,
I apologize for what might be a silly question.
I am interested in doing a one way anova.
This is not too hard in and of itself, either with anova, aov or oneway.test
.
However, I need to
1) get pvalues,
2) do a posthoc analysis with Tukey HSD,
3) and have (sometimes) an unbalanced design.
I just can't seem to put all the pieces together.
Any suggestions?
Thanks in advance,
Dan.
2017 Aug 08
2
how to extract individual values from varcomp?
Hello,
I am trying to use varcomp to decompose the variance across multiple
nested levels on a lme object. I am able to successfully do this and
when I view the varcomp object I can see the individual values /
estimates for the variance at different levels.
However, I want to be able to extract each of them separately, as I
need to build a confidence interval using bootstrapping on the sample
2003 Nov 19
2
repeated measure in GLM
I was recently asked to perform a GLM analysis (the person comes from
the JMP world) on a repeated measures design. I have found some things
using aov but I cannot find anything with glm. In fact, multiple
regressions in general with repeated measures seems to be poorly
covered in documentation. I remember SPSS has separate commands to
handle them.
I have within variables SOA and ec and a
2005 Nov 30
3
Random Effects for One-Way Anova
Hello to All.
I'd want to use a one-way ANOVA. This means that I have only one factor, with,
lets say, 5 levels.
I made a dataframe, called "DATA", with two Columns:
A, that is my response, and it is "class numerical".
B, that defines the different levels of my factor, and it is "class factor".
If I want to use a fixed effect model,
I know that the formula I have
2003 Dec 15
2
help in lme
To anyone who can help,
I have two stupid questions, and one fairly intelligent question
Stupid question (1): is there an R function to calculate a factorial of a number? That is...is there a function g(.) such that g(3) = 6, g(4) = 24, g(6) = 720, etc?
Stupid question (2): how do you extract the estimated covariance matrix of the random effects in an lme object?
Intelligent question