similar to: ctree

I try tu use mob() with my data.frame ('data.frame': 288 obs. of 81 variables; factors, numerics and ordered factors) My response is a binary variable and I should use for modelling a logistic regression (family=binomial). I read in the "MOB" Vignette that I could use a formula like this if I would like to have only partitioning variables apart from the response.

Hi, I am trying to get the terminal nodes of a plot of a ctree object to look nice. Using the iris data I have: library(party) mtree <- ctree(Species ~ ., data=iris) plot(mtree,terminal_panel=node_barplot(mtree)) The terminal nodes don't display the species names because the names are displayed horizontally. ?I would like to reduce the size of the labels and make the terminal nodes

Hi, from what I understand about the canonical correlation function 'cancor', it looks for correlations in two sets of variables, each represented in matrix form. Right? Sounds exactly like what I need. I have tried the following but I am not sure how to interpret the results. AudioPCs <- c(ArTHarF0PCA$x[,2], ArTHarF1PCA$x[,2], ArTHarF2PCA$x[,2], ArTHarF3PCA$x[,2],

Dear all, I didn't find any graphical solution in the package "mva" to plot the canonical scores from a CCA (canonical correlation analysis). Does anybody knows how to plot or has anybody already programmed : - the map of the canonical scores, - the graph of the canonical weights, - the correlation circle i.e. the canonical loadings ? Thank you for help ...

Hi, Given that cdplot() is used to produce the conditional density of a categorical y along a numerical x, it seems natural that it could be used with a date or time x (such as 'POSIXct'). Is this desirable? If so, I've created a patch that would allow this, by coercing the POSIXct x variable to produce the density, but use the original POSIXct x to draw the x axis. Index:

Hello, I wanted to create a histogram, but somehow I got stuck... The interval limits are: x = 1, 2, 3, 3.5, 4.5, 5, 5.5 The interval widths are therefore: 1, 1, 0.5, 1, 0.5, 0.5 Nothing I tried worked... Can anyone help me please? Thanks Tobias -- View this message in context: http://www.nabble.com/Histogram-tf4315900.html#a12288850 Sent from the R help mailing list archive at Nabble.com.

Hello: I am using rpart and would like more background on how the splits are made and how to interpret results - also how to properly use text(.rpart). I have looked through Venables and Ripley and through the rpart help and still have some questions. If there is a source (say, Breiman et al) on decision trees that would clear this all up, please let me know. The questions below pertain to a

Dear all, I am quite new with R and I have a problem with plotting a binomial regression line in a plot. This is what I type in: > model<-glm(Para~Size,binomial) > par(mfrow=c(1,1)) > xv<-seq(3.2,4.5,0.01) > yv<-predict(model,list(area=xv),type="response") > plot(Size,Para) > lines(xv,yv) The error message that I get is: > Error in xy.coords(x, y) :

Hi List, I am unsuccessfully trying to beautify barplot outputs from ctree. For example I would like to rotate x-axis lables and resize/change font/type. mtree <- ctree(ME ~ ., data = mammoexp) plot(mtree,terminal_panel=node_barplot(mtree,col="black",fill=NULL, beside=TRUE, ylines=NULL, widths=1,gap=NULL, reverse=FALSE,id=FALSE))

Hi Mark, I don't know wether you recived a sufficient reply or not, so here are my comments to your problem. Supressing the constant term in a regression model will probably lead to a violation of the classical assumptions for this model. From the OLS normal equations (in matrix notation) (1) (X'X)b=X'y and the definition of the OLS residuals (2) e = y-Xb you get - by

Dear r-help list, I'm trying to use multiple imputation for my MSc thesis. Having good exemples using the Hmisc package, I tried the aregImpute function. But with my own dataset, I have the following error : Erreur dans matxv(X, xcof) : columns in a (51) must be <= length of b (50) De plus : Warning message: In f$xcoef[, 1] * f$xcenter :   la taille d'un objet plus long n'est pas

Hi Prof. Williams, thanks for your suggestion. The updated code is below. It turns out it was a matter of displaying the second column in yval to get the number of N and subtracting it from the n column in the frame to get the number of Y remaining in a binary example. once this is added now the function returns the rules along with Y and N count affected by the resulting rule. I am ccing

I've got the graphic with this: a=c(120,40,75,85,55,75,55,90,90,55,155) plot(ecdf(a)) How can I draw lines to show the lower quatile and upper quartil? I only know there is the argument lines(), but I don't know how to use it -- View this message in context: http://r.789695.n4.nabble.com/How-can-I-draw-a-line-in-empirical-distribution-function-tp3090675p3090675.html Sent from the R

Rers: I can't seem to locate the command to suppress what appears to be a faint internal grid when running the following command to make a filled contour plot of some data I have (x,y,z being the inputs):

Hi all, I am using R1.5.0 under Unix, I have a couple of questions here. 1. My program is running out of memory. I am writing a program to grow a list of trees using rpart() on a subset of a large dataset(5807x693) with a different response for every tree. I saw that after each tree was constucted, 116 MB of data was being added to the Vcells. I have no idea what this data is. My dataset is

> ---------- Forwarded message ---------- > Date: Wed, 5 Oct 2011 21:09:41 +0000 > From: Chris <christopher.a.hane at gmail.com> > To: r-help at stat.math.ethz.ch > Subject: Re: [R] Party extract BinaryTree from cforest? > > I found an internal workaround to this to support printing and plot type > simple, > > tt<-party:::prettytree(cf at ensemble[[1]],

We are trying to make a decision tree using rpart and we are continually running into the following error: > fit_rpart=rpart(ENROLL_YN~MINORITY,method="class") > summary(fit_rpart) Call: rpart(formula = ENROLL_YN ~ MINORITY, method = "class") n= 5725 CP nsplit rel error 1 0 0 1 Error in yval[, 1] : incorrect number of dimensions ENROLL_YN is a

Dear colleagues, following John Fox' advice in this article (http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-cox-regression.pdf), I'm trying to create a new data frame to examine the differential survival curves from a combination of covariates. These are derived from a Cox Proportional Hazards model I fit to data about the diffusion of a particular policy across American

How do I calculate the confidence interval for the value x given by the intersection of two quadratics (i.e. parabolas)? I fit two quadratics of the form: y = C1 + B1*x + A1*x^2 y = C2 + B2*x + A2*x^2 to two sets of points N1 and N2. I test for whether they intersect, if they do then I calculate the roots of: 0 = (C1 - C2) + (B1 - B2)*x + (A1 - A2)*x^2 to determine where they intersect

R-1.7.0 built on NetBSD 1.6, but the validation test suite failed: Machinetype: Intel Pentium III (600 MHz); NetBSD 1.6 (GENERIC) Remote gcc version: gcc (GCC) 3.2.2 Remote g++ version: g++ (GCC) 3.2.2 Configure environment: CC=gcc CXX=g++ LDFLAGS=-Wl,-rpath,/usr/local/lib make[5]: Entering directory `/local/build/R-1.7.0/src/library' >>> Building/Updating