search for: pillai

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 first case (test=c(list)) I got error message which probably means I can only call one test at a time. I thought I could get ride of this by adding "s" to test; in this case (tests=c(list)...

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

...xtable method for the Anova function in CAR, or has anybody used a different set of functions to import Anova results into a table in an Sweave document? If not, any handy hints on how to write a good homebrew based on the output of Anova would be much appreciated - in particular, where can I find Pillai test statistics and the formula for computing approximate F that is used in the CAR package? I have Fox' "R and S-Plus Companion to Regression". I'm aware that this is a very basic question - thank you for your patience! Maria Wolters -- Dr Maria Wolters CSTR / HCRC, Sch...

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

...ure what version he used. Thanks very much, David Booth Kent State University -----Original Message----- From: dvdbooth at cs.com To: kberk at ilstu.edu Sent: Sun, 15 Mar 2009 7:01 pm Subject: Re: MANOVA Data Ken, Did you notice that Wilks, Roy, etc p-values are all the same?? Pillai is almost the SAS result.? Can't figure it out.? I'll submit a bug report. What's Velleman going to talk about?? Thanks for looking at the R. Best, Dave -----Original Message----- From: Ken Berk <kberk at ilstu.edu> To: dvdbooth at cs.com Sent: Sun, 15 Mar 20...

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 higher F-stat value (~60) than the MANOVA fstat (~20), but the other 15 v...

...ter.39 1 Sum of squares and products for the hypothesis: (Intercept) (Intercept) 0.6765062 Sum of squares and products for error: (Intercept) (Intercept) 13.05917 Multivariate Tests: Treatment Df test stat approx F num Df den Df Pr(>F) Pillai 1 0.0492517 0.6216377 1 12 0.44574 Wilks 1 0.9507483 0.6216377 1 12 0.44574 Hotelling-Lawley 1 0.0518031 0.6216377 1 12 0.44574 Roy 1 0.0518031 0.6216377 1 12 0.44574 ------------------------------------------ Term: Hormone...

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,

Hi, I wanted to know whether libvirt-1.2.10 supports HyperV Server 2012 and/or SCVMM Server 2012. There is a requirement for us to deploy a VHD image to the HyperV Server 2012 SP2. We are planning to do this using libvirt-1.2.10. Thanks Paddy

...: LBT'= 0 where B is the matrix of regression coefficients; L is a matrix with rows giving linear combinations of the regressions coefficients; the transformation matrix T has the same meaning as in the anova.mlm function. An example and some bare-bones code is listed below (code depends on the Pillai, Wilks etc. functions defined in src/library/stats/R/mlm.R). Example model: 3 dependents, 1 between-subjects factor with 4 levels set.seed(123) Y <- cbind(rnorm(100), rnorm(100), rnorm(100)) A <- factor(rep(c(1,2,3,4), each=25)) fit <- lm(Y ~ A) Example 1: simple contrast: compare leve...

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

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 deal with this error? T...

...n 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 number of multivariate MANOVAs, including some test statistics (Wilks' Lambda, Pillai's Trace, etc.) on various hypothesis (no CONTEXT effect, no TARGET_SATZ effect, no CONTEXT*TARGET_SATZ effect). Unfortunately, I can't seem to be able to figure out how to reproduce this in R. There's obviously little sense in treating A1II2..A1NN2 as separate responses, so I came u...

...&CP Matrix for TRA E = Error SS&CP Matrix S=2 M=-0.5 N=4.5 Statistic Value F Num DF Den DF Pr > F Wilks' Lambda 0.02042803 32.9813 4 22 0.0001 Pillai's Trace 1.24157026 9.8222 4 24 0.0001 Hotelling-Lawley Trace 35.12690920 87.8173 4 20 0.0001 Roy's Greatest Root 34.75791615 208.5475 2 12 0.0001 ... Manova Test Criteria and Exact F Sta...

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,

Hi, On a few of the hpc cluster nodes, i am seeing a new lustre error that is pasted below. The volumes are working fine and there is nothing on the oss and mds to report. LustreError: 5080:0:(import.c:607:ptlrpc_connect_interpret()) data3-OST0000_UUID at 192.168.2.98@tcp changed handle from 0xfe51139158c64fae to 0xfe511392a35878b3; copying, but this may foreshadow disaster

My system: Suse 9.3 Prof. I downloaded Wine: wine-0.9.7-SuSELinux100.i586.rpm. I get the message libstdc++.so.6 not available. What can I do? Please note: I am not a linux expert. What little I know, I acquired by myself. Regards, suren -- The woods are lovely dark and deep But I have promises to keep And miles to go before I sleep And miles to go before I sleep (Robert Frost)

Hi, In a routed network environment, without the router , we want to use the shorewall as the firewall/router. The ISP has assigned the following set of IP addresses. WAN IP for subnet 1 (DATA) 220.227.202.X/30 ( to be assigned to eth0 of the shorewall) WAN IP for subnet 2 (Voice) 220.227.202.Y/30 ( to be assigned to eth1 of the shorewall) Addresses assigned for Subnet 1 by

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, I am running samba and lprng for 4 hp laserjets. Sometimes a samba print job goes to the wrong printer. Normally everything works just fine. Shown below is the lprng acct log file enries for the concerned print job. Please notice that the print job was supposed to go to hp4050-mail, but ended up in hp8100. Any suggestions are much appreciated. Thanks Balagopal jobstart