When we analyse ordered categorical data, with categories II,
should we expect an equivalent analysis if the ordering is changed to
II ? Not with a cloglog link (see example below). However, I
suspect the cloglog analysis on the reordered categories is equivalent
to a loglog analysis on
After a one-way anova, se.contrast computes the standard error of a
contrast, but not the value of the contrast itself. Wouldn't this be
useful? Am I missing something?
==
I.White
ICAPB, University of Edinburgh
Ashworth Laboratories, West Mains Road
Edinburgh
The REML loglikelihood includes a term -(1/2)logdet(X'WX) where X is the
design matrix for the fixed effects and W is the inverse covariance matrix
for the observations. Under reparametrisation, X becomes XM with M a
non-singular matrix, and the REML loglikelihood changes by logdet(M).
On Tue, 1
Can anyone explain why with latest version of R (2.0.1) on FC3, installed
from R-2.0.1-0.fdr.2.fc3.i386.rpm, update.packages() produces the message
/usr/lib/R/bin/Rcmd exec: INSTALL: not found.
Indeed /usr/lib/R/bin seems to lack various shell scripts (INSTALL,
REMOVE, etc).
Is it possible in R to create a customized family to be used with glm()? I
see no mention of this possibility in the documentation. (S-plus has
something called make.family).
==
I.White
ICAPB, University of Edinburgh
Ashworth Laboratories, West Mains Road
Is there any way to get a numerical summary of the values of a difftime
object? E.g.
TimeToWean - difftime(WeanDate, BirthDate, units = days)
I can repeat
sum(TimeToWean == 20)
with as many other values are needed to build up a frequency table,
but is there a simpler way in R version 1.8.1?