similar to: lme newbie question

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 am using lme to handle repeated measurements. So far i can follow the examples from the book from pinheiro and bates. Now i get the problem , that i have "nested" repeated measuremnts, and i cant find out how to do the grouping part of the lme formula. 1.The simple problem ist that i have different Samples , from which i make repeated measurements (each sample is measured 6 times) and

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

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

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.

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

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

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

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

Dear Madam or Sir, I need to do multiple comparisons following nlme analysis (Compare the effects of different treatments on a response measured repeatedly over time; fixed = response ~ treat*time). On the web I found the notion that one might use the L argument from ANOVA. Do you have an example to show how this works together with nlme? Are there other ways to do a post-hoc analysis in

Hi, Another question about graphics... Suppose I have: par(mfrow = c(4, 6)) which will set up a figure region for 24 plots. But, I actually only need 23 plots. In fact I wrote a little function, foo.plot, which will draw all the 23 plots at once with: foo.plot(x) where x is a data frame. Now, suppose I want to draw several set of these 23 plots, but want to put them all in one postscript

Dear all, I am working with R Version 0.62.2 (July 10, 1998) on a SuSE Linux 5.2. I have troubles in understanding the function image (). Can someone enlighten me? I'd like to draw a square image-plot (with every little part of the image a square as well) and add some (straight) lines to it. Also, I create an empty z matrix which I need to fill step by step. Thus I tried to find out how

Dear R users, I have not found anything on this in the archives. Does anyone know whehther the parameter use= is not functioning in cor or enlighten me what it is supposed to do? My R version is "R version 1.8.1, 2003-11-21" on Windows 2000. I am hoping to be able to update to 1.9.1 as soon as it has appeared (we are not allowed here to install software on our own and thus I am trying

Dear all, I think this question has been raised very recently but I think I didn't see repsonses on the list nor did I find the thread in the archives. (If it is there I'd appreciate directions ...) I am about to install Linux on a Compaq notebook. One important thing I want to do is being able to run R. I have fiddled with different distributions of Linux (LST, SuSE, Debian) and can

Dear all, in a current project I have a pretty huge .Rdata. Thus I was working with R --vsize 100 --nsize 1000000. Today when I tried to restart R I get the following error message: Error: a read error occured Fatal error: unable to restore saved data (remove .RData or increase memory) I increased memory up to --vsize 180 and --nsize 2000000, but the error reoccurrs. Is there a way to know

Dear all, well, today I just seem to run into these problems ... (I am still working on R 0.63.1 on a SuSE 5.2). I would like to have a data.frame that also contains character (and NOT factor) vectors. Two things that I tried (the second worked in R 0.62.2): > test.df <- data.frame (X= character (4)) > test.df$X [1] <- "a" Warning: invalid factor level, NAs generated or:

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