similar to: anova function to test the difference between two coefficients in nlme package

Dear all, This question may be too basic quesition for this list, but if someone has time to answer I will be happy. I have tried to find out, but haven't found a consice answer. As an example I use "Pinheiro, J. C. & Bates, D. M. 2000. Mixed-effects models in S and S-PLUS. Springer, New York." page 225, where rats are fed by 3 different diets over time, which body mass has

HI Dave, My comment was based? on: " >The main question with this test was if the interaction term is significant (i.e. growth rate). However, my question is could I also look at the p-values of the main effects to ?>say if body mass increase significant with body mass?" Here, the result shown were from the summary of the linear model.?? We report the p-values of the main

Dear R-help, In nlme package, anova () can be used to test the difference between two coefficients as shown on page 225 of "Mixed Effects Models in S and S-Plus": anova(fm2BW.lme, L = c("Time:Diet2" = 1, "Time:Diet3" = -1)) Now my question is instead of test the difference between two coefficients, can I use anova to test the significance of the sum of two

  Hola,   Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?   http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5   Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.   Muchas gracias.   Eva [[alternative HTML

Les consulto por lo siguiente: Se me ocurrió usar datos de MCMCglmm en otros procedimientos de R, comparar los resultados y el tiempo como facilidad de uso, para lo cual use LME4, los resultados fueron parecidos, entendiendo que internamente los algoritmos tienen diferencias, es lógico que los resultados no coincidan exactamente, pero aproximadamente tienen que ser compatibles, en este caso

I'm working with a nested model (mixed). I have four factors: Patients, Tissue, sex, and tissue_stage. Totally I have 10 patients, for each patient, there are 2 tissues (Cancer vs. Normal). I think Tissue and sex are fixed. Patient is nested in sex,Tissue is nested in patient, and tissue_stage is nested in Tissue. I tried aov and lme as the following, > aov(gene ~ tissue + gender +

When I try to run the example from Variogram with an lme object, I get an error (although summary works): R : Copyright 2005, The R Foundation for Statistical Computing Version 2.2.1 (2005-12-20 r36812) ISBN 3-900051-07-0 ... > fm1 <- lme(weight ~ Time * Diet, BodyWeight, ~ Time | Rat) Error: couldn't find function "lme" > Variogram(fm1, form = ~ Time | Rat, nint =

I encounter such problem with lmer under FreeBSD, but not under Windows. Anyone knows why? Thanks. > example(lmer) lmer> (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) Error in UseMethod("as.logical") : no applicable method for "as.logical" > traceback() 9: as.logical(EMverbose) 8: as.logical(EMverbose) 7: lmerControl() 6:

I have a simple dataset of a numerical dependent Y, a numerical independent X and a categorial variable Z with three levels. I want to do linear regression Y~X for each level of Z. How can I do this in a single command that is without using lm() applied three isolated times? -- View this message in context: http://r.789695.n4.nabble.com/Linear-regression-on-several-groups-tp2322835p2322835.html

Hi, I am trying to find a way to automate production of page numbers in plots produced using the lattice package. To do this, I started playing around with setHook which works fine with vanilla plot, but the hook 'plot.new' doesn't appear to be relevant to the lattice package. I was wondering if there is a alternative someone can suggest for use with the lattice package. Example

Hi all, I am encountering a strange issue with the Matrix package. I have just built R-devel from source on my macbook pro, and I wonder if others can reproduce this problem. I will give example code to go along: Starting a fresh R session: R version 2.6.0 Under development (unstable) (2007-08-30 r42697) Copyright (C) 2007 The R Foundation for Statistical Computing ISBN 3-900051-07-0 ... >

Dears Sirs During my computational work I encountered unexpected behaviour when calling "ar" function. I want to select the order p of the autoregressive approximation by AIC criterion and sometimes an error occurs. Example: # time series

Hi Dave, Nice driver optimization from Shirley, but requires a new virtio hook. Do you want to take both? I have nothing else overlapping it. Cheers, Rusty.

Hi Dave, Nice driver optimization from Shirley, but requires a new virtio hook. Do you want to take both? I have nothing else overlapping it. Cheers, Rusty.

Hello, I have a data frame consisting of four columns and would like to sort based on the first column and then write the sorted data frame to a file. > df <- read.table("file.txt", sep="\t") where file.txt is simply a tab-delimited file containing 4 columns of data (first 2 numeric, second 2 character). I then do, > df_ordered <- df[order(df$V1), ] OR,

Hi, I'm writing an AGI script to manage outgoing calls. We need to interrogate a database to work out which line a particular user is allowed to use for outgoing calls. However, I cannot find a way for my AGI script to access the SIP username. Does anyone know if this is possible (even if it is just passing a variable from extensions.conf to the script)? Thanks in advance, David

On Fri, 10 Nov 2017 03:41:20 -0500 Bill Shirley <bill at KnoxvilleChristian.org> wrote: > No it isn't shown as a folder.? All folder directories here begin with a dot. > i.e.? .INBOX? .Trash? .Drafts > > Bill No, they don't. me thought that, too. But using the rainloop webmail interface on top of such a config showed the sieve folder in the overview. Sometimes you can

Hello all, Anybody found this problem before? I kept hitting this issue for 2.6.31 guest kernel even with a simple network test. INFO: task kjournal:337 blocked for more than 120 seconds. "echo 0 > /proc/sys/kernel/hung_task_timeout_sec" disables this message. kjournald D 00000041 0 337 2 0x00000000 My test is totally being blocked. Thanks Shirley

Hello all, Anybody found this problem before? I kept hitting this issue for 2.6.31 guest kernel even with a simple network test. INFO: task kjournal:337 blocked for more than 120 seconds. "echo 0 > /proc/sys/kernel/hung_task_timeout_sec" disables this message. kjournald D 00000041 0 337 2 0x00000000 My test is totally being blocked. Thanks Shirley

Dear R help, In limma package, contrasts.fit() function is very useful. I am wondering whether there is a similar function for lme object, which means given a mixed linear model fit, compute estimated coefficients and standard errors for a given set of contrasts. Thanks, Shirley