headers
Christian Hennig | 1 May 2006 01:03
Picon
Picon

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
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>

*** --- ***
Christian Hennig
(Continue reading)

Permalink | Reply | 1 comment |
headers
Gabor Grothendieck | 1 May 2006 01:22
Picon

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
> > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
> >
>
(Continue reading)

Permalink | Reply | 0 comments |
headers
Scott Cunningham | 1 May 2006 03:09
Picon

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
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
(Continue reading)

Permalink | Reply | 7 comments |
headers
Guojun Zhu | 1 May 2006 04:31
Picon
Favicon

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
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Permalink | Reply | 1 comment |
headers
Gabor Grothendieck | 1 May 2006 04:55
Picon

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
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>

______________________________________________
R-help <at> stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
(Continue reading)

Permalink | Reply | 0 comments |
headers
lalitha viswanath | 1 May 2006 04:58
Picon
Favicon

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
(Continue reading)

Permalink | Reply | 4 comments |
headers
Gabor Grothendieck | 1 May 2006 05:15
Picon

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
(Continue reading)

Permalink | Reply | 3 comments |
headers
Gabor Grothendieck | 1 May 2006 05:35
Picon

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
(Continue reading)

Permalink | Reply | 2 comments |
headers
M Senthil Kumar | 1 May 2006 13:55
Picon

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
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Permalink | Reply | 0 comments |
headers
Brian Quinif | 1 May 2006 08:32
Picon

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
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Permalink | Reply | 2 comments |

Gmane