search for: murad

...usters are numbered sequentially starting at 1. one of the clustering programs I use (snob) assigns more or less arbitrary integer ids to clusters starting from 3! (clusters 1 and 2 have special meaning in snob). the modified code fixing both problems is included below, changes are commented. best Murad silhouette.default <- function (x, dist, dmatrix, ...) { cll <- match.call() if (!is.null(cl <- x$clustering)) x <- cl n <- length(x) if (!all(x == round(x))) stop("`x' must only have integer codes") k <- length(clid <- sort...

...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. Department of Biochemistry and Molecular Biophysics College of Physicians and Surgeons of Columbia University 630 West 168th Street. New York, NY 10032 Tel: 212-305-6884 Fax: 212-305-6926

...; of columns, on the left side a bar plot of d[,1] is plotted while on the right side a separate bar plot of d[,2] is plotted. how can I combine the two? actually, while on the subject of histograms. is it possible to plot a 3D-histogram in R (a true 3D bar plot, without using image). many thanks Murad -- Murad Nayal M.D. Ph.D. Department of Biochemistry and Molecular Biophysics College of Physicians and Surgeons of Columbia University 630 West 168th Street. New York, NY 10032 Tel: 212-305-6884 Fax: 212-305-6926

...atement (input/output). Anybody has had to deal with a similar situation. any advice for a workaround? both emacs/ess are relatively recent versions (installed a few months ago). I tried using ess-eval-buffer/region instead of cutting and pasting and the same thing happens for me. many thanks -- Murad Nayal M.D. Ph.D. Department of Biochemistry and Molecular Biophysics College of Physicians and Surgeons of Columbia University 630 West 168th Street. New York, NY 10032 Tel: 212-305-6884 Fax: 212-305-6926

..."class", "matrix"), ...) { ... deleted code ... } <environment: namespace:rpart> > debug(predict.rpart,pos="package:rpart") Error: Object "predict.rpart" not found how can I 'debug' non-visible functions, like predict.rpart? many thanks Murad Nayal

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, for a vector 'data' of about > 2000 elements ranging from -1,1 I do: > > hist(data,probability=TRUE) > > I get a h...

...ake the relevant tukeyHSD output look like a dist object and trick hclust into using it. I am somewhat of an amateur in the field (and R) and I am probably making that obvious. any guidance to the 'right' path to approach this (privately or on the list) is really appreciated. many thanks Murad -- Murad Nayal M.D. Ph.D. Department of Biochemistry and Molecular Biophysics College of Physicians and Surgeons of Columbia University 630 West 168th Street. New York, NY 10032 Tel: 212-305-6884 Fax: 212-305-6926

...o NAs, NULLs), except maybe for the fact that the rows sum to 1. > sum(is.na(cavint.p.r)) [1] 0 > sum(is.nan(cavint.p.r)) [1] 0 > I thought kmeans should select initial centers from the data if not given explicitly! any idea what might be going wrong? I am running R 1.7.0 many thanks Murad

Is there such a thing as a CentOS rescue disk? Google gives many links to such a thing, but when I pursue them I do not get to an URL where I can actually download one. -- Timothy Murphy e-mail: gayleard /at/ eircom.net School of Mathematics, Trinity College, Dublin 2, Ireland

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

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 minimization. # len -- Length of time series. # im -- Use only len/im frequencies. # peri -- Periodogram of data. { m <-...

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

...y- 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 values and want to restore the graphics to some sane state after a long R session. many thanks -- Murad Nayal M.D. Ph.D. Department of Biochemistry and Molecular Biophysics College of Physicians and Surgeons of Columbia University 630 West 168th Street. New York, NY 10032 Tel: 212-305-6884 Fax: 212-305-6926

...y and coverage values). here i get mediocre accuracy but respectable coverage of "I". This is not strictly an R question, but I thought someone on the list might have had recent experience with these types of problems and can offer some comments about such an approach. many thanks -- Murad Nayal M.D. Ph.D. Department of Biochemistry and Molecular Biophysics College of Physicians and Surgeons of Columbia University 630 West 168th Street. New York, NY 10032 Tel: 212-305-6884 Fax: 212-305-6926

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 minimization. # len -- Length of time series. # im -- Use only len/im frequencies. # peri -- Periodogram of data. { m <-...

> -----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...

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