similar to: gtk+1.2.10

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

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

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

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

Hi, I'd like to add lines to the panels of a bwplot, showing (and connecting) the means of the displayed distributions. I didn't find anything in the documentation or the examples that suggests how to accomplish this easily. Has anybody done this already, and willing to share? Thanks, RenE J.V. Bertin College de France/LPPA 11, place Marcelin Berthelot 75005 Paris, France

>From the documentation to the quantile function: ... `quantile(x,p)' as a function of `p' linearly interpolates the points ( (i-1)/(n-1), ox[i] ), where `ox <- order(x)' (the ``order statistics'') and `n <- length(x)'. This gives `quantile(x, p) == (1-f)*ox[i] + f*ox[i+1]', where `r <- 1 + (n-1)*p', `i <- floor(r)', `f

Here's the patch. It looks like it does what I want it to do, but I haven't yet tried it with multiple-line entries and the like. Please let me know if anything turns up that I should look into (I attached it too, in case hotmail wraps). RenE J.V. Bertin College de France/LPPA 11, place Marcelin Berthelot 75005 Paris, France ## A couple of 'diff -c' patches to modify the

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

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

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

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

You could use cut. The key calculation would be: w <- .05; eps <- 1e-5 breakpoints <- seq(min(kk), max(kk), .05) breakpoints <- floor( (breakpoints + (w/2) + eps) / w) * w values <- cut(kk, c(breakpoints, Inf), right = FALSE) values <- ordered(values) If you don't like the labels produced add lab = breakpoints as a cut arg. On 11/10/05, RenE J.V. Bertin

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

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

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

hello, sorry for posting what may be a simple question: i do have a matrix of coordinates (positional judgments, see below) and now want to calculate and plot the corresponding error ellipse. can anyone help me with the exact steps/syntax? > xyDat X Y 1 -0.49 -2.13 2 0.91 0.48 3 0.20 -2.80 4 -0.76 -3.23 5 -0.36 2.50 6 1.38 1.24 7

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

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

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

Negative numbers in (encapsulated) postscript generated by R loose their sign when opened in Adobe Illustrator 8 (french). This was so in 1.5.x and remains in 1.6.0. I suppose that this is a known issue, although it is rather a flaw in Illustrator than in R, if my analysis of it is correct. I sent a lengthy mail containing a solution (ISOLatin1Encoding has to be redefined) to the r-help mailing