similar to: trellis plot problem with R-1.9.0-1

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

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

Hi all, I have a strange problem and rigth now I can't figure out a solution. Trying to calculate an ANOVA with one between subject factor (group) and one within (hemisphere). My dependent variable is source localization (data). My N = 25. My data.frame looks like this: > ML.dist.stack subj group hemisphere data 1 1 tin left 0.7460840 2 2 tin left

Hi everyone, I'm a newbie to R and to linear mixed effects modeling, so please have mercy :-) Just wanted to check, whether what I'm doing is alright. I've collected data concerning tonotopic organization of the auditory cortex in humans, and I have approximately 1400-1800 data/time points per person (13 in total). Observations were made how the focus of neuronal activity changed

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>

Hello, I've got a big problem. I'm using R for geostatistical analyses, especially the field-package. I try to generate plots after the kriging process with help of image.plot(..., col=terrain.colors, ...). Everything works fine, but I want to reverse the color-palettes (heat.colors, topo.colors or gray()) to get darkest colors at highest data-values instead the other way round. Could

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

Dear R community members, I think I am asking a very simple question, but I really looked up in the faqs and manuals and found nothing helpful. I am trying to plot a data frame with the following structure (this is just a small extract): glo conc odor line series X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 X11 X12 X13 1 0 AIR LN1 UP -0.488

Dear R community members, I am very new to R, so I hope I didn´t oversee the answer to my question in the manuals and faqs. I am trying to create a histogram containing several data sets (I would like to see how the dependent variable is distributed for the different groups). This should be possible using the function multhist contained in the plotrix package. (I am using R 2.2.1, in windows XP).

Our Tyan VX50 server, which is otherwise stable, can be crashed by a specific user who starts a Java application. Relevant facts are: x86_64 CentOS 4.3 fully updated (but otherwise default config) 8 Dual-Core Opteron 870 (2 GHz), giving 16 cores 32 GB memory (running at DDR333) The panic is reproducible within a couple of minutes. I tried the 2.6.9-42 largesmp kernel from J.Baron's website

Hi, I would like to report a bug in OpenSSH-1.2.2 (release) under Linux. Under certain conditions the sshd client process segfaults while doing the password authentification. I have observed the behaviour with glibc 2.0.7 (non-PAM), glibc 2.1.1 (PAM) and glibc 2.1.2 (PAM), when 'RhostsRSAAuthentication yes' is chosen in sshd_config. It appears to happen regardless whether PAM is used or

Hi, Here is a patch for release 1.2.2 that allows the use of '=' instead of whitespace when specifying options. For options on the commandline, it can be useful to be able to avoid whitespace in some situations. best regards and thanks for the patch regarding segfaulting with PAM, Stefan ------------------------------------------------------------------- Email: Stefan.Heinrichs at

hi rsync-list, i have an authentication problem with my rsync-server. @ERROR auth failed on module XXX i`ve searched for an solution that fixes this problem in the mailinglist archives of rsync. but no one solves my problem. i`m running rsync version 2.5.6cvs protocol version 26 on debian 3.0 with kernel 2.4.22 my rsyncd.config looks like this: #######################################

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 all! Although the Psycho-Toolbox for Matlab is free software, Matlab isn't. I'm planning to do an experiment where it's essentail to travel to the subjects, not let the subjects come to where the Matlab licences are :-( So I need to use a free software for my experiment if I don't want to by an extra Matlab licence (which I don't want to). Did anyone ever try to do

I seem to have stumbled across another vulnerability in DIP. It appears to allow any user to gain control of arbitrary devices in /dev. For instance, I have successfully stolen keystrokes from a root login as follows... (I could also dump characters to the root console) $ whoami cesaro $ cat < /dev/tty1 <------ root login here bash: /dev/tty1: Permission denied

I have some doubts about the validity of my procedure to estimeate linear contrasts ina a factorial design. For sake of semplicity, let's imagine a one way ANOVA with three levels. I am interested to test the significance of the difference between the first and third level (called here contrast C1) and between the first and the seconda level (called here contrast C2). I used the following