CLASS-L Archives

July 2007

CLASS-L@LISTS.SUNYSB.EDU

Options: Use Proportional Font
Show HTML Part by Default
Condense Mail Headers

Message: [<< First] [< Prev] [Next >] [Last >>]
Topic: [<< First] [< Prev] [Next >] [Last >>]
Author: [<< First] [< Prev] [Next >] [Last >>]

Print Reply
Content-Type:
multipart/alternative; boundary="----_=_NextPart_001_01C7BD6D.62571B6A"
Sender:
"Classification, clustering, and phylogeny estimation" <[log in to unmask]>
Subject:
From:
Peter Flom <[log in to unmask]>
Date:
Tue, 3 Jul 2007 05:26:30 -0700
MIME-Version:
1.0
Comments:
Reply-To:
"Classification, clustering, and phylogeny estimation" <[log in to unmask]>
Parts/Attachments:
text/plain (1958 bytes) , text/html (2587 bytes)
William Shannon  wrote
<<< 
Do you have anyway to reduce the number of covariates before partitioning?  I would be concerned about the curse of dimensionality with 900 variables and 1,000 data points.  It would be very easy to find excellent classifiers based on noise.  Some suggest that a split data set (train on one subset randomly selected from the 1,000 data points and test on the remaining) overcomes this.  However, if X by chance due to the curse of dimensionality discriminates well than it will discriminate well in both the training and test data sets.
>>>>

and suggested the following experiment:
<<<<
1. Simulate a dataset consisting of 1,000 data points and 900 covariates where each covariate value comes from a normal(0,1) (or any other distribution) -- everything independent from each other.

2. Randomly assign the first 500 data points to group 1 and the second 500 data points to group 2

3. Fit your favorite discriminator to predict these two groups and see how well you can with random data.

4. After identifying the best fitting model removes those covariates and redo the analysis.
>>>

I did this.  

Using party, there were no viable splits in the original data.

Here is my code:


<<<
treetestdata <-  as.data.frame(mvrnorm(n = 1000, mu=rep(0,900), Sigma = diag(900)))
treetestdata$GROUP <- rep(c("G1","G2"), each=500)
treetest <- ctree(as.factor(GROUP) ~ ., data = treetestdata)

plot(treetest)
>>>

the result was a single node, with 500 subjects in each of the two groups.

There was, thus, no way to do steps 3 or 4.

Now, it's true that these variables are uncorrelated and mine in real life are correlated.  I can play around with that a little bit, but don't have time to do so right now.  If others are interested in playing around with this structure, I'd appreciate seeing any results.



----------------------------------------------
CLASS-L list.
Instructions: http://www.classification-society.org/csna/lists.html#class-l


ATOM RSS1 RSS2