similar to: anova model refinement/clustering 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

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

Hello all, I am trying to use svm (from the e1071 package) to solve a binary classification problem. The two classes in my particular data set are unequally populated. class 'I' (for important) has about 3000 instances while class "B" (for background) has about 20,000. experimenting with different classifiers I realized that in cases where such an asymmetry exists there is a

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 have sent this message to this list the July, 7th. It was about a problem in MASS package. Until now there is no change in the devel version. As the problem occurs in a package and not in the R-core, I don't know if the message should have been sent here. Anyway, I have added a copy to Pr Ripley. I hope it could have been fixed. Sincerely Marc Le 09/07/2017 ? 16:05, Marc Girondot via

Folks, I'm using stepAIC(MASS) to do some automated, exploratory, model selection for binomial and Poisson glm models in R 1.3. Because I wanted to experiment with the small-sample correction AICc, I dug around in the code for the functions glm.fit stepAIC dropterm.glm addterm.glm extractAIC.glm and came across something I just don't understand. stepAIC() passes dropterm.glm() a

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

Dear R Help, I have noticed some inconsistent behaviour of add1 and drop1 with a weighted linear model, which affects the interpretation of the results. I have these data to fit with a linear model, I want to weight them by the relative size of the geographical areas they represent. _________________________________________________________________________________________ > example

Dear all, I have tried to modify the code of MASS:::stepAIC.default(), dropterm() and addterm() to use AICc instead of AIC for model selection. The code is appended below. Somehow the calculations are still not correct and I would be grateful if anyone could have a look at what might be wrong with this code... Here is a working example: ## require(nlme) model1=lme(distance ~ age + Sex, data =

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

Hi, I encountered a problem when I am trying to write my own function which contains another function. To simplify a problem, I tried the following simplified function, hope someone can idenfity the problem for me. I have a simple data frame called "testdata" as following: >

Hello, I'm currently dealing with a rather strange problem when using the function "stepAIC" ("MASS" package). The setting is the following: From model learning data sets ("learndata"), I want to be able to build prediction functions (in order to save them in a file for further use). This is done by the function "pred.function" (see below). Therein,

This is to announce the new R-package ?ordinal? that implements cumulative link (mixed) models for ordinal (ordered categorical) data (http://www.cran.r-project.org/package=ordinal/). The main features are: - scale (multiplicative) as well as location (additive) effects - nominal effects for a subset of the predictors (denoted partial proportional odds when the link is the logistic) - structured