[R] Lattice bwplot: Adding mean with panel.points fails in vertical arrangement

2012-10-19 Thread Henning Wildhagen
Dear collegues, given a structure of data like this: ## Data ### set.seed(100) a <- rnorm(60,10,3) s <- c(rep("A",20),rep("B",20),rep("C",20)) p <- c(rep("d",6),rep("e",6),rep("f",6),rep("g",6),rep("h",6)) df <- data.frame(a,s,p) i would like to draw a lattice bwplot

[R] ggplot2: Controlling line width of panel borders

2010-12-06 Thread Henning Wildhagen
Dear R-users, i encountered some problems when trying to adjust the line width of the axes and stripes in a plot created with ggplot2. I use the "barley" dataset of the lattice package to illustrate my problem: library(ggplot2) library(lattice) barley[["SD"]] <- 5 limits <- aes(ymax=barley$yie

[R] Error bars in lattice barchart with groups

2010-11-30 Thread Henning Wildhagen
Dear R-users, i want to plot gene regulation data in a lattice barchart. To illustrate the problem i encounter, the following code uses the "barley"dataset: library(lattice) barley[["SD"]] <- 5 PLOT <- barchart(data=barley, yield~variety|site, groups=year,origin=0, as.table=TRUE, scales=list(

[R] Concatenating two vectors into one

2009-05-18 Thread Henning Wildhagen
Dear users, a very simple question: Given two vectors x and y x<-as.character(c("A","B","C","D","E","F")) y<-as.factor(c("1","2","3","4","5","6")) i want to combine them into a single vector z as A1, B2, C3 and so on. z<-x*y is not working, i tried several others function, but did not get to

[R] postscript()-resolution

2009-05-08 Thread Henning Wildhagen
Dear users, another question concerning graphics for publications. My favourite journal wants .eps-graphics, and from older postings i adapted the following code: postscript(file="Figure1.eps", title="Figure 1", width=11.5, height=8, paper="a4",onefile=FALSE) However,  when checking the proper

Re: [R] Tinn-R pdf()

2009-04-08 Thread Henning Wildhagen
nning > > Original-Nachricht > Datum: Wed, 8 Apr 2009 12:29:58 +0200 > Von: "ONKELINX, Thierry" > An: "Henning Wildhagen" , r-help@r-project.org > Betreff: RE: [R] Tinn-R pdf() > > Dear Henning, > > You need to print() lattice plo

[R] Tinn-R pdf()

2009-04-08 Thread Henning Wildhagen
Dear R and Tinn-R users, i recently switched to Tinn-R and sending code to R works fine (R 2.8.1, Tinn-R 2.2.0.2, OS Windows XP). However, i encountered a problem when trying to send plots to pdf files like this: library(lattice) pdf("plot1.pdf") PLOT<-(xyplot, ...) PLOT dev.off() The file "p

[R] Updating packages under R 2.8.1

2009-01-16 Thread Henning Wildhagen
Dear users, i just installed the lastest version of R, 2.8.1 on my computer (OS Windows XP). Then i tried to update the packages copied from my old R version by >update.packages(ask=F) However i get the following warning: "Warning: unable to access index for repository http://cran.ch.r-projec

[R] (no subject)

2009-01-16 Thread Henning Wildhagen
Dear users, i just installed the lastest version of R, 2.8.1 on my computer (OS Windows XP). Then i tried to update the packages copied from my old R version by >update.packages(ask=F) However i get the following warning: "Warning: unable to access index for repository http://cran.ch.r-projec

Re: [R] xyplot strip=function for two conditioning variables

2008-08-05 Thread Henning Wildhagen
change the side of the vertical strips? Cheers, Henning > > Original-Nachricht > Datum: Mon, 4 Aug 2008 12:25:06 -0700 > Von: "Deepayan Sarkar" <[EMAIL PROTECTED]> > An: "Gabor Grothendieck" <[EMAIL PROTECTED]> > C

[R] xyplot strip=function for two conditioning variables

2008-08-04 Thread Henning Wildhagen
Dear list, for a data structure like in df: set.seed(100) Treatment<-rep(c("Nitrogen","Carbon", "Sulfur"),each=9) week<-rep(c(1,5,9),3,each=3) genes<-rep(c("18s", "EF1b", "NR"),9) copies<-rnorm(27, 100,40) df<-data.frame(Treatment,week,genes,copies) i wrote this code for a xyplot: libra

[R] axis.break on Date-x-axis in lattice xyplot

2008-07-28 Thread Henning Wildhagen
Dear list, i am using the following code to produce  a lattice xyplot, but the axis.break-function is seemingly not executed. Date<-seq(as.Date("2006-08-29"), as.Date("2007-08-28"), by="2 weeks") Period<- var1<-rnorm(27, 9000, 3000) var2<-rnorm(27, 5,25000) var3<-rnorm(27, 10

Re: [R] t-test for multiple variables

2008-07-18 Thread Henning Wildhagen
ders <[EMAIL PROTECTED]> > An: Henning Wildhagen <[EMAIL PROTECTED]> > CC: r-help@r-project.org > Betreff: Re: [R] t-test for multiple variables > > Hi Henning, > > Henning Wildhagen wrote: > > Given a dataframe of three columns, where col1 is of

[R] t-test for multiple variables

2008-07-18 Thread Henning Wildhagen
Dear list, i try to do my first statistics with R. Given a dataframe of three columns, where col1 is of type factor, col 2 and col3 are numeric and pairs of observations i would like to perform a paired t-test for each level of col1. I would like to avoid specifying the levels of col1 manually

Re: [R] barchart with bars attached to y=0-line

2008-07-17 Thread Henning Wildhagen
ried with npc, but also did not work)? Thanks a lot for your help, Henning > > Original-Nachricht > Datum: Wed, 16 Jul 2008 09:48:35 -0700 > Von: "Deepayan Sarkar" <[EMAIL PROTECTED]> > An: "Henning Wildhagen" <[EMA

[R] barchart with bars attached to y=0-line

2008-07-16 Thread Henning Wildhagen
Dear R users, i am using the following code to produce barcharts with lattice: Compound<-c("Glutamine", "Arginine", "Glutamate", "Glycine", "Serine", "Glucose", "Fructose", "Raffinose", "Glycerol", "Galacglycerol", "Threitol", "Galactinol", "Galactitol") Class<-c("aminos","aminos","aminos","am

[R] xYplot customizing y-axis scaling

2008-07-10 Thread Henning Wildhagen
Dear list, using the packages Hmisc and lattice i produced some nice xYplots. However, since the data range of the conditioning variable is very big, i need to define more than one y-scale for the plot (in some panels you just see a flat line of data points very close to the x-axis), e.g. diffe

[R] ggplot2: scaling and tick mark of x-axis

2008-07-03 Thread Henning Wildhagen
Dear list, i am trying to plot data in a way like this: Date<-seq(as.Date("2006-08-29"), as.Date("2007-08-28"), by="2 weeks") var<-rnorm(27, 1, 2500) err<-rnorm(27, 3000, 1000) df<-data.frame(cbind(Date, var, err)) library(ggplot2) g <- ggplot(df, aes(x=Date, y=var, min=(var - err), max=(

[R] Problem with "date"-x-axis in lattice

2008-06-20 Thread Henning Wildhagen
Hello list, i am trying to plot a continous variable y against a "date" variable, both in one dataframe named "df", using a code like this library(lattice) plot1<-xyplot(y~date, data=df, type="b") "date" is of class "Date", of course, format="%d.%m.%Y", and spans two calendar years. The proble