similar to: as.dist() too tolerant

The print.dist() method in the mva package currently prints the triangular distance matrices with the ``trick'' of setting the (diagonal and) upper triangular part to NA, and then uses print(mat, na = "") to print that matrix. This is very much undesired if there are true NAs. Reproducible example code: (x <- cbind(c(1,NA,2,3), c(NA,2,NA,1))) (d <- dist(x)) ## does not

Vince Carey found this (thank you!). Since the fix to the problem is not entirely obvious, I post this to R-devel as RFC: help(dist) says: >> Missing values are allowed, and are excluded from all computations >> involving the rows within which they occur. If some columns are >> excluded in calculating a Euclidean, Manhattan or Canberra >> distance, the sum is

I'm diverting this to the more appropriate mailing list, R-help, since heatmap() is standard R function. >>>>> "Johan" == Johan Lindberg <johanl at kiev.biotech.kth.se> >>>>> on Fri, 02 Jan 2004 10:04:14 +0100 writes: Johan> I am trying to plot a matrix of m-values in a heatmap Johan> with "average linkage". The rows

Here is an example: The legends (mainly the factor level names) are cut off on the right. Somehow the internal calculation which computes horizontal space for the legend is not flexible enough. ## Call a new graphics window {with default par()s !}: get(getOption("device"))() par(mfrow = c(2,2)) ## part of example(interaction.plot) _improved_ using with() : data(OrchardSprays)

As many of you have probably realized, the spam filtering at @stat.math.ethz.ch has been dead for since yesterday (09-25) ~16:50 till today ~08:30. The sudden death may have been caused by unrelated installation of some perl modules (spamassassin *is* running on perl) by our IT staff. We are very sorry for this event. On the bright side: You have been able to get a glimpse of what you are

In some contexts, I find the current behavior of rank() very `suboptimal'. We have the argument na.last = {TRUE | FALSE | NA } where the first two cases treating NAs (almost) as if they were == +Inf or == -Inf whereas the 3rd case just drops NAs. For the typical ``Rank Transformation'' that is recommended in EDA in several contexts, I would however want something else, namely keep

>>>>> "daheiser" == daheiser <daheiser@gvn.net> >>>>> on Tue, 6 Apr 2004 04:24:35 +0200 (CEST) writes: daheiser> It is an error in the algorithm. "it" being the behavior reported in bug report PR#1228 --- [too bad you didn't use the whole string "PR#1228" in your subject; if you had, no new report would have

>>>>> "MM" == Martin Maechler <maechler@stat.math.ethz.ch> >>>>> on Thu, 11 Dec 2003 18:29:05 +0100 (CET) writes: >>>>> "Torsten" == Torsten Hothorn <torsten@hothorn.de> >>>>> on Thu, 11 Dec 2003 18:03:07 +0100 (CET) writes: Torsten> On Thu, 11 Dec 2003, Jeffrey Chang wrote:

We have locally collected good experience the last few weeks using procmail-spam filtering, inclduing spamassassin to the procmail filters. The last two spams that went past the majordomo filters (the "traditional ones I've been using) through to R-devel both were not delivered to me, be but caught by the filters. I have now collected a bit more of a week experience using the new scheme

{BCC'ed to three Core groups} Following the recommendation of the mailman developers, I have activated occasional "VERP"ing for our mailing lists. (and will turn it off again, after about a day or so). Here is the mailman "comment-docu" on this : # These variables control the format and frequency of VERP-like delivery for # better bounce detection. VERP is Variable

This is (hopefully) a "once in ten years" event, happening today, from 17--19 CET (= 16--18 UTC = 11--13 US Eastern == 8--10 US Pacific Time) The main `vault' (security, air condition,..) where our mail server hypatia is located will undergo a major `network re-connection' event that will affect a considerable part of ETH -- including all the mailing lists

>>>>> "Jacob" == Jacob van Wyk <jlvw at rau.ac.za> writes: Jacob> Could anybody please help: I simply want to generate Jacob> random samples from a multinomial distribution with Jacob> fixed n and given probability vector p. Jacob> How can I do this? I've recently written the following {for the R package "normix"}. The

During the last 24 hours or so, quite a few e-mails have been delayed for the R-help (and in one case the R-announce) mailing lists. The double posting of "New version of ipred package" to R-announce is my fault, not the posters. I apologize for all inconveniences. The reason is a mail server upgrade and reconfiguration. This should all be past now and the delayed mail hopefully soon

>>>>> "MikeM" == Michael M Meyer <mikem at salter-point.com> >>>>> on Sat, 07 Sep 2002 18:40:01 -0700 writes: MikeM> which(x==max(x)) will return the indices fo all MikeM> elements that are equal to the max. So the first max MikeM> would be which(x==max(x))[1] (At least this works in MikeM> R, I assume it does in

>>>>> "Achim" == Achim Zeileis <zeileis@ci.tuwien.ac.at> writes: Achim> Just a small documentation bug: something went wrong Achim> for the description of the graphical parameter `las' Achim> on the HTML version of the `par' help page: `las' is Achim> missing and instead a line from the description is Achim> used.

I've already fixed the bugs, but as with the last one, this is not critical enough to allow breaking current R-devel's code freeze. I hope I will have corrected it for 1.5.1.. ## Here is code reproducing the problems; ## I use try(.) whenever I know the current versions of R would ## give an error: library(splines) x <- c(1:3,5:6) y <- c(3:1,5:6) (isP <- interpSpline(x,y))#

>>>>> "MikG" == M GRUM <M.GRUM@CGIAR.ORG> writes: MikG> I'm resending this bug report with a new example. As MikG> seen below, cut.dendrogram gives an error message for MikG> some heights, but not for others and with some MikG> datasets adn not others. I can't see why. MikG> Last time I unwittingly sent my message with

Currently, as.data.frame() , i.e., at least its "matrix" method, still coerces logicals to factors. I think this is not desired and not according to the changes to read.table() for R 1.4.0 which allow logicals. Does anybody see a reason against a change which would allow logical matrices to become data.frames with logical variables ? Martin Maechler <maechler@stat.math.ethz.ch>

Executive Summary: This only concerns predict() at new points when the formula contained things like poly() or ns() {in the RHS}. Because I hadn't have time to fix this yet, I submit it. [R version 1.5.1 and at least 1.4.1, probably all versions at all:] For these, R 1.5.x does not differ from 1.4.1 : Inspite of R 1.5.0's "NEWS" entry and the ?Safepredict help page,

Dear R-devel subscribers, I'm feeling quite sheepish at the moment: When I had activated the change from majordomo to mailman this morning, I accidentally didn't activate the spam(+some virus) detection part. I've only activated it half an hour ago. Indeed there were twos message of the following kind sent from R-devel -- I did never seem them personally, because we *do* filtering