similar to: Multivariable model

Dear all, How do I make R show the R-squared (deviance explained by the model) in a logistic regression? Below is how I write my syntax. Basically I want to investigate density-dependence in parasitism of larvae. Note that in the end I perform a F-test because the dispersion factor (residual deviance / residual df) is significantly higher than 1. But how do I make R show the

Dear all I have made some lattice graphs (bwplot(dead ~ treat | group,...). 'group' is a factor with three levels (Artemia larvae, Abalone larvae and Abalone spat), and the result is a graph with three panels corresponding to the levels in 'group', ordered in alphabetical order from bottom to top, as expected. How does one re-order the order in which the levels of

Background: OS: Linux Ubuntu Dapper R : Version 2.3.1 (2006-06-01) Emacs 21.4.1 ESS Colleagues This is neither a bug report, nor a complaint, so I think this is the appropriate list. I am getting an odd font problem with a pdf generated using pdf(file = "../figures/PCNM.models.larvae.fluor.zoo.stability.2001.and.2002.pdf", title="", family="serif",

Dear All, I'm confortable with xyplot(...) and panel.lmline(...) statements (at least I thought I did :). I've used the following code to plot the decline in log-abundance of fish larvae (no.larvae) with age (age.cls, 4 to 27 days-old) for specific dates of sampling (day, 9 dates). I further plotted data with different colors and regression lines for ages 5-14 d and 17-23 d in a 7-by-1

I have taken the liberty of including the R-help mailing list on this reply as that is the appropriate place to discuss lmer results. On 4/5/06, Andreas Svensson <andreas.svensson at bio.ntnu.no> wrote: > Hello again > I have now recieved some helpful hints in this matter and will summarize them but first let me reiterate the problem: > > I had two treatments: 2 types of food

> Background: > OS: Linux Ubuntu Dapper > release: R 2.3.0 > editor: GNU Emacs 21.4.1 > front-end: ESS 5.2.3 > ----------------------------- > Colleagues > I have two spatial datasets (latitude, longitude, fish eggs) and (latitude, longitude, fish larvae) at the same 280 stations (i.e. 280 cases). I want to determine if the 2 datasets are spatially correlated. In other

Hi there, I'm using the lattice package to create an xy plot of abundance vs. depth for 5 stages of barnacle larvae from 5 species. Each panel of the plot represents a different stage, while different loess smoothers within each panel should represent different species. However, I would like depth to be on the y-axis and abundance to be on the x-axis, because this is more intuitive as an

Hi, I have some difficulty in figuring out whether I am doing correct or not. A brief introduction about the work: It is a light/dark choice test conducted in insect larvae.  The response is binary (0- present in dark area, 1-present in light area) and the experiment is run for 15 min, so there are 15 repeated measurements per individual larva at 1 min intervals.  The factors which affect

*Thanks* all those who took the time to help me (even if the "question" was not related to - the use of - R). Now I think I can soundly make my point w/ the referee (can I use your replies? If so I intend to properly cite its use?!?). Regards, Eduardo Esteves ps - Sorry for not explaining the "biological details" of my posting: RNA/DNA is the ratio of RNA content to

Dear All, Help is needed! I have a matrix with frequencies of fish larvae per length class (var. sl) and age-group (var. median.no) obtained with >k<-table(cut(sl,(5:22)),median.no) >k[2:5,1:5] #to ilustrate k 4 5 6 7 (6,7] 3 1 0 0 (7,8] 3 0 1 0 (8,9] 3 4 3 5 (9,10] 3 15 7 13 from this matrix I would like to obtain the mean age per length class i.e.

Dear R list, I have a sequence of weekly observations of number of adults and larvae in various size classes from a butterfly population living in a subtropical area with pronounced wet and dry seasons. Wet and dry seasons are each defined 26 weeks long with fixed start and end dates. The data span 103 weeks (two seasons each of wet and dry) with some missing weeks. What I would like to do is

Hello, I have an unbalanced mixed model design with two fixed effects "site" (2 levels) and "timeOfDay" (4 levels) and two random effects "day" (3 consecutive days) and "trap" (6 unique traps, 3 per site). The dependent variable is the body length ("BL") of insect larvae from 7 to 29 individuals per trap (104 individuals in total). To account

Dear R users, x-values (EI) = Adw, EG1, LA1, Ad1, LA2, LA3...(14 levels, insect stages) y-valus = antpop within the boxplot function x-values are ordered alphabetically Idea: x-values ranked by list order (insect stage: Egg stage 1 is followed by Larvae 1 and not by Egg stage 2 as it would be in an alphabetically order) Problems with the order(tapply()) function: variable lengths

Hi Everyone, I posted a similar question about a week ago, but haven't gotten any replies -- I'm afraid that's because my previous question was too vague. Let me try again with a more specific question, and I hope someone can help. NOTE, I know I should be using the newer lme4 package, I just haven't had a chance to update my version of R yet, so the question below relates

Dear Sebastian Luque (and All R Users) With the following code I managed to plot different characters and regression lines for panels 2 ("Day of year 101") and >4 ("Days of year" 151, 157 and 172): xyplot(log(no.larvae)~age.cls|factor(day),data=mortal, layout=c(7,1),aspect=5/3, xlab="Age class (d)",ylab="Ln(Abundance)",

Hello I am trying to analyse mortality data from fish larvae using lme from the package nlme as well as using lmer in the package lme4 Response is DeathDay Fixed factor is Treatment Random factors are Clucth, Cup Design: Cup nested in Clutch If I do this in lme, I use the syntax: model1 <- lme(DeathDay ~ Treatment, random=~ 1 | Clutch/Cup) summary(model1) I get the first part of the output,

Dear All, (this msg is a statistics/computing question to the list) I'm trying to implement a modern-version of a (classic) "several-step protocol" in Fishery Biology (due to Bhattacharya, 1967): analysis of length-frequency distribution of fish larvae to id cohorts and later estimate growh rates! I've trouble with the 1st step: using kernel density estimation (KDE) function

Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 --------------------------------- Colleagues I have solved the miscoding with subscript and groups in panel.superpose. I was neglecting to pass in the variables correctly. "sardine.egg.T.S.space.2001.and.2002.exp" <- function() { library(lattice) library(grid) ##

Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 --------------------------------- Colleagues I'm plotting fish egg densities in temperature-salinity space (i.e. mean temperature ~ mean salinity | fish egg densities) for two survey years. The bivariate plot is overlaid on a temperature ~ salinity plot (not means) showing the distribution of

I am a Professor of Statistics at Indira Gandhi Krishi Vishwavidyalaya, Raipur, India. While teaching in class about analysis of variance using R, I was doing a one-way analysis for the two data-sets given below in the R-class. I got a typical error in "emmeans" package, please help: Data-set-1: -------------- Medley and Clements (1998) investigated the impact of zinc contamination (and