similar to: Problem with lme4

I?m entirely stumped on this particular issue and really hoping someone has some advice for me. I am running a covariant model in lm I would like to give the standard errors or the confidence intervals for the fitted lines. I?ve been using the dataset OrangeSprays where I want lines for each level of treatment over the covariant ?colpos?. I?ve been able to calculate intercepts and slopes for

I feel bad even asking, but: Rgames> data(OrchardSprays) Rgames> model<-lm(decrease~.,data=OrchardSprays) Rgames> model Call: lm(formula = decrease ~ ., data = OrchardSprays) Coefficients: (Intercept) rowpos colpos treatmentB treatmentC 22.705 -2.784 -1.234 3.000 20.625 treatmentD treatmentE treatmentF treatmentG treatmentH

hello all, I am trying to use the function getCovariateFormula(nlme) in conjunction with the library lme4. When I load both packages I get the following message and the getCovariateFormula function no longer works: library(nlme) library(lme4) Attaching package 'lme4': The following object(s) are masked from package:nlme : contr.SAS getCovariateFormula

Hi All, I'm trying to draw boxplots. I'm having a hard time to get "ticks labels" on multiple panels using 'alternating' option. # R Code: # May not be the best example, please just look into 'scales' option library(lattice) data(OrchardSprays) dta <- subset(OrchardSprays, OrchardSprays$rowpos %in% c(1,2,3)) # Original # This works fine, as you can see

Hi All, Can anyone please help me with getting a "single title" and "legend" for both the plots in the following R code. I'll eventually be using .wmf file. # R code: library(lattice) p1 <- xyplot(decrease ~ treatment, OrchardSprays, groups = rowpos, type = "a", main="Same title / legend", auto.key = list(space = "right",

Here is an example: The legends (mainly the factor level names) are cut off on the right. Somehow the internal calculation which computes horizontal space for the legend is not flexible enough. ## Call a new graphics window {with default par()s !}: get(getOption("device"))() par(mfrow = c(2,2)) ## part of example(interaction.plot) _improved_ using with() : data(OrchardSprays)

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@r-project.org # ###################################################### <<insert bug report here>> Is this a bug in the parser? Specifically, the trace

hello everybody, I want to print two plots in one png file, I tried several options but i didn't succeed the first plot (bwplot) print to the defined position, but the second (ggplot) doesn't Any idea? Thanks a lot Christophe # Example: #------------------------------------- library(ggplot2) library(lattice) library(grid) one <- bwplot(decrease ~ treatment, OrchardSprays, groups

Hi all, I want to retrieve the stats from a 'bwplot' with one factor. I have read the help for 'panel' function and I'm aware of the option 'stats' which defaults to 'boxplot.stats' but I didn't understand it well and therefore I am unable to get what I need. Thanks in advance. Héctor Villalobos Windows XP, R 2.5.0 > bwplot(decrease ~ colpos |

Hi all, I believe I've run into this before, but I seem to have totally forgotten. No headway in the last couple of hours either. How do I make sure that points and margins remain the same absolute size as I change the resolution of a device? (I'm running 2.9.1 patched, on a Win XP-machine) Many thanks in advance, Gustaf Ps: As an afterthought, might it be that this behaviour is

Hi, After installing the latest versions of lme4, Matrix, VR on a Debian Box, I have run into a problem. When I use library(lme4) on R Version 2.1.0 (2005-04-18) I get the following error message: > library(lme4) Loading required package: Matrix Error in autoloader(name = "confint", package = "MASS") : autoloader did not find 'confint' in 'MASS'

Hello all, I wrote a small function to add labels for outliers in a boxplot. This function will only work on a simple boxplot/formula command (e.g: something like boxplot(y~x)). Code + example follows in this e-mail. I'd be happy for any suggestions on how to improve this code, for example: - Handle boxplot.matrix (which shouldn't be too hard to do) - Handle cases of complex

Hi All, I need a file which I can import to MS Word, I'm trying win.metafile(), but it does not seem to support multiple print commands at once (please see below). Is there an alternative to get plots file which can be used in powerpoint/word? # R code: # this does not work; but same thing works with pdf() library(lattice) win.metafile("test.wmf") p1 <- xyplot(decrease ~

Dear Rhelp, I would like to get alternate tick labels for the xyplot: library(lattice) library(grid) xy <-xyplot(decrease ~ treatment, OrchardSprays, ??? main= "Some plot", ??????? groups = rowpos, type = "a", ??? ??????????? page???? = function(n){ ???????????????????????? grid.text(LETTERS[j], ???????????????????????? y = 0.95, ???????????????????????? x = 0.15,

Dear R-ers, I'm not sure if this is a missing feature, a support request, or stupidity on my part, but nevertheless, its a question. Is it possible to add titles to colorkey legends? As far as I can tell, there is a command to do it for normal "key" legends, but not for "colorkeys". eg it works for a normal key, created through auto.key xyplot(decrease ~ treatment,

Dear R-help, I can't plot nls objects for some reason. The following example, taken from help(plot.nls) illustrates the problem: > data(Orthodont) > fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject) > plot(fm1, resid(., type = "p") ~ fitted(.) | Sex, abline = 0) Error in do.call(plotFun, as.list(args)) : couldn't find function

I found that predict.gnls failed with a wierd error message about a non-numeric argument to a binary vector in one of three nearly identical uses. Error in Inh/Ki : non-numeric argument to binary operator (Inh and Ki are arguments to the function used in the formula for the object whose predictions were requested). It turns out that the problem is in getCovariateFormula(). The final line in

I'm sure I'm being stupid so flame away... R2.2.1 on Windoze (boohoo) latest updates of packages. I'm exploring a dataset (land) with three variables looking at an narrowly unbalanced two group (GROUP) ANCOVA of a randomised controlled trial analysing endpoint score (SFQ.LOCF.ENDPOINT) entering the baseline score (SFQ.BASELINE) as covariate and the following work fine: > res.same

I'm looking at Pinheiro & Bates "Mixed-Effects Models in S and S-PLUS" at the moment. Several datasets are used, one of which is called "PBIB" (a partially balanced incomplete block design). All the other datasets can be found somewhere or other in R. However, I cannot locate PBIB, and it does not seem to be mentioned in the latest edition of the R Full Reference

Hi, I guess this is a final plea, and maybe this should go to R-help but here goes. I am writing a set of functions for calibration and prediction, and to calculate standard errors and intervals I need the variance function to be evaluated at new prediction points. So for instance fit<-gnls(Y~SSlogis(foo,Asym,xmid,scal),weights=varPower())