similar to: subset, once more

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

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

I upgraded to gtk+-1.2.10 and try to recompile R. This gtk release has its headers in ${prefix}/include/gtk-1.2/g{t,d}k/ (as can be seen from gtk-config --cflags). The gnome module build doesn't seem to account for this, and thus can't find gtk.h &sui. Is this a known issue (aka bug)? RenE J.V. Bertin College de France/LPPA 11, place Marcelin Berthelot 75005 Paris, France

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

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

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

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

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

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

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

Hello - and happy newyear to all of you! I've got some data that I'm plotting with bwplot, a 3x2x3 design where the observable decreases with the principle independent factor, but at different rates. I'd like to get lattice to impose not a single set of axes ranges identical for all panels, but ranges that are identical for each panel row or each column. Effects will stand out much

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

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

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

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