similar to: Linear Mixed Models with nlme, more than one random effect

Greetings, I have just been given a VPS with a fresh, minimal Centos 4.2 installation. When I ran "yum update" it went on updating some tenths of packages and signalled this error: Updating : perl ####################### [24/88] Updating : rpm-python ####################### [25/88] Updating : MAKEDEV

Dear All, I plan to replace an error disk that is part of an LV. from LVM how-to it could be done with using pvmove to move all PE from old disk to new disk.But the howto also said that pvmove is slow. Anyone has experience using pvmove on 2TB disk? Is it possible to make all PE on the old disk empty so I don't have to do pvmove (assuming that I can make a free space >= 2TB). Thank you in

Hi, I am trying to compile a package i am developing on my home computer where i am running Windows NT 4.0 service pack 6. Everything works fine on my other computer at work running windows 98. I am using the cygwin 1.1 bash shell and tolls with the MingW compiler, version 2.95.2-1 installed seperately and first in the path. I have active perl 6.16 as well. Here is the error I get when i run

Hello, I am using R2.5 under Windows. Looks like the following statement vars <- (obj$sigma^2)*vw in getVarCov.gls method (nlme package) needs to be replaced with: vars <- (obj$sigma*vw)^2 With best regards Andrzej Galecki Douglas Bates wrote: >I'm not sure when the getVarCov.gls method was written or by whom. To >tell the truth I'm not really sure what

Dear all, I try to run the example of lmer and get the following error message. > library(lme4) > example(lmer) lmer> (fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)) [[1]] Error in get(x, envir, mode, inherits) : variable "as.dpoMatrix" was not found This error message is similar to what I get with other models. It looks like the mer class has a slightly

Maybe I do not explain well what I would like to do. I do not want to change the labels of the axis, but the scale. What I want is a general procedure for changing the scale. Its like using a logarithmic scale on a plot. Labels are the same, but the increases of x along the x-axis are defined by a known monotone and continuous function. Florent Bresson ----- Message d'origine ---- De :

Hello R-users, I am using gls function in R to fit a model with certain correlation structure. The medol as: fit.a<-gls(y~1,data=test.data,correlation=corAR1(form=~1|aa),method="ML") mu<-summary(fit.a)$coefficient With the toy data I made to test, the estimate of mu is exactly equal to the overall mean of y which can not be true. But, if I make a toy data with y more than two

On Sep 15, 2004, at 8:40 AM, Henrik Bach wrote: > Hi, > > Has some one committed these patches into CVS? Not yet... checked 10 minutes ago... I'm attaching the cvs diff -u version of them (I've read Reid Spencer post about it) > Paolo, thanks for the patches. In the near future I'll try to compile > them. > > I've installed STLport 6.4.2 with vc6. May I

Hi, I'm try to compute the minimum sample size needed to have at least an 80% of power, with alpha=0.05. The problem is that empirical proportions are really small: 0.00154 in one case and 0.00234. These are the estimated failure proportion of two medical treatments. Thomas and Conlon (1992) suggested Fisher's exact test and proposed a computational method, which according to their table

he copiado los archivos para windows95 en un directorio rnew, donde tengo el ejecutable rjune.exe, al correr este archivo me aparece un menu con las opciones File, workspace, edit, windows y help. La ventana de la consola permacene abierta y no se como se puede iniciar la escritura de un programa ni como compilarlo Yo estoy interezado en este lenguaje ya que trabajo con datos espaciales -

Hi, Has anyone got Gossiptel working with Asterisk? - I am having real problems getting it to register - i'm just getting timeout errors. Thanks --ian

Hello! I have a quite complicated setup. In my network on each interface there is bandwidth limitation for each user. Booth outgoing (on interface itself) and incoming (attached IMQ) traffic. There is main HTB class which limits bandwidth for whole interface and HTB subclasses for each user. Filtering is done with hashing filters. This setup was working correctly. But now in the network I

Hi, I need to use the grant table mechanism, but I can''t find any example regarding how to do it. Could you please provide some examples? Thanks Carlo -- È molto più bello sapere qualcosa di tutto, che sapere tutto di una cosa. Blaise Pascal _______________________________________________ Xen-devel mailing list Xen-devel@lists.xensource.com http://lists.xensource.com/xen-devel

Hi all, I have thousands of strings like these ones: "1159_1; YP_177963; PPE FAMILY PROTEIN" "1100_13; SECRETED L-ALANINE DEHYDROGENASE ALD CAA15575" "1141_24; gi;2894249;emb;CAA17111.1; PROBABLE ISOCITRATE DEHYDROGENASE" and various others.. I'm interested to extract the code for the protein (in this example: YP_177963, CAA15575, CAA17111). I

Hello, How do I get the standard deviations for the random effects out of the lme object? I feel like there's probably a simple way of doing this, but I can't see it. Using the first example from the documentation: > fm1 <- lme(distance ~ age, data = Orthodont) # random is ~ age > fm1 Linear mixed-effects model fit by REML Data: Orthodont Log-restricted-likelihood:

Hi, Has some one committed these patches into CVS? Paolo, thanks for the patches. In the near future I'll try to compile them. I've installed STLport 6.4.2 with vc6. Henrik --- Got Freedom? Software Freedom Day 2004 - 28th of August http://www.softwarefreedomday.org/ --- >From: Paolo Invernizzi <arathorn at fastwebnet.it> >Reply-To: LLVM Developers Mailing List

Hi- I''m trying to implement an auto-complete using Yahoo''s autocomplete server, which basically suggests search completions: http://ff.search.yahoo.com/gossip?output=fxjson&command={your searchterms}'' if you want to see the results in your browser, use xml instead: http://ff.search.yahoo.com/gossip?output=xml&command={your searchterms}'' So, I can

I've just put up (12:00 noon NZ) yet another tmp.zip with the fixes for Peter Dalgaards two bugs. Niels reports that the 3.51 problems are now different but not solved. I'm trying to locate a 3.51 machine and will post another fix. Can I get some feedback on how people want these "patches" delivered? Should I be numbering them and keeping the old patches? I don't want to

Hi you all! this is my strange problem of the day! I use Debian stable, linux 2.6.18-5-686, samba Version 3.0.24. Here's my smb.conf [global] workgroup = NO1KNOWS realm = NO1KNOWS bios name = PBT server string = Rob's Samba dns proxy = no os level = 64 log file = /var/log/samba/log.%m max log size = 50 syslog = 0 panic action = /usr/share/samba/panic-action %d security = user encrypt

Dear All, I have a database of 200 observations named myD. In the dataframe there are a column named code (with codes varying from 1 to 77), a column named "prevalence" with some quantitative measurements are given and an column named Pr_mean, with no values. I would like to set a cycle to compute the average of prevalence values for each different code and store the averages under the