similar to: Analysis of unbalanced data in nlme or car

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 (ranges from negative to positive) and the exploratory variable is a categorical variable (treatment). The experiment was blocked in time (3 blocks / replicates per block) and it is unbalanced - at least 1 replicate per block. I am

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

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

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

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

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

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

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

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

--- -------------------------- CISTI'2016 ----------------------------- 11? Conferencia Ib?rica de Sistemas y Tecnolog?as de Informaci?n 15 a 18 junio de 2016, Isla Gran Canaria, Espa?a http://www.aisti.eu/cisti2016/ ------------------------------------------------------------------- Nos satisface invitar a la comunidad acad?mica y empresarial a presentar trabajos a

--- -------------------------- CISTI'2016 ----------------------------- 11? Conferencia Ib?rica de Sistemas y Tecnolog?as de Informaci?n 15 a 18 junio de 2016, Isla Gran Canaria, Espa?a http://www.aisti.eu/cisti2016/ ------------------------------------------------------------------- Nos satisface invitar a la comunidad acad?mica y empresarial a presentar trabajos a

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

> > 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,

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

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

Hi, I was wondering if there is more flexibility in the output of dendrograms when plotting a hclust object. I can't seem to find information on how to change the default output of a "hanging" style tree with the axis on the right to a left-to-right plot with and axis on the bottom. Example code follows: library(vegan) #loads the "vegan" module that compuptes ANOSIM

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

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