<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"> <html xmlns="http://www.w3.org/1999/xhtml" xml:lang="en" lang="en"> <head> <meta http-equiv="Content-Type" content="text/html;charset=utf-8" /> <title> document_view — bioconductor.org </title> <base href="" /> <meta name="generator" content="Plone - http://plone.org" /> <meta content="sfadmin" name="DC.creator" /> <meta content="2006-10-10 13:20:39" name="DC.date.created" /> <meta content="2006-10-10 13:20:39" name="DC.date.modified" /> <meta content="Document" name="DC.type" /> <meta content="text/plain" name="DC.format" /> <!-- Basic crude style for Netscape4.x - This can be removed if you don't want the special NS4 look - it will still work, just be plain text instead. Column layout for Netscape4.x included --> <link rel="Stylesheet" type="text/css" href="../../../../../ploneNS4.css" /> <!-- Column style sheet. --> <style type="text/css" media="screen"><!-- @import url(http://bioconductor.org/ploneColumns.css); --></style> <!-- Main style sheets for CSS2 capable browsers --> <style type="text/css" media="screen"><!-- @import url(http://bioconductor.org/plone.css); --></style> <!-- Old style sheet from Plone 1.0, remove tal:condition="nothing" if you need to use the old styles. Will be removed in Plone 2.1. --> <!-- Alternate style sheets for the bigger/smaller text switcher --> <link rel="alternate stylesheet" type="text/css" media="screen" href="../../../../../ploneTextSmall.css" title="Small Text" /> <link rel="alternate stylesheet" type="text/css" media="screen" href="../../../../../ploneTextLarge.css" title="Large Text" /> <!-- Style sheet used for printing --> <link rel="stylesheet" type="text/css" media="print" href="../../../../../plonePrint.css" /> <!-- Style sheet used for presentations (Opera is the only browser supporting this at the moment) --> <link rel="stylesheet" type="text/css" media="projection" href="../../../../../plonePresentation.css" /> <!-- Internet Explorer CSS Fixes --> <!--[if IE]> <style type="text/css" media="all">@import url(http://bioconductor.org/ploneIEFixes.css);</style> <![endif]--> <!-- Custom style sheet if available --> <style type="text/css" media="all"><!-- @import url(http://bioconductor.org/ploneCustom.css); --></style> <link rel="shortcut icon" href="../../../../../favicon.ico" type="image/x-icon" /> <link rel="search" href="../../../../../search_form" title="Search this site" /> <link rel="up" href="../MachineLearning.1" title="Up one level" /> <!-- Disable IE6 image toolbar --> <meta http-equiv="imagetoolbar" content="no" /> <!-- A slot where you can insert CSS in the header from a template --> <!-- min-width support for IE via Javascript, required for tableless --> <!--[if IE]> <script type="text/javascript" src="http://bioconductor.org/plone_minwidth.js"></script> <![endif]--> <!-- A slot where you can insert elements in the header from a template --> <!-- A slot where you can insert javascript in the header from a template --> <!-- Common Plone ECMAScripts --> <!-- Pull-down ECMAScript menu, only active if logged in --> <!-- old google analytics here --> <!-- Define dynamic server-side variables for javascripts in this one --> <script type="text/javascript" src="../../../../../plone_javascript_variables.js"> </script> <script type="text/javascript" src="../../../../../plone_javascripts.js"> </script> <!-- Old JS from Plone 1.0, remove tal:condition="nothing" if you need to use the old pop-ups. Will be removed in Plone 2.1 --> </head> <body class="section-workshops"> <div id="visual-portal-wrapper"> <div id="portal-top"> <a href="ML-lab.R#documentContent" class="hiddenStructure">Skip to content.</a> <h1 id="portal-logo"> <a href="../../../../../index.html">bioconductor.org</a> </h1> <div id="portal-slogan"><p>Bioconductor is an open source and open development software project<br /> for the analysis and comprehension of genomic data.</p></div> <h5 class="hiddenStructure">Sections</h5> <ul id="portal-globalnav"><li id="portaltab-index_html" class="plain"><a href="../../../../../index.html" accesskey="t">Home</a></li><li id="portaltab-GettingStarted" class="plain"><a href="../../../../../GettingStarted" accesskey="t">Getting Started</a></li><li id="portaltab-overview" class="plain"><a href="../../../../../overview" accesskey="t">Overview</a></li><li id="portaltab-download" class="plain"><a href="../../../../../download" accesskey="t">Downloads</a></li><li id="portaltab-docs" class="plain"><a href="../../../../../docs" accesskey="t">Documentation</a></li><li id="portaltab-biocpub" class="plain"><a href="../../../../../pub" accesskey="t">Publications</a></li><li id="portaltab-workshops" class="selected"><a href="../../../../../workshops" accesskey="t">Workshops</a></li><li id="portaltab-cabig" class="plain"><a href="http://wiki.fhcrc.org/caBioc" accesskey="t">caBIG</a></li></ul> </div> <div class="visualClear"></div> <!-- The wrapper div. It contains the three columns. --> <div id="portal-columns" class="visualColumnHideNone"> <!-- start of the main and left columns --> <div id="visual-column-wrapper"> <!-- start of main content block --> <div id="portal-column-content" class="topmargin1"> <div id="content" class=""> <div class="documentContent" id="region-content"> <a name="documentContent"></a> <h1 class="documentFirstHeading">ML-lab.R</h1> <div class="documentDescription"></div> <div class="plain"> <p>################################################### ### chunk number 1: loadlib ################################################### library("Biobase") library("ALL")</p> <p>data(ALL, package="ALL")</p> <p>ALLBs = ALL[,grep("^B", as.character(ALL$BT))]</p> <p>ALLBCRNEG = ALLBs[, ALLBs$mol == "BCR/ABL" | ALLBs$mol =="NEG"]</p> <p>ALLBCRNEG$mol.biol = factor(ALLBCRNEG$mol.biol)</p> <p>numBN = length(ALLBCRNEG$mol.biol)</p> <p>ALLBCRALL1 = ALLBs[, ALLBs$mol == "BCR/ABL" | ALLBs$mol == "ALL1/AF4"]</p> <p>ALLBCRALL1$mol.biol = factor(ALLBCRALL1$mol.biol) numBA = length(ALLBCRALL1$mol.biol)</p> <p>################################################### ### chunk number 2: IQRfilter ################################################### lowQ = rowQ(ALLBCRNEG, floor(0.25 <em> numBN)) upQ = rowQ(ALLBCRNEG, ceiling(0.75 </em> numBN)) iqrs = upQ - lowQ ##you can also do: ##iqrs = esApply(ALLBCRNEG, 1, IQR)</p> <p>giqr = iqrs > quantile(iqrs, probs=.75) sum(giqr)</p> <p>BNsub = ALLBCRNEG<a href="ML-lab.R#refgiqr,">[giqr,]</a></p> <p>################################################### ### chunk number 3: distEx ################################################### dSub <- BNsub<a href="#ref1:60,">[1:60,]</a> eucD <- dist(t(exprs(dSub))) eucD@Size eucM <- as.matrix(eucD)</p> <p>################################################### ### chunk number 4: distHM ################################################### library("RColorBrewer") hmcol <- colorRampPalette(brewer.pal(10, "RdBu"))(256) heatmap(eucM, sym=TRUE, col=hmcol, distfun=function(x) as.dist(x))</p> <p>################################################### ### chunk number 5: helperfuns ################################################### closestN = function(distM, label) { loc = match(label, row.names(distM)) names(which.min(distM<a href="ML-lab.R#reflabel,-loc">[label,-loc]</a>)) }</p> <p>closestN(eucM, "03002")</p> <p>################################################### ### chunk number 6: rowttests ################################################### library("genefilter") tt1 = rowttests(BNsub, "mol.biol") numToSel <- 50</p> <p>################################################### ### chunk number 7: topNumToSel ################################################### tt1ord = order(abs(tt1$statistic), decreasing=TRUE)</p> <p>top50 = tt1ord<a href="ML-lab.R#ref1:numToSel">[1:numToSel]</a></p> <p>BNsub1 = BNsub<a href="ML-lab.R#reftop50,">[top50,]</a></p> <p>################################################### ### chunk number 8: rowIQRsFun ################################################### rowIQRs = function(eSet) { if (is(eSet, "exprSet")) numSamp = length(sampleNames(eSet)) else numSamp = ncol(eSet) lowQ = rowQ(eSet, floor(0.25 <em> numSamp)) upQ = rowQ(eSet, ceiling(0.75 </em> numSamp)) upQ - lowQ }</p> <p>################################################### ### chunk number 9: standardize ################################################### standardize = function(x) (x - rowMedians(x)) / rowIQRs(x)</p> <p>exprs(BNsub1) = standardize(exprs(BNsub1))</p> <p>################################################### ### chunk number 10: heatmap ################################################### library("RColorBrewer") hmcol <- colorRampPalette(brewer.pal(10, "RdBu"))(256) spcol <- ifelse(BNsub1$mol.biol=="BCR/ABL", "goldenrod", "skyblue") heatmap(exprs(BNsub1), col=hmcol, ColSideColors=spcol)</p> <p>################################################### ### chunk number 11: knn1 ###################################################</p> <p>library("class")</p> <p>a1 = knn.cv(t(exprs(BNsub1)), BNsub1$mol.biol)</p> <p>ctab1 = table(a1, BNsub1$mol.biol)</p> <p>errrate = (ctab1["BCR/ABL", "NEG"] + ctab1["NEG", "BCR/ABL"])/sum(ctab1)</p> <p>################################################### ### chunk number 12: selectk ################################################### alist = NULL for(i in 1:4) alist[<a href="ML-lab.R#refi">[i]</a>] = knn.cv(t(exprs(BNsub1)), BNsub1$mol.biol, k=i)</p> <p>sapply(alist, function(x) { ct1 = table(x, BNsub1$mol.biol) (ct1["BCR/ABL", "NEG"] + ct1["NEG", "BCR/ABL"]) / sum(ct1) })</p> <p>################################################### ### chunk number 13: MLInterfaces ################################################### library("MLInterfaces") </p> <p>################################################### ### chunk number 14: knn ################################################### knnResult <- MLearn("mol.biol", BNsub1, "knn", 1:50) knnResult</p> <p>################################################### ### chunk number 15: knn-confusion ################################################### confuMat(knnResult)</p> <p>################################################### ### chunk number 16: xval ################################################### knnXval <- xval(BNsub1, "mol.biol", knnB, xvalMethod="LOO")</p> <p>################################################### ### chunk number 17: xval-table ################################################### table(given=BNsub1$mol.biol,predicted=knnXval)</p> <p>################################################### ### chunk number 18: xval-feature-selection-function ################################################### t.fun<-function(data, fac) { (abs(rowttests(data,data[<a href="#reffac">[fac]</a>], tstatOnly=FALSE)$statistic)) }</p> <p>################################################### ### chunk number 19: realXval ################################################### BNx = BNsub exprs(BNx) = standardize(exprs(BNx))</p> <p>lk3f <- xval(BNx, "mol.biol", knnB, xvalMethod="LOO", 0:0, fsFun=t.fun, fsNum=50) table(given=BNx$mol.biol, predicted=lk3f$out)</p> <p>################################################### ### chunk number 20: ################################################### sessionInfo()</p> </div> <div class="discussion"> </div> </div> </div> </div> <!-- end of main content block --> <!-- start of the left (by default at least) column --> <div id="portal-column-one"> <div class="visualPadding"> <br><br><br> <!-- disabled left slot image <img tal:replace="structure nocall:here/pict.jpg" /> --> <div class="portlet" id="portlet-navigation-tree"> <div> <h5>Navigation</h5> <div class="portletBody"> <div class="portletContent odd"> <a href="../../../../../GettingStarted/index.html" accesskey="n" class="navItem navLevel1" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">Getting Started</span> </a> <a href="../../../../../overview/index.html" accesskey="n" class="navItem navLevel1" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">Overview</span> </a> <a href="../../../../../download/index.html" accesskey="n" class="navItem navLevel1" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">Downloads</span> </a> <a href="../../../../../docs/index.html" accesskey="n" class="navItem navLevel1" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">Documentation</span> </a> <a href="../../../../../pub/index.html" accesskey="n" class="navItem navLevel1" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">Publications</span> </a> <a href="../../../../index.html" accesskey="n" class="navItem navLevel1" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">Workshops</span> </a> <a href="../../../../2010/index.html" accesskey="n" class="navItem navLevel2" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2010</span> </a> <a href="../../../../2009/index.html" accesskey="n" class="navItem navLevel2" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2009</span> </a> <a href="../../../../2008/index.html" accesskey="n" class="navItem navLevel2" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2008</span> </a> <a href="../../../../2007/index.html" accesskey="n" class="navItem navLevel2" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2007</span> </a> <a href="../../../index.html" accesskey="n" class="navItem navLevel2" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2006</span> </a> <a href="../../index.html" accesskey="n" class="navItem navLevel3" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">R/Bioc Intro Course - October</span> </a> <a href="../index.html" accesskey="n" class="navItem navLevel4" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">labs</span> </a> <a href="../cDNA/index.html" accesskey="n" class="navItem navLevel5" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">cDNA (Two Color) Microarray</span> </a> <a href="index.html" accesskey="n" class="navItem navLevel5 currentNavItem" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">MachineLearning</span> </a> <a href="../../../BioC2006/index.html" accesskey="n" class="navItem navLevel3" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">BioC2006</span> </a> <a href="../../../aucklandwksp/index.html" accesskey="n" class="navItem navLevel3" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">Workshops on S, R and Bioconductor</span> </a> <a href="../../../biocintro_april/index.html" accesskey="n" class="navItem navLevel3" title="April Intro R/Bioc Course at FHCRC"> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">R/Bioc Intro Course</span> </a> <a href="../../../rforbioinformatics/index.html" accesskey="n" class="navItem navLevel3" title="Course materials. The course took place at FHCRC Seattle in January 2006"> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">R for Bioinformatics</span> </a> <a href="../../../../2005/index.html" accesskey="n" class="navItem navLevel2" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2005</span> </a> <a href="../../../../2004/index.html" accesskey="n" class="navItem navLevel2" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2004</span> </a> <a href="../../../../2003/index.html" accesskey="n" class="navItem navLevel2" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2003</span> </a> <a href="../../../../2002/index.html" accesskey="n" class="navItem navLevel2" title="Workshops and courses 2002"> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">2002</span> </a> <a href="../../../../../developers/index.html" accesskey="n" class="navItem navLevel1" title=""> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">Developers</span> </a> <a href="../../../../../News/index.html" accesskey="n" class="navItem navLevel1" title="BioC Project News"> <!-- Disable Folder icon <img src="" height="16" width="16" alt="" class="navIcon" tal:condition="not:isAnon" tal:attributes="src python:portal_url+'/'+sibling.getIcon(1); title sibling/Type;" /> --> <span class="navItemText">News</span> </a> </div> </div> </div> </div> </div> </div> <!-- end of the left (by default at least) column --> </div> <!-- end of the main and left columns --> <!-- start of right (by default at least) column --> <div id="portal-column-two"> <div class="visualPadding"> <div> <!-- The Related Items box --> <script type="text/javascript" src="plonesearchbox_utils.js"> </script> <div class="portlet"> <!-- <h5 i18n:translate="searchbox_title">Search Box</h5> --> <div class="portletBody"> <form name="searchbox_form" id="searchbox_form" method="post" onsubmit="return process()" action="http://bioconductor.org/workshops/2006/biocintro_oct/labs/MachineLearning/redirectToUrl" target="_blank"> <div class="portletContent odd"> <div id="searchbox_fields"> <input type="text" name="toSearch" id="toSearch" size="22" /> </div> <select name="choice" id="choice" onchange="return generateForm();" style="margin: 3px 0px 0px 0px;"> <option value="http://bioconductor.org/workshops/2006/biocintro_oct/labs/MachineLearning/ML-lab.R/search?SearchableText=" id="search_portal">In this site</option> <option value="http://bioconductor.org/workshops/2006/biocintro_oct/labs/MachineLearning/ML-lab.R/search?path=/Plones/rgentlem/bioconductor/workshops/2006/biocintro_oct/labs/MachineLearning/ML-lab.R&SearchableText=" id="search_rubric">In this folder</option> <option value="http://google.com/search?sitesearch=www.bioconductor.org&q=%(text)s">Google this site</option> </select> <input class="context searchButton" type="submit" value="Search" style="margin: 3px;" /> </div> <!-- <div class="portletContent odd"> <input class="context searchButton" type="submit" value="Search" style="margin: 3px;" i18n:attributes="value"/> </div> <div class="portletContent even"> <a href="" tal:attributes="href string:${portal_url}/search_form" i18n:translate="advanced_search_link">Advanced search</a> </div> --> </form> <script type="text/javascript"> <!-- /*in case the first choice of the combobox is an url with several fields*/ clearForm(); generateForm(); --> </script> </div> </div> </div> <div class="portlet" id="portlet-news"> <h5>News</h5> <div class="portletBody"> <div class="portletContent odd"> <div class="portletDetails"> <a href="../../../../../News/2009-10-28-Release-2.5" class="date" title="BioC 2.5 Released">2009-10-26</a> <p> BioC 2.5, consisting of 352 packages and designed to work with R 2.10.z, was released today. </p> </div> </div> <div class="portletContent even"> <div class="portletDetails"> <a href="../../../../../News/RinNYT" class="date" title="Print Media Recognition">2009-01-07</a> <p> R, the open source platform used by Bioconductor, featured in a series of articles in the New York Times. </p> </div> </div> <div class="portletContent odd"> <a href="../../../../../news" class="portletMore"> More... </a> </div> </div> </div> </div> </div> <!-- end of the right (by default at least) column --> </div> <!-- end column wrapper --> <div class="visualClear"></div> <hr class="netscape4" /> <div id="portal-footer"> © 2003-2009 BioConductor. All Rights Reserved. <script type="text/javascript"> var gaJsHost = (("https:" == document.location.protocol) ? "https://ssl." : "http://www."); document.write(unescape("%3Cscript src='" + gaJsHost + "google-analytics.com/ga.js' type='text/javascript'%3E%3C/script%3E")); </script> <script type="text/javascript"> try{ var pageTracker = _gat._getTracker("UA-357281-1"); pageTracker._trackPageview(); } catch(err) {} </script> </div> </div> </body> </html>