similar to: MANOVA proportion of variance explained

Hi, Suppose I have: > summary(manova(plank.man)) Df Pillai approx F num Df den Df Pr(>F) plankton.new[, 1] 1 0.5267 9.8316 6 53 2.849e-07 *** Residuals 58 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 My understanding is the MANOVA summary returns a list.

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

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 Message----- From: dvdbooth at cs.com To: kberk at

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

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

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

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

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

Hello everybody ! I try to obtain in R eigenvectors and canonical analysis on MANOVA results, but I don't find how to process? In particular, I would be interesting to obtain "standardized canonical coefficients" of the canonical variates. There analysis give some information on the correlation between response variates. My data are organised in 2 terms (one is continu, one is a

Hi I need to compare the performance of two sludge inertization methods. For that i want make a manova Wilks test. Description of the experiment: After the calcination at different temperatures my calcinated sludge are submeted to the lixiviation test. In my tables i show the concentration of the some elements in the extract phase. The results: Method A (calcination at 1100 C) Chromium

Dear All, I would be very appreciative of your help with the following 1). I am running multivariate multiple regression through the manova() function (kindly suggested by Professor Venables) and getting two different answers for test=c("Wilks","Roy","Pillai") and tests=c("Wilks","Roy",'"Pillai") as shown below. In the

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

I just finished a MANOVA test and got the following output: > summary(M, test="Pillai") Df Pillai approx F num Df den Df Pr(>F) as.factor(X) 3 1.1922 6.5948 36 360 < 2.2e-16 *** Residuals 129 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' '

Dear R fans, I have got a difficult sounding problem. For fitting a linear model using continuous response and then for re-fitting the model after excluding every single variable, the following functions can be used. library(MASS) model = lm(perf ~ syct + mmin + mmax + cach + chmin + chmax, data = cpus) dropterm(model, test = "F") But I am not sure whether any similar functions is

Hi, sorry to bother you again, but I can't figure it out myself and I also can't find any in-depth documentation about it... Consider the following SAS code (A1II2... contain the measurements for 40 subjects): proc glm; model A1II2 A1IN2 A1NI2 A1NN2 = /nouni; repeated CONTEXT 2, TARGET_SATZ 2; title "A1 500-900 ms"; This produces not only the univariate ANOVAs, but also a

Dear r-helpers, I have a query regarding use of contrasts in MANOVA. summary(manova(model)) gives me only result of test for overall difference. Would you be so kind and give me a hint how to get the same test statistics (e.g.Pillai's) and P values for the predefined contrasts? Best regards Ondrej Mikula -- Institute of Animal Physiology and Genetics Academy of Sciences of the Czech

Dear all, is there anyone who can help me extracting p-values from an Anova object from the car library? I can't seem to locate the p-values using str(result) or str(summary(result)) in the example below > A <- factor( rep(1:2,each=3) ) > B <- factor( rep(1:3,times=2) ) > idata <- data.frame(A,B) > fit <- lm( cbind(a1_b1,a1_b2,a1_b3,a2_b1,a2_b2,a2_b3) ? sex,

Dear friends, I'm stuck with the following problem: I would like to do a multiple regression with muliple dependent variables, and i don't know how to write my formula in the model. Someone in previous messages offered the cbind method, but the result is just as many regression as the number of dependent variables, it is just a time saving method; i'm looking for a method closer to a

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

I would like to test the hypothesis that the difference between pairs, for several variables, is zero. This is easily done separately for each variable with: lm(Y ~ rep(0, nrow(Y))) where Y is a matrix whose columns are the differences for each variable between pair members. However, I would like to get an overall probability across all variables from a Wilks or Pillai-Bartlett statistic as in