similar to: GLMM - Am I trying the impossible?

Hello, I have a problem concerning estimation of GLMM. I used methods from 3 different packages (see program). I would expect similar results for glmm and glmmML. The result differ in the estimated standard errors, however. I compared the results to MASS, 4th ed., p. 297. The results from glmmML resemble the given result for 'Numerical integration', but glmm output differs. For the

Dear R users, I'm hoping that more experienced users will be able to assist me in examining the model fit of a mixed generalised linear model. The example using the data 'bacteria' within the MASS package will hopefully illustrate what I would like to acheive; library(MASS) library(nlme) attach(bacteria) # y being output and the trt - treatment group being an explanatory variable.

Dear R users I encounter a problem when i perform a generalized linear mixed model (binary data) with the lmer function (package lme4) with R 2.1.0 on windows XP and the latest version of package "lme4" (0.96-1) and "matrix" (0.96-2) both options "PQL" and "Laplace" for the method argument in lmer function gave me the same results (random and fixed effects

I have a data frame that looks like this: > gctablechromonly[1:5,] refseq geometry gccontent X60_origin X60_terminus length kingdom 1 NC_009484 cir 0.6799 1790000 773000 3389227 Bacteria 2 NC_009484 cir 0.6799 1790000 773000 3389227 Bacteria 3 NC_009484 cir 0.6799 1790000 773000 3389227 Bacteria 4 NC_009484 cir 0.6799

Is the labeling/naming of levels in the documentation for the predict.glmmPQL function "backwards"? The documentation states "Level values increase from outermost to innermost grouping, with level zero corresponding to the population predictions". Taking the sample in the documentation: fit <- glmmPQL(y ~ trt + I(week > 2), random = ~1 | ID, family =

Hi, I was recently misfortunate enough to have to use regular expressions to sort out some data in R. I'm working on a data file which contains taxonomical data of bacteria in hierarchical order. A sample of this file can be generated using: tax.data <- read.table(header=F, con <- textConnection(' G9SS7BA01D15EC Bacteria(100) Cyanobacteria(84) unclassified G9SS7BA01C9UIR

Sorry about this one being long, and I apologise beforehand if there is something obvious here that I have missed. I am new to creating my own functions in R, and I am uncertain of how they work. I have a data set that I have read into a data frame: > gctable[1:5,] refseq geometry X60_origin X60_terminus length kingdom 1 NC_009484 cir 1790000 773000 3389227 Bacteria 2

Hello Everybody, I have an eps figure an awesome bacteria and a plot (generated using R) also in eps format. Now it looks like there is space for only one figure and I have to insert the picture of the bacteria into the plot. Is there a way to insert figures (eps/png/jpg) in to plots (may be control over placement of figures in the plot as well?) ? By plots I mean data represented using axes and

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 am interested in correctly testing effects of continuous environmental variables and ordered factors on bacterial abundance. Bacterial abundance is derived from counts and expressed as percentage. My problem is that the abundance data contain many zero values: Bacteria <-

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 Rusers, I have used R,S-PLUS and SAS to analyze the sample data "bacteria" in MASS package. Their results are listed below. I have three questions, anybody can give me possible answers? Q1:From the results, we see that R get 'NAs'for AIC,BIC and logLik, while S-PLUS8.0 gave the exact values for them. Why? I had thought that R should give the same results as SPLUS here.

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

Dear R community, I would like to compare the degree of aggregation (or dispersion) of bacteria isolated from plant material. My data are discrete counts from leaf washes. While I do have xy coordinates for each plant, it is aggregation in the sense of the concentration of bacteria in high density patches that I am interested in. My attempt to analyze this was to fit negative binomial

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

Hi ALL, I have very simple question regarding pattern matching. Could anyone tell me how to I can use R to retrieve string pattern from text file. for example my file contain following information SpeciesCommon=(Human);SpeciesScientific=(Homo sapiens);ReactiveCentres=(N,C,C,C,+ H,O,C,C,C,C,O,H);BondInvolved=(C-H);EzCatDBID=(S00343);BondFormed=(O-H,O-H);Bond+

> library(MASS) > attach(bacteria) > table(y) y n y 43 177 > y<-1*(y=="y") > table(y,trt) trt y placebo drug drug+ 0 12 18 13 1 84 44 49 > library(lme4) > model1<-lmer(y~trt+(week|ID),family=binomial,method="PQL") Error in match.arg(method, c("Laplace", "AGQ")) : 'arg' should be one of

I am new to R, and I couldn't find the answers to my question in a faq. This could however be because I didn't know what to look for...:) I have three classes of data, data for bacteria, archaea and eukaryotes. I wish to display these in a histogram where all of the values are used to calculate each column. But, I want each column split in three, where the size of each coloured area