similar to: Corrupt trees

Dear All When drawing a classification tree with plot(mytree) text(mytree) the conditions are written just before the nodes branch. My question is: can one be certain that those conditions refer to the left-side branches? (The R documentation surprisingly lacks the information that I am asking for.) Thanks in advance, Paul

I have a number of x, y observations (Time, Conc) for a number of Subjects (with subject number Subj) and Doses. I can plot the individual points with xyplot fine: xyplot(Conc ~ Time | Subj, Groups=Dose, data=myData, panel = function(x,y) { panel.xyplot(x, y) panel.superpose(???) # Needs more here } ) I also like to plot on

Dear all, I don't understand the following behaviour: Running compar.gee (in library ape ) with and without the option 'data', it give me different results Example: .... Start R .... > load("eiber.RData") > ls() [1] "gee.na" "mydata" "mytree" > library(ape) > # runnig with the option data= mydata > compar.gee(alt ~ R,

Hi all, I used rpart to fit a model, where the covariates in the model are categorical variables. Then I plotted the tree (mytree) and used the command "text" to add labels to the tree. In the nodes of the tree, the values of the covariates are represented with a, b or c. Is there a way to show the real value(s) of the variable in the nodes instead of a, b or c ? I found that the

Hello all, I have a problem with the way attributes are dealt with in the function xmlToList(), and I haven't been able to figure it out for days now. Say I have a document (produced by nmap) like this: > mydoc <- '<host starttime="1365204834" endtime="1365205860"><status state="up" reason="echo-reply"

Dear users, I'll appreciate your help with this (hopefully) simple problem. I have a model object which was fitted to inputs X1, X2, X3. Now, I'd like to use this object to make predictions on a new data set where only X1 and X2 are available (just use the estimated coefficients for these variables in making predictions and ignoring the coefficient on X3). Here's my attempt but, of

Hello All, Am I misreading the documentation? The text.rpart documentation says: "label a column name of x$frame; values of this will label the nodes. For the "class" method, label="yval" results in the factor levels being used, "yprob" results in the probability of the winning factor level being used, and ?specific yval level? results in the probability of

Hello I am trying to store output from a loop into an empty matrix. The current code I am using is: > M<-mat.or.vec(11,89) > for (j in list(3,91)) + {M[,(j-2)]<-pic(datain[,j], mytree)} datain is a matrix (11,91). I only want to use the pic() function on the columns 3:91. When I use this code output is stored into the 1st column of M and the 89th column of M but no other columns. Any

I am having trouble displaying a dendrogram of evolutionary relationships (a phylogram imported from the ape package) as the vertical component of a heatmap, but keeping the hierarchical clustering of the horizontal component. The relationships of the vertical component in the generated heatmap are not that of the dendrogram, although the ordering is. In more detail, I am attempting to generate

Hi, I am using rpart decision trees to analyze customer churn. I am finding that the decision trees created are not effective because they are not able to recognize factors that influence churn. I have created an example situation below. What do I need to do to for rpart to build a tree with the variable experience? My guess is that this would happen if rpart used the loss matrix while creating

Hi Everyone, Is there a way to get predict.rpart() to return the nodes reached by the new examples in addition to the predicted probabilities it already returns? In other words, I would like to know the leaf node in the tree object that each new example data drops down to. Thanks in advance for your help. Osei

Is there way on import or add to exclude directories? (Obvious candidate being log or tmp dirs.) I''m imagining something like this: $ svn import mytree file:///usr/local/svn/newrepos/some/project --exclude /tmp --exclude /log \ I have scoured "Version Control with Subversion" without luck but maybe I don''t understand some terminology somewhere... -- Posted via

Dear R Help, I read through the archives pretty extensively before sending this email, as it seemed there were several threads on using predict with GLM. However, while my issue is similar to previous posts (cannot get it to predict using new data), none of the suggested fixes are working. The important bits of my code: set.seed(644) n0=200 #number of observations

Hello, I have estimated a simultaneous equation model (similar to Klein's model) in R using the system.fit package. I have an identity equation, along with three other equations. Do you know how to explicitly identify the identity equation in R? I am also trying to forecast the dependent variables in the simultaneous equation model, while incorporating the identity equation in the

Dear all, I'm using the "ape" package in R and want to draw a phylogenetic tree with not only the tip labels but also some labels for the edges. e.g. Mark the edge AB as "m" in the tree ABC. Couldn't find a way to do that. Can someone help? Thanks, Hua

Hi, I created a logistic regression model using the glmnet package. This model is of class "glmnet" or "lognet". I wanted to plot a calibration curve for this model using the calibrate() function from rms package, but the objects used are different, rms requires a fit from lrm(). Is there another function for getting the calibration plot for this glmnet object, or can anyone

Hi, Suppose i have generated an object using the following : fit <- rpart(Kyphosis ~ Age + Number + Start, data=kyphosis) And when i print fit, i get the following : n= 81 node), split, n, loss, yval, (yprob) * denotes terminal node 1) root 81 17 absent (0.7901235 0.2098765) 2) Start>=8.5 62 6 absent (0.9032258 0.0967742) 4) Start>=14.5 29 0 absent (1.0000000

Hello, Please let me know if this is a bug or a feature. I used to be able to do the following with 2.5.5. On a remote host I have /mytree and I mirrored at my local host as a regular user doing rsync -avz --delete kent@host:/mytree /. with 2.5.6 I get mkstemp .file.aFhK3N failed: Permission denied probably because rsync2.5.6 is attempting to create a temporary file on the root directory,

Dear all, I would like to model the relationship between y and x. y is binary variable, and x is a count variable which may be possion-distribution. I think it is better to divide x into intervals and change it to a factor before calling glm(y~x,data=dat,family=binomail). I try to use rpart. As y is binary, I use "class" method and get the following result. >

Hi Frances, I have not touched the system.fit package for quite some time, but to solve your problem the following two pointers might be helpful: 1) Recast your model in the revised form, i.e., include your identity directly into your reaction functions, if possible. 2) For solving your model, you can employ the Gau?-Seidel method (see https://en.wikipedia.org/wiki/Gauss%E2%80%93Seidel_method).