similar to: Two apparent bugs in aov(y~ *** -1 + Error(***)), with (PR#6523)

I think there are two bugs in aov() that shows up when the right hand side of `formula' contains both `-1' and an Error() term, e.g., aov(y ~ a + b - 1 + Error(c), ...). Without `-1' or `Error()' there is no problem. I've included and example, and the source of aov() with suggested fixes below. The first bug (labeled BUG 1 below) creates an extra, empty stratum inside

I believe you are right, but can you please explain why anyone would want to fit this model? It differs only in the coding from aov(y ~ a + b + Error(c), data=test.df) and merely lumps together the top two strata. There is a much simpler fix: in the line if(intercept) nmstrata <- c("(Intercept)", nmstrata) remove the condition (and drop the empty stratum later if you

When compiling R 1.9.1 on a Linux system with the environment variable LANG=3Dno_NO, I get a similar error as reported in PR#6958: ## FACTS: $ locale LANG=3Dno_NO LC_CTYPE=3D"no_NO" LC_NUMERIC=3D"no_NO" LC_TIME=3D"no_NO" LC_COLLATE=3D"no_NO" LC_MONETARY=3D"no_NO" LC_MESSAGES=3D"no_NO" LC_PAPER=3D"no_NO"

Hi, I'm interested in doing a repeated measures anova using aov. The procedure is nicely described in section 6.7.1, pp. 24-27 of Baron and Li's "Notes on the use of R for psychology experiments and questionnaires," and I've reproduced their example exactly. My own problem is almost identical to theirs: rawdat<-c(1.6530074e+001, 1.2124254e+001, 1.0040371e+001,

As far as I can tell, there are no facilities in the "base" R for image analysis. Is there any extension or package for R that implements image analysis? I'm thinking of features such as o Reading and writing image files, o displaying indexed or "true color" images, o resizing, cropping, etc, of images, o filters for manipulating images (smoothing, edge detection, etc) o

I'm using R for calculations and plots in a statistical paper. Is there a canonical reference to use in the bibliography? (Or is it not common to acknowledge computer programs in journal papers?) -- Bj?rn-Helge Mevik <bhm at math.uio.no> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi, I have a problem with aov(). I used it many times, but now I have new data, tried to use it in the same way. I get I strange Error message that i can't understand (see below). I guess it is caused by incompatible data types. The dataframe also looks ok to me. I use: spk

Bj??rn-Helge Mevik writes: > > tools to make it easy to convert R packages to S-PLUS. > > Not the other way around as well? Actually, we'll be discussing tools to make packages in general that can work in both S-PLUS and R, and also how to make some S-PLUS-only libraries available in an open-source environment as well as working with package authors to port R packages to work in

Ahhrg List, there was a hitch in the previous code. Updated and functional version follows after my initial question: I am trying to randomise a dataframe with mixed factors/numeric variables and return a new (randomised) dataframe with the same columns (as factors/numeric). simgap<-function(x) { ma<-max(x) mi<-min(x) Xout<-runif(length(x),min=mi,max=ma) return(Xout) }

Hello I can see several postings on this mailing-list that people have problem with mounting ext3 partition with mode data=journal. See URL's: https://www.redhat.com/archives/ext3-users/2004-March/msg00000.html https://www.redhat.com/archives/ext3-users/2004-March/msg00050.html We are going to use ext3 on a Compact Flash disk in true IDE mode. We need this filesystem to be as safe and

On Fri, Sep 4, 2015, 9:24 AM Bj?rn-Helge Mevik <b.h.mevik at usit.uio.no> wrote: arnaud gaboury <arnaud.gaboury at gmail.com> writes: > After a few days of reading and headache, I finally gave a try at > building R from source with Intel MKL and ICC. Documentation and posts > on this topic are rather incomplete, sometime fantasist et do not give > much explanations about

Hello I try again. Can anybody of you acknowledge or not if mode data=journal in ext3 is safe to use in Linux kernel 2.6.x? Wee need to have a very consistent and integrity for our filesystem, and it would then be desired to journal both data and metadata. But if this mode can corrupt the filesystem as both Phil White and Nicolas Kowalski has experienced, it may be more advised to use mode

On Mon, Aug 24, 2015 at 11:29 AM, Bj?rn-Helge Mevik <b.h.mevik at usit.uio.no> wrote: > arnaud gaboury <arnaud.gaboury at gmail.com> writes: > >> - Intel MKL: this is part of Intel Parallel Studio and is a paid >> software. Now, there is the MKL package distributed by >> Revolutionanalytics, but I am not certain how this can be distributed >> for free. Is

Hi I have made a R function 'convolve2' for convolution of two real valued vectors based on Rs 'convolve' with option type="open" - see below. (exp.length and irf.length are variables set in another part of the program) I wish to implement the function convolve2 in C and use it in a function used from R with .Call - e.g. I need to call fft in C. All I can find in the

Hi all, I want to label an axis with exponents, but can't get it done with expression. Any hints would be very welcome! # simulated data, somewhat similarly distributed to my real data: set.seed(12); d <- rbeta(1e6, 0.2,2)*150 ; d <- d[d>1e-8] hist( d? , breaks=100) # now on a logarithmically scaled axis: hist(log10(d), breaks=100, xaxt="n") abline(v=

Looks like I can't edit posts. Also of possible interest, I was using v1.0 of the unofficial patch by scient, available here (http://www.weplayciv.com/forums/downloads.php?do=file&id=240). However, I don't suspect that made a difference. If I get round to more thorough testing I'll add a version and some more notes to the appDB. Comments welcome.

Hi, I'm studying math, and i have to make an analysys using Randomised Block Design. I have two factors, i know how to do this in Statistica, but how to do this in R, i read some manuals but the only thing that i have found was 2 factor ANOVA. Please could someone help me, or give some usefull links ?? Krzytsztof Suwada

[sorry if this arrives in duplo: it doesn't show up in the archives and it seems that the address I posted this from originally is no longer functional] Hello, I think I could do with some suggestions concerning the following problem. I have data from a set of experiments on motion sickness where for each subject, I have 1) personal data like age and gender 2) a subjective rating of

Hi, OK, I will describe it in more detail.In the meantime I installed the newest debian package of r-base etc. and the error message changed from:Error in "names<-.default"(`*tmp*`, value = nmstrata) : names attribute must be the same length as the vectorto:Error model is singular in: aov(rati ~ cs * spk + Error(sub/(cs * spk)), data = rating.df) The model is a repeated

Hello I'm very much a beginner on meta-analysis, so apologies if this is a trivial posting. I've been sent a set data from separate experimental studies, Treatment and Control, but no measure of the variance of effect sizes, numbers of replicates etc. Instead, for each study, all I have is the mean value for the treatment and control (but not the SD). As far as I can tell, this forces