similar to: lattice question: independent per-row or per-column scaling?

Hello, I have some data in a data.frame where for each of a number of subjects, I have scores for all of a number of symptoms. Subjects are subdivided in a number of groups, which have unequal sizes. I'd like to plot between-group differences in the scores on the various symptoms. Ideally, that would be in a form as would be produced by > bwplot( Score~Symptom ) but I'm not sure

RenE J.V. Bertin wrote: > Hello, and (somewhat belated) best wishes for 2005. > > Can one order not to draw outliers in bwplot, or at least exclude them from the vertical axis scaling? If so, how (or what doc do I need to consult)? > The options that have this effect in boxplot() do not appear to have any effect with bwplot (although outline=FALSE in boxplot does *not* change the

On Friday 17 September 2004 13:52, RenE J.V. Bertin wrote: > Hello again, > > I am doing regressions (using panel.lmline() (and panel.abline( > rlm(...))) ) inside a panel method which I pass to bwplot(). > > What I would like to do is create a boxplot of categorised data > (binned on the independent variable), and superpose a regression line > which is calculated using the

>From: Patrick Connolly <p.connolly@hortresearch.co.nz> >To: "RenE J.V. Bertin" <rjvbertin@hotmail.com> >Subject: Re: [R] changing (core) function argument defaults? >Date: Thu, 18 Nov 2004 11:43:10 +1300 > >On Wed, 20-Oct-2004 at 07:48PM +0200, RenE J.V. Bertin wrote: > >|> Hello,

>From: Patrick Connolly <p.connolly@hortresearch.co.nz> >To: "RenE J.V. Bertin" <rjvbertin@hotmail.com> >Subject: Re: [R] changing (core) function argument defaults? >Date: Thu, 18 Nov 2004 11:43:10 +1300 > >On Wed, 20-Oct-2004 at 07:48PM +0200, RenE J.V. Bertin wrote: > >|> Hello,

Let's say I have a within-subject experiment with 2 observables, obs1 and ob2 and 2 independent factors, fac1 and fac2. I can do summary( aov( obs1~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) summary( aov( obs2~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) to test the 2 observables separately. > summary( fit<-manova( cbind(obs1,obs2)~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) gives

On 6/23/05, RenE J.V. Bertin <rjvbertin at gmail.com> wrote: > Hello, > > I was just having a look at the aov function source code, and see that when the model used does not have an Error term, Helmert contrasts are imposed: > > if (is.null(indError)) { > ... > } > else { > opcons <- options("contrasts") >

This may be of interest: http://www.ijs.si/software/snprintf/ It looks like this version may well be more complete and conform the standards than the version that comes with R. BTW: I am currently patching unix/sys_std.c and modules/gtkconsole.c such that they don't store lines in the history that are identical to the previous line. Is there any interest in posting those patches here? RenE

Concerns: R 1.5.0 gui version, Windows (downloaded binary) and Linux (installed from sources). # Load the data from the attached file: kk<-read.table("__filename__", header=1) # attach the data: attach(kk) Snr<-factor(Snr) # fool around with a call to anova.glm(): anova.glm( aov( nFD~Type+size+Modality+Error(Snr/(Type+size+Modality)) ) ) # Error: object nFD not found # Well, I

Rene, Look at ?format. Sean On Sep 17, 2004, at 9:21 AM, RenE J.V. Bertin wrote: > Hello, > > I can't seem to find the way to modify the precision with which > paste() prints its floating point numbers, more precisely the number > of decimal digits printed. This is apparently not controlled by > options( digits= ), and there is no appropriate argument to paste() >

It doesn't print anything: the summary.aov (or summary.aovlist) print method does. ?summary.aov tells you the structure of the objects they return. On Tue, 15 Feb 2005, RenE J.V. Bertin wrote: > I'd like to annotate a plot with the output of summary(aov(model)), > ideally just with the significant effects. I don't find a means to > redirect what that command prints into

Dear expeRts, How can I get ... (1) different y-axis scales for each row (2) while having the same y-axis scales for different columns? I coulnd't manage to do this with relation="free" [which gives (1) but not (2)]. I also tried relation="sliced", but it did not give the same y-axis scales within each row (see the fourth row). Further, it "separates" the

Take a look at package mgcv. Hope this helps. --Matt -----Original Message----- From: RenE J.V. Bertin [mailto:rjvbertin at despammed.com] Sent: Thursday, February 27, 2003 1:39 PM To: r-help at stat.math.ethz.ch Subject: [R] multidimensional function fitting Hello, I have been looking around for how to perform a multidimensional, arbitrary function fit (in any case non-linear; more below),

Well... this may be a red herring after all, but it's an interesting one! It turns out I have got DYLD_LIBRARY_PATH set. This gets exported to R (or rather, affects the loader that loads the R process?), and somehow mixes up the way symbols are exported. Something that applies to launching R from the command line (/usr/bin/R) AND via LaunchServices (open -a R). Once I noticed that invoking

Suppose I have a factor size with levels "small", "medium" and "large". Then, when I subset this factor: >ss<-size[size!="medium"] to get at the extremes, >levels(ss) .... Levels: large medium small The same happens with >subset( size, size!="medium") I understand that the resulting factor inherits the possible levels from its

Why are my coefficients getting appended with a 1? It borks a match I do later against the original list that doesn't have the random 1 added to the end. > linearModel[[1]] Call: lm(formula = modelSource ~ +UNITBUILD + UNITDB + ITBUILD + ITDB + UATBUILD + UATDB + HOGANCODE + RCF + ReleaseST1 + ReleaseST2 + ReleaseBLA + Small.Bank.Acquisitions + HLY.NewYear + HLY.MLK + HLY.PRES +

Dear List, I know this has been covered before, but I don't seem to be able to get it right. I am constructing a boxplot in lattice and can't get the scales in the correct alphebetical order. I have already read that this is due to the way factors are treated, and I have to redefine the levels of the factors. However, I have failed. As a simple example: library(lattice)

I have a Lattice plot in which I want to adjust the number of tick marks used, and I want to have the drawn grid reflect that change. Here is what I'm doing: bwplot(var1 ~ var2, data=df, scales=list(tick.number=10), panel=function(...) { panel.grid(h=0,v=-1,...); panel.stripplot(col="gray40", pch="|", cex=2, ...); panel.bwplot(...); })

Dear All, There is no advantage of GUI over CLI, IMO. The real issue is the answer to the questions: "What should I do next?" or "What am I allowed to do here?" A "nice" interface, not necessarily GUI, will offer friendly answers: "I was expecting you to do _this_" or "In this situation you are allowed to do _these things_" You see, it's all

There is an incompatibility between R's configure/build process, and gcc 3.2.0 (and maybe earlier). Gcc 3.2.0 includes -I/usr/local/include by default, on my system (default configuration), and cpp0 issues a warning when this flag is issued additionally by the user (= R's configure). This warning ends up in the dependency sections of the Makefile that are generated