similar to: What is my replication unit? Lmer for binary longitudinal data with blocks and two treaments.

I am trying to run an ANCOVA with defined error terms. Thus I have to use AOV and not lm. my response variable is proportion of mice paw prints on track plates. These plates were placed on plots that had vegetation and fruit manipulated to two levels each (present or absent), and were sampled monthly for 14 months (repeated measures). The fully crossed factor design was also blocked. My sample

I have a command which returns a data fram if I am not mistaken: sc <- split(x, list(x$Category, x$SubCategory), drop=TRUE) Now I wish to get the Category and SubCategory that the data was split on. So my first attempt would be: sc[[1]]$Category[1] But that yields [1] (Unknown) 46 Levels: (Unknown) 10" Plates 7" Plates (Dessert) ... WOMEN Is all I want is '(Unknown)'.

Dear helpers, I tried these models to run in the package nlme, but allways got the same error message... I have a correlation in 5 sessions within a field (n=12) with ten traps in one field. res2a <- lme(response~x+y+z+treatment),correlation = corARMA(form = ~ session|trapfield, p = 1, q = 0), random=~1|field, na.action=na.omit, data=plates, method="ML") res2a <-

Dear All, I'm rather a beginner on nested ANOVAs, so here goes with my 2 questions; Qu 1: I'm modelling the number of galls on a leaf (the response variable) as a function of; the tree on which I find the leaf, the branch on which I find the leaf. Then, the tree and the branch are both random factors, and I'm quite happy that I should write; aov(galls~tree/branch +

I am trying to fit a formula to my data, but I just can't find the right way to do it. My experiment consists of manipulating FRUITS and VEGETATION to two levels each(intact or removed) on 12 experimental plots. This leaves me with 4 treatment combinations Fruit intact Vegetation removed Fruit int. Veget int. Fruit rem. Veget rem. Fruit rem. Veget. intac those treatements are distributed

I was wondering if anyone could help me with this, simple problem. I am essentially following the example on Hadley's webpage (http://had.co.nz/ggplot2/geom_errorbar.html), but it still doesn't make any sense to me. df <- data.frame(trt = factor(c("intact", "intact", "removed", "removed")), coon = c(0.093, 0.06, 0.057, 0.09), group =

Hi everyone. I am trying to add a formula to my data using the groupedData function. My experiment consists of randomized block design using fruits, vegetation and time as factors. The idea is to see if fruits, vegetation and time explain the abundance of mice. I am using tree density as a covariate. So I tried to fit the following structure to my data. >

I have to provision several dozen * users to a seperate building on our campus in the same subnet. Ordinarily, I'd just run a gigabit cat6 cable to another switch if it doesn't violate the 100 metre rule, but this building is several hundred metres away from my backbone. My only option for cabling to the remote building is copper. My plan is to provision them with a Linux bridge with 4

Hi, I am newbie of R. I a currently using multdrc object to generate fitting curve and IC50. My 384 well format raw data contains multi dose response curves. My script goes through set of data then produce curve and ic50. Here is my sudo code: For (plateid in platelist) { Input data (plateid) as matrix Curve fitting model4logistic <- multdrc(rdata ~ ld, logDose=10) }

Is there an easy way, for MCLUST in particular, to force its graphical output to all black and not use colors to show the clusters/groups? It shows the groups quite nicely by using triangles, crosses, etc. Unfortunately, my printer doesn't convert reds and yellows too nicely and the result is a indistinct figure on the page. (Yes, I have access to a color printer but my intended publication

Hello, I'm having a problem concerning r-binding datasets. I have six datasets, from six different plates, and two different days. I want to combine these datasets together for analysis. Datasets from day 2, have all the same columns than datasets from day 1. However in addition, there are few columns more in day 2. Thus, using rbind for this, results a error, because the objects are not the

dear R experts---I am struggling with the requirements to prepare my files for my printers. I am printing in 2/2 format, which means cyan and black for me, which they take from my color-separated pdf files. R comes into play, because it produces all the figures that are embedded in my book (pdflatex). now, TeX has no problems producing CMYK files. However, R produces RGB files (for

Hello, I'm trunking with an ITSP that, when treating an outbound to an unknown destination, either: - send a SIP error code (I can't be more explicit, at the moment), - or cast a pre-recorded audio message using Early Media. At the same time, I'm also trunking with Contact Center solution which doesn't support Early Media. Beside asking my ITSP to treat calls consistently or

Hi Karl, Thanks - works perfectly now, I think I need to keep the "L" plates on for now :) Thanks again, Jon ----- Original Message ----- From: "Karl Heyes" <karl@xiph.org> To: "icecast" <icecast@xiph.org> Sent: Thursday, May 13, 2004 8:36 PM Subject: Re: [icecast] FreeBSD v5.2.1-RELEASE and Ices 2.0.0... <p>> On Thu, 2004-05-13 at 19:26, Jon

On Wednesday 12 November 2008 00:33, Bill Wendling wrote: > > Please take a look and send feedback toi the list. I'd like to get > > the > > process moving early in the 2.4 cycle. > > Hi Dave, > > Here are my opinions: > > I like the idea of regular validation tagging. However, I think that > it should be as automated as possible. I'm worried that

Hohoho, it''s Xmas time, kids! We have a beta full of goodies for everyone who''s been nice this year. Rails 4.1 is packed to the gills with more marvelous real-world feature extractions, bug fixes, and the tireless polish only a community full of Rails elves could bestow it with. While this is just a beta release, it''s arguably a lot better tested and ready than

On Nov 10, 2008, at 12:59 PM, David Greene wrote: > Back during the LLVM developer's meeting, I talked with some of you > about a > proposal to "validate" llvm. Now that 2.4 is almost out the door, > it seems a > good time to start that discussion. > > I've written up a detailed proposal and attached it to this > message. The goal > is to ease

Hohoho, it''s Xmas time, kids! We have a beta full of goodies for everyone who''s been nice this year. Rails 4.1 is packed to the gills with more marvelous real-world feature extractions, bug fixes, and the tireless polish only a community full of Rails elves could bestow it with. While this is just a beta release, it''s arguably a lot better tested and ready than

Hi, I'm have some data on a grid (specifically high throughput assay plates) and am interested in evaluating measures of spatial autocorrelation to flag plates for corrections. I have been using moran.test and geary.test from the spdep package. My approach is as follows: ## plate is a matrix of data coords <- expand.grid(1:32, 1:48) x <- as.numeric(plate) moran.test(x,

Sender: centos-bounces at centos.org On-Behalf-Of: m.roth at 5-cent.us Subject: Re: [CentOS] TOT - Cemtos 7 : Systemd alternatives ? Message-Id: <a6d79e0d18fdd5ac7d1c714e45d9ed3d.squirrel at host290.hostmonster.com> Recipient: ibudiman at overstock.com -------------- next part -------------- An embedded message was scrubbed... From: <m.roth at 5-cent.us> Subject: Re: [CentOS] TOT -