search for: sspe

...B2 3 A1 B3 4 A1 B4 5 A2 B1 6 A2 B2 7 A2 B3 8 A2 B4 9 A3 B1 10 A3 B2 11 A3 B3 12 A3 B4 > attach(poke.idata) > pokeAnova =Anova(multmodel,idata=poke.idata,idesign=~Afac*Bfac,type="III") Error in linear.hypothesis.mlm(mod, hyp.matrix, SSPE = SSPE, idata = idata, : The error SSP matrix is apparently of deficient rank = 4 < 6 Thanks for any help or advice. And thanks for the 'car' package, which is a great asset to my course. I'm just stuck on this one example. colleen moore http://www.childrenandpollution.org...

...products for error: [,1] [,2] [,3] [,4] [1,] 4.4957632 1.084781 2.196134 -0.7621832 [2,] 1.0847809 173.494183 20.425963 -7.7303319 [3,] 2.1961336 20.425963 8.584237 -3.8182268 [4,] -0.7621832 -7.730332 -3.818227 2.3045286 1-pf(diag(linhyp2$SSPH)/diag(linhyp2$SSPE),1,79) 1-pt(sqrt(diag(linhyp2$SSPH)/diag(linhyp2$SSPE)),79) ## The last two lines of code are what I am stuck with; I thought using diag() on the SSPH and SSPE matrices should give me the sums of squares for hypothesis and error; so division should yield the F values? I would be most grateful fo...

I'm trying to determine the lack of fit for regression on the following: data <- data.frame(ref=c(0,50,100,0,50,100), actual=c(.01,50.9,100.2,.02,49.9,100.1), level=gl(3,1)) fit <- lm(actual~ref,data) fit.aov <- aov(actual~ref+Error(level),data) According to the information I have, the lack of fit for this regression is the

...ated with a particular set of species and their abundances. I have 20 rows (samples) split into 6 categories, and 17 columns (species). I am getting the following error message, which I don't understand: > can<-candisc(mod, data=canIN) Error in linear.hypothesis.mlm(mod, hyp.matrix.2, SSPE = SSPE, V = V, : The error SSP matrix is apparently of deficient rank = 16 < 17 Does anyone has any experience in candisc? -- Tanya. [[alternative HTML version deleted]]

Hi, I have encountered this error when attempting a One-way Repeated-measure ANOVA with my data. I have read the "Anova in car: SSPE apparently deficient rank" thread by I'm not sure the within-subject interaction has more degrees of freedom than subjects in my case. I have prepared the following testing script: rm(list = ls()) suppressMessages(require("plyr")) suppressMessages(require(ez))...

Hello all, I am getting the following error: Error in linear.hypothesis.mlm(mod, hyp.matrix.1, SSPE = SSPE, V = V, : The error SSP matrix is apparently of deficient rank = 7 < 11 After running: mod.ok <- lm(as.matrix(dat[,-1]) ~ DC, data=dat) (av.ok <- Anova(mod.ok, idata=idata, idesign=~week)) Although if I jitter the data in "dat", the function seems to work. What shou...

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

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,

...1:B2" .. .. ..$ : chr [1:2] "A1:B1" "A1:B2" ..$ sex:A:B : num [1:2, 1:2] 26.4 2.3 2.3 0.2 .. ..- attr(*, "dimnames")=List of 2 .. .. ..$ : chr [1:2] "A1:B1" "A1:B2" .. .. ..$ : chr [1:2] "A1:B1" "A1:B2" $ SSPE :List of 8 ..$ (Intercept): num [1, 1] 34459 .. ..- attr(*, "dimnames")=List of 2 .. .. ..$ : chr "(Intercept)" .. .. ..$ : chr "(Intercept)" ..$ sex : num [1, 1] 34459 .. ..- attr(*, "dimnames")=List of 2 .. .. ..$ : chr &quo...