similar to: fromulation of random effects (nlme)

Doug, I have been attempting to learn a little bit about nlme without too much documentation except the online help. The Latex file in the nlme directory looks interesting but uses packages that I do not have so that I have not been able to read it. I have run the example from gnls to compare it with the results I get from my libraries (code below - I have not included output as it is rather

Full_Name: Madeleine Yeh Version: 1.9.1 OS: AIX 5.2 Submission from: (NULL) (151.121.225.1) After compiling R-1.9.1 on AIX 5.2 using the IBM cc compiler, I ran the checks. One of them failed. Here is the output from running the check solo. root@svweb:/fsapps/test/build/R/1.9.1/R-1.9.1/tests/Examples: ># ../../bin/R --vanilla < stats-Ex.R R : Copyright 2004, The R

In Win95/NT you can use xgobi (program and R interface) as well, provided you will install an x-server on your machine (a free one is available from http://www.microimages.com/freestuf/mix/). Regards! Zdenek Skala ***************************************************** Uli Flenker wrote: I'm not familiar with S/S+, but I think the "xgobi" software package can do what you want

Dear R-experts, I have a problem to formulate non linear mixed effects when more than one explanatory variable is present. I'm using R 1.1.1 under Linux. The version number of the nlme library is 3.1-7. "help(nlme)" says that among several possibilities it is possible to code fixed effects according to fixed = a1 + a2 + ... + an ~ X1 + X2 + ... + Xm where the left hand side

Hello everybody, since I installed R-0.90.0 (Linux 2.0.25) I have problems with generating printable output with 'help()'. Default output to the terminal ('echo $PAGER' gives 'less') works o.k.. When invoking 'help(topic,offline=T)' all the usual procedures seem to start: 'top' reports invocation of latex as well as of dvips. But after finishing there is

Hello everybody, I observed some strange behaviour of R (0.64.1/Linux 2.0.25) when trying to use FORTRAN-libraries. I first downloaded "ALGORITHM AS 304" from StatLib. All I changed at the code, was to substitute "REAL"-routines by "DOUBLE PRECISION". Compiling went o.k.: > g77 -fpic -O2 -shared -o rtest.o fisher.f ("g77 --version" gives 2.7.2.1)

Dear R-helpers, I frequently run into problems when I modify elements of factors (R 1.0.0, Linux 2.0.25). For example, after splitting a data frame accordinng to whatever criterion, it might well happen that not all levels of a factor are present in all new data frames. However, R doesn't seem to care about. When performing some kind of analysis seperately on the new frames, each time the

Hello everybody, I observed some strange behavior of R (0.64.1/Linux 2.0.25) when trying to use FORTRAN-libraries. I downloaded "ALGORITHM AS 304" from StatLib. All I changed at the code, was to substitute "REAL"-routines by "DOUBLE PRECISION". Compiling went o.k. ... > g77 -fpic -O2 -shared -o rtest.o fisher.f ("g77 --version" gives 2.7.2.1)

Dear R-helpers, for the first time since R-0.6x I have to face installation problem. I switched to a new Linux-box (SuSE 9.0) and installed as many libraries as possible. R-1.9.0.tgz went to /usr/local/lib and was extracted without problems. After running properly for a while, 'configure' (called without any switches) gave "configure: creating ./config.status" which was

We had problems with our mail server. I have no idea whether this message has reached the list or whether there were any replies. So I retry ... ------------------------ Date: Wed, 5 May 1999 16:42:55 +0200 (CEST) From: "Uli Flenker; Raum 704" <uli at pcinternet.biochem.dshs-koeln.de> To: R-Mailing-List <r-help at stat.math.ethz.ch> Subject: ANOVA "ex post"

Hi all, I am having issues comparing models with lmer. As an example, when I run the code below the model summaries (AIC, BIC, loglik) differ between the summary() and anova() commands. Can anyone clear up what's wrong? Thank you! Darren Ward library(lme4) data(sleepstudy) fm1<-lmer(Reaction ~ Days + (1|Subject), sleepstudy) summary(fm1) fm2<-lmer(Reaction ~ Days +

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

Hello everybody, I frequently have to produce barplots from variables the natural variation of which is far away from zero. A typical range would be, say from -16 to -28. So, the following example should give an "honest" presentation: x<-c(-20,-22,-21,-28) barplot(x,beside=T,ylim=c(-18,-30)) But using anything else but "0" for the first element of ylim leads to strange

Dear R community, I'm trying to understand this behavior of TukeyHSD. My goal is to obtain defensible, labelled multiple comparisons of an interaction term. Firstly, if I plot the TukeyHSD from the model that calculates its own interactions, then the y-axis labels appear to be reflected on their median when compared to the text output of the TukeyHSD statement. The labels are integers.

Dear all, I noticed the following in the call of lme using msVerbose. fm1 <- lme(distance ~ age, data = Orthodont, control = lmeControl(msVerbose=T)) 9 318.073: -0.567886 0.152479 1.98021 10 318.073: -0.567191 0.152472 1.98009 11 318.073: -0.567208 0.152473 1.98010 fm2 <- lme(distance ~ age, random =~age, data = Orthodont,

Hello everybody, shame on me if I have overlooked something (CRAN, StatLib), but I think I've searched carefully. Maybe it's just too obvious to see for me (happens frequently). Is a function available that extracts multiple comparison of means from objects produced by "aov()" ? Sorry if the term is not correct, I translated from german word by word ("Multiple

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

The subject says it all really. Question 1. Here is some code created to illustrate my problem, can anyone spot where I'm going wrong? Question 2. The reason I'm following a manual specification of models relates to the fact that in reality I am using mgcv::gam, and I'm not aware that dredge is able to separate individual smooth terms out of say s(a,b). Hence an additional request,

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 =