#mathematical modeling

[cauchyfeelings]

My future husband doesn’t need to worry about getting me a ring he needs to worry about getting me the maxed out Mac Studio so I can run all my analyses and mathematical/computational/theoretical models seamlessly.

[mathprofessorquotes]

Rabbits aren’t immortal.

Mathematical modelling lecturer

[bpod-bpod]

Use the Force

For a tumour to spread and grow, its cells must pull against the surrounding microenvironment of supporting molecules like elastic collagen. This jostling for space generates tension inside the cell – a struggle that drugs might exploit to limit cancerous growth. To investigate these forces, here researchers swap real tumours for lab-grown tumour ‘spheroids’ (black) which strain against the surrounding collagen under a microscope, creating starburst patterns of ‘deformations’ over a 12 hour period (top left to bottom right). Mathematical models simulate these patterns with the aim of predicting the contractility of different spheroids – how much their internal structure contracts when straining. Similar techniques could be used to predict the response of tumours growing in different tissue environments and provide a model for testing drugs aiming to push back on cancer.

Written by John Ankers

Image from work by Christophe Mark, Thomas J Grundy, Pamela L Strissel and David Böhringer, and colleagues

Department of Physics, Friedrich-Alexander University Erlangen-Nürnberg, Erlangen, Germany

Image originally published with a Creative Commons Attribution 4.0 International (CC BY 4.0)

Published in eLife, April 2020

You can also follow BPoD on Instagram, Twitter and Facebook

[luxxeonthings]

Linked Star Ball For 3d Printing | Blender 2.8 | 3d Tutorial

This tutorial will feature some of the latest hotkeys and shortcuts in the new Blender 2.8 rc2 or higher, and we model this starball shape using real-world units (millimeters) which could then be exported for 3d printing.

[rebrobindoesmath]

How I Attempted to Mathematically Model the Spread of Memes on the Internet, pt 3 (The Final Chapter)

Welcome back! If you have not read the earlier parts of this series, you can find them here:

Part 1 (Intro to Markov Chains)

Part 2 (Setting Up The Model)

This will be the last part of this series and will primarily involve ~actually~ using the model. Finally, amirite?

Things I’m Going To Touch On:

More About I(t)

Fitting I(t) To Data

Predicting Probabilities

For the Future

Here we go!

More About I(t)

First, let’s recall the equations we obtained in the last post. 

We can figure out the long term behavior of this I(t) by taking the limit as t approaches $\infty$. Let’s note that since $\alpha$, $\beta$ and $\gamma$ are less than 1, the exponential terms will go to zero when taking the limit. So, 

What does this mean in context?

This $\frac{\gamma}{\beta+\gamma}$ tells us the probability of someone being in the infected state in the long term (that is, after the initial “outbreak” of the meme). We can also verify this from our steady state distribution, but I’ll leave that as an exercise. 

Additionally, I think it’s interesting to note that this long term value does not depend whatsoever on $\alpha$. This means that the probability for moving from susceptible to infected has no effect on the long term probability of being in the infected state. 

Fitting I(t) to Data

I obtained my data from Google Trends. Here’s what I did:

First, I looked up a particular meme as a search term on Google Trends and got a spreadsheet of popularity values. Google plots the data so that the search term’s peak popularity is given a value of 100 and then the rest is scaled in relation to that value. Since I’m concerned with probabilities, I can’t use the values given by Google. So, I scaled the values down so that they were between 0 and 1.

Now that I have my adjusted values, I want to try and fit my I(t) equation to them. I did this by performing regression analysis using SageMath to find the $\alpha$, $\beta$, and $\gamma$ values that best fit I(t) to the data (Will provide my code upon request!). 

Below is the Google Trends graph for “dat boi”. I do want to note that each data point obtained from this graph is for a one week time span. 

And here’s a graph of the best fit I(t) (*drumroll*):

This is awesome! Look how closely I(t) matches the actual data. In this case, the $\alpha$, $\beta$, and $\gamma$ that best fit I(t) are:

$\alpha \approx$ 0.5354

$\beta \approx$ = 0.3308

$\gamma \approx$ = 0.03285

We can also see that $\frac{\gamma}{\beta+\gamma} = 0.09$.

I performed this analysis with 10 different memes/search terms. I’ll present the results in the table below, just listing the $\alpha$, $\beta$, and $\gamma$ values. 

Quick note: those terms with an asterisk had data points that were taken daily, those with a cross were monthly, and those with nothing were taken weekly.

I included “Flint, Michigan” because although it isn’t a meme, it does follow the same viral pattern that the memes do. This is due to the increased media coverage due to the Flint Water Crisis. I do want to stress that I am in no way trying to diminish the significance of the Flint Water Crisis or the struggle that those involved are facing, but rather I want to show that this model can be applied to more things than just memes. 

It’s neat that we can find these probabilities with memes that have already died down from their “outbreaks,” but is there a way we can try and predict the $\alpha$, $\beta$, and $\gamma$ values for up and coming memes? This brings us to the next section of this post. 

Predicting Probabilities

I want to try to predict the $\alpha$, $\beta$, and $\gamma$ values for a particular meme given an initial set of data points. Let’s assume that we have the data points up to and including the peak popularity. We’ll talk about why we need to assume this in a little bit. 

Let’s start by attempting to predict $\alpha$. Briefly recall that $\alpha$ is the probability of moving from the susceptible state to the infected state. I know that $\alpha$ has to do with how quickly the meme rises in popularity. Because of this, I approximate $\alpha$ using slopes in the following process:

Take the set of points up to and including the peak point.

Find all of the slopes between consecutive pairs of points.

Take the maximum over all these slopes. This is the estimated $\alpha$. 

