similar to: summary() dispatch puzzle

>>>>> Ben Bolker <bbolker at gmail.com> >>>>> on Fri, 15 Jul 2016 16:45:50 -0400 writes: > I'm sorry I haven't boiled this down to a more-minimal example yet, > but ... > I'm working on an S3 method (tidy.merMod, in the 'broom' package). It > normally handles 'merMod' objects from the lme4 package,

Hi, The R CMD check is successful in R 3.0.1 but fails to install package lmerTest under R-devel r63910, Here is what I get: ** preparing package for lazy loading Error in reconcilePropertiesAndPrototype(name, slots, prototype, superClasses, : no definition was found for superclass "merMod" in the specification of class "merModLmerTest" In DESCRIPTION file I have:

(This time with the r-help in the recipients...) Be careful when mixing lme4 and lmerTest together -- lmerTest extends and changes the behavior of various lme4 functions. >From the help page for lme4-anova (?lme4::anova.merMod) > ?anova?: returns the sequential decomposition of the contributions > of fixed-effects terms or, for multiple arguments, model >

Please reread my point #1: the tests of the (individual) coefficients in the model summary are not the same as the ANOVA tests. There is a certain correspondence between the two (i.e. between the coding of your categorical variables and the type of sum of squares; and for a model with a single predictor, F=t^2), but they are not the same in general. The t-test in the model coefficients is simply

I have no problem setting up my mixed model, or performing anova or lsmeans on my model?s outputs. However, performing lsd mean separation is giving me fits. So I do not have a problem when using two-way anova model. When using the code: fit.yield.add <- lm(data = ryzup, Yield ~ Rep + Nitrogen + Treatment) LSD.test(fit.yield.add, trt = "Nitrogen", alpha = 0.1, console = TRUE)

> On 30 Sep 2017, at 14:22 , Robert Baer <rbaer at atsu.edu> wrote: > > > > On 9/29/2017 3:37 PM, Rolf Turner wrote: >> On 30/09/17 07:45, JLucke at ria.buffalo.edu wrote: >> >> <SNIP> >> >>> >>> The conceptual paradigm for R is only marginally commensurate with >>> that of >>> standard statistical

Dear All, I'm trying to use waldtest to test poolability (parameter stability) between two logistic regressions. Because I need to use robust standard errors (using sandwich), I cannot use anova. anova has no problems running the test, but waldtest does, indipendently of specifying vcov or not. waldtest does not appear to see that my models are nested. H0 in my case is the the vector of

> On Oct 2, 2017, at 2:05 AM, David <dasolexa at hotmail.com> wrote: > > Dear list, > > > I am just starting on analysis of count data in R 3.4.0. My dataset was obtained from counting particles on a surface before andd after a cleaning process. The sampling positions on the surface are pre-defined and are the same before and after cleaning. I have ~20% of 0's. I

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 On 15-03-23 12:55 PM, Thierry Onkelinx wrote: > Dear Ben, > > Last week I was struggling with incorporating lme4 into a package. > I traced the problem and made a reproducible example ( > https://github.com/ThierryO/testlme4). It looks very simular to > the problem you describe. > > The 'tests' directory contains

hi, i am running R with JGR and Deducer under Ubuntu Mint Lucid. since i updated to R 2.12 i can no longer make graphs - R just crashes. the full error report was posted to launchpad (https://bugs.launchpad.net/ubuntu/+source/openjdk-6/+bug/675905) without any success so far. anyone with a similar problem or helpful suggestions? thanks! kat

Diverted from ESS-help. It is really about R's LANGUAGE, locales and then a wish for sessionInfo() : >>>>> Sparapani, Rodney <rsparapa at mcw.edu> >>>>> on Thu, 20 Feb 2014 19:39:33 +0000 writes: >> Just guessing: perhaps this is something that is set when >> R is initialized and not queried every time something is >>

Hi R community! I've fitted three mixed-effects regression models to a thousand bootstrap samples (case-resampling regression) using the lme4 package in a custom-built for-loop. The only output I saved were the inferential statistics for my fixed and random effects. I did not save any output related to the performance to the machine learning algorithm used to fit the models (REML=FALSE).

The attached code produces corrupted pdfs (test2.pdf, test4.pdf and test5.pdf). The resulting pdf depends on how many cores are available on the machine. I don't see why there should be any difference between the pdfs (exept for the timestamp). Doing many operations involving mclapply can increase the size of the resulting pdf by ten times! Thank you for checking this. require(multicore)

Dear David, thanks ever so much for your answer. Do you mean predicting the original values based on the fitted model and then comparing observed vs. predicted by, for example, a scatterplot? Thanks, David ________________________________ De: David Winsemius <dwinsemius at comcast.net> Enviado: lunes, 2 de octubre de 2017 18:18:36 Para: David Cc: R-help Asunto: Re: [R] Help on

Hi, Let's say one has an object with multiple classes, and a generic function to apply to it has associated S3 methods for more than one of those classes. Further, the method it chooses (I do not know how; some order in the class vector?) is not the suitable one and it produces an error. Would there be some way to make the generic function to choice the correct method, or in case that for

dear members, I am using arfima() from forecast package to model a time series. The following is the code: > LYGH[[202]] [1] 45.40 3.25 6.50 2.15 > arfima(LYGH[[202]]) Error in .fdcov(x, fdf$d, h, nar = nar, nma = nma, hess = hess, fdf.work = fdf$w) : NA/NaN/Inf in foreign function call (arg 5) I tried viewing .fdcov() with the following code:

On 9/29/2017 3:37 PM, Rolf Turner wrote: > On 30/09/17 07:45, JLucke at ria.buffalo.edu wrote: > > <SNIP> > >> >> The conceptual paradigm for R is only marginally commensurate with >> that of >> standard statistical software. >> You must immerse yourself in R to become proficient. > > Fortune nomination. For newer list members wondering what

Hi I haven't found a way to produce a tabulation from factor data with NA values using xtabs. Please find a minimal example below, it's also on R-pubs [1]. Tested with R 3.1.2 and R-devel r67720. It doesn't seem to be documented explicitly that it's not supported. From reading the code [2] it looks like the relevant call to table() doesn't set the "useNA"

Hello, I've problems installing several packages in my R on Fedora 27 64 bit. I found out that it has to do something with a missing compiler (libgfortran.so.3, see below). It works if I downgrade the current version of libgfortran to the specified version by downloading libgfortran-6.2.1-2.fc25.x86_64.rpm and manually installing it. However, I don't want to mess up my system,

>>>>> akshay kulkarni >>>>> on Wed, 31 May 2023 20:55:33 +0000 writes: > dear members, > I am using arfima() from forecast package to model a time > series. The following is the code: >> LYGH[[202]] > [1] 45.40 3.25 6.50 2.15 >> arfima(LYGH[[202]]) > Error in .fdcov(x, fdf$d, h, nar = nar, nma = nma,