similar to: as.double( factor(something) )??

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

RenE J.V. Bertin wrote: > I find myself doing lots of tests like > > >>subset( data, symptoms=='a' | symptoms=='e' | symptoms=='z' .... ) > > > with symptoms one of the factors contained in the data frame. > > and I wonder if there is not an existing operator or function which implements this sort of repeated conditional in a more

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

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() >

On Sun, 10 Oct 2004, RenE J.V. Bertin wrote: > Could I ask some hints/help in interpreting the following ANOVA results, > please? This concerns an experiment where I study the incidence and > severity of motion sickness. I have Sickness.norm, a subjective > discomfort/sickness estimate, normalised to 0..1, the session time T > (normalised to 0..1 and binned in 0.2 wide bins) and a

>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,

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

Default arguments are evaluated in the function frame, not in the calling environment (nor in the same place as explicit arguments). > Which to me reads that a with statement as above is equivalent to > > > attach(data) ; aov.SS1(y=Obs) ; detach(data) > > Or is that just wishful thinking?? The latter. On Thu, 7 Oct 2004, RenE J.V. Bertin wrote: > Hello, > >

On Thu, 16 Sep 2004, RenE J.V. Bertin wrote: > Dear all, > > A few quick questions about interpreting and evaluating the results of > linear regressions, to which I hope equally quick answers are possible. > > 1) The summary.lm method prints the R and R^2 correlation coefficients > (something reviewers like to see). It works on glm objects and (after > tweaking it to

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") >

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),

Starting html help in the current version of R has a very annoying side-effect. It indiscriminantly removes $HOME/.R, and replaces it with a virgin copy. I discovered that when all of a sudden I got complaints about my startup "library" not being found. Below is a modified version of the script that doesn't do this. It is not perfect yet (it shouldn't try to recreate links

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

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

New to R, I had the bad idea to send a bug report about '[' not knowing it had a drop= argument. Now, I wonder about the absence of this argument in subset... In both availabe methods (see below), there is a ... argument, but this argument is not used in either. Rather, subset.data.frame explitictly passes drop=F in 1 instance. Before I start patching (for my own use): what is the

Hello, This is a usage question for others with experience of R under the Aqua Mac OS X interface. Basically, I don't succeed in importing PDF files (created with the 'Save As' menu to Quartz device windows) into Illustrator. Versions up to (and including?) 10 loose paths (lines/polygons) and or fill them in black. Version CS (11) on the Mac imports almost correctly, but, "to

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

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

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