This process is why we have to assume the peak is included in our data set. If we didn’t have the peak, we wouldn’t know if there was a slope that would be higher than the maximum we found over the points we have. 

Using this method over all 10 examples gives the following results for $\alpha$.

This seems promising at first, but with a closer look, the differences for some of the $\alpha$ values are pretty large. Is there a way to make it an even better approximation?

Spoiler alert: yes. 

Using SPSS, I performed various regressions with this estimated $\alpha$ as the independent variable and the actual $\alpha$ value obtained from the best fit I(t) as the dependent variable. See the table below:

Let’s quickly note that for every regression performed, the p-value is < 0.05, and so the relationship between $\alpha$ and our approximation of $\alpha$ is statistically significant (via F-Test of Overall Significance). We want to pick the equation that gives us the highest R-Square value. The quadratic and the cubic equations have the same R-Square value, so I’ll pick the quadratic just so I have one less term to deal with. 

Thus,

where $\alpha_{est}$ is that largest slope value. 

Using this new way to estimate $\alpha$, we find our new approximate $\alpha$ values by plugging in the maximum slope the equation above:

We can see that the differences between the estimated and the actual $\alpha$ values are much smaller than in the previous estimate! The bold values are those that have decreased significantly.

Now we have a pretty nice approximation for $\alpha$!

Can we do the same for $\beta$?

Let’s try! Since $\beta$ is the probability of moving from the infected state to the recovered state, we can’t really use the slopes as an approximation. Instead, we’ll perform a regression analysis with $\alpha$ as the independent variable and $\beta$ as the dependent variable to see if there’s a connection. 

In this case, we will pick the cubic equation. Thus,

Let’s look at how this approximation compares to the actual $\beta$ values obtained from this model. 

The average difference between the expected and the actual $\beta$ with this approximation is 0.06214. I believe it’s safe to say that provided we have an accurate $\alpha$ value, we can obtain an estimated $\beta$ value that is somewhat close to the actual $\beta$. 

Of course, there is probably a method that will better predict $\beta$, but for now, this is the best we’ve got.

What about $\gamma$?

Well, here’s where we have problems. Unfortunately, there is no significant relationship between $\alpha$ and $\gamma$ or $\beta$ and $\gamma$, as determined by performing various regressions. Because of this, we can’t use the same methods of approximation that we used for $\alpha$ and $\beta$. This brings me to the final section of this post. 

For The Future

We’ve done a lot with this already, but there is always more to be done! There are several things that my model does not account for:

If a meme has two peaks instead of just one (i.e. it died and then became popular again)

Extremely sharp increase followed by an extremely sharp decrease (say, in the matter of a couple days)

If $\beta + \gamma = \alpha$ or if $\beta + \gamma > 1$. This is because of issues with our I(t) equation. $\beta + \gamma  = \alpha$ gives a denominator of zero and $\beta + \gamma > 1$ makes $1-\beta-\gamma$ negative, which causes issues when raising it to t. 

Problem 3 is an issue because initially our Markov chain has no restrictions on what $\alpha$, $\beta$, or $\gamma$ could be, yet when I come up with the I(t) equation, we somehow have these problems. It’d be nice to have a way to circumvent this issue! 

Additionally, we still cannot predict what $\gamma$ will be, and until we can do so, we can’t create an equation given just an initial data set. Finding a way to $\gamma$ is the first thing on the list for future work. 

These are all things to consider more in the future, and I’d (of course) love to hear any ideas on how to combat these issues. 

Thanks for reading!

As always, if there are any questions/comments/concerns or even suggestions for another post, please feel free to send an ask!

Hope all is well and (as always) stay positive!

[neurosciencenews]

Mathematical Modeling Uncovers Mysteries of HIV Infection in the Brain

After uncovering the progression of HIV infection in the brain thanks to a new mathematical model developed by a UAlberta research team, clinicians and researchers are developing a nasal spray to administer drugs more effectively.

The research is in Journal of Neurovirology. (full access paywall)

[cauchyfeelings]

while alive:

study_math() think_existentially() write_theory() suffer_in_code() repeat()

[cauchyfeelings]

When I pursued math it was in the pursuit of truth scrutinized by rigor. There is just something I find so beautiful when it remains true against all odds. When it connects seamlessly across each of its subdomains. The invariances essentially. And tbh I fell in love with it. The truth I craved from life.

Coding is nice because it provided a way to simulate a phenomenon and visualize something that would remain too abstract for public understanding (and maybe my own understanding as well). Building something that only existed in your mind is honestly pretty cool (despite how shitty and ugly it might be lmao). And diving into how it’s built to understand that process itself? The meta? Also very cool.

But at the end of the day these are all just perspectives. I climbed these mountains because I wanted to see the world. And the world looks so beautiful when you look at it from different perspectives. However there comes a point when you must ask yourself, am I losing sight of the world because all I see now is this particular perspective? And by all means if the perspective itself is your world then that’s absolutely wonderful and I’m genuinely happy for you. But for me I feel like I’ve found myself further and further away from my world. My home. The more I go in pursuit of trying to understand my home the more I feel like I’ve strayed further from that home. And it kinda hurts because you feel lost and homeless in a sense. If you had just stayed home in the first place maybe you wouldn’t have felt so lost. But at the same time how can you regret not having stayed? All the experiences and knowledge you’ve gained was inspired by your home. All you wanted was to understand your home, the world, in every way possible. And it made you want to better yourself, expand yourself, so you can see it clearly through all perspectives. So you can understand it as much as possible. So how can you regret not having stayed?

But also it’s like what The Alchemist taught us. Maybe the treasure has been home all this time. You just needed to cross your own mountains and deserts to know and come back home. To what remains true to you. To what is invariant.