# Russell's Blog

## A desirable extinction

Posted by Russell on March 25, 2010 at 3:19 p.m.
Some weeks ago, Buzz (my cat) escaped out my front door while I was carrying my bicycle into the apartment. For ten or twenty minutes, he romped through the ivy and bushes around my apartment while I followed him around rattling a bag of cat treats. Eventually, he let me pick him up and take him back inside. Naturally, he picked up a few fleas. Naturally, they have multiplied.

Oddly, the fleas don't seem to like Neil very much, nor do they like me. It's just poor Buzz that's beset by the nasty little critters.

### Figure 1: A flea.

As it happens, I've been thinking about endogenous metrics for estimating the sampling quality of an environmental shotgun sequencing dataset, and Buzz's little problem presented an opportunity to play with a simplified problem. So, I have decided to make Buzz, or rather his fleas, into a small experiment in ecology. I am going to try to see if I can drive them into extinction.

Now, this is normally what a pet owner does when they discover their pet has contracted some sort of annoying parasite, but I decided to take a more quantitative approach.

### Figure 2: A cat.

It's simple enough to count fleas on a cat, if the cat is willing to cooperate. Buzz loves the flea comb, and will gleefully hop onto the coffee table and wait to be combed if you show it to him. So, in the interest of science, I convinced my roommate to count the number of passes I made with the flea comb and how many fleas I captured (posterity will remember your efforts, Mehdi). Using his tally, I plotted the cumulative number of passes verses the cumulative number of fleas.

### Figure 3: Fleas captured

As expected, it became somewhat more difficult to capture the next flea as more fleas were captured, suggesting a depletion curve. The value of the asymptote should be the actual number of fleas on Buzz at the time, and reaching that number would imply local extinction for the fleas. Of course, there are probably other fleas lurking about that would recolonize Buzz. In principle, if I were to repeat the exercise frequently enough, Buzz would become a sink for fleas, and their migration to his fur would gradually deplete them from the environment.

There are a couple of different ways to model the impact of the combing on the flea population, with various advantages and disadvantages. All we really want to do here is to estimate the value of the asymptote, and so a simple model is probably sufficient. I showed this data to my fried Sharon Shewmake, an economics graduate student. Sharon, after editorializing on the endeavor ("Ew."), suggested this very simple model.

Assume that Buzz is not going to sit still long enough for the fleas to reproduce, for more fleas to migrate to his fur, and that the fleas already on his fur are going to stay put unless captured. Thus, there is a fixed initial population which only changes as a result of capturing fleas. Next, we assume that any given flea is equally likely to be captured on a single pass of the comb. So, the expectation value for number of fleas captured on a single pass is the product of the current population and the probability of capturing a flea.

$N_{i+1}=N_i-pN_i$

where N is the population of fleas and p is the probability of any particular flea being captured on a single pass. One could tart this up a bit by modeling it as a stochastic process and executing a bunch of Monte Carlo trials until the outcomes converge, but that seems like overkill for a simple single variable problem like this. We will put up with the intellectual inconvenience of capturing fractional fleas.

This is a little easier to see if we let N represent the number of fleas remaining on the cat, rather than the number of fleas captured.

$N_{i+1}=qN_i$

If we stretch our credulity far enough to imagine this as a continuous function, we can express it as a differential equation.

$\frac{d N}{d t}=\alpha N$

Sorry if this bothers you. Not only are we extracting fractional fleas, but we are now modeling the combing process as a sort of flea-killing-combine continuously mowing its way through the fur. This is a model, so you shouldn't be surprised to find massless rope and spherical cows. Anyway, it has a nice easy solution.

$N(t)=e^{\alpha t}$

Well, what the heck. This is a decaying function, so let's pluck a minus sign out of the exponential factor, and maybe tack on a scale factor for the initial population.

$N(t)=N_0 e^{-\alpha t}$

While we're at it, why don't we go back to letting the function stand for the number of fleas captured, rather than the fleas on the cat.

$N(t)=N_0(1-e^{-\alpha t})$

