similar to: add confidence intervales to xyplot for ANCOVA and extracting info

Hi, what would be a good way of enhancing the ancova function from the HH package when using a 2 factor ANCOVA? The current behaviour for the "ancova" function from package HH is: ---------------------------------------------- | P1 || P1 || P3 || P4 | | PS | <- the lattice strip ---------------------------------------------- | x|| x ||x ||x | |

Dear R-Users, I am resending this message just to reminder my question regarding the calculation of a bootstrap confidence intervals for a GAM plot. I am trying to apply a bootstrap to a GAM in order to calculate the 95% confidence intervals for a smooth curve obtained by the ?plot.gam? function of the mgcv package. Nonetheless, I am getting some difficulties in transposing the results for

Hi All, I have a data set which I need to plot and show the values of one of the variables as a second x-axis. library(lattice) year<-c(2001,2002,2003,2004,2005,2006) fac<-c("arts","arts","arts","sci","sci","sci") staff<-c(95,98,99,32,31,36) part<-c(32,31,33,15,16,14) df1<-data.frame(year,fac,staff,part)

Hello, I'm making a simple plot using xYplot in the Hmisc library and having problems with labeling the values on the x-axis. Using the reproducible example below, how can I have the text (jan, feb,mar, etc.) in place of 1:12. Thanks, AB x <- c(seq(0,0.5,by=0.1),seq(0.5,0,by=-0.1)) ci <- rnorm(12,0,sd=0.1) xupper <- x + ci xlower <- x - ci mo.fac <- c("jan",

Hi, is there a grayscale setting for lattice plots? I like the default color settings. I also like the settings that are available for setting black and white with something like this: --8<---------------cut here---------------start------------->8--- ltheme <- canonical.theme(color = FALSE) ## in-built B&W theme ltheme$strip.background$col <- "transparent" ##

Dear List, I have been writing a Lattice function to draw what we call stratigraphic diagrams, these are diagrams with a panel for each species showing a time series of abundance, but drawn vertically to represent time passing from bottom of plot towards to top. I am most of the way there with this now, but I want to do away with the strip on each panel and instead draw a custom top axis with a

Hello, is there any function in the package snow to check for a really running cluster? The function checkCluster only checks the variable cl. And the variable is still available after stopping the cluster! ( a simple solution would be deleting the cluster variable cl in the function stopCluster) > library(snow) > cl <- makeCluster(5) 5 slaves are spawned successfully. 0

Dear list, this may not be related to R but rather to my OS, but I do not understand the issue of compiling R packages deeply enough to figure out the exact cause of the problem. I am trying to install a R package from source as it is not yet available under Cran (Rssa, downloaded here: https://github.com/asl/rssa). Running sudo R CMD INSTALL asl-rssa-6f458e4.tar.gz from the console

Dear list, in R under wheezy the following error occurs: > library(cairoDevice) > Cairo() Error in Cairo() : Graphics API version mismatch It seems to be Debian-specific, as I tried the same under OpenSUSE 11.3 and it worked. I tried installing cairoDevice via `apt-get install r-cran-cairodevice' and with `install.packages("cairoDevice")' in R (without any

Dear Developers, I'm updating my R installation via svn (currently: Revision 44626). For some time now (sorry don't know the exact revision number) make for R devel yields the following error make[4]: *** Keine Regel vorhanden, um das Target ?../../../src/include/R_ext/GraphicsBase.h?, ben?tigt von ?g_fontdb.o?, zu erstellen. Schluss. (translation: no rule to make target ...

Dear userRs, "playing around" with combinations of replicate() and random number generating functions inside a self-defined "wrapper" function I encounterd a puzzling behaviour. The following are intentionally simple (and rather nonsense-) examples to isolate the relevant aspects. Please, note the seemingly "inconsistent" behaviour for the second call of

With newest R devel #sessionInfo() R Under development (unstable) (2012-01-10 r58085) Platform: amd64-portbld-freebsd10.0 (64-bit) locale: [1] de_DE.ISO8859-15/de_DE.ISO8859-15/de_DE.ISO8859-15/C/de_DE.ISO8859-15/de_DE.ISO8859-15 attached base packages: [1] stats graphics grDevices utils datasets methods base I get the following error when I try to build and install lubridate from

Dear R list, I get a strange error in R: *** caught segfault *** address 0x6951c20, cause 'memory not mapped' Traceback: 1: .C("spline_eval", z$method, nu = as.integer(n), x = as.double(xout), y = double(n), z$n, z$x, z$y, z$b, z$c, z$d, PACKAGE = "stats") 2: spline(gam.data$x[, col.data], gam.smooths.all$fit[, m], xout = gam.results.global[m, ,

Dear UseRs, here is an example scenario presenting my problem: Multiplying a dsyMatrix with a numeric vector results in an error (unfortunately in German due to my locale): > (M1 <- Matrix( c( 1, 2, 2, 2, 1, 2, 2, 2, 1), nrow = 3)) 3 x 3 Matrix of class "dsyMatrix" [,1] [,2] [,3] [1,] 1 2 2 [2,] 2 1 2 [3,] 2 2 1 > M1 %*% 1:3 Fehler in x %*%

I have an application that uses I18n, but in the first deployment, I don''t want to enable different languages per user/request because there are very few users who all speak german). Therefore I have config.i18n.default_locale = :de in config/ application.rb. I have no code whatsoever to set/change I18n.locale in my app. In development environment, it works as expected: the

Hello list, here is what I stumbled upon: 1> test <- data.frame(time=as.POSIXct((1:2),origin="2000-1-1")) 1> test time 1 2000-01-01 00:00:01 2 2000-01-01 00:00:02 1> rbind(test,b=1:2) Fehler in as.POSIXct.numeric(value) : 'origin' muss angegeben werden When I try to attach an additional row to a dataframe with a row containing dates I get the

Hi there is it possible that pdfs generated using the pdf() function with default settings leads to loss of information? I was plotting copy number changes from Agilent 180k data in form of rectangles (rect()) while each rectangle represents one region of copy number change. When plotting into a pdf I noticed that some very small rectangles do not appear (even after extensive zooming) in the pdf

Dear list, I have a balance connected to the serial port of a windows machine ("COM1") and I read the text output of the balance with scan("COM1", what="character", sep="\n", n=1) after calling the previous line I press the print key on the balance which triggers sending one line of text to the serial connection and with R 2.9.2 I get something like Read

According to the instructions given at https://www.asterisksounds.org/de I converted and installed German prompts successfully and for numbers, I can successfully listen to a German female voice counting or telling the date/time. But unlikily, somehow the voicemail prompt is still English, although my general language settings are "de". I use pjsip.conf, not sip.conf. In

Hi, I am analyzing a data set with greater than 1000 independent cases (collected in an unrestricted manner), where each case has 3 variables associated with it: one, a factor variable with 0/1 levels (called XX), another factor variable with 8 levels (X) and a third response variable with two levels (Y: 0/1). I am trying to see if X1 has an effect on the relationship between X2 and the