similar to: Covarianz matrix from an aov model

The original subject line is false. It *does* reposition, but `origin' was being ignored, so your request was to reposition to the beginning of the file, and that is what happened. This was already fixed in R-patched: from NEWS there: o seek() on file connections was ignoring the `origin' argument. It's always worth looking at the current patched version. On Tue, 29 Jan 2002

Thanks for the advice, folks! Shravan Vasishth <vasishth at ling.ohio-state.edu> Try pico. There's no substitute for/comparison with emacs, though, IMHO. but just to be sure: you don't have to have ESS to use (x)emacs for editing R scripts. Sven Garbade <garbade at psy.uni-muenchen.de> ESS can be installed via rpm on Mandrake. baron at cattell.psych.upenn.edu (Jonathan

Dear Martin, I think that your plan makes sense. It's too bad that aov() behaved differently in this respect from lm(), and thus created more work, but it's not be a bad thing that the difference is now explicit and documented. I expect that that other problems like this will surface, particularly with contributed packages (and I know that you're aware that this has already happened

Hi, I'm confused about how to specify random and fixed factors in an aov() term. I tried to reproduce a textbook example: One fixed factor (Game, 4 levels) and one random factor (Store, 12 levels), response is Points. The random factor Store is nested in Game. I tried > str(kh.df) `data.frame': 48 obs. of 4 variables: $ Subj : Factor w/ 48 levels

How about -- d is a vector as you described it, then define d1 <- d[1:(length(d)-1)] d2 <- d[2:length(d)] d3 <- d1!=0 & d2!=0 & (d1-d2)!=0 I think d3 is the final vector you wanted. I didn't test it, so I'm not sure my logic is right, but I think this is definitely a path you can take... -----Original Message----- From: owner-r-help at stat.math.ethz.ch

Hi all, I computed a two way anova with two within-subjects fixed effects: diffi (2 levels) and cond (4 levels) and one random effect subject, and one dependent variable accuracy: summary(aov(accuracy ~ diffi*cond + Error(sub/(diffi+cond)), data=phd.df)) Error: sub Df Sum Sq Mean Sq F value Pr(>F) Residuals 4 136.549 34.137 Error: sub:diffi Df Sum Sq Mean Sq F value

Hi, I am seeking help with an error when running the example from R Documentation for TukeyHSD. The error occurs with any example I run, from any text book or website. thank you... > plot(TukeyHSD(fm1, "tension")). Error in plot(confint(as.glht(x)), ylim = c(0.5, n.contrasts + 0.5), ...) : error in evaluating the argument 'x' in selecting a method for function

I'm baffled. When I run make check after installing from source, I get a Error 2. From my understanding of how these things work, it would appear to be coming from this (as at the end of base-Ex.Rout.fail: > has.VR <- require(MASS, quietly = TRUE) Attaching package 'MASS': The following object(s) are masked from package:base : confint confint.lm nclass.FD nclass.scott

Full_Name: Sven Sandin Version: 2.0.0 OS: SuSE Linux 9.0 Submission from: (NULL) (81.227.17.135) Have just compiled and installed R-2.0.0.tar.gz running SuSE9.0. The function vcov do not accept "coxph" object as input any longer. The same R-program running R1.9.1 do work. R-program attached below. Exporting the coxph object from R2.0.0 to R1.9.1 I get vcov ouput in R1.9.1. Exporting

Hi all, how do I change the magnification of xaxis annotation? I need a smaller axis text size in some of my plots. I tried cex.axix=0.5 in plot(), but this doesn't work. Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hello everyone, I've conducted a Type III repeated-measures ANOVA using Anova() from the car package, based on the suggestions at http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/(option 3) and http://languagescience.umd.edu/wiki/EEG#ERP_ANOVA_in_R. My ANOVA has two factors: Condition (3 levels) and Region (6 levels) and their interaction. Below is code to run the Anova

Hi there, some statistical programs (e.g. SPSS) calculate a correction of the degrees of freedom in a repeated measure analysis of variance (see Greenhouse-Geisser (1958) or Huynh-Feld (1976)) by a factor epsilon. This factor is used to correct the deg. of freedom to get a corrected f-test. Is this also possible with R? Thanks, Sven P.S.: I read in the lm help page: singular.ok logical,

Hi list, can I extract the cov.unscaled ("the unscaled covariance matrix") from a gls fit (package nlme), like with summary.lm? Background: In a fixed effect meta analysis regression the standard errors of the coefficients can be computed as sqrt(diag(cov.unscaled)) where cov.unscaled is (X'WX). I try do do this with a gls-fit. Thanks, Sven

Achim Zeileis writes: > Sven Garbade wrote: > > > > Hi, > > > > can I extract the p-values after a pairwise comparision with > > pairwise.t.test()? The object of class "pairwise.htest". > > This is a list with one component being "p.value". You can extract it as > usual with > pairwise.t.test()$p.value Sorry, I

Some days before I asked about a recommended linux environment for R. Looks like, in principle it runs on all major linux distributions (debian red hat, suse). There was a report of a problem with dynamically loading Fortran on Red Hat 6.2. Someone encountered trouble compiling library mclust with f77 and succeded with f2c, but other people have succeeded compiling mclust with f77, so there is

Hi, I am a method developer in drug discovery. I have developed a similarity searching method for small (drug-like) molecules. I want to compare the performance of my method to the performance of other methods. Similarity searching methods are commonly assessed by their ability to (re)discover a set of molecules that is avtive versus a given target, given one randomly selected query molecule.

Dear R Users, I want to write a function that applies to the dataframe and variables that were used in a previous call to lm or aov. In order to do this, I need to write a function that applies to the output of lm or aov, and yields the names of the dataframe and variables that were used in the lm or aov analysis. For example, suppose that I give the command: aov.out <- aov( Rt ~

<ripley at stats.ox.ac.uk> wrote: > Hint 2: in the absence of balance, ...., and lme can do that Would it be possible to make aov like wrappers to the various special lm variants allowing for a uniform syntax for anova? aov(resp~f1*f2+Error(S/(f1*f2))) ## uses lm aov.lme(resp~f1*f2+Error(S/(f1*f2))) ## uses lme aov.rlm(resp~f1*f2+Error(S/(f1*f2))) ## uses rlm ... I'd do it, but I

Hi all, I don't understand the behavior of is.integer(): > x <- integer() > is.integer(x) [1] TRUE > x <- 10 > is.integer(x) [1] FALSE > x <- 1:10 > is.integer(x) [1] TRUE Why is.interger() returns FALSE if x has only one element? And how can someone check if x is an integer but contains only one value? (R 1.5.0 on Linux i386) Thanks, Sven

>>>>> Fox, John <jfox at mcmaster.ca> >>>>> on Thu, 14 Sep 2017 13:46:44 +0000 writes: > Dear Martin, I made three points which likely got lost > because of the way I presented them: > (1) Singularity is an unusual situation and should be made > more prominent. It typically reflects a problem with the > data or the