Re: [R] possible bug in flexclust

From: <Friedrich.Leisch_at_r-project.org>
Date: Thu, 05 Jun 2008 13:47:53 +0200

>>>>> On Wed, 4 Jun 2008 13:40:58 -0400, >>>>> Gary Lewis (GL) wrote:

> Hi - Writing to see if someone can suggest whether a problem warrants a bug
> report. It concerns the use of stepFlexclust in the flexclust package. The
> problem concerns the size of clusters returned.

> Versions: R-2.7.0 on Windows XP; RODBC_1.2-3

> code snippet:
> r8 <- stepFlexclust(df,8,nrep=100,FUN=kcca, family=kccaFamily("kmedians"))
> summary(r8) ## returns cluster sizes of 51, 115, 218, 73, 118, 140, 311,
> and 118
> trev <- data.frame(r8_at_cluster)
> sqlSave(ch, trev, rownames="unitid")

> in MySQL:
> select r8cluster, count(*) from trev group by 1 order by 1 ## returns
> cluster sizes of 51, 116, 217, 73, 117, 140, 312, 118

> problem:
> 4 of the 8 clusters differ in sizes. They're close but I'd think they should
> be exact.

Strange. What numbers do you get from

         table(r8_at_cluster)

? Which version of flexclust?

Best,
Fritz

-- 
-----------------------------------------------------------------------
Prof. Dr. Friedrich Leisch 

Institut für Statistik                          Tel: (+49 89) 2180 3165
Ludwig-Maximilians-Universität                  Fax: (+49 89) 2180 5308
Ludwigstraße 33
D-80539 München                     http://www.statistik.lmu.de/~leisch
-----------------------------------------------------------------------
   Journal Computational Statistics --- http://www.springer.com/180 
          Münchner R Kurse --- http://www.statistik.lmu.de/R

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Thu 05 Jun 2008 - 13:33:15 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 05 Jun 2008 - 14:30:37 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive