similar to: Forced inclusion of varaibles in validate command as well as step

It looks from str(SA) that Response IPS1 is a data.frame of class "anova", which probably cannot be coerced to vector. Maybe you can use unlist() instead of as.vector() Or something like SA[["Response IPS1"]]["as.factor(WSD)",] ## to select the first row only, even maybe with unlist() Without a better REPRODUCIBLE example, I cannot tell more (maybe some others

Dear prof. Harrell, I'm not able to use the force option with fastbw, here an example of the error I've got (dataset stagec rpart package): > fitstc <- cph(Surv(stagec$pgtime,stagec$pgstat) ~ age + eet + g2 + grade + gleason + ploidy, data=stagec) > fbwstc <- fastbw(fitstc,rule="aic",type="individual") > fbwstc Deleted Chi-Sq d.f. P Residual d.f.

Greetings guRus -- If a variable, e.g., 'varname', is a character string, e.g. varname <- "datavector", and I want to apply a function, such as table(), to datavector, what syntax or method will do so using only the variable varname? This seems similar to indirect addressing, but I have not seen a method for it in the R manuals. Is there a general name for such indirect

Hello, I have simulated 30 observations from a binomial(5,0.1) distribution. Now I need to make frequency table( that means I need to tally how many 0's , 1's 2's....... 5's) I know that the simple R function table() will do this, but I am afraid that some times I may get zero frequency for some particular values (for example in the above there are 5-0's 10-1's ,

I would like to produce a bar plot with varying-width bars. Here is an example to illustrate: ww <- c(417,153,0.0216,0.0065,556,256,0.0162,0.0117, + 726,379,0.0358,0.0501,786,502,0.0496,0.0837, + 892,591,0.0785,0.0795) yy<-t(t(array(ww,c(2,10)))) barplot(yy[,2*1:5],las=1,space=c(.1,.5),beside=T) produces a barplot of 5 pairs of bars that are of equal width