similar to: Suggestions for enhanced routines for "mlm" models.

Full_Name: Yves Rosseel Version: 2.2.1 OS: i686-pc-linux-gnu Submission from: (NULL) (157.193.116.152) Dear developers, For the single-model case, the anova.mlm() function does not seem to handle multi-parameter predictors (eg factors) correctly. A toy example illustrates the problem: Y <- cbind(rnorm(100),rnorm(100),rnorm(100)) A <- factor(rep(c(1,2,3,4), each=25)) fit <- lm(Y ~ A)

Hi, Am newish to stats and R, so I certainly appreciate any help. Basically I have 50 inidividuals whom I have 6 photos each of their optic nerve head. I want to check that the orientation of the nerve head is consistent, ie the 6 replicates show minimal or preferably no rotation differences. I'll draw an arbitrary line between some blood vessels (same reference in each set of replicates) and

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

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

Hello all, I have a question on manova in R: I'm using the function "manova()" from the stats package. Is there anything like a stepwise (backward or forward) manova in R (like there is for regression and anova). When I enter: step(Model1, data=Mydata) R returns the message: Error in drop1.mlm(fit, scope$drop, scale = scale, trace = trace, k = k, : no 'drop1'

Discover "The Secret to Supercharge your MLM!" A must read... Hi there, The secret is out! Here's the mail we've all been waiting for! Read carefully and take immediate action on it! I've discovered an amazing formula that will give your MLM an enormous enrollers explosion. You'll benefit hugely if you use it with YOUR primary MLM! -THE ULTIMATE RECRUITMENT

I had a look at some influence measures, and it seems to me that currently several methods handle multiple lm (mlm) objects wrongly in R. In some cases there are separate "mlm" methods, but usually "mlm" objects are handled by the same methods as univariate "lm" methods, and in some cases this fails. There are two general patterns of problems in influence measures:

>>>>> steven pav >>>>> on Thu, 19 Jul 2018 21:51:07 -0700 writes: > It seems that confint.default returns an empty data.frame > for objects of class mlm. For example: > It seems that confint.default returns an empty data.frame for objects of > class mlm. Not quite: Note that 'mlm' objects are also 'lm' objects, and so it is

Full_Name: Yves Rosseel Version: 2.2.1 OS: Submission from: (NULL) (157.193.116.152) The code below sketches a possible implementation of a function 'mlh.mlm' which I think would be a good complement to the 'anova.mlm' function in the stats package. It tests a single linear hypothesis of the form H_0: LBT'= 0 where B is the matrix of regression coefficients; L is a matrix

I would like to perform multiple comparisons or post-hoc testing on the independent variable in an Anova.mlm object generated by the Anova function of the car package. I have defined a multivariate linear model and subsequently performed a repeated measures ANOVA as per the instructions in section #3 of the following comprehensive tutorial on the subject from the Gribble lab at UWO:

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

Dear list, I have been using multivariate multiple regression (MMR) in the form lm(Y~X) where Y and X are matrices of response and predictor variables. I know that summary(mlm.object) would give the usual lm statistics for each response variable separately and that anova.mlm(mlm.object) will give the analysis of variance table of the mlm object. However, anova.mlm (also manova(mlm.object))

I'm an experimental psychologist and when I run ANOVA analysis in SPSS, I normally ask for a test of non-sphericity (Box's M-test). I also ask for output of the corrections for non-sphericity, such as Greenhouse-Geisser and Huhn-Feldt. These tests and correction factors are commonly used in the journals for experimental and other psychology reports. I have been switching from SPSS to R

[hope this doesn't get posted twice, my first mail bounced] Hi, again, this is about the call bug in anova.mlm. I originally reported it as PR#7898 and I suggested a fix at PR#7904. (None of these message were forwarded to R-Devel, instead I received a bounce telling me that my provider's SPF settings are incorrect. They are, though, so there seems to be a problem with the

When I fit a multivariate linear model, and the formula is defined outside the call to lm(), the method summary.mlm() fails. This works well: > y <- matrix(rnorm(20),nrow=10) > x <- matrix(rnorm(10)) > mod1 <- lm(y~x) > summary(mod1) ... But this does not: > f <- y~x > mod2 <- lm(f) > summary(mod2) Error en object$call$formula[[2L]] <- object$terms[[2L]]

I'm trying to extend some work in the car and heplots packages that requires getting a table of multivariate means for one (or later, more) terms in an mlm object. I can do this for concrete examples, using aggregate(), but can't figure out how to generalize it. I want to return a result that has the factor-level combinations as rownames, and the means as the body of the table

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

Hello, I've just fitted a model with multi-responses, and I get an object of class "lm" "mlm". My problem is that as soon as I invoke the predict method for a dataframe "newdata", the methods runs and give me back prediction at the fitting points but not for newdata. Does someone has an explanation for this behavior, and some ideas to make predict.mlm work

predict.mlm fails with the error: Object "X" not found. CAUSE: looking at the predict.mlm function, we have: x <- model.matrix(object, newdata) two lines later, we have: pred <- X[, piv, drop=FALSE] %*% object$coefficients[piv,] SOLUTION(?): change the case of one of the x's, since it seems like what is meant in the second line is lowercase x... platform:

On Sat, 12 Aug 2006, takahashi kohske wrote: > Dear list members: > > I'd like to one-way repeated measured anova by using mlm. > I'm using R-2.3.1 and my code is: > > dat<-matrix( c(9,7,8,8,12,11,8,13, 6,5,6,3,6,7,10,9, > 10,13,8,13,12,14,14,16, 9,11,13,14,16,12,15,14), > ncol=4, dimname=list(s=1:8, c=1:4)) >