similar to: Summary: proc mixed vs. lme

Dear All, Comparing linear mixed effect models in SAS and R, I found the following discrepancy: SAS R random statement random subj(program); random = ~ 1 | Subj -2*loglik 1420.8 1439.363 random effects variance(Intercept) 9.6033 9.604662

Dear all, I'm working with Cox-regression, because data could be censored. But in this particular case not. Now I have a simple example: PRO and PRE are (0,1) coded. The response is not normal distributed. We are interested in a model which could describe interaction. But my results are depending strongly in the choose of the contrast option. It is clear that there is some dependence in

Dear Experts, May be the problem is still solved, however I tried to find the answer in the archives: I use: > R.version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 7.1 year 2003 month 06 day 16

Hi, I'm using: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 5.1 year 2002 month 06 day 17 language R and I would like to apply: > update.packages() trying URL

Dear R-helpers, First of all I have nothing against the French language! But now my problem, yesterday I installed R 2.1.1 and I had to experience that my RGui is speaking French. My windows locals is French (Switzerland). I'm used to English and I want to reset my RGui to English. I was seeking for the solution in the archives, however not successfully. By the way the searchable archives

Hello, I'm using: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 8.1 year 2003 month 11 day 21 language R and I get the following error with: library(Design) df <- list(pre=c(0,, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1,

Dear all, I have a question about linear models. Is there something comparable to the estimate statement in SAS? E.g. proc mixed data=test; class x1 x2; model y=x1|x2 /s; estimate 'x1' x1 -1 1 x1*x2 -0.5 -0.5 0.5 0.5; estimate 'x2' x2 -1 1 x1*x2 -0.5 0.5 -0.5 0.5; run; Thanks for the help, Dominik Dominik Grathwohl Biostatistician Nestl? Research

Dear R-helpers: My question is how do I efficient and valid correct for multiple tests in a repeated measurement design: Suppose we measure at two distinct visits with repeated subjects a treatment difference on the same variable. The treatment differences are assessed with a mixed model and adjusted by two methods for multiple tests: # 1. Method: Adjustment with library(multcomp)

Hello, ? I believe that in R, it is not possible to analyze mixed effect-models when the distribucion is not gaussian (p.e. binomial or poisson), isn't? ? Somebody can suggest me alternative? ? thanks ? xavi ? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Dear r-helpers, I have a question about multiple testing. Here an example that puzzles me: All matrixes and contrast vectors are presented in treatment contrasts. 1. example: library(multcomp) n<-60; sigma<-20 # n = sample size per group # sigma standard deviation of the residuals cov1 <- matrix(c(3/4,-1/2,-1/2,-1/2,1,0,-1/2,0,1), nrow = 3, ncol=3, byrow=TRUE, dimnames =

Dear Prof. Ripley, you mention the theory of perceptrons. Could you please point me to an introduction paper or book? Thanks in previous, Dominik > -----Original Message----- > From: ripley at stats.ox.ac.uk [mailto:ripley at stats.ox.ac.uk] > Sent: dimanche, 10. novembre 2002 18:55 > To: Ben Liblit > Cc: r-help at stat.math.ethz.ch > Subject: Re: [R] building a formula for

I assigned a class the first problem in Pinheiro & Bates, which uses the data set PBIB from the SASmixed package. I have recently downloaded 2.0.1 and its associated packages. On trying library(SASmixed) data(PBIB) library(nlme) plot(PBIB) I get a warning message Warning message: replacing previous import: coef in: namespaceImportFrom(self, asNamespace(ns)) after library(nlme) and a

Hi, I want to free memory used by R. The usage of rm and gc give no result. I'm running an algorithm consuming a huge memory and I need to recover the memory used by R between 2 call of my algorithm. Thank you in advance for your help. e-mail: sofiane.lariani at rdls.nestle.com Sofiane Lariani

Based on the stepAIC help, I have assumed that it only was for lm, aov, and glm models. I gather from the following correspondence that it also works with lme models. Thomas Lumley 07:40 a.m. 28/04/03 -0700 4 Re: [R] stepAIC/lme problem (1.7.0 only) Prof Brian Ripley 04:19 p.m. 28/04/03 +0100 6 Re: [R] stepAIC/lme problem (1.7.0 only) Prof Brian Ripley 06:09 p.m. 29/04/03 +0100 6 Re: [R]

Dear list, I ran a mixed effect model using R 1.4.1 and SAS 8.0 on the SIMS data found in the SASmixed package and found that the degrees of freedoms for fixed effects are very different. From R, df = n - v -1 where n is total # of observations, v is the # of levels for the grouping factor. From SAS df = v -1. Am I wrong about this or can somebody explain which is correct and why? Thanks a

I suppose that multcomp in R and multicomp in S-Plus are related and it appears that it is possible to use multicomp with lme in S-Plus given the following correspondence on s-news sally.rodriguez at philips.com 12:57 p.m. 24/04/03 -0400 7 [S] LME summary and multicomp.default() Is it possible to use multicomp with lme in R and if so what is the syntax from a simple readily available

Dear all, In the SASmixed package there is an example of an analysis of a split-plot experiment. The model is fm1Semi <- lme( resistance ~ ET * position, data = Semiconductor, random = ~ 1 | Grp) where Grp in the Semiconductor dataset is defined as ET*Wafer. Is it possible to specify the grouping directly some way, e.g. like fm1Semi <- lme( resistance ~ ET * position, data =

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

Dear all, I want to get a series of boxplots (grouped by two factors) and I want to overlay the original observations and the following code does almost what I want: library(ggplot) ddf <- data.frame(x=factor(rep(LETTERS[1:4], each=30)), y = runif(120,0,10), grp = factor(rep(rep(1:3, 10), 4))) ggplot(ddf, aes(x, y, colour=grp)) + geom_boxplot() + geom_point() Yet the position of the points

Dear list members, I am new to R and to the R-help list. I am trying to perform a mixed-model analysis using the lmer() function. I have a problem with the output anova table when using the anova() function on the lmer output object: I only get the numerator d.f., the sum of squares and the mean squares, but not the denominator d.f., F statistics and P values. Below is a sample output, following