Search Archives

Hi Michael, > I have a nested anova, with random factor "lakefac" nested within > factor "fishfac" (fixed), with an additional fixed factor > "Habfac". If I consider everything as fixed effects, it's addmittedly > not the correct model, but I can at...

> But I recently realized something. Most of the variables that I've > tested as fixed effects are properties of the subject (e.g. Race, > Gender, etc.). Is it correct to be using a random effect > Subject that is nested within...

Dear Emmanuel, Thank you! Stupid me ... A nice weekend to all of you! Lorenz > Am 12.09.2014 um 17:55 schrieb "Emmanuel Curis" <emmanuel.curis at parisdescartes.fr>: > > Double check your results, you will see that there is agreement...

Dear R users, my question is about groupwise modeling of spatial autocorrelation in mixed effect models: In a soil study, five different sites are to be compared respect to specific soil properties. A regular sampling grid with 32 points in...

Hi Ignacio, following this link: http://www.ne.su.se/polopoly_fs/1.216115.1426234213!/menu/standard/file/clustering1.pdf you can download a documentation for Arai's cl-function that is mentioned in the link in your email. I...

Hi Chris, > I am trying to fit the following linear model to logged per capita > fecundity data (ie number of babies per female) for a mouse: > > RsNRlS <- glm(formula = ln.fecundity ~ summer.rainfall + N + > lagged.rainfall + season, ....) > > I am using...

The two models that you compare differ in their fixed effects. As far as lme had been concerned, the anova gave a warning because one should not compare two models differing in their fixed effects if they had been estimated...

Dear Michael, This is possibly not what you should do. As has been recently suggested to me on this list, you could use bootMer with a function that computes the sum(s) that you are interested in. Then use boot...

Dear all, in a current project I have a pretty huge .Rdata. Thus I was working with R --vsize 100 --nsize 1000000. Today when I tried to restart R I get the following error message: Error: a read error occured...

Thanks for the details, I suspected something like that. I think that begs the question: what is the meaning of quantile regression through the origin? If the tau=.5 line does not pass through 1/2 the data how do...

> I have been going over and over the examples in MASS > and the Pinheiro and Bates example, but cannot get my model > to run correctly with either aov or lme. > > Could someone give me a hand with the correct model...

Dear Ben, Many thanks for your input. [... snip] > > Now, apart from this aspect, can confint be tweaked to calculate not > > only the confidence interval of the 'raw' parameters but also for > > some function of the parameters? If not, do I...

Dear Martin, Many thanks for this explanation which, of course, is very reasonable ;-) But - and I may be real slow on this - why is the same seemingly not true for the random effects as well (summary and confint give the...

Dear Edzer Pebesma, thank your for your advice - you were right that the semivariance values are standardized by default in the output of the nlme-Variogram function. So, with the adequate "resType"-option in that function, I could reproduce nearly...

Dear all, based on Bens suggestion (and it has come up again in another post) I made a first try in using bootMer and boot.ci. > [... snip] > > > > Now, apart from this aspect, can confint be tweaked to calculate not > > > only...

Dear all, I am currently evaluating different models regarding haemodynamic reactions in the brain. Thus the realtive concentration of O2Hb and HHb are my response variables. These concentrations are recorded over relatively short stimulation periods (75 sec) and are presently...

Dear all, I am trying to use ace on a 446 x 14 Matrix, with two categorical variables in the first two columns. If I try to run it I get: > aceiH <- ace (data.mat, resp, cat= 1:length(x...

> 1. The first I used to do in SPSS and I would like to be able > to do it in R as well. > This is the hierarchical model I would like to use: a continuous > variable explained by factor A...

All, So that suggests that .GlobalEnv[["X"]] is more efficient than get("X", pos=1L). What about .GlobalEnv[["X"]] <- value, compared to assign("X", value)? Dave On Wed, Dec 3, 2014 at 3:30 PM, Peter Haverty <haverty.peter at...

Dear all, In a current project (2 x 2 x 2 factorial design) we are interested in calculating p-values for binary outcomes (we are aware that such an approach is not unequivocal but circumstances are such that results will...