search for: mlmlist

...it just like this from the console, it works just fine. If, however, I call it from within a function, e.g. using fct <- function(data.k) { # same code with data.n replaced by data.k } fct(data.n) It gives me > fct(data.n) Call: lm(formula = data.k ~ 1) Coefficients: ... Error in anova.mlmlist(object = mlmfit, mlmfit0, test = "Spherical") : Object "mlmfit" not found What causes mlmfit to "disappear" between two lines? I haven't got the slightest clue where to look for an answer... Thank you very much for your help. Bela Bauer -- Weitersagen:...

...ly. (PR#7824) o select.list(multiple=TRUE) now detects and tries again for invalid text input. o add1.[g]lm could give strange results with interaction terms when the model and the upper scope had different orders for the main effects. (PR#7842) o A bug had sneaked into the anova.mlmlist() code, affecting the Greenhouse-Geisser epsilon. Code wrongly assumed a matrix to be symmetric. (Thanks to Bela Bauer.) o anova.mlmlist() and mauchley.test() are now more tolerant to rank deficiency in the M and X matrices (also when they are implicitly generated via model.matrix())....

...ly. (PR#7824) o select.list(multiple=TRUE) now detects and tries again for invalid text input. o add1.[g]lm could give strange results with interaction terms when the model and the upper scope had different orders for the main effects. (PR#7842) o A bug had sneaked into the anova.mlmlist() code, affecting the Greenhouse-Geisser epsilon. Code wrongly assumed a matrix to be symmetric. (Thanks to Bela Bauer.) o anova.mlmlist() and mauchley.test() are now more tolerant to rank deficiency in the M and X matrices (also when they are implicitly generated via model.matrix())....

...call for the fit objects, but the same problems occurs for idata. Hence, there needs to be an additional change made to anova.mlm. I basically changed anova.mlm to the following: #... if(length(list(object, ...)) > 1){ cl <- match.call() cl[[1]] <- as.name("anova.mlmlist") cl[[2]] <- list(object, ...) if (!is.null(cl$idata)) { cl$idata <- idata } return(eval(cl)) } # ... Additionally, anova.mlmist has to be changed to accept a list as parameters. I've attached a patch of my current code against R-2.1.0...

Dear WizaRds, I am sorry to bother you with a newbie question, but although I tried to solve my problem using the various .pdf files (Introduction, help pages etc.), I have come to a complete stop. Please be so kind as to guide me a little bit along my way of exploring multivariate analysis in R. I want to test wether the means-vector mu1 of X, consisting of the means per column of that matrix

...d at least provide the estimated residual covariance matrix (or SSD structure). vcov currently inherits from "lm" leading to a completely arbitrarily scaled matrix. The correct matrix is a Kronecker product of the unscaled covariance matrix and estVar. anova.mlm, anova.mlmlist, drop1.mlm, add1.mlm These can (seemingly...) be obtained by relatively small modifications of their lm counterparts. The actual test calculations need to be excised from summary.manova (generalize? e.g. mvlin.test(SSD1, SSD2, method="Pillai")). It should be possible...

...ng directly. o summary.manova() uses a stabler computation of the test statistics, and gains a 'tol' argument to allow highly correlated responses to be explored (with probable loss of accuracy). Similar changes have been made to anova.mlm() and anova.mlmlist(). o Sweave() now writes concordance information inside a \Sconcordance LaTeX macro, which allows it to be inserted into PDF output. o system.time() now uses lazy evaluation rather than eval/substitute, which results in more natural scoping. (PR#11169) o...

...ng directly. o summary.manova() uses a stabler computation of the test statistics, and gains a 'tol' argument to allow highly correlated responses to be explored (with probable loss of accuracy). Similar changes have been made to anova.mlm() and anova.mlmlist(). o Sweave() now writes concordance information inside a \Sconcordance LaTeX macro, which allows it to be inserted into PDF output. o system.time() now uses lazy evaluation rather than eval/substitute, which results in more natural scoping. (PR#11169) o...