search for: manova

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: pcor.manova&...

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 th...

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 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-va...

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...

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 an...

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 a...

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...

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. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tu...

Dear All, I am trying to perform MANOVA. I have table with 504 columns(species) and 36 rows) with two grouping (season and location) Zx <- Z[c(4:504)] Zxm <- as.matrix(Z) m<- manova(Zxm~Season*location, data=Z) when I do summary.aov, I get respond for each species but summary.manova summary.manova(m) :" residuals have ra...

...ender", "antidepressant","PR", "dBP", "QRS")) ami$gender <- as.factor(ami$gender) ami$TCAD <- ami$TCAD/1000 ami$drug <- ami$drug/1000 library(car) fit.lm <- lm(cbind(TCAD, drug) ~ gender + antidepressant + PR + dBP + QRS, data = ami) fit.manova <- Manova(fit.lm) fit1.lm <- update(fit.lm, .~ . - PR - dBP - QRS) fit1.manova <- Manova(fit1.lm) Is there an easy way to find out whether the reduced model is adequate? I am thinking of something similar to the anova() function, I guess? Many thanks and best wishes, Ranjan

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. -----------------------------...

Hi, Does anyone know what the error message mean? > Boot2.Pillai <- function(x, ind) { + x <- as.matrix(x[,2:ncol(x)]) + boot.x <- as.factor(x[ind, 1]) + boot.man <- manova(x ~ boot.x) + summary(manova(boot.man))[[4]][[3]] + } > > man.res <- manova(as.matrix(pl.nosite) ~ + as.factor(plankton.new[,1]))$residuals > boot2.plank <- cbind(plankton.new[, 1], man.res) > boot.sep <- boot(boot2.plank, Boot2.Pillai, R = 1000, +...

Dear friends, Sorry for this somewhat generically titled posting but I had a question with using contrasts in a manova context. So here is my question: Suppose I am interested in doing inference on \beta in the case of the model given by: Y = X %*% \beta + e where Y is a n x p matrix of observations, X is a n x m design matrix, \beta is m x p matrix of parameters, and e is a normally-distributed random matrix wi...

Dear R-users; Previously I posted a question about the problem of rank deficiency in summary.manova. As somebody suggested, I'm attaching a small part of the data set. #*************************************************** "test" <- structure(.Data = list(structure(.Data = c(rep(1,3),rep(2,18),rep(3,10)), levels = c("1", "2", "3"), class = "fac...

Hello all, I have a question on manova in R: I'm using the function "manova()" from the stats package. Is there anything like a stepwise (backward or forward) manova in R (like there is for regression and anova). When I enter: step(Model1, data=Mydata) R returns the message: Error in drop1.mlm(fit, scope$drop, scale...

Hi.? There appears to be a bug in R function manova.? My friend and I both ran it the same way as shown below (his run) with the shown data set. His results are shown below. we both got the same results.? I was running with R 2.3.1. I'm not sure what version he used. Thanks very much, David Booth Kent State University -----Original Mess...

Hi, I have a nested unbalanced data set of four correlated variables. When I do univariate analyses, my factor of interest is significant or marginally significant with all of the variables. Small effect size but always in the same direction. If I do a MANOVA instead (because the variables are not independent!) then my factor is far from being significant. How does that come about? I have found a mention of a so-called Rao's paradox, which seems to deal with exactly this phenomenon. Does anyone know more about it, e.g. a reference? The next st...

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...

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