[R] update.formula has 512 char buffer?

2012-02-28 Thread Chris Hane
Hello, I am trying to "paste" together a formula to use in the mob function of party. This means the formula will be of the form y ~ x1+ ...+xM | z1+..zN. I am doing some preliminary fits of y ~ x1+ ...+xM, then want to add the conditional part of the equation using update(). Here's the test cod

[R] Party extract BinaryTree from cforest?

2011-10-05 Thread Chris Hane
Hello, I want to plot one of the trees from a cforest object: >data(iris) >cf <- cforest(Species, data=iris) >From the docs, cf@ensemble contains a list of BinaryTrees: ensemble: Object of class "list", each element being an object of class " BinaryTree <../../party/help/BinaryTree%2dclass>". So

[R] Design Survplot performance

2011-07-21 Thread Chris Hane
Hi, I have a Cox PH model that's large for my server, 120K rows, ~300 factors with 3 levels each, so about 1000 columns. The 300 factors all pass a preliminary test of association with the outcome. Solving this with cph from Design takes about 3 hours. I have created the fit with x=T, y=T to save

[R] convert rpart tree to phylo (in ape)

2010-03-05 Thread Chris Hane
Hello, Has anyone written a conversion from rpart tree class to phylo or hclust trees? The plot.phylo method allows much more readable layouts. Thanks! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.eth

[R] cut and re-factor data

2009-09-22 Thread Chris Hane
Hello R-users, I have a data frame with a factor of ages in 5 year increments, and various count data for each age group. I only have this summary information in R at the moment. I want to create a new factor that aggregates the age factors if the existing factors have insufficient counts. Then I