1 May 2006 01:03

Re: Number of Clusters

```On Sun, 30 Apr 2006, kumar zaman wrote:

> Dear R users, I am interested in clustering in R. In SAS we have some criteria for determining the number of
clusters using the PROC CLUSTER procedure, which are "CCC" cubic clustering criterion (Sarl 1981),
Psuedo F (PSF), and Psuedo T square (PST). My question is do thsese criterion exists in R, I tried to search
and got one hit (BIC) in Mclust, which I am aware of, any input is greatly appreciated.

Function cluster.stats in package fpc gives you some more.
(See Gordon book reference on the help page.)

Christian

>
>  Ahmed Albatineh, PhD
>  Nova Southeastern University
>  FL, USA
>
>
> ---------------------------------
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help <at> stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
>

*** --- ***
Christian Hennig
```

1 May 2006 01:22

Re: Duration labels on plot axes

```That should have been:

library(chron)
plot(times(0:23/24), 0:23)

On 4/30/06, Gabor Grothendieck <ggrothendieck <at> gmail.com> wrote:
> The times class in chron will give hours and minutes. e.g.
>
> library(times)
> plot(times(0:23/23), 0:23)
>
> and you could modify chron:::axis.times for the others.
>
> On 4/30/06, Duncan Murdoch <murdoch <at> stats.uwo.ca> wrote:
> > I have a variable containing a measurement of a duration in seconds
> > which I would like to use in a plot, with the axes labelled in a format
> > like %H:%M:%S (or possibly %Hh%Mm) if the duration is more than an hour,
> > or just %M:%S if more than a minute, or just decimal seconds if short
> > enough.  Is there a class with an axis.* method defined that has
> > behaviour something like this?
> >
> > Duncan Murdoch
> >
> > ______________________________________________
> > R-help <at> stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> >
>

```

1 May 2006 03:09

R-GUi

```I'm trying to install R-GUI on my Mac OSX, but for some reason it is
not working.  Anyone with a Mac able to answer why I would be getting
the following error?  I am currently running Mac OS X 10.4.6, which
it seems to not be recognizing.

***

R.app GUI 1.15 for R 2.3.0
for Mac OS X 10.4.4 and higher

This package contains universal Cocoa-based R GUI for R 2.3.x for Mac
OS X 10.4.4 (Tiger) or higher.

Requirements:
- Mac OS X 10.4.4 or higher
- R framework 2.3.x installed in its default location /Library/
Frameworks

The application called R will be installed by default in your
Applications folder, but you can change the location during the
installation.

This is a universal binary for both PowerPC and Intel based Macs.

______________________________________________
R-help <at> stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help

```

1 May 2006 04:31

Yet another help needed

```I have a big data.frame with abou 20 column and 60,000
rows for analyze.  Let us say I had a column "a".  I
want to generate a new column which value should be
the average of the 60 a before the current column.
Let us say very row is time t_i.  I need to calculate
a(t_(i-60))+a(t(i-59)+...+a(t(i-1)).  Also I want to
get another number by run regression of
a(t(i-60)),..., a(t(i-1)) on b(t(i-60)),...,b(t(i-1)).
Is there any simple densed code for this?  Thanks.

______________________________________________
R-help <at> stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help

```
1 May 2006 04:55

Re: Yet another help needed

```Look at
?filter
?embed
rollmean in the zoo package
running in the gtools package
runmean in the caTools package

The last one is probably the fastest.

On 4/30/06, Guojun Zhu <shmilylemon <at> yahoo.com> wrote:
> I have a big data.frame with abou 20 column and 60,000
> rows for analyze.  Let us say I had a column "a".  I
> want to generate a new column which value should be
> the average of the 60 a before the current column.
> Let us say very row is time t_i.  I need to calculate
> a(t_(i-60))+a(t(i-59)+...+a(t(i-1)).  Also I want to
> get another number by run regression of
> a(t(i-60)),..., a(t(i-1)) on b(t(i-60)),...,b(t(i-1)).
>  Is there any simple densed code for this?  Thanks.
>
> ______________________________________________
> R-help <at> stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
>

______________________________________________
R-help <at> stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
```

1 May 2006 04:58

table of means/medians across bins used for a histogram

```Hi
I am trying to get a table of means of parameter 1
across BINS of parameter 2.

I am working in proteomics and a sample of my data is
as follows

cluster-age clock-rate(evolutionary rate) scopclass
0.002       10                              A
0.045       0.1                             B
0.13        15                              A
0.15        34                              D
....
....
....
....

Scop class has only 9 distinct categories (A-I)
Whereas cluster-age and clock-rate are discrete
variables greater than 0.

I am trying to do two things with this kind of data,
out of which I managed to accomplish one thanks to the
documentation and pre-existing queries on the mailing
lists.
1. Plot a histogram of the age distribution with scop
class category superimposed on each bin. I managed to
do this with barplot2.
2. Now I am trying to plot a scatter plot of the age
v/s the clock-rate. However to eliminate possible
```

1 May 2006 05:15

Re: table of means/medians across bins used for a histogram

```My understanding is that you want to replace each rate with its average
over the associated bin and then plot age against that.  In that
case try this:

> DF  # test data
age rate bin
1 0.002 10.0   A
2 0.045  0.1   B
3 0.130 15.0   A
4 0.150 34.0   D
> with(DF, plot(ave(rate, bin), age))

Assuming they
are stored in vectors
the columns are age, rate, bin we would have

plot(ave(clock, bin), age)

On 4/30/06, lalitha viswanath <lalithaviswanath <at> yahoo.com> wrote:
> Hi
> I am trying to get a table of means of parameter 1
> across BINS of parameter 2.
>
> I am working in proteomics and a sample of my data is
> as follows
>
> cluster-age clock-rate(evolutionary rate) scopclass
> 0.002       10                              A
> 0.045       0.1                             B
> 0.13        15                              A
```

1 May 2006 05:35

Re: table of means/medians across bins used for a histogram

```Or perhaps a bit simpler:

plot(age ~ ave(clock, bin), DF)

On 4/30/06, Gabor Grothendieck <ggrothendieck <at> gmail.com> wrote:
> My understanding is that you want to replace each rate with its average
> over the associated bin and then plot age against that.  In that
> case try this:
>
> > DF  # test data
>    age rate bin
> 1 0.002 10.0   A
> 2 0.045  0.1   B
> 3 0.130 15.0   A
> 4 0.150 34.0   D
> > with(DF, plot(ave(rate, bin), age))
>
> Assuming they
> are stored in vectors
> the columns are age, rate, bin we would have
>
> plot(ave(clock, bin), age)
>
> On 4/30/06, lalitha viswanath <lalithaviswanath <at> yahoo.com> wrote:
> > Hi
> > I am trying to get a table of means of parameter 1
> > across BINS of parameter 2.
> >
> > I am working in proteomics and a sample of my data is
> > as follows
```

1 May 2006 13:55

Re: R-GUi

```
<Message snippe for brevity>
Hi,
|Requirements:
|- Mac OS X 10.4.4 or higher

This is ok.

|- R framework 2.3.x installed in its default location /Library/
|Frameworks

But is the above one installed in your system ?

HTH,

Senthil

______________________________________________
R-help <at> stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help

```
1 May 2006 08:32

problem installing Econometrics view

```When I try to install hte Econometrics view I get the following error:

CRAN task view Econometrics not available in: install.views("Econometrics")

I have already install the ctv package and loaded it before trying to
install the above...

Any ideas as to what's going on?

______________________________________________
R-help <at> stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help