Re: [R] survreg 3-way interaction

2011-01-28 Thread tdenes
The fire devoured my laptop, but the PC is still working... :) Sorry for the lazyness, you are totally right. Here goes a reproducible example, which resembles the main features of our dataset: # three experimental factors, 30 subjects expgrid <- expand.grid(F1=0:1,F2=0:1,F3=0:1,id=1:30) # set

[R] survreg 3-way interaction

2011-01-28 Thread tdenes
Dear All, I was wondering why survreg (in survival package) can not handle three-way interactions. I have an AFT (accelerated failure time) model with a frailty term. If the model is restricted to two-way interactions, survreg and gamlss.cens (with random() term) give very similar results, except

Re: [R] unicode&pdf font problem RESOLVED

2011-01-13 Thread tdenes
Dear David, Thank you for your efforts. Inspired by your remarks, I started a new google-search and found this: http://stackoverflow.com/questions/3434349/sweave-not-printing-localized-characters SO HERE COMES THE SOLUTION (it works on both OSs): pdf.options(encoding = "CP1250") pdf() plot(1,ty

Re: [R] unicode&pdf font problem

2011-01-13 Thread tdenes
Hi! Sorry for the missing specs, here they are: > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 12.1 year 2010 month 12 day16 svn rev5385

[R] unicode&pdf font problem

2011-01-12 Thread tdenes
Dear List, I would like to print a plot into pdf. The problem is that the character \U0171 is replaced by a simple 'u' (i.e. without accents) in the pdf file. Example: # this works fine plot(1,type="n") text(1,1,"print \U0171") # this fails pdf("trial.pdf") plot(1,type="n") text(1,1,"print \U01

[R] gamlss & multcomp

2011-01-10 Thread tdenes
Dear Experts, We conducted an experiment with 2 within-subject factors. We have very good reason to use gamlss, which works fine for our dataset, but unfortunately in the final model the vcov matrix can not be produced. This possibility is documented in the gamlss manual, but no hint is given how

[R] multistate multiple failure times

2010-09-08 Thread tdenes
Dear List, I would greatly appreciate if you could give us advice how to analyze multistate multiple failure times in R. We treated event type as an independent variable in the survreg model, but this is clearly not the best solution... If you are more interested, see the details below: First I