search for: derosario

Hi all, several packages print messages during loading. How do I avoid to see them when the packages are in the defaultPackages? Here is an example. With this in ~/.Rprofile ,----[ ~/.Rprofile ] | old <- getOption("defaultPackages") | options(defaultPackages = c(old, "filehash")) | rm(old) `---- I get as last line when starting R: ,---- | filehash: Simple key-value

I have a doubt about the calculation of tests for fixed effects in mixed-effects models. I have read that, except in well-balanced designs, the F statistic that is usually calculated for ANOVA tables may be far from being distributed as an exact F distribution, and that's the reason why the anova method on "mer" objects (calculated by lmer) do not calculate the denominator df nor a

I have encountered a problem with faceted bar plots. I have tried to create something like the example explained in the ggplot2 book (see pp. 126-128): library(ggplot2) mpg4 <- subset(mpg, manufacturer %in% c("audi", "volkswagen", "jeep")) mpg4$manufacturer <- as.character(mpg4$manufacturer) mpg4$model <- as.character(mpg4$model) base <-

.../chippeakanno-package-getAllPeakSequence-problem-tp3830284p3836422.html > Sent from the R help mailing list archive at Nabble.com. > > > > ------------------------------ > > Message: 2 > Date: Fri, 23 Sep 2011 13:07:57 +0200 > From: "Helios de Rosario"<helios.derosario at ibv.upv.es> > To:<jcbouette at gmail.com> > Cc: r-help at r-project.org > Subject: Re: [R] How make a x,y dataset from a formula based entry > Message-ID:<4E7C84AD0200000C0000A9F8 at mailhost.biomec.upv.es> > Content-Type: text/plain; charset=UTF-8 > > To sepa...

A (late) update to this question: On Fri Aug 17 07:33:29, Henrik Singmann wrote: > Hi Diego, > > I am struggeling with this question also for some time and there does > not seem to be an easy and general solution to this problem. At least I > haven't found one. > However, if you have just one repeated-measures factor, use the solution > describe by me here:

Dear R users, I want to announce an update of the package "phia", version 0.2-0, now on CRAN: <http://cran.r-project.org/web/packages/phia/> "phia" (Post-Hoc Interaction Analysis) is aimed at the analysis of the expected values and other terms of in linear, generalized, and mixed linear models, on the basis of multiple comparisons of factor contrasts, and is specially

Dear R users, I want to announce an update of the package "phia", version 0.2-0, now on CRAN: <http://cran.r-project.org/web/packages/phia/> "phia" (Post-Hoc Interaction Analysis) is aimed at the analysis of the expected values and other terms of in linear, generalized, and mixed linear models, on the basis of multiple comparisons of factor contrasts, and is specially

When I fit a multivariate linear model, and the formula is defined outside the call to lm(), the method summary.mlm() fails. This works well: > y <- matrix(rnorm(20),nrow=10) > x <- matrix(rnorm(10)) > mod1 <- lm(y~x) > summary(mod1) ... But this does not: > f <- y~x > mod2 <- lm(f) > summary(mod2) Error en object$call$formula[[2L]] <- object$terms[[2L]]

Hi, I'm afraid that the function read.octave from package "foreign" has some problems with the ASCII data format exported by new versions of Octave (later than 3.2.X). It fails even for a simple case as: [Octave code:] octave:1> x=1; octave:2> save -ascii testdata.mat x [Now in R:] > octavedata <- read.octave('testdata.mat') Mensajes de aviso perdidos In

Hi, I know that the type of model described in the subject line violates the principle of marginality and it is rare in practice, but there may be some circumstances where it has sense. Let's take this imaginary example (not homework, just a silly made-up case for illustrating the rare situation): I'm measuring the energy absorption of sports footwear in jumping. I have three models (S1,

For some time I have been looking for a convenient way of performing post-hoc analysis to Repeated Measures ANOVA, that would be acceptable if sphericity is violated (i.e. leaving aside post-hoc to lme models). The best solution I found was John Fox's proposal to similar requests in R-help: http://tolstoy.newcastle.edu.au/R/e2/help/07/09/26518.html