similar to: Survival aareg problem

It's a bug in summary.aareg which no one found until now. What's wrong: If dfbeta=TRUE then there is a second estimate of variance calculated, labeled as test.var2. If maxtime is set, then both estimates of variance need to be recalculated by the summary routine. An incorrect if-then-else flow led it to look for test.var2 when it wasn't relevant. My test cases with maxtime also

Hi all, I've ventured into the world of nonparametric survival and I would like to use the "maxtime" option for printing and plotting my aareg fit. However, my fit does not have "test.var2" and this stops the print and plot when adding a maxtime. My code is as follows: Response<-Surv(Time,Event) Model<-aareg(Response~Factor1*Factor2)

Dear all, I resend this mail because it was blocked: I submitted a bug from the r-bug webpage and hypatia seems to block mail that is send from a different IP than that usually associated with the email. Looks like it is currently impossible to correctly submit bugs from the website. However, here is the original bug report: (PR#8161) Dear all, The following shows cases where accessing

Hi, I'm using windows xp and R 2.8.0 - I wonder what is the command to put in a script that has the same effect as when in a plot you choose menu "History" and "Recording". Best wishes Troels -- Troels Ring - - Department of nephrology - - Aalborg Hospital 9100 Aalborg, Denmark - - +45 99326629 - - tring at gvdnet.dk

On 7/28/2009 6:30 PM, smckinney at bccrc.ca wrote: > Hi all, > > When I install a new package, and don't yet know any function names, > I have to play the "poor man's game" to get to the standard > help system "Index" page for the package: You could complain to the package author or maintainer, who should have created a help alias with the package name

Hi all, When I install a new package, and don't yet know any function names, I have to play the "poor man's game" to get to the standard help system "Index" page for the package: Poor Man's Game -=20 Load new package;=20 issue search() command; find position (say N) of loaded package;=20 issue objects(pos =3D N) command; get name of a random function (san

move or copy the directories mass, nnet and class to the library directory - then execute link.html.help() Now execute library(MASS) and data(petrol) Should work. "Troels Ring"

> -----Original Message----- > From: Duncan Murdoch [mailto:murdoch at stats.uwo.ca] > Sent: Tuesday, August 04, 2009 8:03 AM > To: Steven McKinney > Cc: r-devel at stat.math.ethz.ch; R-bugs at r-project.org > Subject: Re: [Rd] Wishlist: Navigate to "Index" page of help when no > topic specified (PR#13860) >=20 > On 7/28/2009 6:30 PM, smckinney at bccrc.ca

Hi Troels, Its in the incoming directory: ftp://ftp.ci.tuwien.ac.at/incoming Try this for a nice demo: > library(timeslab) > arma(alpha, beta, x, 2, 4, 3, 0.2,500, 30,0) Bernd On 01-Jan-99 Troels Ring wrote: > Hi Johannes Wuebben and a happy new year. > Where did you actually put the timeslab ? > Yours > > Troels Ring, MD > Department of Nephrology

> -----Original Message----- > From: Duncan Murdoch [mailto:murdoch at stats.uwo.ca] > Sent: Tuesday, August 04, 2009 5:06 PM > To: Steven McKinney > Cc: 'r-devel at stat.math.ethz.ch'; 'R-bugs at r-project.org' > Subject: Re: [Rd] Wishlist: Navigate to "Index" page of help when no > topic specified (PR#13860) >=20 > On 04/08/2009 7:33 PM,

I am getting an error that I don't understand from R CMD check on my current instance of the survival code. R2.7.1 on Linux. Here is the last of the log * checking line endings in Makefiles ... OK * checking for portable use of $BLAS_LIBS ... OK * creating survival-Ex.R ... OK * checking examples ... OK * checking tests ... make[1]: Entering directory

I was looking at how the model.frame method for lm works and comparing it to my own for coxph. The big difference is that I try to retain xlevels and predvars information for a new model frame, and lm does not. I use a call to model.frame in predict.coxph, which is why I went that route, but never noted the difference till now (preparing for my course in Nashville). Could someone shed light

On 04/08/2009 7:33 PM, Steven McKinney wrote: >> -----Original Message----- >> From: Duncan Murdoch [mailto:murdoch at stats.uwo.ca] >> Sent: Tuesday, August 04, 2009 8:03 AM >> To: Steven McKinney >> Cc: r-devel at stat.math.ethz.ch; R-bugs at r-project.org >> Subject: Re: [Rd] Wishlist: Navigate to "Index" page of help when no >> topic

Dear Douglas Bates. I just downloaded the compiled version (I'm a poor Windows devil, not yet having found the time to move to a more advanced platform...) from NT- the files are dated 30.5-1999 so they are not old - and the problem persisted....wonder what I did wrong ? R : Copyright 1999, The R Development Core Team Version 0.64.0 Patched (May 3, 1999) R is free software and comes with

Dear friends, I have struggled to do what is likely rather simple but cannot get it working. I have a dataframe with 428 obs. of 5 variables on 26 patients. Two variables are responses, one an index of patients (1:26), one a period-indicator :1 :max 3, and one total time from min 1 to max 36 weeks (max 21 entries). A complete set would be 26*21=546 long but here are only 428 entries. Now,

Dear friends - I'm on windows, R 2.7.0 I try again asking if anyone can explain why a single pig of 16 makes so wild swings. Warnings are issued, and they are 1: pseudoinverse used at 482.1 2: neighborhood radius 242.1 3: reciprocal condition number 0 4: at 360 5: radius 14400 6: all data on boundary of neighborhood. make span bigger 7: There are other near singularities as well. 14400 8:

dear friends - I have an ordered factor, ID, in a data.frame, labs, and make a table, aa <- with(labs,table(val >150,ID)) - and now the order in the table is according to the alphabetic order in ID, and not as I want it. I have tried adding [order(unique(labs$ID)] but it doesn't help. I have R 2.7.1 on WinXP. Best wishes Troels -- Troels Ring - - Department of nephrology - -

Dear friends, I'm examining the characteristics of two models that both fit the sodium concentration in 16 pigs quite well under treatment or control conditions. The more complicated model is by anova better than the less complicated model. To take it further I have generated replicate data using the independent variables and parameter estimates under the two models. A clinically

Dear friends - please give me a hand. I have a dataset of 40 patients in two groups observed on three occasions. I only want to plot a line for each patient in the two groups. I use the ltext function to put the patient number but fail to make lattice understand the numbers as unique since apparently it starts all over with 1:20 for each panel instead of respecting my desires to have 1:20

Dear friends - I'm on winXP, R 2.71 - I have with some help dveloped this multivariate normal model, which gives very plausible results in WinBUGS even without any initial values specified. However, when I then try to run the same model via the bugs function in R2WinBUGS with inits specified as inits=NULL the program stops in a dead end. So I have tried to make inits for the bugs function