similar to: SVM question

Hello all, since we're on the topic of R-editors. I am using emacs/ess on a unix workstation (to interact with R and have been having a little problem. I usually write the R commands I need to run in a separate buffer then copy and paste them into the *R* buffer for evaluation. The problem is, if any command is spread over multiple lines emacs/R hangs when I paste it in the R buffer for

Hello, I was wondering whether a Poisson mixture modeler/cluster analysis package is available for R. I scanned CRAN packages and couldn't find anything but I thought I'd ask. If not could anyone recommend a non-R open source package. I have found 'snob' but this program seems a bit hard to use in an automated, non interactive fashion. regards, Murad -- Murad Nayal M.D. Ph.D.

Hello all, This might have been fixed in later versions (I am using R1.7.0), r-help archive contains messages reporting similar problems but no reports of codes fixes. I have encountered a couple of problems using the silhouette function. one occurs when the clustering contains clusters composed of 1 element (Martin Maechler posted code few months ago that fixes a similar problem that occurs

Hello, I am trying to compare two histograms using barplot. the idea is to plot the histograms as pairs of columns side by side for each x value. I was able to do it using barplot before but I can't remember now for the life of me now how I did it in the past: > d [,1] [,2] -37.5 0.0000000000 2.789396e-05 -32.5 0.0001394700 5.578801e-05 -27.5 0.0019804742

Hi, I am trying to refine models of a continuous response variable and a number of categorical predictor variables. I know of some model refinement tools available in R that help in the selection of model terms like dropterm and addterm from MASS etc. However, I would also like to try to refine the model by 'coalescing' some levels of some of the predictor factors. Is there a standard

Hello, Is there a specific command to clear the graphics window. On occasion I need to construct plots using commands that don't clear the graphics window (like text, lines and points etc.) -only- and hence need to clear the graphics completely before hand. also, is there a way to restore the graphics parameters to default values, say in these cases where you forgot to save the original

Hi, I would like to step-through a non-visible function. but apparently I don't know enough about namespaces to get that to work: > methods(predict) ... deleted lines ... [27] predict.rpart* predict.smooth.spline* [31] predict.survreg.penal* Non-visible functions are asterisked > debug(predict.rpart) Error: Object "predict.rpart" not found >

> -----Original Message----- > From: Murad Nayal [mailto:mn216 at columbia.edu] > Sent: Tuesday, 24 June 2003 12:02 PM > Cc: r-help at stat.math.ethz.ch > Subject: [R] help on R programming. > > > > > Hello all, > > I am looking for books to help me gain a firmer grasp on the S/R > programming language , programing / data structures etc. it > seems

Hello! Some CS-guys (the type who knows what Church formalism is) keep asking me questions about formal definitions of R-language that I can not answer (or even understand). Is there some freely available papers which I can throw at them where it would be explained is R functional/OOP/procedural language, does it use weak/strong, dynamic/static typization, does it use lazy or ...(do not know

just realized that the bin value is actually the relative frequency divided by the bin width. sorry for consuming band width. Alas, is there anyway to make hist() calculate relative frequencies irrespective of bin width? thanks Murad Nayal wrote: > > Hello, > > I am rather new to R. in trying to use the hist() command I get behavior > that is somewhat puzzling me, in short,

Hi, I want to use (similar to) las options in lattice (bwplot) plot. Actually I want to have x-axis labels as vertical instead of default horizontal. Thanks in advance for your help. Mahbub

In Splus the code test.lm <- lm(y ~ x, data = test.data) plot(test.lm) generates a graphics window that contains multiple graph sheets that one may choose from via the "page" tabs at the bottom of the window. Is there a way to do this sort of thing in R? As another example, I have some repeated measures data with continuous outcomes and have been working with the nlme library in

Hello, I have been getting the following intermittent error from kmeans: >str(cavint.p.r) num [1:1967, 1:13] 0.691 0.123 0.388 0.268 0.485 ... - attr(*, "dimnames")=List of 2 ..$ : chr [1:1967] "6" "49" "87" "102" ... ..$ : chr [1:13] "HYD" "NEG" "POS" "OXY" ... > set.seed(34) >

Hello, I have an SVM model previously calibrated using libsvm R implementation from the e1071 package. I would like to use this SVM to predict values, from a Java program. I first tried to use jlibsvm and the "standard" java implementation of libsvm, without success. Thus, I am now considering writing data in files from my Java code, calling an R program to predict values, then gather

I am accustomed to the LIBSVM package, which provides cross-validation on training with the -v option % svm-train -v 5 ... This does 5 fold cross validation while building the model and avoids over-fitting. But I don't see how to accomplish that in the e1071 package. (I learned that svm(... cross=5 ...) only _tests_ using cross-validation -- it doesn't affect the training.) Can

Hello all, I'm working with the svm (libsvm) implementation from library(e1071). Currently I'm trying to extend recursive feature elimination (R-SMV) to work with multiclass classification. My problem is that if I run svm for a 3 class problem I get a 2-D vector back from model$coefs, can someone explain me what this values are? I understand them in the 2-class problem where this is a

Dear R-users, I have sent another mail some hour ago about a matlab Code I was trying to translate in R. Actually I have found a simpler code originally written in S-PLUS for the same function. Author's page -> http://math.bu.edu/people/murad/methods/locwhitt/ ============================================================= rfunc_function(h, len, im, peri) # h -- Starting H value for

I am starting to use svm from e1071 and I wonder how exactly crossvalidation is implemented. Whenever I run > svm.model <- svm(y ~ ., data = trainset, cross = 3) on my data I get dirrerent values for svm.model$MSE e.g. [1] 0.9517001 1.7069627 0.6108726 [1] 0.3634670 0.9165497 1.4606322 This suggests to me that data are scrambled each time - the last time I looked at libsvm python

Hi I'm using Libsvm, I made a function to save models, confusion matrix of different models but when I want to use saved model, I can't let suppose output of function is : models=list(model1,model2,....) but how can I use these models? using predict(models[[1]],y) causes error message: "Error in UseMethod("predict") : no applicable method for 'predict'

Dear David, Dear Friends, After any running svm I receive different results of Error estimation of 'svm' using 10-fold cross validation. What is the reason ? It is caused by the algorithm, libsvm , e1071 or something els? Which value can be optimal one ? How much run can reach to the optimality.And finally, what is difference between Error estimation of svm using 10-fold cross validation