similar to: bug in nlme package function predict.lmList (PR#13788)

Dear all, I've got a problem in including square variables in lme and lmlist functions. I've tried to work on Oxboys data of Pinheiro and Bates'book, which consist of the heights of 26 boys, each mesured on nine different occasions : > Oxboys Grouped Data: height ~ age | Subject Subject age height Occasion 1 1 -1.0000 140.50 1 2 1 -0.7479 143.40

Dear Martin, Have you try to create a new variable for age squared, say agesq? If you fit the model using this new variable you should get the coefficients. So your new model is something like height~age+agesq I hope this helps, Saghir > -----Original Message----- > From: MARTIN Ludovic [SMTP:martinl@mathinfo.ens.univ-reims.fr] > Sent: Wednesday, 07 May, 2003 2:02 PM > To:

Full_Name: Vaidotas Zemlys Version: 2.1.1 OS: Ubuntu 05.10 Submission from: (NULL) (213.197.173.50) Hi, This is really a feature request, not a bug. I wrote the mail to R-devel, but nobody answered it. I use Gnome on my computer and sometimes I use its default text editor gedit. It uses gtksourceview library for syntax highlighting. I decided that it would be nice if gedit supported R

Hi, Sweave for some reason does not respect encoding setting in pdf.options. Everything is OK with ps.options. Thus when using non-default encoding for graphs, resulting graphics eps file is correct, but pdf file is not. Workaround is very simple, just use latex->dvips->ps2pdf chain to get the intended pdf document. I tried to find the explanation for this, but I did not find anything. If

On 17 May 2017 at 08:46, Vaidotas Zemlys wrote: | Hi, | | > Le 17 mai 2017 ? 00:42, Dirk Eddelbuettel <edd at debian.org> a ?crit : | > | > | > On 8 May 2017 at 15:39, Vaidotas Zemlys wrote: | > | Hi, | > | | > | Dirk Eddelbuettel advised me to write here. Here is my original letter to him: | > | | > | I would like to enquire about package r-cran-rjava on

Hi, > Le 17 mai 2017 ? 00:42, Dirk Eddelbuettel <edd at debian.org> a ?crit : > > > On 8 May 2017 at 15:39, Vaidotas Zemlys wrote: > | Hi, > | > | Dirk Eddelbuettel advised me to write here. Here is my original letter to him: > | > | I would like to enquire about package r-cran-rjava on Debian jesse. It seems that if default-jre package is not installed, but

Hi, Today I came upon unexpected R behaviour. I did some modelling and the result was R object, about 28MB size (nested list, with matrixes as list elements). When I was saving the session with save.image, the resulting .RData file was 300MB. There were no other large objects: >

Dear list members, A new package, called *midasr* (version 0.1), is now available on CRAN: http://cran.r-project.org/package=midasr This package provides econometric methods for working with mixed frequency data. The package provides tools for estimating the time series MIDAS regression, where the response and explanatory variables are of different frequency, e.g. quarterly vs monthly. The

Hi, I was trying to repeat the estimation of threshold GARCH models from the book "Analysis of Financial Time Series" by Ruey S. Tsay, and I was succesfull, but I had to use "for" loop, which is quite slow. The loop is necessary, since you need to calculate recursive sequence. Is there a faster way to do this in R, without using loops? The model is such: r_t = \mu + \alpha_2

Hi, Is it possible to make R start not from terminal with option --gui=GNOME? Now if in the launcher preferences checkbox run from terminal is not checked, R with option --gui=GNOME doesn't start. If started with checked checkbox, i.e. from terminal, besides R, terminal window also appears, but it is of no use. Yet if it is closed, the R window closes together. Is this a feature, or is

Hi, Dirk Eddelbuettel advised me to write here. Here is my original letter to him: I would like to enquire about package r-cran-rjava on Debian jesse. It seems that if default-jre package is not installed, but openjdk-7-jre is installed, then library(rJava) in R fails. I?ve been bitten by this today and I wonder whether this an issue of mine, or is this a possible bug. My server admin used

Hello R-Helpers, I was trying to use the lmList function to get the lmList graphic similar to Pinheiro and Bates (pg 33). I did not have a problem creating the graphic when I used the Orthodont data frame or 2 other data sets when there are no missing values. My data has missing values. Do I need to remove the missing values before the lmList function will work? for a small example: > a

Dear colleges, I would like to use the lmList function. I have installed the lme4 library but when I try to use the lmList function I get this error message. Error: could not find function "lmList" Here you can see what kind of messages I am getting when installing lme4 package. Does anyone know how to solve this problem? Thanks Rosario > install.packages() --- Please select a

Dear all, I have a question on handling of missing values in lmList. My data set have continuous predictor and response, x and y, and a grouping variable group.id. All these variables have NAs and the data set also has several other variables that also contains NAs. To create the lmList-object seems to work fine: y.list <- lmList(y ~ x | group.id, data=mydata, na.action=na.omit) However,

In R, we currently have the functions anovalist.lm and anova.glmlist S / S-plus has anova.lmlist anova.glmlist On the other hand, the [n]lme package (library) of Doug Bates and Jose Pinheiro has an "lmList" class and an anova.lmList(.) method for that. We are considering to use anovalist.lm and anovalist.glm instead of the S/S-plus names mentioned above. These functions

Hello, in R-2.1.0 I'm trying to prodice trellis plots from an lmList object as described in the help for plot.lmList. I can generate the plots from the help, but on my own data plotting fails with an error message that I cannot interpret (please see below). Any hints are greatly appreciapted. kind regards, Arne > dim(d) [1] 575 4 > d[1:3,] Level_of_Expression SSPos1 SSPos19

Dear community, I'm trying to examine my grouped data following page 6 http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-mixed-models.pdf. I'm trying lmList this way: model.list <- lmList(log(v.dep) ~ log(v2) +log(v4) + v3 | v5, subset = v6=0, data=data) and obtain this error message: In Ops.factor(v4, v4) : | not meaningful for factors My original model is

Hi, Does R have any functions implementing such multinomial regression: (S_t^A,S_t^B)~MN(N_t-Y_{t-1},P_t^A,P_t^B) where MN(n,p_1,p_2) is multinomial distribution with parameters n, p_1, p_2. Here P_t^A and P_t^B are nonlinear functions from predictor variables and parameters which need to be estimated. Here A and B are used for notation, they are not parameters. My second question is about

I'd like to fit a GLM to each of a number of subsets of some data. The `family' argument to `lmList' (in lme4) has given me cause for optimism, but so far I've only been able to achieve linear model fits. For example > df <- data.frame(gp = gp.temp <- factor(rep(1:3, each = 100)), x = x.temp <- rnorm(300), y = rpois(300, exp((-1:1)[gp.temp] + x.temp))) Then a call to

Hello, I installed the newest version of R (3.0.2) as well as the newest version of nlme (3.1-113) in order to use summary.lmList and other nlme functions. Once loading the new library, lmList and summary.lm can be found, but a number of additional functions cannot be found via command. Versions installed and loaded are correct. Any suggestions on how to allow for these functions of the new