similar to: MANOVA extension of paired t-test?

Hello there, I have two data sets with 14 variables each, and wish to do the test for equality of their covariance matrices and mean vectors. Normally these tests should be done by chi square test (box provided) and Hotelling's T square test respectively. Which R functions could do this kind of test? I just find some functions could do for one dimension, but no for multidimension. Some one

Hi everybody ... I have to perform Hotteling's T^2 Test - more generally a MANOVA - on a set of data. Is there a (simple?) possibility to do it in R? The somewhat obvious way to do it would be > summary(aov(Y ~ x1 + x2 ... + xn)) where Y would be a two-column matrix. But this does not work in any possible combination of matrix or factor dimensions! Is it principally not possible

Dear helpeRs, the following data set comes from Johnson/Wichern: Applied Multivariate Statistical Analysis, 6th ed, pp. 304-306. /X <- structure(c(9, 6, 9, 3, 2, 7), .Dim = as.integer(c(3, 2))) Y <- structure(c(0, 2, 4, 0), .Dim = as.integer(c(2, 2))) Z <- structure(c(3, 1, 2, 8, 9, 7), .Dim = as.integer(c(3, 2)))/ I would like to compute Wilk's Lambda in R, which I know is 0.0385.

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

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

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

I would like to conduct an extra sum-of -squares test that compares a full MANOVA model (with all 1st order interactions) to a reduced model (no interactions) to determine if I can drop all interactions at the same time. This is analagous to an extra sum-of-squares F-test in ANOVA, but instead using MANOVA. Is there a command in R that does this? If not, is there a command that calculates

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

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 rank" 24<501. What can it be the reason for this error

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

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.

Hi R-users I'm trying to do multivariate analysis of variance of a experiment with 3 treatments, 2 variables and 5 replicates. The procedure adopted in SAS is as follow, but I'm having difficulty in to implement the contrasts for comparison of all treatments in R. I have already read manuals and other materials about manova in R, but nothing about specific contrasts were found in them,

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

Hi, I've got a dataset with 7 variables for 8 different species. I'd like to test the null hypothesis of no difference among species for these variables. MANOVA seems like the appropriate test, but since I'm unsure of how well the data fit the assumptions of equal variance/covariance and multivariate normality, I want to use a permutation test. I've been through CRAN looking at

Dear All I would need to perform a MANOVA with both fixed (group, sex, group*sex) and random (brood) effects. I wonder if this is at all possible and if R does that. At the moment, I only know that I can run a classic MANOVA with R. Thank you David ______________________________________________ David Costantini, PhD http://www.davidcostantini.it NERC Postdoctoral research associate

Hi All, I have questions about MANOVA which I am still not sure if appropriately I should use it. For example I have a data set like this: BloodPressure (BP) Weight Height 120 115 165 125 145 198 156 99 176 I know that BloodPressure is correlated with both Weight and Height, however colinearity exists between Weight and Height. When I use BP = Weight + Height

Dear colleagues, has anyone an idea how to carry out a nonparametric manova for comparing K groups? has anyone a good package of non-parametric stats? Thank you for your help. -- José David Gómez Químico Farmacéutico Universidad Nacional [[alternative HTML version deleted]]

Hi I am using R for Windows. I was curious if there was a multivariate anova function or something like the manova in SPLUS for Windows 3.2 package? If there is, I am clearly looking in the wrong places... Thanks in advance, Mark Hall -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html

Hello everybody After doing a MANOVA on a bunch of data, I want to be able to make some comment on the amount of variation in the data that is explained by the factor of interest. I want to say this in the following way: XX% of the data is explained by A. I can acheive something like what I want by doing the following: X <- structure(c(9, 6, 9, 3, 2, 7), .Dim = as.integer(c(3,

Dear All, I am running a simulation to obtain coverage probability of Wald type confidence intervals for my parameter d in a function of two parameters (mu,d). I am optimizing it using "optim" method "L-BFGS-B" to obtain MLE. As, I want to invert the Hessian matrix to get Standard errors of the two parameter estimates. However, my Hessian matrix at times becomes