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. How can I do that? I tried /U <- rbind(X,Y,Z) m <- manova(U~rep(1:3, c(3, 2, 3))) summary(m,test="Wilks")/ which gives / Df Wilks approx F num Df den Df Pr(>F) rep(1:3, c(3, 2, 3)) 1 0.162 12.930 2 5 0.01057 * Residuals 6 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1/ I suppose the argument rep(1:3, c(3, 2, 3)) in manova() is not appropriate. Any help is very much appreciated. Dietrich -- Dietrich Trenkler c/o Universitaet Osnabrueck Rolandstr. 8; D-49069 Osnabrueck, Germany email: [EMAIL PROTECTED] ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.