similar to: Validating a mixed-effects model

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

Hello there, I am using a Gigabyte H55M-UD2H (Intel H55-Chipset) together with a Core-i5-650. According to Gigabyte Support the Mainboard does support VT-d. Well, actually support says, that VT-d does work on Windows-7 in XP-Mode. They do not test and support Linux though. Did anybody here have success using VT-d on any H55-Board? Perhaps it is some misconfiguration on my side. I tried a

Dear all, I have a performance problem in terms of computing time. I estimate mixed models on a fairly large number of subgroups (10000) using lme(.) within the by(.) function and it takes hours to do the calculation on a fast notebook under Windows. I suspect by(.) to be a poor implementation for doing individual analysis on subgroups. Is there an alternative and more efficient way for doing

Dear mailinglist, I own a Gigabyte motherboard GA 970A UD3 with IOMMU support. Since the update XSA-36 (also part of the latest debian wheezy pkg), the IO-Virtualisation does not work any more as discussed on this mailinglist [0] and [1]. I like to ask, if there is an "official" solution in sight. I''m not sure about my alternatives. How "dangerous" is the

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 :

Hi Trying to get WinXP HVM to install as domU on AMD. Got latest Xen Kernel and Xen package (NOT unstable or debug). When start, a window pops up for few seconds.. then it dies and another one comes up... and it disappears too. (Xen is trying to restart and giving up because the restart happens in 2 seconds..) It seems DomU is not able to read/understand something from Cirrus VGABIOS.

Hi Trying to get WinXP HVM to install as domU on AMD. Got latest Xen Kernel and Xen package (NOT unstable or debug). When start, a window pops up for few seconds.. then it dies and another one comes up... and it disappears too. (Xen is trying to restart and giving up because the restart happens in 2 seconds..) It seems DomU is not able to read/understand something from Cirrus VGABIOS.

I am wondering if it would be possible (and safe) to develop online banking or other security-related applications with Ruby on Rails, instead of using J2EE frameworks or .NET? Are there any particular implications needed to write e-banking apps with Rails? Security problems? So... is anyone developing these types of applications with RoR? Any advantage or disadvantage? Personal experiences are

Hi all, Does anyone have an example of an autoregressive (AR) time-series model specified as a state space model in R? That is, I want to go beyond the locally linear (constant) model, and fit the following Gaussian AR state process model: Xt = a + (1+b)*Xt-1 + epsilon ,where the model for the observation process is Yt = Xt + tau I have information of the tau's (observation variance)

Does anyone know how to make a plot like Fig. 4 of Paik, et al., New England Journal of Medicine, Dec. 30, 2004? Given survival data and a covariate, they plot a curve giving "Rate of Distant Recurrence at 10 Yr (% of patients)" on the y-axis versus the covariate on the x-axis. They also plot curves giving a 95% confidence interval. Thanks very much. -Ben The information in this

Hello all I would like to visualize the hazard ratios of a cox proportional hazards model. I have seen some good examples in the New England Journal of Medicine: http://content.nejm.org/cgi/content/short/353/26/2747/F3 http://content.nejm.org/cgi/content/short/350/26/2654/F2 (I hope these are open access, but I am not sure) And something similar in the book 'R Graphics' by Paul

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,]

Hey all, I am just born today in the r language and I discover that it the language that I need for my research. I don't too much yet about this language. I want to know if I can find in this language and install the package data set related to " The Framingham Heart Study", see references here: http://content.nejm.org/cgi/content/full/357/4/370

Hi list I tried to install package lattice on FreeBSD, but install.package() does not seem to find it in the repositories, even trying different mirrors. Could it be that lattice package is not available for FreeBSD? -- Armin Goralczyk, M.D. -- Universit?tsmedizin G?ttingen Abteilung Allgemein- und Viszeralchirurgie Rudolf-Koch-Str. 40 39099 G?ttingen -- Dept. of General Surgery University of

Hi list Maybe someone can help with the following problem (thanks in advance): In a function I have a plot and want to add symbols/text only when indicated by a logical vector (which was generated by the function before, not manually like in the following example): plot(1:10, 1:10) lv <- c(T,T,T,F,F,F,T,T,T,F) text(1:10, 1:10, '*', pos = 3) In this example the asterisk is plotted at

Hello Is there an easy way, i.e. a function in a package, to calculate the area under the curve (AUC) for drug serum levels? Thanks for any advice -- Armin Goralczyk, M.D. -- Universit?tsmedizin G?ttingen Abteilung Allgemein- und Viszeralchirurgie Rudolf-Koch-Str. 40 39099 G?ttingen -- Dept. of General Surgery University of G?ttingen G?ttingen, Germany -- http://www.chirurgie-goettingen.de

Basically all platforms have a proper "mkdir -p". Autoconf has install-sh as fallback when using AC_PROG_MKDIR_P, which is also hardened against races if it is recent enough. Update that script and use MKDIR_P whereever possible. Installs fine in my tests and avoids the "mkdir: cannot create directory '?': File exists" messages along the way. Greetings, Eike --

hmm, the mailing list engine seems to have eaten my attachmen (and my signature....) trying again ----- Forwarded message from Michael Adam <ma@sernet.de> ----- Date: Mon, 4 Aug 2008 22:43:55 +0200 From: Michael Adam <ma@sernet.de> To: Dag Ager <dag.ager@enders-ing.de> Subject: Re: [Samba] unable to build 3.2 on HP-UX Hi,