similar to: Compiling R for OS X 10.2.6 (Darwin 6.6)

It may be my lack of unerstanding, but round() seems to me to give inconsistent results when rounding 5s as in the following examples? > round(1.45, 1) [1] 1.4 # OK > round(2.45, 1) [1] 2.5 # shouldn't this be 2.4? > round(1.05, 1) [1] 1.1 # 1.0 ? and signif(): > signif(2.445, 3) [1] 2.44 # OK > signif(3.445, 3) [1]

Sorry for the off-topic post, but there is such a wealth of knowledge here that I thought someone might have suggestions. Geoff MacLachlan and David Peel have a program - EMMIX - available for mixture modelling. It is only available as FORTRAN source code and is "downloadable" from Geoff's site. The code is reputed to run on Unix machines and is compilable according to

Are there any S-Plus or R libraries/packages that do Finite Mixture Analysis following the algorithms similar to those implemented in Geoffrey MacLachlan's EMMIX program? Thanks. Dr. Marc R. Feldesman email: feldesmanm at pdx.edu email: feldesman at attglobal.net fax: 503-725-3905 "Don't know where I'm going. Don't like where I've been. There may be no exit. But

I just updated rpart to the latest version (3.1.0). There are a number of changes between this and previous versions, and some of the code I've been using with earlier versions (e.g. 3.0.2) no longer work. Here is a simple illustration of a problem I'm having with xpred.rpart. iris.test.rpart<-rpart(iris$Species~., data=iris[,1:4], parms=list(prior=c(0.5,0.25, 0.25))) + ) >

Probably wrong group for this, but a quick question. I've just switched from emacs to Xemacs. In reinstalling ESS 5.1.19 I keep getting the following error when loading Xemacs: "Error in init file: Symbol's function definition is void: w32-using-nt" I've debugged the ess-site.el file, which is where the error originates from. The line causing the difficulty is:

Is there any effective way to get distinct geometric plotting symbols and colors for plots involving more than 25 groups? Thanks. Dr. Marc R. Feldesman Professor and Chairman Emeritus Anthropology Department - Portland State University email: feldesmanm at pdx.edu email: feldesman at attglobal.net fax: 503-725-3905 "Don't knock on my door if you don't know my

I've been using the package rpart with R 1.3.0 for Windows to produce simple classification trees for some measurement data from paleontological specimens. Both the rpart documentation and the output confirm that the program produces splits on continuous data that leave "holes" in the data. It is probably of little practical importance, but is there a reason why the binary

At the risk of starting a holy war, I'd like opinions on the "best" current Linux distribution to run "out of box" with current versions of R. I note rpms for various versions of RedHat, SUSE, Mandrake, and Debian. I've had experience with RH 7.1 & 7.2 and it hasn't been entirely pretty. I understand that 7.3 is a bit "better", but have also

I've created a bunch of postscript trees (post.rpart), which I subsequently edit in Adobe Acrobat 4.05 (I need additional labelling that isn't easy to add with R). After editing the labels and annotations using Acrobat, I then export the file to Adobe's .eps format so I can insert the picture into a Microsoft Word 2000 document. The .pdf file distilled from the .ps file is

I ran across this error the other day while using lm(). I confess that I haven't used it for awhile and haven't been tracking the changes between versions of R. However, the piece of code below is a tiny modification of the example in the help file for "lm". I just separated the commands apart for clarity. A similar piece comes from the help file for

What is the generally accepted way to put certain math symbols into plain text email? For example, how does one convey the "is exactly equal to", "is approximately equal to", and "delta" symbols in a standard text-based (e.g. Courier) font? Dr. Marc R. Feldesman email: feldesmanm at pdx.edu email: feldesman at attglobal.net fax: 503-725-3905

I just downloaded version 1.0.0 and several binary libraries (VR, rpart, norm, stataread) - WinNT version. I then converted a file from Stata 6.0 to R format by using the stataread library. The file converts perfectly and I was able to use the VR function lda on the dataframe without difficulty. I then tried to use the same dataframe with RPART. The model statement:

Forgive what may seem to be a trivial question/problem. Below is some simple R 1.2.1(Windows) code with output. > summary(mammals, digits=10) Name Body.Weight Brain.Weight Red Fox :1 Min. : 3.0 Min. : 26.0 Pig :1 1st Qu.: 35.5 1st Qu.: 138.5 Man :1 Median : 100.0 Median : 406.0 Kangaroo:1 Mean : 761.2 Mean :1000.0 Jaguar :1 3rd

Hello, Does logistic regression really provide better results than lda or qda ? (my purpose is not classification but highlighting of discriminant variables). If this is true, where could I get an R script performing stepwise logistic regression ? Thanks -- Daniel AMORESE Lab. M2C "Morphodynamique Continentale et C?ti?re" UMR CNRS 6143 Caen/Rouen Centre de G?omorphologie UCBN

I'm running a four-group, 15-variable discriminant analysis using the MASS function lda(). Discriminant functions 1 and 2 explain about 93% of the variance. I would like to plot the group scatters in a single scatterplot, identify each group's centroid, and encircle the group scatter with a 90% bivariate confidence ellipse. I can do this for each group on a separate plot, but I

All Is there an R command to produce error bars around a plotted point. Crawley's book uses the command error.bar() but my version of R rejects it. Must be an S+ command(?). Thanks REX [[alternative HTML version deleted]]

At the risk of being flamed for mortal stupidity, I'm trying to figure out what could have possibly changed (at my end almost certainly) that make Brian Ripley's posts (in particular) not show up on my mail server (at least not on the one where I receive my email). All my email comes to a central university server from which it is forwarded to several different places where I can get

In my research area, researchers are very stingy about sharing data. In fact, they are frequently downright secretive about these morsels. Typically, we get fed a diet of summary statistics and the assurance that the data are *normal* without any necessary documentation that they really are. With that as background and assuming that the data really are normal, is there any way in R (or any

Before I reinvent the wheel, I have need for a relatively straightforward crosstabulation (2 x n) function. I know that R has table(), ftable(), xtabs(), and summary(xtabs()), but none of these produce a fully "tricked" out cross-tabulation with marginal totals, expected cell frequencies, and an array of statistics about the contingency table. Is there a more complete (something

Hi Are there any functions that de-convolute data into a given number of clusters, rather like the NPMLE GLIM macros from Murray Aitkin and Brian Francis? Basically I would like to code into R the same approach but include the possiblility of some data being censored. In principle the formulae are the same (just replace the likelihood function) but I haven't managed to get my head round the