similar to: Additional Documentation for rpart?

Dear r-help, I am trying to build a new function (to process rpart objects) that will output matrix that has a row for each node and a column for each feature. With each entry in the table is a numerical property at that node for that feature (e.g. surrogate split agreement, improvement). My current trouble is that the only clue to the identity of the feature is stored as the *name* of the

Dear R, Short verse of the question: Is there R code which will calculate surrogate splits and/or delta impurity for decision trees at each node? Long Version: I have local, legacy code which I use to calculate my decision trees. I would like to switch to R, but as I understand it surrogate splits are not implemented. Surrogate splits and feature ranking are described in Breiman et al

Dear r-help, Is there any available description of the components of lqs objects found in the package "lqs"? > names(slts) [1] "crit" "sing" "coefficients" "bestone" [5] "fitted.values" "residuals" "scale" "terms" [9] "call"

Dear r-devel, I am mutating rpart to do calculations on trees. I am trying to extract information from the tree. However, I got a seg. fault. This is the offending line in "rpmatrix.c": deltaI[0][0] = spl->improve; (Commenting it out cures the seg fault) I would like some advice on how to debug this. I have allocated memory with calloc and deltaI[0][0] should be

Hi All, When I tried to do logistic regression and I got the following messages: > SampledW.glm.ALL <- glm (V1 ~ ., family = binomial, data = SampledW) Warning messages: 1: Algorithm did not converge in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, 2: fitted probabilities numerically 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y =

Do those who have had trouble finding stuff in the R-archive used the search engine? http://www.scirus.com/ And as far as threading goes. My mailler "mutt" threads the mailing list quite nicely. --Clayton -- Clayton Springer, Ph. D. Sandia National Laboratories csprin at ca.sandia.gov Biosystems Research Department (925) 294-2143 P.O. Box

Dear R users, I want to do discriminant analysis on my data. I have successfully followed the discriminant analysis in V & R on the iris data: > ir <- rbind (iris3[,,1],iris3[,,2],iris3[,,3]) > ir.species <- c(rep("s",50),rep("c",50),rep("v",50)) > a <- lda(log(ir),ir.species) > a$svd^2/sum(a$svd^2) [1] 0.996498601 0.003501399 > a.x <-

Dear R users, I am a novice user working through V & R. My first question stems from p20 of the second edition.: I perform the following actions: > mydata <- c(2.9,3.4,3.4,3.7,3.7,2.8,2.8,2.5,2.4,2.4) > names(mydata) <- c('c','j','b','e','i','h','g','d','f','a') > mydata c j b e i h

Dear All, I was trying to follow some of the examples in Venables and Ripley "Modern applied ... with S-plus" I have downloaded a copy of the iris data set and loaded into R. : however I cannot use the apply command (from p47): > apply (iris, 2 ,mean) Error in sum(..., na.rm = na.rm) : invalid "mode" of argument > apply (iris, c(2) ,mean) Error in sum(..., na.rm =

I am doing some testing in Samba 3.0. I am using security=ADS. I am wondering if Samba has any support for cross-realm trust. For example, I have one AD forest SANDIA.GOV that has trust with SANDIA2.GOV. I have the Samba server on linux.sandia2.gov. I have a local user account on linux.sandia2.gov called user. When I log into a win2k client as SANDIA2.GOV/user and connect, it works fine.

I need to calculate some metrics such as Mean Time Between Failure (MTBF), etc (see http://www.cs.sandia.gov/~jrstear/ras for a more complete list). I have observations like start end state 1 2005-11-11 09:05:00 2005-11-11 12:20:00 Scheduled Downtime 2 2005-11-12 13:42:00 2005-11-12 14:45:00 Unscheduled Downtime where each row describes

Prof. Harrell, My name is Armida Carbajal, I'm a graduate student intern at Sandia National Laboratories (SNL) and am conducting some research for my thesis project at the University of New Mexico in Statistics for SNL. My project entails a logistic regression and I wanted to create a classification table like the one found in SAS using the function CTABLE. I was running out of ideas on

Hi.. I am using R as a user, another group built everything I am using. But I need to install some packages and my only internet access is via a web-proxy. I have hunted though the docs and I'm not finding anything on setting a proxy in the package install command. Any help would be really appreciated. Thanks Sue Susan Caskey International Biological Threat Reduction Sandia National

This question concerns rpart's facility for user-defined functions that accomplish splitting. I was interested in modifying the code so that in each terminal node, a linear regression is fit to the data. It seems that from the allowable inputs in the user-defined functions, that this may not be possible, since they have the form: function(y, wt, parms) (in the case of the

I have created plots of rpart objects with the fancy option for text.rpart ("fancy" creates ellipses and rectangles and labels branches with splitting criteria). The ellipses and rectangles are supposed to "interrupt" the tree lines (as seen in Therneau and Atkinson 1997, page 48, Fig. 18, http://www.mayo.edu/hsr/techrpt/61.pdf), but they do not, even when I use Therneau and

Hi! I?m trying to calculate de C index (concordance probability) through the somers2 function (library Hmisc). I?m interesting on including the sampling effort as a weight factor for the evaluation of model predictions with real data. I?ve some questions about that: first of all I?m not really sure if I can include sampling effort as a weight factor. Since the weight factor should be a numeric

Hi all, How can the label of the x-axes in the plot() of a stl.object be adapted? e.g., When plotting: plot(stl(nottem, "per")) In the labels of the x-axes is “time”. How can this be changed to e.g., “Time (dekade) “? It does not work with xlab or others anymore… Thanks, Jan _______________________________________________________________________ Ir. Jan Verbesselt Research

Dear R list, From the lrm() binary logistic model we derived the G2 value or the likelihood-ratio chi-squared statistic given as L.R. model, in the output of the lrm(). How can this value be penalized for non-linearity (we used splines in the lrm function)? lrm.iRVI <- lrm(arson ~ rcs(iRVI,5), penalty=list(simple=10,nonlinear=100,nonlinear.interaction=4)) This didn’t work

Hello, We are trying to figure out how to preserve the low-frequency content when we try to convert a .wav file to an .opus file. What settings in Opus that can accomplish this? Sincerely, Vinson Go cid:image001.png at 01C94EE1.4DA09EF0 TBSI - SIMPLICITY in Neuro-Solutions <http://www.tbsi.biz/> Tethered and Wireless Solutions for Neural Recording and Stimulation Vinson

Does SAMBA 3.0 use the PAC information available within a Microsoft Kerberos ticket? Thanks. -dan -------------------------------------- Daniel Wachdorf drwachd@sandia.gov Sandia National Laboratories System Security Research and Integration 505-284-8060