search for: bela_b

Hi, I'm encountering a very odd problem with calls to anova.mlm() from within a function. Consider the following code (data.n is a matrix of numeric values): mlmfit <- lm(data.n ~ 1) mlmfit0 <- lm(data.n ~ 0) print(mlmfit) anova(mlmfit,mlmfit0,test="Spherical") If I run it just like this from the console, it works just fine. If, however, I call it from within a function,

....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. I would suggest PR#7904 to be deleted; I meant it to be a comment on PR#7898, but by mistake added an extra ":" to the subject line. Regards, Bela -- Bela Bauer - bela_b@gmx.net PGP 0x97529F5C http://www.aitrob.de

Hi, I'm using anova.mlm sphericity tests/corrections, and I'm getting different values than SAS. In order to be able to use these values for publications, I'd need to know more about the SAS bug that is mentioned in the Reference Manual. - What exactly causes the different values? - Is it just a slight difference, or can I expect significant differences in H-F/G-G epsilons and

Hi, yet another anova.mlm problem - it doesn't seem to end. This time, I have a setup with a few within-subject factors and a between-subject factor (SGROUP). Consider the most simple case with only one within-factor (apo): > mlmfit0 <- lm(data.n ~ 0 + SGROUP) > mlmfit1 <- lm(data.n ~ 1 + SGROUP) > anova(mlmfit1,mlmfit0,test="Spherical",M=~hemi,X=~1) Analysis of

Hi, I'm currently working on porting some SAS scripts to R, and hence need to do the same calculation (and get the same results) as SAS in order to make the transition easier for users of the script. In the script, I'm dealing with a two-factorial repeated-measures anova. I'll try to give you a short overview of the setup: - two between-cell factors: facBetweenROI (numbering

Hi, I'm still trying to figure out that GLM procedure in SAS. Let's start with the simple example: PROC GLM; MODEL col1 col3 col5 col7 col9 col11 col13 col15 col17 col19 col21 col23 =/nouni; repeated roi 6, ord 2/nom mean; TITLE 'ABDERUS lat ACC 300-500'; That's the same setup that I had in my last email. I have three factors: facSubj,facCond and facRoi. I had this pretty

Hi, sorry to bother you again, but I can't figure it out myself and I also can't find any in-depth documentation about it... Consider the following SAS code (A1II2... contain 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

Hi, I'm still not quite there with my H-F (G-G) correction code. I have it working for the main effects, but I just can't figure out how to do it for the effect interactions. The thing I really don't know (and can't find anything about) is how to calculate the covariance matrix for the interaction between the two (or even n) main factors. I've looked through some books

Hi, I'm still not quite there with my H-F (G-G) correction code. I have it working for the main effects, but I just can't figure out how to do it for the effect interactions. The thing I really don't know (and can't find anything about) is how to calculate the covariance matrix for the interaction between the two (or even n) main factors. I've looked through some books