<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>TheSpicyChicken &#187; Science</title>
	<atom:link href="http://thespicychicken.com/category/science/feed/" rel="self" type="application/rss+xml" />
	<link>http://thespicychicken.com</link>
	<description>n'th generation research</description>
	<lastBuildDate>Thu, 10 Jun 2010 22:14:36 +0000</lastBuildDate>
	<generator>http://wordpress.org/?v=2.9.2</generator>
	<language>en</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
			<item>
		<title>Submitting to &#8220;The Journal&#8221;</title>
		<link>http://thespicychicken.com/2010/05/03/submitting-to-the-journal/</link>
		<comments>http://thespicychicken.com/2010/05/03/submitting-to-the-journal/#comments</comments>
		<pubDate>Mon, 03 May 2010 18:17:16 +0000</pubDate>
		<dc:creator>Spicy</dc:creator>
				<category><![CDATA[Science]]></category>
		<category><![CDATA[nejm]]></category>

		<guid isPermaLink="false">http://thespicychicken.com/?p=203</guid>
		<description><![CDATA[Science it stupid-fun. I mean, honestly, as a graduate student in informatics I get to sit in front of my computer all day, musing about problems, programming scripts, and testing hypotheses. It&#8217;s the greatest job in the world and every so often you hit on something really cool and get to publish a paper. Sometimes [...]]]></description>
			<content:encoded><![CDATA[<p>Science it stupid-fun. I mean, honestly, as a graduate student in informatics I get to sit in front of my computer all day, musing about problems, programming scripts, and testing hypotheses. It&#8217;s the greatest job in the world and every so often you hit on something really cool and get to publish a paper. Sometimes you figure out a new way to analyze the data or you apply an old technique to a new field. Either way, sharing your discovery is both fun and gratifying.</p>
<p>The papers you publish are the lifeblood of a scientist&#8217;s career and as a graduate student will form the basis of my PhD thesis. I&#8217;ve been lucky so far to have worked with some very talented people and publish a couple of papers in my first two years.  However, recently I stumbled upon a discovery which could get me a paper in, by most metrics, the most prestigious journal in the world, The New England Journal of Medicine. Just typing this sentence blows my mind. Anyway, because I don&#8217;t think about much else, I&#8217;m going to blog about my experience submitting a paper to NEJM and will share all the ups and downs&#8230;. cross your fingers for ups!</p>
<hr/><br/></p>
<h3>Update: April 28, 2010</h3>
<p>The goal was to submit the paper today. But right before my advisor and I were about to submit the final draft of the paper, and I realized that I left out a whole set of analyses! Doh, now I have to go back to the database to extract more patient records and analyze another set. Looks like, it won&#8217;t be today.</p>
<hr/><br/></p>
<h3>Update: April 29, 2010</h3>
<p>Okay, made a whole bunch of new figures today, modified the text where needed and sat down with my advisor again. We made some final edits (some of which we should have caught in the previous 16 iterations of the manuscript), but alas we feel great about the paper. Personally, I think it&#8217;s a work of art and couldn&#8217;t imagine it getting rejected <img src='http://thespicychicken.com/wp-includes/images/smilies/icon_biggrin.gif' alt=':D' class='wp-smiley' /> .  Don&#8217;t worry, I&#8217;m not allowing my hopes to get too high.  My advisor let me submit the paper and click on all the &#8220;double-check-your-submission&#8221; links. In his words &#8220;you&#8217;re not going to be doing this that often, so you might as well click around.&#8221; It was fun and it&#8217;s amazing I even get to submit to the NEJM, really.</p>
<p>Paper status is now &#8220;Submitted&#8221; and we got a confirmation email that it will be forwarded to the editors. Time to order some flowers and wine for our collaborators!</p>
<hr/><br/></p>
<h3>Update: April 30, 2010</h3>
<p>I have been neurotically checking the author website at NEJM and today our paper was &#8220;Assigned to an Editor&#8221;, I have no idea exactly what that means, but it&#8217;s movement!</p>
<hr/><br/></p>
<h3>Update: May 2nd, 2010</h3>
<p>I went on the website today (yes, and yesterday too), not expecting any updates over the weekend, but much to my surprise the status of the paper has changed.  It looks like it made passed the editor and they are now looking for peer reviewers.  According to a close friend getting passed the editors is a big deal and a lot of papers don&#8217;t make it passed that point.  Very exciting!  &#8212; I just read on the NEJM website about the editorial process and it seems that it has been read by the Editor-in-chief and also another expert editor and has passed both of their filters.  Here&#8217;s to the reviewers liking it! [crosses-fingers]</p>
<p>Paper status is now &#8220;Searching for Reviewers.&#8221;</p>
<hr/><br/></p>
<h3>Update: May 10th, 2010</h3>
<p>It&#8217;s been one week and one day since the last update, and I was starting to think that maybe the next update would come when the reviews came back. That is not the case, however. Looks like they found some reviewers for the paper! Not sure how long they have to review the paper, probably a few weeks at least I imagine. I<del datetime="2010-05-10T17:52:55+00:00">&#8216;ll check with my PI and add that info in here.</del> Looks like its 2 to 3 weeks.</p>
<p>Paper status is now &#8220;Out for Review.&#8221;</p>
<hr/><br/></p>
<h3>Update: May 20th, 2010</h3>
<p>The reviews of the paper have been returned to the Journal, but a decision has not yet been made. The next step is for the editorial board to discuss the paper and make a decision. I&#8217;m not sure how long this process takes, but every other time the status included the word &#8220;editor&#8221; it only took a few days to get another update. Perhaps I will hear something by weeks end. YIKES! This is getting real.</p>
<p>Paper status is now &#8220;With the Editor.&#8221;</p>
<hr/><br/></p>
<h3>Update: May 25th, 2010</h3>
<p>It&#8217;s been a serious roller-coaster ride and I have had quite the experience submitting to the world&#8217;s most prestigious journal.  The reviews came back a few days ago and we just got a chance to read them today. As far as reviews for The Journal go, they are quite benevolent and really quite positive. They made some great suggestions on how to improve the paper and I have already made most of the changes they suggested (they were quite minor). Unfortunately, for us, the editors were not as favorable on us as the reviewers were and they rejected our paper today. Perhaps the paper is too bold for the New England Journal, perhaps it&#8217;s ahead of its time.  Perhaps, they just don&#8217;t like me or my parentage &#8212; I simply can&#8217;t be sure.  But, what I do know is that this is some of the best work I have ever done and there are many other Journals out there. I will get the good word out about this discovery one way or another. </p>
<p>Paper status is now &#8220;Rejected.&#8221;</p>
]]></content:encoded>
			<wfw:commentRss>http://thespicychicken.com/2010/05/03/submitting-to-the-journal/feed/</wfw:commentRss>
		<slash:comments>8</slash:comments>
		</item>
		<item>
		<title>R Tip: Fitting Sigmoidal Data</title>
		<link>http://thespicychicken.com/2010/02/17/r-tip-fitting-sigmoidal-data/</link>
		<comments>http://thespicychicken.com/2010/02/17/r-tip-fitting-sigmoidal-data/#comments</comments>
		<pubDate>Wed, 17 Feb 2010 07:12:54 +0000</pubDate>
		<dc:creator>Spicy</dc:creator>
				<category><![CDATA[Development]]></category>
		<category><![CDATA[Science]]></category>
		<category><![CDATA[math]]></category>
		<category><![CDATA[R]]></category>
		<category><![CDATA[sigmoid]]></category>

		<guid isPermaLink="false">http://thespicychicken.com/?p=175</guid>
		<description><![CDATA[Using R to fit a sigmoid function to experimental data.]]></description>
			<content:encoded><![CDATA[<p>
  Sigmoid functions are our friends and sometimes you have data which you would like to fit with a sigmoid function. We can use R to find such a fit.  First let us look at a sigmoid function.
</p>
<div class="code">
  y = 1 / (1 + exp( -a*x + b) )
</div>
<p>
  Now let&#8217;s say you are given a vector x and y, say:
</p>
<div class="code">
  x = c(0.00,0.02,0.04,0.06,0.08,0.10,0.12,0.14,0.16,0.18,0.20,0.24,0.26,<br/><br />
0.28,0.30,0.34,0.40,0.42,0.48,0.54,0.56,0.64,1.00)<br/><br />
  y = c(0.409742,0.319277,0.530120,0.377778,0.357143,0.608696,0.315789,<br/><br />
0.692308,0.642857,0.636364,0.750000,0.000000,0.833333,1.000000,<br/><br />
0.000000,1.000000,1.000000,1.000000,1.000000,1.000000,1.000000,<br/><br />
1.000000,1.000000)
</div>
<p>
  In this case the y values here represent probabilities and one thing you’ll notice is that we have probs of 1 and 0. Both of which are bad. So we apply a little “laplace smoothing” to them:
</p>
<div class="code">
  y[y==0] = 0.001<br/><br />
  y[y==1] = 0.999
</div>
<p>
Now let&#8217;s look at what the data looks like.
</p>
<div class="code">
plot(x, y)
</div>
<p><a href="http://thespicychicken.com/wp-content/uploads/2010/02/RTip_Sigmoid_Fig1.png"><img src="http://thespicychicken.com/wp-content/uploads/2010/02/RTip_Sigmoid_Fig1-300x279.png" alt="" title="RTip_Sigmoid_Fig1" width="300" height="279" class="alignnone size-medium wp-image-186" /></a></p>
<p>
Well, it may be sigmoidal, maybe not. For now let’s assume we think it is. Which we do for the most part.
</p>
<p>
Okay, now let’s solve for a line in our sigmoid function:
</p>
<div class="code">
y = 1 / (1 + exp( a*x + b) )<br/><br />
1 + exp( a*x + b) = 1/ y<br/><br />
a*x + b = log ( (1/ y) &#8211; 1 )<br/>
</div>
<p>
Now the left hand side of the equation is a line and the right hand side is some logarithm of the y data. We can plot x versus this right hand side:
</p>
<div class="code">
new_y = log( 1 / y &#8211; 1 )<br/><br />
plot(x, new_y)
</div>
<p><a href="http://thespicychicken.com/wp-content/uploads/2010/02/RTip_Sigmoid_Fig2.png"><img src="http://thespicychicken.com/wp-content/uploads/2010/02/RTip_Sigmoid_Fig2-300x272.png" alt="" title="RTip_Sigmoid_Fig2" width="300" height="272" class="alignnone size-medium wp-image-187" /></a></p>
<p>
Looks pretty interesting and hopefully at this point it also looks kinda linear, which it kinda does.
</p>
<p>
Now let’s fit it with a line:
</p>
<div class="code">
lm.res <- lm( new_y ~ x )<br/><br />
lm.res
</div>
<p>Which produces this output:</p>
<div class="code">
Coefficients:<br/><br />
(Intercept)            x<br/><br />
1.122      -11.647<br/>
</div>
<p>We can also test the significance of the fit with an ANOVA.</p>
<div class="code">
anova(lm.res)
</div>
<p>Which produces this output:</p>
<div class="code">
Analysis of Variance Table<br/><br />
Response: new_y<br/><br />
Df Sum Sq Mean Sq F value    Pr(>F)<br/><br />
x          1 172.80 172.802  14.641 0.0009834 ***<br/>
</div>
<p>
And we can plot the resulting fit in linear space:
</p>
<p><a href="http://thespicychicken.com/wp-content/uploads/2010/02/RTip_Sigmoid_Fig3.png"><img src="http://thespicychicken.com/wp-content/uploads/2010/02/RTip_Sigmoid_Fig3-300x276.png" alt="" title="RTip_Sigmoid_Fig3" width="300" height="276" class="alignnone size-medium wp-image-188" /></a></p>
<p>
Now let’s see how our fit looks back in normal space using our formula with our derived a and b values.
</p>
<div class="code">
a = -11.647<br/><br />
b = 1.122<br/><br />
plot(x, y)<br/><br />
sim_x = (1:101-1)/100<br/><br />
points(sim_x, 1/(1+exp(a*sim_x+b)), type=&#8221;l&#8221;)<br/>
</div>
<p><a href="http://thespicychicken.com/wp-content/uploads/2010/02/RTip_Sigmoid_Fig4.png"><img src="http://thespicychicken.com/wp-content/uploads/2010/02/RTip_Sigmoid_Fig4-300x271.png" alt="" title="RTip_Sigmoid_Fig4" width="300" height="271" class="alignnone size-medium wp-image-189" /></a></p>
<p>
Voila! We have fit a sigmoid function to our data.</p>
]]></content:encoded>
			<wfw:commentRss>http://thespicychicken.com/2010/02/17/r-tip-fitting-sigmoidal-data/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Science Jobs&#8230; and the people that are looking for them.</title>
		<link>http://thespicychicken.com/2009/07/29/science-jobs-and-the-people-that-are-looking-for-them/</link>
		<comments>http://thespicychicken.com/2009/07/29/science-jobs-and-the-people-that-are-looking-for-them/#comments</comments>
		<pubDate>Wed, 29 Jul 2009 16:47:16 +0000</pubDate>
		<dc:creator>Spicy</dc:creator>
				<category><![CDATA[General]]></category>
		<category><![CDATA[Science]]></category>

		<guid isPermaLink="false">http://thespicychicken.com/?p=140</guid>
		<description><![CDATA[A few friends of mine have started a new Science Jobs blog were they muse about what it&#8217;s like to look for and work in science related jobs in this economic climate.  Check it out at SearchScienceJobs.com
]]></description>
			<content:encoded><![CDATA[<p>A few friends of mine have started a new Science Jobs blog were they muse about what it&#8217;s like to look for and work in science related jobs in this economic climate.  Check it out at <a href="http://searchsciencejobs.com">SearchScienceJobs.com</a></p>
]]></content:encoded>
			<wfw:commentRss>http://thespicychicken.com/2009/07/29/science-jobs-and-the-people-that-are-looking-for-them/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>The Science of Hearing &#8211; We need you!</title>
		<link>http://thespicychicken.com/2009/05/26/the-science-of-hearing-we-need-you/</link>
		<comments>http://thespicychicken.com/2009/05/26/the-science-of-hearing-we-need-you/#comments</comments>
		<pubDate>Tue, 26 May 2009 18:42:00 +0000</pubDate>
		<dc:creator>Spicy</dc:creator>
				<category><![CDATA[Science]]></category>
		<category><![CDATA[hearing]]></category>
		<category><![CDATA[mosquito]]></category>

		<guid isPermaLink="false">http://thespicychicken.com/?p=135</guid>
		<description><![CDATA[A friend of mine, and neuroscience graduate student at Case Western Reserve University are conducting a study into the genetic basis of hearing the following sound:
http://tinyurl.com/2wjoxy
NPR recently did a story that teenagers have been using the above sound as a ring tone to evade the prying eyes ears of adults.
So we started a little investigation [...]]]></description>
			<content:encoded><![CDATA[<p>A friend of mine, and neuroscience graduate student at Case Western Reserve University are conducting a study into the genetic basis of hearing the following sound:</p>
<p><a href="http://tinyurl.com/2wjoxy">http://tinyurl.com/2wjoxy</a></p>
<p>NPR recently did a story that teenagers have been using the above sound as a ring tone to evade the prying <span style="text-decoration: line-through;">eyes</span> ears of adults.</p>
<p>So we started a little investigation of our own, and we need your help!  Just click on that sound file, play it, and respond in the comments with your age, gender, and if you can hear it or not.</p>
<p>SCIENCE!</p>
]]></content:encoded>
			<wfw:commentRss>http://thespicychicken.com/2009/05/26/the-science-of-hearing-we-need-you/feed/</wfw:commentRss>
		<slash:comments>12</slash:comments>
<enclosure url="http://tinyurl.com/2wjoxy" length="160625" type="audio/mpeg" />
		</item>
		<item>
		<title>Finally&#8230; back to Science.</title>
		<link>http://thespicychicken.com/2008/11/07/finally_back_to_science/</link>
		<comments>http://thespicychicken.com/2008/11/07/finally_back_to_science/#comments</comments>
		<pubDate>Fri, 07 Nov 2008 18:39:01 +0000</pubDate>
		<dc:creator>Spicy</dc:creator>
				<category><![CDATA[Science]]></category>
		<category><![CDATA[algorithm]]></category>
		<category><![CDATA[CUDA]]></category>
		<category><![CDATA[GPU]]></category>
		<category><![CDATA[k-means]]></category>

		<guid isPermaLink="false">http://thespicychicken.com/?p=52</guid>
		<description><![CDATA[With the election behind and a collective sigh of relief from the science community we can get back to the good stuff&#8211; throwing some fuel into the engine we call SCIENCE!  For my part, I got some much needed time to learn about and implement some simple algorithms utilizing CUDA.  CUDA is a c-extension from [...]]]></description>
			<content:encoded><![CDATA[<p>With the election behind and a collective sigh of relief from the science community we can get back to the good stuff&#8211; throwing some fuel into the engine we call SCIENCE!  For my part, I got some much needed time to learn about and implement some simple algorithms utilizing CUDA.  CUDA is a c-extension from Nvidia for leveraging the massive computational power in your GPU.  In my case my GeForce 8600M graphics card in my MacBook Pro.</p>
<p>I have to give some serious props to the devs of CUDA because it is very straight forward and makes programming for the GPU rather intuitive.  I have been able to implement the Lloyd&#8217;s Algorithm for the <a href="http://en.wikipedia.org/wiki/K-means_algorithm">K-Means problem</a> in just a few days.  The speed increase is graphically represented in the following graph.</p>
<p><a href="http://thespicychicken.com/wp-content/uploads/2008/11/gpuvscpu_kmeans.png"><img class="alignnone size-full wp-image-55" title="gpuvscpu_kmeans" src="http://thespicychicken.com/wp-content/uploads/2008/11/gpuvscpu_kmeans.png" alt="" width="499" height="351" /></a></p>
<p>In green is the time in seconds required to compute the kmeans solution using the CPU and in red is the time required for the GPU.  The variance in the times comes from the nature of Lloyd&#8217;s algorithm which, depding on the input, can take a variable amount of time.  The GPU times are also varied but it is difficult to see in the figure.</p>
<p>These results are just from running the application on my MBP.  If you have a machine with many graphics processors (or many machines each with their own GPU) then you can parallelize this even more and massively boost your computational power.  In my opinion the GPU architecture is the future of scientific computing as so many of our computational challenges can be addressed using parallelization.</p>
<p>For those that are interested I will be writing up some CUDA tutorials covering Installation/Set Up, A &#8220;Hello World&#8221; example, and the full K-Means implementation over at <a href="http://macresearch.org">MacResearch.org</a>.</p>
]]></content:encoded>
			<wfw:commentRss>http://thespicychicken.com/2008/11/07/finally_back_to_science/feed/</wfw:commentRss>
		<slash:comments>4</slash:comments>
		</item>
		<item>
		<title>Chem(o)informatics</title>
		<link>http://thespicychicken.com/2008/10/29/chemoinformatics/</link>
		<comments>http://thespicychicken.com/2008/10/29/chemoinformatics/#comments</comments>
		<pubDate>Wed, 29 Oct 2008 18:02:57 +0000</pubDate>
		<dc:creator>Spicy</dc:creator>
				<category><![CDATA[Science]]></category>
		<category><![CDATA[cheminformatics]]></category>
		<category><![CDATA[chemoinformatics]]></category>

		<guid isPermaLink="false">http://thespicychicken.com/?p=45</guid>
		<description><![CDATA[As a PhD student you start out by doing rotations, quarter-long research projects, in various Professor&#8217;s lab.  In my rotation I am attempting to learn all there is to know about Chem(o)informatics.  I write it that way because there is an apparent dispute as to what exactly to call this field of study (Cheminformatics versus [...]]]></description>
			<content:encoded><![CDATA[<p>As a PhD student you start out by doing rotations, quarter-long research projects, in various Professor&#8217;s lab.  In my rotation I am attempting to learn all there is to know about Chem(o)informatics.  I write it that way because there is an apparent dispute as to what exactly to call this field of study (Cheminformatics versus Chemoinformatics).  As I do with many problems, I turned to Google for the answer. A quick search revealed that the term Cheminformatics is more prominent. &#8220;Ah ha!&#8221; I exclaimed.  For the all knowing Google could not be wrong and so it must be.  However, an itch in the back of my mind made me do a search using Google Scholar.  Alas, to my dismay, it showed the opposite.  In an attempt to break the tie and resolve this conflict I attempted another search in PubMed.  Surely that would settle things.  The results of my studies are below.</p>
<p><a href="http://thespicychicken.com/wp-content/uploads/2008/10/chemo-chem.png"><img class="alignnone size-full wp-image-46" title="chemo-chem" src="http://thespicychicken.com/wp-content/uploads/2008/10/chemo-chem.png" alt="" width="305" height="222" /></a></p>
<p>(Perhaps definatively, my spellcheck says there both wrong.)</p>
]]></content:encoded>
			<wfw:commentRss>http://thespicychicken.com/2008/10/29/chemoinformatics/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
	</channel>
</rss>
