similar to: Error in anova(): objects must inherit from classes

Hello: Will someone PLEASE help me with this problem. This is the third time I've posted it. When I appply anova() to two equations estimated using glmmPQL, I get a complaint, > anova(fm1, fm2) Error in anova.lme(fm1, fm2) : Objects must inherit from classes "gls", "gnls" "lm","lmList", "lme","nlme","nlsList", or

Hello: I am really stuck on this problem. Why do I get an error message with anova() when I compare these two equations? Hope someone can help. ANDREW ____________________________ > fm1 <- glmmPQL(choice ~ day + stereotypy, + random = ~ 1 | bear, data = learning, family = binomial) > fm2 <- glmmPQL(choice ~ day + envir + stereotypy, + random = ~ 1 |

Hello: I am getting an error message when appplying anova() to two equations estimated using glmmPQL. I did look through the archives but didn't finding anything relevant to my problem. The R-code and results follow. Hope someone can help. ANDREW ____________________________ > fm1 <- glmmPQL(choice ~ day + stereotypy, + random = ~ 1 | bear, data = learning, family =

Hi, I'm having some problems regarding the packages lme4 and nlme, more specifically in the denominator degrees of freedom. I used data Orthodont for the two packages. The commands used are below. require(nlme) data(Orthodont) fm1<-lme(distance~age+ Sex, data=Orthodont,random=~1|Subject, method="REML") anova(fm1) numDF DenDF F-value p-value (Intercept) 1

Dear listers, Here we have a strange result we can hardly cope with. We want to compare a null mixed model with a mixed model with one independent variable. > lmmedt1<-lme(mediane~1, random=~1|site, na.action=na.omit, data=bdd2) > lmmedt9<-lme(mediane~log(0.0001+transat), random=~1|site, na.action=na.omit, data=bdd2) Using the Akaike Criterion and selMod of the package pgirmess

My question concerns the logic behind hypothesis tests for fixed-effect terms in models fitted with lme. Suppose the levels of Subj indicate a grouping structure (k subjects) and Trt is a two-level factor (two treatments) for which there are several (n) responses y from each treatment and subject combination. If one suspects a subject by treatment interaction, either of the following models seem

Hi I am working on a Nested one-way ANOVA. I don't know how to implement R code to test the significance of the random factor My R code so far can only test the fixed factor : anova(lme(PCB~Area,random=~1|Sites, data = PCBdata)) numDF denDF F-value p-value (Intercept) 1 12 1841.7845 <.0001 Area 1 4 4.9846 0.0894 Here is my data and my hand

Dear R users I have run an regression and want to extract the p value of the F statistics, but I can find a way to do that. x<-summary(lm(log(RV2)~log(IV.m),data=b)) Call: lm(formula = log(RV2) ~ log(IV.m), data = b[[11]]) Residuals: Min 1Q Median 3Q Max -0.26511 -0.09718 -0.01326 0.11095 0.29777 Coefficients: Estimate Std. Error t value Pr(>|t|)

Hello, I ran two lme analyses and got expected results. However, I saw something suspicious regarding p-level for fixed effect. Models are the same, only experimental designs differ and, of course, subjects. I am aware that I could done nesting Subjects within Experiments, but it is expected to have much slower RT (reaction time) in the second experiment, since the task is more complex, so it

Hi, sorry by this off. I'm still try to understand nested design. I have the follow example (fiction): I have 12 plots in 4 sizes in 3 replicates (4*3 = 12) In each plot I put 2 species (A and B) to reproduce. After a period I make samples in each board and count the number of individuals total (tot) and individuals A and B (nsp). Others individuals excepts A and B are in total of

Dear all, Using the data set and code below, I am interested in modelling how egg temperature (egg.temp) is related to energy expenditure (kjday) and clutch size (treat) in incubating birds using the lme-function. I wish to generate the F-distribution for my model, and have tried to do so using the anova()-function. However, in the resulting anova-table, the parameter kjday has gone from being

Dear R-help list, I would like to perform multiple comparisons for lme. Can you report to me if my way to is correct or not? Please, note that I am not nor a statistician nor a mathematician, so, some understandings are sometimes quite hard for me. According to the previous helps on the topic in R-help list May 2003 (please, see Torsten Hothorn advices) and books such as Venables &

Dear All, I am writing to ask a clarification on what R, and in particular lme, is doing. I have an experiment where fly wing area was measured in 4 selection lines, measured at 18 and 25 degrees. I am using a lme model because I have three replicated per line (coded 1:12 so I need not use getGroups to creat an orederd factor). The lines are called: "18"; "25";

Hi, I have been testing sieve in my setup with qmail-ldap and deliver on a LAN with an artificial domain name. Everything seems to be working as expected, except in cases when autoreplies (vacation, reject messages) need to be tested. The domain name is vmint, and dawnone is the hostname on which mail server is setup, so a users have address like cot at vmint, cute at vmint and dove at vmint *

Dear All, I've come across a problem in predict.lme. Assigning a model formula to a variable and then using this variable in lme (instead of typing the formula into the formula part of lme) works as expect. However, when performing a predict on the fitted model I gan an error messag - predict.lme (but not predictlm) seems to expect a 'properly' typed in formula and a cannot extract

tl;dr anova.lme() claims to provide sums of squares, but it doesn't. And some names are misspelled in ?lme. I can submit all this stuff as a bug report if that's preferred. ?anova.lme says: When only one fitted model object is present, a data frame with the sums of squares, numerator degrees of freedom, denominator degrees of freedom, F-values, and P-values The output of fm1

Hi all I tried to reproduce an example with lme and used the Orthodont dataset. library(nlme) fm2a.1 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1 | Subject) anova(fm2a.1) > numDF denDF F-value p-value > (Intercept) 1 80 4123.156 <.0001 > age 1 80 114.838 <.0001 > Sex 1 25 9.292 0.0054 or alternatively

Dear all, I am starting to use the lme package (and plan to teach a course based on it next semester...). To understand what lme is doing precisely, I used balanced datasets described in Pinheiro and Bates and tried to compare the lme outputs to that of aov. Here is what I obtained: > data(Machines) > summary(aov(score~Machine+Error(Worker/Machine),data=Machines)) Error: Worker

Hello, I am trying to run a mixed effects nested ANOVA but none of my codes are giving me any meaningful results and I am not sure what I am doing wrong. I am a new user on R and would appreciate some help. The experimental design is that I have some frogs that have been exposed to three acoustic Treatments and I am measuring neural activity (egr), in 12 brain regions. Some frogs also called

Greetings, i am testing: - dovecot-1.2-sieve-0.1.11 - dovecot 1.2.3 - dovecot-1.2-managesieve-0.11.8 (with patch) I am trying to use sieve_after to execute some redirection of mails. I have put in the after script: require "include"; require "fileinto"; if header :contains ["To","Cc"] "userA at acasalud.com.ar" { keep; redirect