similar to: Sv: lme problem ?

Dear friends. I tried the session below with 10 MB in both vsize and nsize but didn't get the example work. Is this a problem in LME or in me or both or somewhere else or undefined ? R : Copyright 1999, The R Development Core Team Version 0.64.0 Patched (May 3, 1999) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type

Thank you to Peter Dalgaard, who provided help for some of this solution. It is still not ideal since the dots which are equidistant in the x-direction ideally should respect the reciprocal scale in the y-direction. op <- par("cex"= 0.7) x<-y<-NA plot(x,y,xlim=c(0,10),ylim=1/c(2000,70), axes=F) box() axis(1)

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,

Doug, I thought perhaps that you might be interested in the comparison of lme to the results for the same models fitted by Richard Jones' carma (I just wrote the R interface to his Fortran code). The code to run the example from the lme help and for the equivalent with carma is in the file below. The two main differences in results are 1. the random coefficients covariance matrix is quite

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

Hi everyone, I have a quick and probably easy question about lme for this list. Say, for instance you want to model growth in pituitary distance as a function of age in the Orthodont dataset. fm1 = lme(distance ~ I(age-8), random = ~ 1 + I(age-8) | Subject, data = Orthodont) You notice that there is substantial variability in the intercepts (initial distance) for people at 8 years, and that

Dear Peter. I guess you know that Jim Lindseys code include nordr and ordglm in library gnlm - I attach the htmls which do various linear and nonlinear ordinal regressions - exemplified with just the data mentioned, McCullagh (1980) JRSS B42, 109-142. I had it work very fine. -----Oprindelig meddelelse----- Fra: Peter Malewski <p.malewski at tu-bs.de> Til: Troels Ring <tring at

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

Dear friends. When using in windows version the menu function to change directory, what actions are actually performed, i.e. what commands could I write to a file to make a CD ? I couldn't find it mentioned anywhere - undoubtedly an oversigth on my part. Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk

Dear R developers, first of all let me join the chorus of congratulations for the release of R 1.0.0. Well, done! Unfortunately, I find it necessary to e-mail in a bug report regarding the `nlme' package. On my office machine I experience the following trouble: bossiaea:/opt/R$ R CMD check -c nlme Checking package `nlme' ... Massaging examples into `nlme-Ex.R' ... Running

Dear friends. I've forgotten how to simulate data with a known correlation, e.g. only two columns. I'd be most pleased to be told. Thanks in advance Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Full_Name: Per Broberg Version: 1.00 OS: Windows 98 Submission from: (NULL) (62.20.231.229) I tested my installation with the following: > library(lme) Loading required package: nls Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library "C:\PROGRAM\R\RW1000/library/nls/libs/nls.dll": LoadLibrary failure > data(Orthodont) > fm1

Dear friends. I would like to link an object file made in C. Dyn.load rather wants to link a whole library and I only wanted to link this small object file. I followed instructions given for S - by Clive Loader in his software for changepoint handling and it is his object file I made from http://cm.bell-labs.com/cm/ms/departments/sia/clive/smcp.shar Any dirty tricks or easy oversights ? Yours

Dear friends. I just downloaded tseries dated 27.10.99 (from Aalborg, Denmark) and it seems it lacked something, library ts was not available ? What to do ? > library(tseries) Loading required package: ts Warning: There is no package called `ts' > > Best wishes Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk

Dear friends. I wanted to reproduce a sheet of paper useful for plotting the reciprocal of some biochemical variable (creatinine) over time (Bleyer Am J Kidney Dis 34:576-578,1999). Time on X axis with units years and months. On Y axis, 1/creatinine, but with the original creatinine value indicated - to make it easier for people to use it. When using the graph, the linear decrease in 1/crea over

Using formulas constructed from strings only partially works for me in lme: library(nlme) data(Orthodont) fm2<-lme(as.formula("distance~age"),data=Orthodont,random=~1|Subject) summary(fm2) # works augPred(fm2) # fails #Error in inherits(object, "formula") : #Argument "object" is missing, with no default I assume that my use of as.formula is wrong, but

Dear friends. On the Bayesian averaging homepage http://www.research.att.com/~volinsky/bma.html I found some S code some of which perhaps may run in R. There was a call to an algorithm possibly within S but not supported by R 64.1: "leaps and bounds". I guess it is a minimization step. Can anyone clarify the algorithm and perhaps even give a pointer to some code ? I guess this may be

I am trying to estimate models with subsets using the NLME package. However, I am getting an error in the case below (among others): > subset <- c(rep(TRUE, 107), FALSE) > fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1, subset=subset) Error in xj[i] : invalid subscript type 'closure' > fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1,

Dear friends. Can anyone explain what is happening here ? I am told that no data entry editor is available in this version of R ? The example is not far from that shown in help(data.entry) ? Besides, could an empty spreadsheet entry be an option to start inputting data ? R : Copyright 2000, The R Development Core Team Version 1.0.0 (February 29, 2000) R is free software and comes with

Dear friends. Below is ouput directly from the help on qqnorm.lme.html It do not seem to work as expected - on win98, R 1.1, although the regression is undertaken, so qqnorm misunderstands, or what ? library(nlme) data(Orthodont) fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject) # normal plot of standardized residuals by gender qqnorm(fm1, ~ resid(., type = "p") | Sex,