similar to: Performance problem

Hi, I usually use lme from the nlme library. Now I have read an article about lmer in Rnews and lmer seemed to me more comfortable to use. Unfortunately, I didn't find out how to use covariance structures (e. g. corSymm(), corAR1()). Is there a way to use them similarly as in lme ? Is it implemented ? If somebody knows, please let me know. Thank you very much in advance, Stephan

Claus-Juergen Neumann/BASF-AG/B ASF An r-help at r-project.org 23.09.2008 12:49 Kopie

Full_Name: Liqun Xing Version: 2.6.2 OS: Solaris 9 Submission from: (NULL) (76.182.91.82) When first time launch X11(), it was fine. However, when closed the X11 windon by clicking the close button, or closed it off by command Dev.off() and launched the X11() again, got segfault, and asking for exiting R and other options.

Hi I constructed a mixed-effects model from longitudinal repeated measurements of lab values in 22 patients seperated into two groups with the groups as fixed effect using lme. I thought about using the jackknife procedure, i. e., removing any one subject and calculating the fixed effect, to assess the stability of the fixed effect and thereby validate the model. I suppose this has been done in

Dear "R"-Team, we had some problem to install R-2.7.2 on HP-UX 11.11 on PA-Risc platform. 1.) i load down the Software and made a extract of the tar-file in a specific Directory. 2.) i follow the instructions on the INSTALL file: ./configure 3.) during "configure" we got a error message see below: checking for history_truncate_file... no configure: error:

I tried following commands: amp~time|subject/trial #this was the grouping structure of the data plot(dip,inner=~condition,layout=c(2,2)) after the plot command I obtained this error message: Error in if(!any(cond.max.level - cond.current.level <0)&&(row-1)* : missing value where TRUE/FALSE needed This error only occured in compination with "layout". It was no

Trying following example from Pinheiro and Bates in order to fit an ARMA(1,1) model: library(nlme) fm1Ovary.lme<-lme(follicles~sin(2*pi*Time)+cos(*pi*Time),data=Ovary,random=p dDiag(~sin(2*pi*Time))) fm5Ovary.lme<-update(fm1Ovary.lme,corr=corARMA(p=1,q=1)) I get follwing error message: Error in "coef<-.corARMA"(`*tmp*`, value = c(62.3428455941166, 62.3428517930051 :

I have a dataset that contains numbers between -10 and 0. E.g. x = c(-9.23, -9.56, -1.40, ...) If I no do a > qplot(x, geom="histogram") I get the error: Error: position_stack requires non-overlapping x intervals Strangely, the following both work: > qplot(x * -1, geom="histogram") > qplot(x+100, geom="histogram") Has anyone else encountered this? Is this a

I have recently started a new job at an NHS hospital in Scotland. Since I took up this post 6 months ago I have had an ongoing dispute with the IT secutiry dept. who refuse to install R on my computer. I previously worked in another branch of the NHS where R was widely used and yet there is nothing I can say which will persuade the IT dept here to even visit the website! With some help from our

Fellow R-users, I have a longitudinal data set with missing values in it. I would like to produce a residual plot for each time using panel.xyplot function but I get an error message. Here's a simple example, library(nlme) set.seed(1544) longdata <- data.frame(ID=gl(10,1,50), y=rnorm(50), time = as.numeric(gl(5,10,50)), x = rnorm(50)) longdata$y[5] <- NA longdata$x[35] <- NA

Dear R-helpers, I am trying to construct a confidence interval on a prediction of a gam fit. I have the Wood (2006) book, and section 5.2.7 seems relevant but I am not able to apply that to this, different, problem. Any help is appreciated! Basically I have a function Y = f(X) for two different treatments A and B. I am interested in the treatment ratios : Y(treatment = B) / Y(treatment = A) as

Valued R-community. Please let me make you aware of a job opening with focus on Statistics at BASF in Shanghai: We are the world's leading chemical company because we offer intelligent solutions for our customers and for a sustainable future. We link and develop people with diverse talents all over the world. For you, this means a variety of ways to advance. Not only your performance but

Hi the list, I build a package. They was a mistake in it, but R CMD check did not find it. Is that normal ? Here is what Kurt gets (which is right, I did this mistake): --- 8< ---------------- * checking for code/documentation mismatches ... WARNING S4 class codoc mismatches from documentation object 'LongData-class': Slots for class 'LongData' Code: id other time traj

Hi the list, I am the maintainer of the package kml. I quite often receive some email about a specific bug. The answer is always the same, I tell people to remove the package Matrix, and then kml works fine. I wonder what my code is not compatible with the Matrix package. Here is a (simplified) version of the kml code : --- 8< ------------ library(longitudinalData)

Hi, Is it possible to create Rose Diagrams of wind data (speed & direction) with R?? Best regards, Lars Peters ----- Lars Peters University of Konstanz Limnological Institute D-78457 Konstanz Germany phone: +49 (0)7531 88-2930 fax: +49 (0)7531 88-3533 e-mail: Lars.Peters@Uni-Konstanz.de web: Lars Peters <http://www.uni-konstanz.de/sfb454/tp_eng/A1/doc/peters/peters.html>

I got a set of data which has seasonal trend in form of sinx, cosx, I don't have any idea on how to deal with it. Can you give me a starting point? Thanks, Terry

We need to know a great deal more information, and asking on R-devel would be a lot kinder than blaming this on a bug in R (and please use a meaningful subject line). What version of R is this? Precisely what OS (Solaris 9 comes for several architectures)? Assuming you compiled R yourself, what compilers did you use? (And where did you get those from if they are not Sun compilers?) If got a

Hi, I got a computational problem with lme (nlme library R 1.9.1) using corSymm(). Here is the data: [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.19639793 0.09127954 0.11733288 0.07598273 0.06545106 0.06211532 [2,] 0.22773467 0.10981912 0.16052847 0.38101187 0.18353474 0.24072918 [3,] 0.46743388 0.45733836 0.32191178 0.43356107 0.39159746 0.53984221 [4,]

Dear all, I've got a big problem. I try to analyse my data using R with a mixed model ANOVA without useful results and success. My data are as follows: 3 factors (Treatment, Site, Subsite) with 'Subsite' as random factor and nested into 'Site'. I want to analyse the effects of the three main effects (factorial design to a specified degrees (2)) with the interactions between

I forgot to mention that the problem has nothing to do with the kernel being tainted (due to VMware server). I installed VMware server after being faced with the crashes, to give the Java user a "sandbox" which s/he could crash without taking down the real server. Only then I discovered that VMware server has a memory limitation to 3600MB (I guess because it's not a 64bit