Displaying 20 results from an estimated 3000 matches similar to: "Manova for repeated measures"
2003 Nov 10
8
Memory issues..
Hi dear R-listers, I'm trying to fit a 3-level model using lme in R. My
sample size is about 2965 and 3 factors:
year (5 levels), ssize (4 levels), condition (2 levels).
When I issue the following command:
>
lme(var~year*ssize*condition,random=~ssize+condition|subject,data=smp,method
="ML")
I got the following error:
Error in logLik.lmeStructInt(lmeSt, lmePars) :
2004 Jan 14
3
summary() within a function
I have the following function
f <- function {
...
model <- lm(rttx[,1] ~ rttx[,2] + 0);
summary(model);
...
}
while summary(model) shows the summary if I execute the function line by line
in the Command Line Interface, if I call f() summary is silent
how to solve it? or is there workaround?
thanks in advance
Ulisses
Debian GNU/Linux: a dream come true
2004 Feb 11
6
AGREP
Hi all, I have two questions
1 - I have the version 1.4.1 of R, and it doesn't have the 'agrep'
function in the base library. Is there a way to make this funcion
avaliable in R 1.4.1? I mean, how to 'copy' it from R 1.8.1 and 'paste'
it in R 1.4.1?
2 - The AGREP function doesn't give me the Levenshtein distance (edit
distance). Is there a function in R that does
2004 Feb 23
2
Error in multiple xyplots
Dear R-listers,
I got an error when I try to plot two grouped data into a single
win.metafile device:
library(lattice)
trellis.device(device="win.metafile",color=F,filename="Profiles-Var1.wmf")
par(mfrow=c(1,2))
# First plot
xyplot(log(v)~t|id,data=Con.20,
main="Group A: Control, Var-1",
xlab="ocasion",ylab="var1",
panel=function(x,y) {
2006 Nov 09
2
Repeated Measures MANOVA in R
Can R do a repeated measures MANOVA and tell what dimensionality the statistical variance occupies?
I have been using MATLAB and SPSS to do my statistics. MATLAB can do ANOVAs and MANOVAs. When it performs a MANOVA, it returns a
parameter d that estimates the dimensionality in which the means lie. It also returns a vector of p-values, where each p_n tests
the null hypothesis that the mean
2010 Oct 29
1
Repeated Measures MANOVA
Hello all,
Is there an r function that exists that will perform repeated measures MANOVAs? For example, let's say I have 3 DVs, one between-subjects IV, and one within-subjects IV. Based on the documentation for the manova command, a function like that below is not appropriate because it cannot take Error arguments.
manova(cbind(DV1,DV2,DV3) ~ BetweenSubjectsIV * WithinSubjectsIV +
2007 Apr 05
1
MANOVA with repeated measurements
Hello,
I have a question regarding performing manova. I have an experiment where I want to measure 10 output variables for 3 different measurement methods. Since each of the methods requires some user interaction, I would also like to include repeated measures for each of the output variables to include intraobserver variability in the design.
How can I perform such a repeated-measures
2009 Aug 11
1
nested repeated measures MANOVA using adonis
I am trying to apply a permuation-based MANOVA (Anderson 2001) to a set of
morphological data from three ecomorphs of fish reared under two different
conditions and measured at two points during ontogeny. I will supply a
distance matrix based on Procrustes distances calculated outside of vegan.
I have not found an example of a design such as this for adonis. However, I
have designed my factors
2011 May 04
1
Str info. Thanks for helping
It looks from str(SA) that Response IPS1 is a data.frame of class "anova", which probably cannot be coerced to vector.
Maybe you can use unlist() instead of as.vector()
Or something like
SA[["Response IPS1"]]["as.factor(WSD)",] ## to select the first row only, even maybe with unlist()
Without a better REPRODUCIBLE example, I cannot tell more (maybe some others
2006 Feb 16
2
MANOVA: how do I read off within and between Sum-of-Squares info from the manova result?
Hi all,
I am experimenting the function "manova" in R.
I tried it on a few data sets, but I did not understand the result:
I used "summary(manova_result)"
and "summary(manova_result, test='Wilks')"
and they gave a bunch of numbers...
But I need the Sum-of-Squares of BETWEEN and WITHIN matrices...
How do I read off from the R's manova results?
Any
2004 May 24
2
Manova and specifying the model
Hi,
I would like to conduct a MANOVA. I know that there 's the manova() funciton and the summary.manova() function to get the appropriate summary of test statistics.
I just don't manage to specify my model in the manova() call. How to specify a model with multiple responses and one explanatory factor?
If I type:
2007 Feb 22
1
MANOVA usage
Hello,
I had a couple questions about manova modeling in R.
I have calculated a manova model, and generated a summary.manova output
using both the Wilks test and Pillai test.
The output is essentially the same, except that the Wilks lambda = 1 -
Pillai. Is this normal? (The output from both is appended below.)
My other question is about the use of MANOVA. If I have one variable which
has a
2003 Nov 22
3
summary.manova and rank deficiency
Hi all,
I have received the following error from summary.manova:
Error in summary.manova(manova.test, test = "Pillai") :
residuals have rank 36 < 64
The data is simulated data for 64 variables. The design is a 2*2 factorial with 10 replicates per treatment. Looking at the code for summary.manova, the error involves a problem with qr(). Does anyone have a suggestion as to how to
2001 Nov 08
2
Manova in R vs. SAS
While I was helping a SAS-using friend with an analysis I noticed some
differences in the multivariate test statistics, approximate F statistics,
and p-values in the manova function using R and proc GLM using SAS. The
univariate coefficients are identical. Is there a reason to expect R and SAS
to give different results?
Thanks,
Bill Kristan.
2004 Feb 15
1
manova() with a data frame
I'm trying to learn to use manova(), and don't understand why none of
the following work:
> data(iris)
> fit <- manova(~ Species, data=iris)
Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
incompatible dimensions
> fit <- manova(iris[,1:4] ~ Species, data=iris)
Error in model.frame(formula, rownames, variables, varnames, extras,
2006 Mar 30
2
Unbalanced Manova
Dear all,
I need to do a Manova but I have an unbalanced design. I have
morphological measurements similar to the iris dataset, but I don't have
the same number of measurements for all species. Does anyone know a
procedure to do Manova with this kind of input in R?
Thank you very much,
Naiara.
--------------------------------------------
Naiara S. Pinto
Ecology, Evolution and Behavior
1
2008 Jul 15
2
extracting elements from print object of Manova()
Hi there,
Does anyone know how to extract elements from the table returned by Manova()?
Using the univariate equivalent, Anova(), it's easy:
a.an<-Anova(lm(y~x1*x2))
a.an$F
This will return a vector of the F-values in order of the terms of the model.
However, a similar application using Manova():
m.an<-Manova(lm(Y~x1~x2))
m.an$F
Returns NULL. So does any attempt at calling the
2012 Aug 25
2
Standard deviation from MANOVA??
Hi,
I have problem getting the standard deviation from the manova output.
I have used the manova function: myfit <- manova(cbind(y1, y2) ~ x1
+ x2 + x3, data=mydata) .
I tried to get the predicted values and their standard deviation by using:
predict(myfit, type="response", se.fit=TRUE)
But the problem is that I don't get the standard deviation values, I only
2012 Mar 19
1
car/MANOVA question
Dear colleagues,
I had a question wrt the car package. How do I evaluate whether a
simpler multivariate regression model is adequate?
For instance, I do the following:
ami <- read.table(file =
"http://www.public.iastate.edu/~maitra/stat501/datasets/amitriptyline.dat",
col.names=c("TCAD", "drug", "gender", "antidepressant","PR",
2003 Nov 20
4
p value in MANOVA
Dear R users,
Can anyone tell me how to get the p value out of the output from
summary.manova?
I tried all the methods I can think of, but failed.
Many thanks
Yu-Kang
_________________________________________________________________
¥ß§Y¥Ó½Ð MSN Mobile ªA°È¡G¦b±zªº¤â¾÷¤W¦¬µo MSN Hotmail
http://msn.com.tw/msnmobile