similar to: grouping for lme with nested repeated measurements

Hello All, I have a more statistical question, and how this is implemented in R. The problem is the following: We have 2 different solutions (samples), which are filtered and then the concentration of the filtrate is measured. We want to evaluate how the filter proces and the concentration measurement influences the detection of the difference of the two solutions and which step has which

Hello, Is there like the kruskal wallis test in relation to ANOVA (no restrictions on normallity and variance homogenity) something (in R) for MANOVA? thanks -- Dr.Nicolaas Busscher Universit?t GH Kassel Nordbahnhofstrasse: 1a, D-37213 Witzenhausen Phone: 0049-(0)5542-98-1715, Fax: 0049-(0)5542-98-1713

enclosed a simple R script (and a data file) , with calls lda similar to the example with the iris data in the documentation. it is not working and i dont understand the error message. can anybody help me? i am using R 1.5.1 (2002.06.17) on debian woody stable. thanks -- Dr.Nicolaas Busscher Universit?t GH Kassel Nordbahnhofstrasse: 1a, D-37213 Witzenhausen Phone: 0049-(0)5542-98-1715, Fax:

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

I don't think, this has been answered: > I'm trying to run a 3-way within-subject anova in lme with 3 > fixed factors (Trust, Sex, and Freq), but get stuck with handling > the random effects. As I want to include all the possible random > effects in the model, it would be something more or less > equivalent to using aov > > > fit.aov <- aov(Beta ~ >

enclosed a simple R script (and a data file, and the output) , with calls lda similar to the example with the iris data in the documentation. it is not working and i dont understand the error message. can anybody help me? i am using R 1.5.1 (2002.06.17) on debian woody stable. I would like to avoid updating now, because i want to keep the system in "stable". can somebody test, if the

Does someone out there have an example of R-code for a split-split plot ANOVA using aov or another function? The design is not balanced. I never set up one in R before and it would be nice to see an example before I tackle a very complex design I have to model. Thanks, Mike Mike Saunders Research Assistant Forest Ecosystem Research Program Department of Forest Ecosystem Sciences University of

Dear all, is this a bug or do a miss something (this concerns version 0.63.1 on Linux)? > boxplot (c (1, 2, 3, 4, 5), ylim= c (-5, 5)) Warning: parameter "ylim" couldn't be set in high-level plot() function > boxplot (c (1, 2, 3, 4, 5), pars= list (ylim= c (-5, 5))) Warning: parameter "ylim" couldn't be set in high-level plot() function > and actually in

I am trying to perform a model 3 ANOVA for a 2 factor (say factor A and factor B) anova in which factor A is fixed and factor B is random. Therefore, the error term for the test of factor A should be the A:B interaction term and the error terms for B and A:B should be the model residual (within) term. I have tried to work out how to specify such error strata using aov, however, I have had

Dear all, I do not seem to grasp how contrasts are set for ordered factors. Perhaps someone can elighten me? When I work with ordered factors, I would often like to be able to reduce the used polynomial to a simpler one (where possible). Thus, I would like to explicetly code the polynomial but ideally, the intial model (thus, the full polynomial) would be identical to one with an ordered factor.

Hi there Please excuse this elementary question, but I have been fumbling with this for hours and can't seem to get it right. I have a nested anova, with random factor "lakefac" nested within factor "fishfac" (fixed), with an additional fixed factor "Habfac". If I consider everything as fixed effects, it's addmittedly not the correct model, but I can at

Hi I try to implement a simple 2-factorial repeated-measure anova in the lme framework and would be grateful for a short feedback -my dependent var is a reaction-time (rt), -as dependent var I have -the age-group (0/1) the subject belongs to (so this is a between-subject factor), and -two WITHIN experimental conditions, one (angle) having 5, the other 3 (hands) factor-levels;

Dear all, Many thanks to all who always so quickly respond with helpful hints! Currently I am working with R 0.63.1 on SuSE 5.2 which I installed this morning within about half an hour: my deepest compliments to the developpers! (I had still 0.62.2 before.) This time I am trying the following: > xy <- list (a= c (1, 2, NA), b= c (2, 3)) > list.name <- "xy" > get

I do not think anyone has answered this. > I'm trying to run a simple one-way ANCOVA with the lmer > function in R package lme4, but have encountered some > conceptual problem. The data file MyData.txt is like this: > > Group Subj Cov Resp > A 1 3.90 4.05 > A 2 4.05 4.25 > A 3 4.25 3.60 > A 4 3.60 4.20 > A 5 4.20 4.05 > A 6 4.05 3.85

Hi ! I try to compare some models obtained from glmmPQL. model1 <- glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4 +I(freq8_4^2), random=~1|num, binomial); model2 <- glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4 , random=~1|num, binomial); anova(model1, model2) here is the answer : Erreur dans anova.lme(model1, model2) : Objects must

En innebygd og tegnsett-uspesifisert tekst ble skilt ut... Navn: ikke tilgjengelig Nettadresse: https://stat.ethz.ch/pipermail/r-help/attachments/20080129/5f8a6ac4/attachment.pl

Hello. I am fitting a two-level mixed model which assumes equality of variance in the lowest-level residuals across groups. The call is: fit3<-lme(CLnNAR~CLnRGR,data=meta.analysis, + na.action="na.omit",random=~1+CLnRGR|study.code) I want to test the assumption of equality of variances across groups at the lowest level. Can someone tell me how to do this? I know that one

Hi everybody, I've found a very interesting problem in R: the if-else statement doesn't work in a main program. Sounds crazy, but true. I tried this very easy example, and I got syntax error at the "else" line. Example: ------- a <- 1 if ( a == 1 ) print("yes") else print("no") -------- I tried on Windows and on Linux, but none of them works.

I am trying to run a GLMM on some binomial data. My fixed factors include 2 dichotomous variables, day, and distance. When I run the model: modelA<-glmmPQL(Leaving~Trial*Day*Dist,random=~1|Indiv,family="binomial") I get the error: iteration 1 Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 >From looking at previous help