This gives us a nice function to use for a linear regression. A little help from scipy, and we find that the initial population is estimated at 39.7 fleas, and the decay factor is 0.011.

### Figure 4: Flea population

I captured 34 fleas, so that means I missed about five or six. In order to be reasonably confident that I'd captured all 39 fleas, I would have had to continued for about 400 passes with the comb, instead of 173. Buzz is a patient cat, but he started to loose interest around 120 passes, and had to be fetched back onto the coffee table a few time times during the last 50 passes. My guess is that 400 passes would require some kind of sedative. On the other hand, he does seem to like Guinness, so there may be something to that.

Science has been served. I'm going to the pet store to buy some flea collars.

r4 on May 13, 2013 at 2:47 p.m.

1 x carte r4 et la cle USB pour y inserer une carte Micro SD.
La carte mémoire Micro SD n'est pas incluse.abercrombie

kredyt na start on May 13, 2013 at 4:35 p.m.

I was so lucky to visit this site I really want to recommend this to everyone cause theres no doubt about its content I really love this one thanks for sharing this!I hope that I can see new interesting topics Im looking forward to thatkredyt na start

Small Business Directory on May 14, 2013 at 5:16 p.m.

Eventually, he let me pick him up and take him back inside. Naturally, he picked up a few fleas. Naturally, they have multiplied.

Adult Drivers Texas on May 15, 2013 at 8:45 p.m.

I conceive feeling to have the accumulation that you are card.

Free Samples Online on May 16, 2013 at 3:29 p.m.

Wonderful site you got here! Please maintain updating, I will def read more. Itll be in my bookmarks so far better update! Free Samples Online

mechanic schools hq on May 16, 2013 at 4:29 p.m.

Explanation why don’t you need a whole lot more in maximize your web page, not wearing running shoes can be safer to fall all-around your blog looking car. Nice telling these sorts of motivating knowledge in you. Important net, good job.!
mechanic schools info | mechanic schools info | mechanic schools guide | mechanic schools guide | mechanic schools blog

r4 pour ds on May 18, 2013 at 6:01 a.m.

I conceive feeling to have the accumulation that you are card.

coach outlet on May 19, 2013 at 9:26 p.m.

A blog is something that can Coach Factory only benefit you and not hurt you, most of Coach Factory Online air max the time that is. Once you establish a presence online via louis vuitton outlet blog you then create more potential followers for you and your business that you didn’t previously have. You see there is much to gain when you create a blog, so read through this article and see how blogging can help you.

Choose a domain name that immediately tells potential readers what your blog is about. coach outlet It’s not likely that you louboutin are going to ralph lauren be able to procure a name like but, your blog is more likely to be about some particular aspect Hogan of your subject. Incorporate that aspect into the domain name along with your louis vuitton Coach Factory Outlet overall focus.

Don’t be afraid to stretch out your louis vuitton hand and ask your reader for a donation. Your loyal readers, in particular, will be likely to donate some to coach outlet your cause. If your blog is valuable enough, people will realize it. They will also realize that, longchamp not Coach Factory louis vuitton outlet Online only does Coach Outlet it cost you money to produce your blog, your time is valuable.

Make sure that you have a blog mailing list started early. The sooner you get this started, the more time you coach outlet will have gucci to make that list larger. Once your blog is more established, this list will be used to bring in money, and you will be thankful that you already took care of this.

R4i on May 31, 2013 at 6:03 a.m.

I conceive feeling to have the accumulation that you are card.

springbed murah on June 12, 2013 at 12:15 p.m.

You're so cool! I don't think I've read anything like this before. So good to find somebody with some original thoughts on this subject. Thanks for starting this up. tava tea

lyza on June 13, 2013 at 5:52 a.m.

That is a wonderful post. Thank you so much for sharing. I will surely check other information posted on this site.
marble tiles

gregory on June 14, 2013 at 12:01 p.m.

Impressive articles and information are posted on this site. I'm going to share this to my friends and colleague. Sarasota Real Estate Information

TDS meter on June 22, 2013 at 7:54 a.m.

Hey, people utilized to write great posts, however the previous several blogposts happen to be somewhat boring... My partner and i miss your current great writing. Earlier several blogposts are merely a bit beyond observe!

