similar to: Analysis of poorly replicated array data

The Folding at Home client seems to be a wonderful stressor of the CentOS server systems I've been burning in. I've been running two copies of the client on physical dual-core and single die dual-core machines for a few days now. Nary a whimper from the machines and only a slight elevation in temperature (they're 1RU rackmount systems with aggressive/loud fans). In contrast, I

Hello, I was using R (v.2.11.1, 32 bits) and I did the upgrade to R (v.2.12.0, 64 bits). I followed the instructions in R´s FAQ (What´s the best way to upgrade, question 2.8) and updated my packages. However, now, I can´t use the library "vars". When I call it, there is an error message concerning the package "MASS" which couldn´t be updated because it seems to be no more

Dear list, I am looking for an implementation of random rotation matrix generation in R to do a rotation test: I want to use the matrices to create random multivariate normal matrices with common covariance structure and mean based on an observed data matrix. The rRotationMatrix-function in the mixAK-package is an option, but as far as I can tell I need to draw rotation matrices with determinant

Hello, I am looking for some help with this question: how could I test structural breaks in a instrumental variables´s model? For example, I was trying to do something with my model with three time series. tax_ivreg <- ivreg(l_y ~ l_x2 + l_x1+ dl_y | lag(l_x2, -1)+lag(l_x2, -2)+ lag(l_x1, -1)+lag(l_x1, -2)+lag(l_y, -1)+lag(l_y, -2), data=tax1) summary(tax_ivreg) ## after estimating it,

I need to rewrite a software in R, that runs on a cluster. I thought Rmpi will be good for me but I can't find any help other than the Rmpi-manual that actually only describe the functions in the Rmpi package. Can someone point me to some usefull guide? For example, I would like to run a for-statement on several processors (a subset of the statement on each processor) but I can't figure

Hi folks, is it possible to read more things about R behavior in multiprocessor / multihost environment ? Is there any distributed computation project associated to it ? Laurent

According to the documentation at http://docs.puppetlabs.com/puppet/2.7/reference/lang_conditional.html#if-statements, you can test for truth with syntax like this: if $is_virtual == ''true'' { AFAICT that''s not the case. The quotes around ''true'' make the comparison always fail. Here''s my test case class users {

Dear all, I would like to use the values in vegaggr.BLMCMR02$colony str(vegaggr.BLMCMR02) `data.frame': 1678 obs. of 3 variables: $ vegtype : Factor w/ 27 levels "2010","2020",..: 3 4 5 19 4 5 19 5 $ colony : Factor w/ 406 levels "0","1","10","100",..: 1 1 1 1 2 2 2 $ Totvegproparea: num 0.00055 0.03956 0.95705

I just discovered the great piece of software that is available with the survey package. Many thanks and 'Hats off' to Thomas Lumley. While package survey covers analysis of features of objects sampled (in clusters, strata) I could not find analysis of features of repeated relationsships between sampled objects (in clusters, strata). My understanding is that it is not adequate to treat

Hello, I have a couple questions regarding generalized linear mixed models specifically around fitting the random effects terms correctly to account for any pseudo-replication. I am reading through and trying to follow examples from Zuur et al. Mixed Effects Models and Extensions in Ecology with R, but am still at bit unsure if I am specifying the models correctly. Background information: Our

Hi All, I'm working with some SAS code to analyze an experiment set up as follows: 66 subjects (colonies) treated with a random treatment (1-8) and measured at three time points. The data structure looks like: input colony tmt y1 y2 y3; y=y1; date=*1*; output; y=y2; date=*2*; output; y=y3; date=*3*; output; datalines; 1

I have a data frame containing the Id, Mother, Father and Sex from about 10,000 animals in our colony. I am interested in graphing simple family trees for a given subject or small number of subjects. The basic idea is: start with data frame from entire colony and list of index animals. I need to identify all immediate relatives of these index animals and plot the pedigree for them. We're

Hi, I have a data set with 999 observations, for each of them I have data on four variables: site, colony, gender (quite a few NA values), and cohort. This is how the data set looks like: > str(dispersal) 'data.frame': 999 obs. of 4 variables: $ site : Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 1 2 2 ... $ gender: Factor w/ 2 levels "0","1":

I'm not quite sure whether this is a plain browsing problem, or a wins problems, or ... I had wins running on one machine (Samba(2.0.5a)/SCO Unix), but my main file server on another (Samba (2.2.1a)/Linux). I moved the wins support to the Linux server, and at the same time I made the Linux server the Domain Master Browser, but now some devices have dissapeared from the network. The devices

Hi, I have been trying to re-order several items in a trellised barchart display in lattice, but can't seem to figure it out. ###sample code, Stage and Colony have 2 and 3 levels respectively. barchart(Activity ~ Percent | Stage + Colony, data = Percent.df, horizontal = TRUE, layout = c(2,3), xlab = "Percent Time Engaged in Activity", ylab = "Activity") This

I've downloaded Samba 2.2.6 (binaries) from SCOs Skunkware site, but when I try and run smbclient I get a message : dynamic linker : smbclient : error opening /usr/local/lib/libncurses.so.4 Killed I've checked and there is no libncurses.so.4 (or libncurses.<anything> for that matter) on the system. I've also failed to find anything relevant in either the list archive or SCOs

Hello, I am Mudit Raj Gupta, fourth year student of M.S. (Hons.) Chemistry and B.E. (Hons.) Electronics and Instrumentation at BITS-Pilani ( http://www.bits-pilani.ac.in/). I am interested in *Machine Learning and Computaional Intelligence*. I have an interest in implementing various existing algorithms and developing modified/new algorithms related to machine learning and computational

Hello, I have been attempting to run a poisson glmm using lme4 for some time now and have had a lot of trouble. I would say 9 times out of 10 I receive the following warning: CHOLMOD warning: %h Error in mer_finalize(ans) : Cholmod error `not positive definite' at file:../Cholesky/t_cholmod_rowfac.c, line 432 My data are counts of microbe colony forming units (CFUs) collected from

Hi everybody, If I am correct, you can compare a model with random effect with the same model without the random effect by using the nlme function, like this: no.random.model <- gls(Richness ~ NAP * fExp, method = "REML", data = RIKZ) random.model <- lme(Richness ~NAP * fExp, data = RIKZ, random = ~1 | fBeach, method = "REML")

Dear All, I wonder if anyone can please offer any advice on a model including 2 fixed effects and 1 random effect, as well as a covariate? The experimental design is as follows: I have a two by two factor design, where the two factors, Age (A) and Group size (G), both have 2 levels (old or young, and 1 or 3 respectively), and I am interested in the effect of these factors upon a continuous