similar to: Label individual points in lattice by and ID

All I have been looking into the books on performing statistics in R, in particular I am interested in General Linear Mixed Models, for Randomized Complete Block Design Experiments The list I have come away with so far is: Mixed Effects Models in S and S-plus by Pinheiro (2002)

All I want to overlay two variables on the same plot following their appropriate grouping. I have attempted to use subscripting in panel with panel.xyplot, but I can't get the grouping to follow into the panel...here is an example... dat<-data.frame( y= log(1:10), y2=10:19, x=1:10, grp = as.factor(1) ) dat2<-data.frame( y= log(10:19), y2= 20:29, x=1:10, grp = as.factor(c(2)) )

I am having trouble with variable scoping inside/outside functions. I want to use get() to grab a named and quoted variable as an input to a function, however the function can't find the variable when it is entered into the function call, only when it is named in the main environment. I obviously am not so clear on variable scoping, and ?get really doesn't clear up my confusion. I am

R-Help I have a data set which uses a DateTime string as follows : "2009-06-30 18:14:59" While I have been able to convert to DateTime properly time <- strptime(as.character(dat$Time),format='%d %b %Y %T') #Convert to dateTime string I would like to use the time of day "hour" as a *factor* level. I have found that I can convert the date time to a factor time

I am attempting to replicate some of my experience from SAS in R and assume there are best methods for using a combination of summary(), subset, and which() to produce a subset of mean values by categorical or ordinal factors. within sas I would write proc means mean data=dataset; class factor1 factor2 var variable1 variable2; RUN; producing an output with means for each variable by factor

I want to use a strip.custom as with useOuterStrips for three conditioning variables. useOuterStrips restricts this to only two conditioning variables, and I cannot figure out how to write strip.custom properly to do this. library(lattice) mtcars$HP <- equal.count(mtcars$hp) #with two factors x2<-xyplot(mpg ~ disp | HP + factor(cyl), mtcars) useOuterStrips(x2) #with three factors

Hi, I have a lattice lot conditioned on two variables. Example code is: library(lattice) x <- data.frame(d=runif(100), f1=sample(c('yes', 'no'),100,replace=TRUE), f2=c(rep('Run1',30),rep('Run2',30),rep('Run3',40))) histogram(~d | f1 + f2, x) In the plot, for a given value of f2, there are two panels, one for 'n'

Hi, I am trying to add text to the bottom of a lattice bwplot with multiple panels. I would like to add a label below each boxplot, but the labels do not come from the data. I've tried the following, code: f1 <- c(rep(c(rep("a", 3), rep("b", 3), rep("c", 3)), 2)) f2 <- c(rep("A", 9), rep("B", 9)) dv <- c(0.9, 0.8, 0.85, 0.6, 0.65,

Dear R- Community, to learn i reanalysed some data provided and analysed by Zuur et. al. in their book "Mixed effect models and Extensions in Ecology with R". When i run the last command i get a warning message i dont understand. Loyn<- read.table(file = "loyn.txt",header = TRUE) Loyn$L.AREA<- log10(Loyn$AREA) fGRAZE <-factor(Loyn$GRAZE) M0<- lm(ABUND~ L.AREA

Both vector query's can select the values from the data.frame as written, however in the first form assigning a value to said selected numbers fails. Can you explain the reason this fails? dat <- data.frame(index = 1:10, Value = c(1:4, NA, 6, NA, 8:10)) dat$Value[dat$Value == "NA"] <- 1 #Why does this fails to work, dat$Value[dat$Value %in% NA] <- 1 #While this does

Dear R-list, I would like to use the lattice library to show several groups on the same graph. Here's my example : ## the data f1 <- factor(c("mod1","mod2","mod3"),levels=c("mod1","mod2","mod3")) f1 <- rep(f1,3) f2 <-

Hi, I use a Debian Sarge with samba 3.0.23c (from testing) and cups 1.1.23. For unexpected reason, some printers are sometime stopped (pause), and I need to start then again. There is no trace of a "stop" in the log of cups. (and I'm the only one who can pause them, anyway). here is an example : cat 10.lpr.log | grep imp-cao | grep -E "(stop|start)" Nov 10 14:15:41

Hello all, I have a problem with starting the winbind daemon. I want to connect to a win2003 server domain. I get the following error message when i start winbindd with: winbind -d10 -i winbindd version 3.2.7 started. Copyright Andrew Tridgell and the Samba Team 1992-2008 initialize_winbindd_cache: clearing cache and re-creating with version number 1 Could not fetch our SID - did we join?

Hello, I'd like to produce a ggplot where the order of factors within facets is based on the average of another variable. Here's a reproducible example. My problem is that the factors are ordered similarly in both facets. I would like to have, within each facet of `f1', boxplots for 'x' within each factor `f2', where the boxplots are ordered based on the average of x

Hi, Is there any way to suppress plotting of panels that don't actually contain any information? I have tried using 'drop.unused.levels=TRUE', but there doesn't seem to be any effect. Here is an example: library(lattice) # some fake data: d <- data.frame(x=runif(20), x.class=rep(letters[1:5], each=4), f1=rep(letters[1:2], each=10), f2=rep(letters[10:19], each=2) ) # plot

Dear list, I have a two character vector with two different values in them (two each, that is). Naturally, when I use these vectors as grouping factors in a lattice plot, I get four panels. Now, one of the possible four combinations will never have data, so one panel will always be empty. So now my question is: can I force R / Lattice not to print it, and give me only the 3 filled panels? I have

Hi List, I am running R 2.8.0 on a Windows XP machine, running ggplot2 version 0.8.1 I want to label the bars in a faceted grid barplot. Reproducible R code is given below: #### reproducible facet barplot ##### library(ggplot2) # Dataset from which to create the barplot ml <- rep(1:10,2) vals <- rnorm(20,mean = 10, sd=1) type <- c(rep("MAPE",10),rep("AIC",10))

I am suddenly no longer able to see any existing files on one of my shares, from my Windows clients. I have a Linux server (Fedora Core 8) which used to be my firewall and router as well as file and print server. There are two Samba shares on it I have mounted on various Windows clients: my home directory and a multimedia directory with my music and videos, etc. I mount the former as H: on all my

==================================================================== "What I just said is the fundamental, end-all, final, not-subject-to-opinion absolute truth, depending on where you're standing." Steve Martin ==================================================================== Release Announcements ===================== This is is the latest stable release of Samba

==================================================================== "What I just said is the fundamental, end-all, final, not-subject-to-opinion absolute truth, depending on where you're standing." Steve Martin ==================================================================== Release Announcements ===================== This is is the latest stable release of Samba