similar to: Wide confidence intervals or Error message in a mixed effects model (nlme)

I am analyzing a dataset on the effects of six pesticides on population growth rate of a predatory mite. The response variable is the population growth rate of the mite expressed as ln(Nfinal/Nstarting) of the mite, where N final the population of the mite at the end of the experiment and N starting the population of the mite at the beginning of the experiment. Each of the six treatments was ran

I am running a binomial glm with response variable the no of mites of two species y->cbind(mitea,miteb) against two continuous variables (temperature and predatory mites) - see below. My model shows overdispersion as the residual deviance is 48.81 on 5 degrees of freedom. If I use quasibinomial to account for overdispersion the dispersion parameter estimate is 2501139, which seems

I am running an lme model with the main effects of four fixed variables (3 continuous and one categorical – see below) and one random variable. The data describe the densities of a mite species – awsm – in relation to four variables: adh31 (temperature related), apsm (another plant feeding mite) awpm (a predatory mite), and orien (sampling location within plant – north or south). I have read

Hello, I am struggling to figure out how to analyse a dataset I have inherited (please note this was conducted some time ago, so the data is as it is, and I know it isn't perfect!). A brief description of the experiment follows: Pots of grass were grown in 1l pots of standad potting medium for 1 month with a regular light and watering regime. At this point they were randomly given 1l of one

I am running a two-way anova with Type III sums of squares and would like to be able to understand what the different SS mean when I use different contrasts, e.g. treatment contrasts vs helmert contrasts. I have read John Fox's "An R and S-Plus Companion to Applied Regression" approach -p. 140- suggesting that treatment contrasts do not usually result in meaningful results with Type

Hi all, I have been trying to calculate Type III SS in R for an unbalanced two-way anova. However, the Type III SS are lower for the first factor compared to type I but higher for the second factor (see below). I have the impression that Type III are always lower than Type I - is that right? And a clarification about how to fit Type III SS. Fitting model<-aov(y~a*b) in the base package and

Does anyone know how to extract residuals in lmer? Here's the error I get: > crop.lme=lmer(response~variety*irrigation*pesticide+(1|rep)+(1|rep: pesticide)+(1|rep:pesticide:irrigation), crop.data) > qqnorm(crop.lme) Error in qqnorm.default(crop.lme) : y is empty or has only NAs > resid(crop.lme) NULL Thanks! --Jake

Full_Name: Menelaos Stavrinides Version: 1.3. 1 OS: Windows 98 Submission from: (NULL) (193.129.76.90) When model simplification is used in glm (binomial errors) and anova is used two compare two competitive models one can use either an "F" or a "Chi" test. R always performs an F test (Although when test="Chi" the test is labeled as Chi, there isn't any

Hi every body: I need some help with kendall.global. The example in the manual seems not working well, and cannot used with my data, always the same error. data(mite) > mite.hel <- decostand(mite, "hel") > > # Reproduce the results shown in Table 2 of Legendre (2005), a single group > mite.small <- mite.hel[c(4,9,14,22,31,34,45,53,61,69),c(13:15,23)] >

This is not a request for help, but I thought the readers of this mailing list might appreciate the following, anyway. I have been working with US EPA's (that's the United States Environmental Protection Agency) Office of Pesticides on a preliminary assessment for the effects of exposure to multiple organophosphate pesticides on human health (as required by the Food Quality Protection

Good morning, I wrote a little code in R which has to show two graphs but I can get only one. How can I adress the graphs in two files? Second, I'd like, always in the same code, to add a legend to a graph. Better, I'd like to put in such a legend a new item whose color could remind the colour ol the columns it refers to in the plot. I wrote: leg.txt<-c("control

Hello, I am trying to obtain f-values for response (independent) variables from a CCA performed in vegan package, to see which ones of them have significative influence in my dependent variables (like the manual selection in canoco), but I can't find any function (or package) that do such a thing. The dependents variables are species data, and the independents are ambiental data. Than you.

Dear all i?m working with a program i?ve made in R (using functions that others created) to run my program i need a sample. if i generate the sample using for example, rnorm(n, mu, sigma) i have no problem but if i obtain a sample from a column in excel and i copy it, the program says that there is a mistake: it says "Error en `[.data.frame`(data, indices) : undefined columns

Is there a soft phone for PocketPC or iPaq? If not, is someone working on it? I will be more than willing to contribute my mite if needed. Thanks, -- sudhir

I am trying to mirror only select directories from one server to another with rsync through its daemon mode. Server A /export /home /A-do /A-not /A-copy /das /htdocs /docs /em /htdocs /docs /psb /htdocs /docs /X-do /X-not /X-copy Server B /export /home /extranet

Ok here's the deal. I set up Icecast and Oddcast, got them to communicate properly, but cannot access the live audio remotely. On other computers in the LAN it works perfetly. Currently I operate wirelessly behind a D-Link 514 router (a mite outdated by now, yes?) and have opened all necessary ports, turned my firewalls off, set this computer as a DMZ, pretty much everything I can think of -

I have a table with some data structured something like this: book ---- title (varchar) pages (int) completed_on (date) While I can find the shortest book completed in each year with something like this: select min( pages ), year from book group by year I would also like to get the title of the shortest book in each year. Is there any way to do this by modifying the previous query to pull

> > Message: 5 > Date: Sat, 28 Apr 2007 11:37:51 -0700 > From: "J.H." <warthog19 at eaglescrag.net> > Subject: Re: [syslinux] USB stick meets ISO system > To: Geert Stappers <stappers at stappers.nl> > Cc: syslinux at zytor.com > Message-ID: <1177785471.25140.29.camel at localhost.localdomain> > Content-Type: text/plain > > On Sat,

Hi all, I've installed CentOS 4.4 in a box which previously had dual boot (XP SP2 and Ubuntu) I can run CentOS without problem but I can't boot XP SP2. When I try it the system simply hangs... [root at seth ~]# fdisk -l Disco /dev/hdb: 40.0 GB, 40020664320 bytes 255 cabezas, 63 sectores/pista, 4865 cilindros Unidades = cilindros de 16065 * 512 = 8225280 bytes Disposit. Boot Start

Aside from the replay attacks discussed, there are some other attack vectors on the cookie_session store. I appreciate (and admire!) Jeremy''s good humor on all of this: > Planting the seed here led to quick ripening and plenty of pesticide. > Thanks for the fish, all. > > jeremy Anyway, here''s what we came up with: 1. Brute Force SHA512 can be computed _very_ fast.