similar to: Mixed model ANOVA with a nested design

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

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

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

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

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 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: #######################################

Hey, I have an double loop like this: chunk <- list(1:10, 11:20, 21:30) for(k in 1:length(chunk)){ print(chunk[k]) DummyCatcher <- NULL for(i in chunk[k]){ print("i load something") dummy <- 1 print("i do something") dummy <- dummy + 1 print("i do put it together") DummyCatcher = rbind(DummyCatcher, dummy) } print("i save a chunk

On Fri, Jul 29, 2005 at 10:08:13AM +0100, james wrote: > SVN root: svn://svn.xapian.org/xapian > Changes by: james > Revision: 6355 > Date: 2005-07-29 10:08:13 +0100 (Fri, 29 Jul 2005) > > Log message (6 lines): > omindex.cc: add --preserve-nonduplicates / -p option to not delete any > documents that aren't updated, in replace duplicates mode

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

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 :

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).

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

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

For a while libguestfs.org/builder has redirected to builder.libguestfs.org. This change just makes virt-builder go direct to the subsite instead of via the redirect. --- builder/libguestfs.conf.in | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/builder/libguestfs.conf.in b/builder/libguestfs.conf.in index 633a0abaf..3b9e3854f 100644 --- a/builder/libguestfs.conf.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

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

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

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