search for: airct

Hello, I've been using ctree and have developed a 55 node - 28 terminal solution. As can be imagined, the plot is difficult to travel down each of the major branches. I've read the help files for ctree I saw where terminal nodes can be color coded. plot(airct, type = "simple") > plot(airct, terminal_panel = node_boxplot(airct, col = "blue", + fill = hsv(2/3, 0.5, 1))) Here is my question: Since my model has 55 nodes and 28 terminal nodes,(ie many branches) is it feasible to color code the each of the major branches and track t...

Hello All, I am using library "party" and I have found a curious/strange behaviour when trying to save the output of a ctree in a file via jpeg/png command. If you use: ################ library(party) airq <- subset(airquality, !is.na(Ozone)) airct <- ctree(Ozone ~ ., data = airq) plot(airct, terminal_panel = node_boxplot, drop_terminal = FALSE) ############### you get a perfect and nice graph, but when you try the equivalent but trying to put the output in a file: ################# library(party) airq <- subset(airquality, !is.na(Ozo...

I have run the ctree function, and my dependent variable is broken into 3 categories: low cost, moderate cost and high cost. When i plot the results (eg. using plot(test.ct)), the plot shows, at the very bottom of each node, the probability of falling into each cost category. When i print the actual results (eg. using print(test.ct)), i get all of the backup

Hello All, I am using library "party" and I have found a curious/strange behaviour when trying to save the output of a ctree in a file via jpeg/png command. If you use: ################ library(party) airq <- subset(airquality, !is.na(Ozone)) airct <- ctree(Ozone ~ ., data = airq) plot(airct, terminal_panel = node_boxplot, drop_terminal = FALSE) ############### you get a perfect and nice graph, but when you try the equivalent but trying to put the output in a file: ################# library(party) airq <- subset(airquality, !is.na(Ozo...

I'm trying to understand how to set graphical parameters for trees created with the party package.  For example take the following code:   library(party)     data(airquality)     airq <- subset(airquality, !is.na(Ozone))     airct <- ctree(Ozone ~ ., data = airq,                    controls = ctree_control(maxsurrogate = 3))     plot(airct)   My problem is, I've got a ctree that has 14 terminal nodes, and as a result of the default graphical paramters of plot(ctree), all of the text is completely undreadable (decison...

Dear R-helpers, I am looking for a segmentation package that gives some "tree identifier" as output for every observation in the data set (my response variable is binary). I have skimmed through "rpart", "ada" and "adabag": The output "trees" gives you the formula, but I have to run several thousand segmentations on different data sets and it

Dear R.Users, using ctree() (from "party" library) on a data.frame, I want to append a column with the references for the groups/segments detected. While these nodes are easy readable in output, I need a vector for my obs. Hints? Cheers -- Daniele Medri

...glmres) # Berserk script # this shows that even for very few cases, my intended way of doing TS C-C is quite time-consuming. # #case<-c(5,0,2,1) #airq<-c(10,2,3,3) # #nctrl<-0*case #for (i in 1:length(case)) #{ # nctrl[i]<-sum(case[-i]) #} # #airca<-rep(airq,case) #airct<-rep(airq,nctrl) # #cases<-rep(1,sum(case)) #ctrls<-rep(0,sum(nctrl)) # #mi<-c(cases,ctrls) #airq<-c(airca,airct) #a<-Sys.time() #air.cl<-clogit(mi~airq) #b<-Sys.time() #b-a #a<-Sys.time() #air.cl<-clogit(mi~airq, method="approximate") #b<-Sys...

...m(qres) fivenum(glmres) # Berserk script # this shows that even for very few cases, my intended way of doing TS C-C is quite time-consuming. # #case<-c(5,0,2,1) #airq<-c(10,2,3,3) # #nctrl<-0*case #for (i in 1:length(case)) #{ #? nctrl[i]<-sum(case[-i]) #} # #airca<-rep(airq,case) #airct<-rep(airq,nctrl) # #cases<-rep(1,sum(case)) #ctrls<-rep(0,sum(nctrl)) # #mi<-c(cases,ctrls) #airq<-c(airca,airct) #a<-Sys.time() #air.cl<-clogit(mi~airq) #b<-Sys.time() #b-a #a<-Sys.time() #air.cl<-clogit(mi~airq, method="approximate") #b<-Sys.time() #b-a