[R] Clustering for social/complex networks

2006-02-03 Thread Kel Lam
Hi group,

Here is another clustering question.  Is there
anything available in R that I can cluster groups with
overlaps?  It looks like a Venn Diagram with two or
more circles overlapping one another.  

I realize that Hierarchical Clustering (hclust in
Stats) only group individuals into distinct groups. 
My scenario is I have to group Physicians that each of
them may belong to more than one group.  That reflects
the complexity of relationhip, something analogous to
Social Network.  

Thanks!

kel

__
R-help@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] Clustering Question

2006-01-26 Thread Kel Lam
Hi group,

My case has N physicians with each seeing M patients. 
One physician could have seen a group of patients, or,
a patient could have been seen by multiple number of
physicians.  In order words, there are overlaps.  Now,
I have the following NxM matrix

  Patient#1 Patient#2 Patient#3 ...
Patient#m
Physician#1   10  1 ...   
 0
Physician#2   11  1 ...   
 1
Physician#3   01  0 ...   
 1
. ..  . ...   
 .   
. ..  . ...   
 .   
Physician#n   11  0 ...   
 0

1 indicates previous encouter and 0 otherwise.  My
aim is to identify physician group practice based on
the common patients they see.  Any suggestion on which
R package would best serve this purpose?  Thank you so
much!

Regards,
Kelvin

__
R-help@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] Combine vector of different length

2005-09-15 Thread Kel Lam
Hi group,  

Is there a quick way to cbind vector of different
length?  I should have checked out the previous
postings but somehow I can't access the list from CRAN
website.  Thanks for your help.

Regards,
Kelvin

__
R-help@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] Nested Logistic Regression

2005-05-10 Thread Kel Lam
Hi group,

I did a quick search in the archive but still couldn't
find a function that performs nested logistic
regression.
Please point me in the right direction.  Thanks a
million!

Regards,
Kel

__
R-help@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


RE: [R] gbm (Kelvin Lam)

2005-03-07 Thread Kel Lam
I have got some help from Greg during the course of my
boosting procedure.  I maybe able to help out.  What
exactly are you looking for?

Message: 46
Date: Fri, 4 Mar 2005 12:07:06 -0600
From: WeiWei Shi [EMAIL PROTECTED]
Subject: [R] gbm
To: R-help@stat.math.ethz.ch
Message-ID:
[EMAIL PROTECTED]
Content-Type: text/plain; charset=US-ASCII

Hi, there:
Is there anyone who read the codes for gbm package
before? Before i
sent this email, I also sent an email to ask for help
from the author,
Greg. But still I am wondering if someone here can
share some
understanding like the roadmap or document on the
implementation
too.

Thanks,

Ed.

__
R-help@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] gamm with correlation structure question

2005-01-21 Thread Kel Lam
Dear group,

I am trying to use gamm() in mgcv.  Here's the
scenario.  The data frame has approx. 110K
observations with information on paediatric
readmission binary outcome (Y/N) and total volume of
their most responsible physician as the covariate. 
Since any physician can have multiple patients, the
data contains clustering structure which I am trying
to account for.  My original formula is

a - gamm(readmission~s(volume,fx=F),
correlation=corCompSymm(0.1,form=~1|physician),
family=binomial)

During the first iteration, I received the warning 
Incompatible formulas for groups in random and
correlation.  Frankly I don't understand what it
means.  So I looked up the archive and saw one posted
by Prof. Thomas Lumley on lme as follows:

http://tolstoy.newcastle.edu.au/R/help/00a/0913.html

So I changed my formula as this:

a - gamm(readmission~s(volume,fx=F),
random=list(physician=~1),family=binomial)

However I still get an error message saying Lapack
routine dgesv : system is exactly singular.  

Can anyone suggest a better way to deal with my
problem.  Please let me know if you need further
details.  Thank you so much.

Regards,
Kel

__
R-help@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