similar to: summary(aov(...)) into a string?

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

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

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

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

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

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

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

I'm trying to write a function to automate doing a variance analysis, part of which involves doing some further calculations. The method I've been using isn't very robust, if variable names change then it stops working. For this dummy data > dput(assayvar,"") structure(list(Run = structure(c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 4L), .rk.invalid.fields =

Hi, I don't quite understand what are the return values of aov. I know that it has 'coefficients'. But I need to know what all the other return values are. Can somebody let me know how to figure them? Value: An object of class 'c("aov", "lm")' or for multiple responses of class 'c("maov", "aov", "mlm",

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

Hi all, I search the archive for finding a simple solution for using TukeyHSD with a multistratum aov result (a repeated emasure anova). The Question have been asked but I've found no clear answer. res<-aov(y~Fa*Fb+Error(Subject/(Fa*Fb)) ) I think that the problem is that res is an aovlist object instead of the "aov" object required by TukeyHSD. Is there an easy solution to

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

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

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

Hello, I am trying to print out ANOVA summaries --- as returned by summary (aov(...))--- in tab-delimited format, either to a file or to the console, but so far I haven't been able to figure out a solution. My particular ANOVA has repeated measures, so it is an instance of class aovlist. Can anyone offer any suggestions? I am relatively new to R, so please forgive me if I am missing

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

Dear all, Type III SS time again. This case trying to reproduce some SPSS (type III) data in R for a repeated measures anova with a betwSS factor included. As I understand this list etc, if I want type III then I can do library(car) Anova(lm.obj, type="III") But for the repeated measures anova, I need to include an Error-term in the aov() call (Psychology-guide from Jonathan Baron)