similar to: GLMM random effects

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

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

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

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

Trying to plot multiple lines from a simple matrix with headers (eight observations per column). I will be doing a number of these, with varying numbers of columns, and do not want to enter the header names for each one (I got frustrated and just wrote them out, which did work). Data reads fine, first plot is fine, but when i use the code at the bottom for a for i loop it tells me that x and y

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

Dear R-help, First of all, thank you VERY much for any help you have time to offer. I greatly appreciate it. I would like to write a function that, given an arbitrary number of factors from a data frame, tabulates the number of occurrences of each unique combination of the factors. Cleary, this works: > table(horse,date,surface) <SNIP> , , surface = TURF

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

Our software company is looking for a technical support/customer service agent This person should be knowledgeable with FreeBSD (installation and troubleshooting), Windows, hardware support (diagnosing hardware problems, replacing memory, cards, etc...), building out hardware (servers/workstations) for customers. Experience with JMS, and/or FoxPro desirable but not required. Some light book

Greetings, I have "inherited" a cDNA macroarray dataset that is structured as follows. Three different stressors were tested. For each stressor, there are two treatments (control and stressed). For each treatment, two biological replicates exist, and these are paired (i.e., there is a stressed array for colony A and a control array from this same colony). For one of these samples,

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

Members of this mailing list may be interested in the following new book: Beginner's Guide to GLM and GLMM with R. - A frequentist and Bayesian perspective for ecologists - Zuur AF, Hilbe JM and Ieno EN This book is only available from: http://www.highstat.com/BGGLM.htm This book presents Generalized Linear Models (GLM) and Generalized Linear Mixed Models (GLMM) based on both

Here is a patch for clustered xAPIC and Unisys ES7000 support in Xen. To turn the features on add CONFIG_X86_CLUSTERED_APIC and CONFIG_ES7000 defines to xen/include/asm-x86/config.h Signed off by Aravindh Puthiyaparambil <aravindh.puthiyaparambil@unisys.com> _______________________________________________ Xen-devel mailing list Xen-devel@lists.xensource.com

Hello all I'm trying to do a simulation that involves identifying the minimum point between two peaks of a (usually) bimodal distribution. I can do this easily if there are only two peaks: CnBdens<-density(Ys/Xs) #probability density function for ratio of Ys to Xs for(p in 1:512) ifelse(CnBdens$y[p]>CnBdens$y[p-1],peak1<-p,break) #identifies first peak in probability

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

How do I know which edible ink cartridges to get for my printer? I'm a cake deco hobbyist and I want to start printing photo cakes for friends and family. The edible photo systems online look pretty expensive, but I hear that you can buy a cheap Canon printer and then just buy the edible inks and paper. But how will I know which edible ink cartridges match which Canon printer??? Please offer

Greeting R helpers J I am not familiar with R but I have to use it to analyze data set that I have (30,000 20,000) I want to change the structure of the dataset and I am wondering how that might be possible in R A main data looks like this: some entities are empty Age No. Age No. Age No. Center1 5 2 8 7