similar to: nested anova not giving expected results

Dear all, I am trying to repeat an example from Sokal and Rohlfs "Biometry" -- Box 11.4, example of a randomized-complete-blocks experiment. The data is fairly simple: series genotype weight 1 pp 0.958 1 pb 0.985 1 bb 0.925 2 pp 0.971 2 pb 1.051 2 bb 0.952 3 pp 0.927 3 pb 0.891 3 bb 0.892 4

Dear All, I am studying some process measurement time series in R and trying to identify time delays using cross-correlation function ccf. The results have however been bit confusing. I found a couple of years old message about this issue but unfortunately wasn't able to find it again for a reference. For example, an obvious time shift is observed between the measurements y1 and y2 when the

Hi, I am performing a repeated measures 2-way ANOVA to assess the influence of plant and leaf on aphid fecundity. Fecundity is measured for each aphid on a single leaf. Here is what I typed. wingless <- reshape(Wingless, varying =

Dear list, I have data on insect survival in different cages; these have the following structure: deathtime status id cage S F G L S 1.5 1 1 C1 8 2 1 1 1 1.5 1 2 C1 8 2 1 1 1 11.5 1 3 C1 8 2 1 1 1 11.5 1 4 C1 8 2 1 1 1 There are 81 cages and

Dear list members, In re-running with GLMM() from the lme4 package a generalized-linear mixed model that I had previously fit with glmmPQL() from MASS, I'm getting a warning of a convergence failure, even when I set the method argument of GLMM() to "PQL": > bang.mod.1 <- glmmPQL(contraception ~ as.factor(children) + cage + urban, + random=~as.factor(children) + cage +

Lately R has been behaving strange on my Linux (Ubuntu 7.10) machine, with occasional segfaults. Today something else and reproducible happened: If I type the code below (meant for calibrating data), I get the error message that "the C stack usage is too close to the limit". calcurve <- cbind(1:2e4, 1:2e4, 1:2e3); #dummy curve, real one is more complex caldist <-

Dear all, I wrote the following batch script on a iMac, and ran it on a linux mosix cluster. tu <- read.table("cage.mm5.tags.rna_lib.CAA-CAJ.tu-reshape.table") tu_reshaped <- t(reshape(tu[1:50,], direction="wide", timevar="tu", idvar=c("rna","lib"))) write.table(tu_reshaped, "cage.mm5.tags.rna_lib.CAA-CAJ.tu-reshaped.table")

Dear All, I would like to ask you how to accomplish a little tricky data manipulation. I have a large dataset, looking something like: temp line cage number 18 18 1 6678.63 18 18 1 7774.458 18 18 1 7845.902 18 18 1 9483.578 18 18 1 8983.555 18 18 1 9181.052 18 18 1 9458.696 18 18 1 8138.616 18 18 1 7981.994 18 18 1 7556.491 18 18 1 7672.137 18 18 1 6607.776 18 18 1 8383.65 18 18 1 7129.852 18 18

Hi, I'm trying to figure out how anova works in R by translating the examples in Sokal And Rohlf's (1995 3rd edition) Biometry. I've hit a snag with planned comparisons, their box 9.4 and section 9.6. It's a basic anova design: treatment <- factor(rep(c("control", "glucose", "fructose", "gluc+fruct",

Say we have a tab delimited file called bug.txt Part Rep Cage Hb pupae 1 1 S 32 1 M 34 L 42 2 S 36 M 28 L 36 read.delim("bug.txt") Part Rep Cage Hb.pupae 1 1 1 S 32 2 NA 1 M 34 3 NA NA L 42 4 NA NA NA 5 NA NA NA 6 NA 2 S 36 7 NA NA M 28 8 NA NA L 36 >

Hi all, we have discovered a segfault in nsd-patch when renaming slave zone in nsd config file if some data for this zone still exists in the IXFR diff database. In my case, the zone "black" was renamed to "blackinwhite": > root at ggd115:/cage/nsd/var/nsd/zones#nsd-patch -c > /cage/nsd/etc/nsd-dns-slave.conf > reading database > reading updates to database >

On Nov 11, 2020, at 6:37 PM, Valeri Galtsev <galtsev at kicp.uchicago.edu> wrote: > > how do you map failed software RAID drive to physical port of, say, SAS-attached enclosure. With ZFS, you set a partition label on the whole-drive partition pool member, then mount the pool with something like ?zpool mount -d /dev/disk/by-partlabel?, which then shows the logical disk names in

Hi I am estimating the following coxph function with stratification and frailty?where each person had multiple events. m<-coxph(Surv(dtime1,status1)~gender+cage+uplf+strata(enum)+frailty(id),xmodel) ? > head(xmodel) id enum dtime status gender cage uplf 1 1008666 1 2259.1412037 1 MA 0.000 0 2 1008666 2 36.7495023 1 MA 2259.141 0 3 1008666

Dear All, I would like to ask a couple of questions on a LME model. I tested 4 selection lines at 4 food concentrations against a standard competitor stock. I had 3 replicate cages per selection line. In each cage I have 10 vials. I counted the number of wild type flies and competitor stock emerging in each vial. My main question is: is there any difference between selection lines? I did fit

Dear all, I'm puzzled by the behavior of stringHeight in the grid package. Consider the following test, library(grid) test <- function(lab="dog", ...){ g1 <- textGrob(lab) g2 <- rectGrob(height=grobHeight(g1), width=grobWidth(g1)) gg <- gTree(children=gList(g1,g2), ...) print(c("height:", convertUnit(stringHeight(lab), "mm",

Just installed stellarium-0.10.2-1.el6.rf(i686) last night with the Add/Remove Software tool. I'm getting video corruption and this message in the terminal: StelAppGraphicsScene: drawBackground needs a QGLWidget to be set as viewport on the graphics view A bit of Google work turned up this on the Ubuntu forum: "The new version of Stellarium 0.10.3 resolves my display problems.

Lots of opinions out there and I have read and read, so I have a home server that I have finally setup with a hot swap hdd cage and what I am planning on doing is copying my data drive every other day and rotating them offsite, haven't figured out how often though. So I did my first test last night and the backup drive was formatted with ext3 but looking to try to speed things up and was

I frequently want to test for differences between animal size frequency distributions. The obvious test (I think) to use is the Kolmogorov-Smirnov two sample test (provided in R as the function ks.test in package ctest). The KS test is for continuous variables and this obviously includes length, weight etc. However, limitations in measuring (e.g length to the nearest cm/mm, weight to the nearest

Dear All, I am attempting to use forward and/or backward selection to determine the best model for the variables I have. Unfortunately, because I am dealing with patients and every patient is receiving treatment I need to force the variable for treatment into the model. Is there a way to do this using R? (Additionally, the model is stratified by randomisation period). I know that SAS can be

Hi to all, I usually use R for simple 2D graphs so now I cannot solve the following problem (which I suppose very simple too). I'm using scatterplot3d. I want to add dotted lines parallel to the 3 axes, in order to create a sort of cage in the 3D box (so to represent subdomains in my dataspace). I saw that points3d can be used to draw points or lines to the existing plot but actually