Hello everyone, Does anyone know how to implement the argument "unit" in R's plclust function ? I used to use Splus where this argument exists but it has not been implemented in R's plclust. The reason why I switched from Splus to R is that Ward's method is not implemented for S's hclust whereas it is implemented for R's hclust. What I would need is S's plclust and R's hclust together, either in R or in S. Many thanks to the ones who will help me solve that trick... Arnaud ************************* Arnaud DOWKIW Department of Primary Industries J. Bjelke-Petersen Research Station KINGAROY, QLD 4610 Australia T : + 61 7 41 600 700 T : + 61 7 41 600 728 (direct) F : + 61 7 41 600 760 ************************** ********************************DISCLAIMER**************************** The information contained in the above e-mail message or messages (which includes any attachments) is confidential and may be legally privileged. It is intended only for the use of the person or entity to which it is addressed. If you are not the addressee any form of disclosure, copying, modification, distribution or any action taken or omitted in reliance on the information is unauthorised. Opinions contained in the message(s) do not necessarily reflect the opinions of the Queensland Government and its authorities. If you received this communication in error, please notify the sender immediately and delete it from your computer system network.
On Tue, 6 May 2003, Dowkiw, Arnaud wrote:> Does anyone know how to implement the argument "unit" in R's plclust > function ? I used to use Splus where this argument exists but it has not > been implemented in R's plclust. The reason why I switched from Splus to > R is that Ward's method is not implemented for S's hclust whereas it is > implemented for R's hclust. What I would need is S's plclust and R's > hclust together, either in R or in S. > Many thanks to the ones who will help me solve that trick...All `unit' does is to rescale the `height' component. Follow this example: data(USArrests) hc <- hclust(dist(USArrests), "ward") plot(hc) hc1 <- hc hc1$height <- rank(hc$height) plot(hc1) -- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595