(PR#1964) The attached function working fine with R 1.3.0 but giving problem with R 1.5.1 (PR#1964)

[Yannis Yatracos]

     Dear Professor Ripley:

I hope you can allow me to interfere between you and Kaushik.
I think I hav misled him.

1) The function cluproj is used in the function pp1.
2) I called "newiris" just the numerical data from "iris".
3) Writing: clus.proj(as.matrix(newiris),rep(1,4))
gave the following:
Error in (seq(n - 1) * mean(XS) - cumsum(XS[seq(n - 1)])) * diff(XS)/((n -  :
        non-conformable arrays
4) A colleague made clusproj work by writing as.vector in front of the numerator
and denominator:
 as.vector((seq(n -1)*mean(XS)-cumsum(XS[seq(n-1)]))*diff(XS))/
as.vector(((n-1)*var(XS)))
5) Here is the result: 
> clusproj(as.matrix(newiris),rep(1,4))
[1] 12.43743 

6) Note that the denominator is a scalar.

I hope things are more clear now.

Best regards,
Yannis Yatracos


-----Original Message-----
From: ripley@stats.ox.ac.uk [mailto:ripley@stats.ox.ac.uk]
Sent: Friday, August 30, 2002 5:54 PM
To: Kaushik Bhattacharyya
Cc: r-devel@stat.math.ethz.ch; R-bugs@biostat.ku.dk; Yannis Yatracos
Subject: RE: (PR#1964) The attached function working fine with R 1.3.0
but giving problem with R 1.5.1 (PR#1964)


1) That is not where you said the error was.

2) iris is a data frame with five columns, one a factor, and that example
should not have worked in 1.3.0 (and I doubt it did).

Please do read the section on BUGS in the FAQ.

On Fri, 30 Aug 2002, Kaushik Bhattacharyya wrote:

> Yes. It passed the "make check". The problem is coming in Tru64 Unix also.
>
> Here is the error message.
>
> For example:
>
> > clusproj(newiris,rep(1,4))
> Error in X %*% a : requires numeric matrix/vector arguments
>
> Here newiris is the iris dataset.
>
> Thanks,
>
> Kaushik
>
> -----Original Message-----
> From: ripley@stats.ox.ac.uk [mailto:ripley@stats.ox.ac.uk]
> Sent: Friday, August 30, 2002 4:41 PM
> To: Kaushik Bhattacharyya
> Cc: r-devel@stat.math.ethz.ch; R-bugs@biostat.ku.dk; Yannis Yatracos
> Subject: RE: (PR#1964) The attached function working fine with R 1.3.0 but giving problem with R 1.5.1 (PR#1964)
>
>
> And the error message is?
>
> Since you mention this is Solaris, did the the installation pass `make check'?  There is a known problem with current gcc compilers on Solaris generating incorrect code.
>
> On Fri, 30 Aug 2002, Kaushik Bhattacharyya wrote:
>
> > The division part of the following code is not executing -
> >
> > (seq(n-1)*mean(XS)-cumsum(XS[seq(n-1)]))*diff(XS)/((n-1)*var(XS))
>
>
> --
> Brian D. Ripley,                  ripley@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 272860 (secr)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595
>
>

-- 
Brian D. Ripley,                  ripley@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 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._