malima on June 25, 2013 at 3:38 p.m.

Many people would ask you to correct the "balance" of the piece by correcting the factual errors it contains. I am not of this view. If the Republican party wishes to use dishonest portrayals of events in its campaigns, it is welcome to do so. wv-accident-lawyer.com
wv-car-accident-lawyer.com
wv-employment-lawyer.com
wv-medical-negligence-lawyer.com
wv-personal-injury-lawyer.com
wv-truck-accident-lawyer.com
yellowcatstone.com
yourelectronicparalegal.com
yourocklegal.com
zeusportclub.com
zeusportclubs.com

zhului on July 13, 2013 at 9:09 a.m.

The value of the asymptote should be the actual number of http://www.test-insides.com/000-780-exam.html fleas on Buzz at the time, and reaching that number would imply local extinction for the fleas. Of course, there are probably other fleas lurking about that would recolonize Buzz. In principle, testinside 000-780 if I were to repeat the exercise frequently enough, Buzz would become a sink for fleas, and their migration to his fur would gradually deplete them from the environment.

coque galaxy s3 on July 22, 2013 at 6:55 a.m.

If we stretch our credulity far enough to imagine this as a continuous function, we can express it as a differential equation.

coquehtc on August 01, 2013 at 12:24 p.m.

coquehtc http://www.coquehtc.org/

S. Nibhaana Saathiya on August 16, 2013 at 12:07 a.m.

I read your post and like it so much.Such a great information here.

jtry on September 18, 2013 at 5:53 a.m.

Anti-Aging Therapy is the combination of medical science, Buy Steroids Online biotechnology and life conditioning protocols, one of the most innovative treatments available today is Hormone Replacement Therapy (HRT), Nandrolone Phenylpropionate this is the process by which vital hormones are replaced in a prescribed schedule to Nandrolone Decanoate mimic the body’s natural process and at more youthful levels.

Midsummer on September 30, 2013 at 6:52 a.m.

This web site is really a walk-through for all of the info you wanted about this and didn’t know who to ask. Glimpse here, and you’ll definitely discover it…. Midsummer Midsummer Day

fitflop on October 02, 2013 at 6:48 a.m.

This type of message always inspiring and I prefer to read quality content, so happy to find good place to many here in the post, the writing is just great, thanks for the post.

Thomas on October 02, 2013 at 8:19 p.m.

I have not much time to visit many websites. But today I found this site accidently. I explored much information which is useful to my life. Thanks! futbol de risa winning minesweeper gillette soccer plaisir betting secrets japanese soccer begin attacker sports en bets that guarantee soccer on internet futbol betting tips soccer wagering trends

laonodeh on October 02, 2013 at 8:47 p.m.

You may post on the professional credentials for the blog owner. You could express it's outstanding. Your blog experience can springboard your click through. futbol de risa
sports betting secrets
sports betting secrets
park league football
staked nfl solutions
sports betting secrets
soccer leagues
free predictions online

crocs murah on November 04, 2013 at 10:23 a.m.

Great info and right to the point. I don't know if this is really the best place to ask but do you folks have any ideea where to employ some professional writers? Thx

Electrolizer on November 07, 2013 at 7:06 a.m.

Outstanding Publish. I have not been viewing the website lately. Took a visit again and there were some terrific feedback on the website. Outstanding post. Keep up the great work.

Florist Jogja on November 14, 2013 at 10:49 a.m.

This excellent article assisted me very much! Saved your site, extremely great categories everywhere that I read here! I really appreciate the info.

miami on November 19, 2013 at 1:24 a.m.

Impressive information you have here. I am glad to find this page, I will certainly check this site regularly. Read daily miami real estate news article directory.

Florist Jogja on November 19, 2013 at 11:23 a.m.

Thank you for the excellent content I did experienced studying it, I will be sure to save your website and definitely will come returning from again. I want to motivate that you proceed your excellent job, have a excellent day.

 Name: Email: optional; will not be displayed URL: optional Comment: URLs auto-link and some tags are allowed:

.