[R] polr (MASS) link functions

2005-07-11 Thread I M S White
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

[R] se.contrast

2004-10-27 Thread I M S White
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

Re: [R] polynomials REML and ML in nlme

2005-02-02 Thread I M S White
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

[R] Failure of update.packages()

2005-02-10 Thread I M S White
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).

[R] own family for glm()

2004-06-11 Thread I M S White
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

[R] summary(difftime(...))

2004-04-15 Thread I M S White
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?