similar to: GLM contrasting question

When syncing one backup system (running dirvish, but that's not really relevant) to a second long-term backup system which uses btrfs snapshots, I get the following output: rsync: link "/parasite.oob.telegraaf.net/vservers-var-www/tree/elo/var/www/lokaalx/openx-htdocs/etc/xmlcache/.svn/prop-base/cache_etc-changes_schema-tables-core-12934.bin.svn-base" (in backup-aquagirl) =>

I have a host-parasite association matrix in which parasite species are rows and host species columns and cells contain the frequency of interactions. Some parasites are associated with many hosts, and some hosts harbor several parasites, and I want to repeatedly select only one single representative host per "generalized" (multi-host) parasite to create a new matrix in which no hosts

Dear all, I have tried to calculate a GLMM fit with lmer (lme4) and glmmPQL (MASS), I also used glm for comparison. I am getting very different results from different functions, and I suspect that the problem is with our dataset rather than the functions, but I would appreciate help in deciding whether my suspicions are right. If indeed we are attempting the wrong type of analysis, some

Hi everyone, I determined the presence of three types parasites in a passerine bird over two years. I would like to create a bar chart that shows the proportion infected on the y and year/parasite on the x such that each type of parasite is grouped together (single label) and a bar for each year . This would show if there have been changes in the prevalence of a the parasite over two years.

Hello all, I am hoping to use survival analysis to examine whether parasite attack increases nest death in a species of social wasp. I therefore have data for 1. Whether the nest "died" in the 6 week census period ("Status", where 1=died, 0=survived) 2. The day number of death/last recorded day it was observed alive. 3. Whether the nest was attacked by the parasite (0/1 as

Dear All, Apologies for sending this email to both list, but at this point I'm not sure which one could help me the most. I have 4 sets of data, 1 test and 3 different sets of controls. The measurements are binary, with a matrix of 0 and 1 I'm measuring across time (rows, ~815) the behaviour of organelles in the cell by microscopy in response to different stimuli (several measurements

Hi, I have a problem concerning discrepances between R (which I use) and Statistica (which uses my supervisor). I can't say what is the origin of these differences but unfortunately my supervisor doesn't know that either. Our response variable is number (or presence/absence) of parasites in rodents and explanatory variables are presence/absence of several alleles. The rodents were

Announcing the the 0.2.0 release of parasite and the camping_generator. Parasite (from the project page -- http://parasite.rubyforge.org/): Camping app developers no longer have any reason to envy their Ruby on Rails friends: Parasite brings generators, environments, and other Rails-y goodness to the world of Camping app development. Parasite is currently at version 0.2. The Parasite package is

Dear all, I'm building binomial mixed-model using lme4 package. I'm able to obtain outputs properly except when I include two particular variables: date (from 23 to 34; 1 being to first sampling day) and Latitude (UTM/100 000, from 55.42 to 56.53). No "NA" is any of those variables. In those cases, I get the warning message: "nlminb failed to converge" I tried to

Hi R people! I am looking for some suggestions writing an if-else function. The idea is to characterize different plots containing counts of variables (here parasites). If a plot has a count equal or higher than 4 for any parasite the function should return a 1 else a 0. Later I can loop the function over all plots. Here I have a little subset of my data: VariablePAR Plot1

Hi, I have found quite a few posts on normality checking of response variables, but I am still in doubt about that. As it is easy to understand I'm not a statistician so be patient please. I want to estimate the possible effects of some predictors on my response variable that is nº of males and nº of females (cbind(males,females)), so, it would be:

Hi, I am trying to fit a non-linear model for a parasite dataset. Initially, I tried log-transforming the data and conducting a 2-way ANCOVA, and found that the equal variance of populations and normality assumptions were violated. Gaba et al. (2005) suggests that the Weibull Distribution is best for highly aggregated parasite distributions, and performs better (lower type 1 and 2 error rates)

I am analysing parasite egg count data and am having trouble with glm with a negative binomial family. In my first data set, 55% of the 3000 cases have a zero count, and the non-zero counts range from 94 to 145,781. Eventually, I want to run bic.glm, so I need to be able to use glm(family= neg.bin(theta)). But first I ran glm.nb to get an estimate of theta: > hook.nb<- glm.nb(fh,

Is there a cable length limit to use with a fxsmod (TDM400B) ? I noticed that if I use a very long cable 148 feet, I can hear some parasite : I hear the radio ! If I use a short cable, there is no problem... Thanks for help, St?phane HENRY

Hi, Default register allocator tries to reuse the same registers over and over again even if register file have a plenty of registers to use. This creates parasite false dependencies and makes scheduling less effective. How to instruct allocator (may be override some virtual function in mine backend?) that it is profitable to use maximum number of available registers with minimal dependencies?

Hi, I am running a mixed-effect model with a nested-random effect. I am interested in gut parasites in moose. I has three different type of treatment that I applied to moose which are from different "families". My response variable is gut parasites and the factors are moose families which is nested within treatment. My data is balanced. To answer this question, I used the lme function

All, I see this item on the merb-0.0.8 release notes Added rails session parasite mode. I am assuming merbs can piggy back off the rails current session. How does one achieve that ? Thanks Fernand -------------- next part -------------- An HTML attachment was scrubbed... URL: http://rubyforge.org/pipermail/merb-devel/attachments/20070112/0803f46e/attachment.html

Dear R help, I am having a problem with the Design package and my problem is detailed here. I fit a cox model to my data and validate the Somer's Dxy using the Design package. (Because of computation time problem, i only try 10 bootstrap samples for the time being) This is the model without stratification: > library(Design) >

Hi all, I have run a (glm) analysis where the dependent variable is the gender (family=binomial) and the predictors are percentages. I get a warning saying "fitted probabilities numerically 0 or 1 occurred" that is indicating that quasi-separation or separation is occurring. This makes sense given that one of these predictors have a very influential effect that is depending on a

Hi, I'm totally new to R so I apologise for the basic request. I am looking at the incidence of a disease over two time periods 1990-1995 and 2003-2008. I have counts for each year, subdivided into three disease categories and by males/females. I understand that I need to analyse the data using poisson regression and have managed to use the pois.daly function to get age-sex adjusted rates and