Message Boards Message Boards

Help needed --- Every Japanese to be called SATO by 2531"

Hello!

I'm a teacher of journalism at the international EDJ school in Nice FR

My students have access to Mathematica 14 (but they mainly use chat-driven notebooks as the wolfram language is sort of complicated to them).

Anyway. As part of the "Data Journalism" course I want to give them the assignment to read this article https://www.dailymail.co.uk/news/article-13266707/Everyone-Japan-called-Sato-year-2531-countrys-marriage-laws.html

It tells the story of a Tohoku University economics professor that concluded that, given how married couples in jpn drop one of their family names in favour of both individuals taking the same one everyone will be called "SATO" by 2531.

I'd like the students to perform the same simulation for their own country (after finding out the candidate family name -- for Italy that would probably be "Rossi")

Now, after some digging I've found in this PDF https://think-name.jp/assets/pdf/Sato_estimation_yoshida_hiroshi.pdf

the method used to perform the calculation for Japan:

Handling of Past Data

First, we obtained the number of people with the Sato surname in Japan from the data provided and published by "Myoji-yurai.net" (https://myoji-yurai.net/), which covers more than 99.04% of Japanese surnames.

Next, we divided the number of people with the Sato surname by the total population of Japan for each year (estimated by the Ministry of Internal Affairs and Communications) * 99.04% to obtain the "ratio of the Sato surname in a given year t": x(t).

From the change in the Sato surname between the latest years of 2022 and 2023, we calculated the one-year growth rate [Rho] in the Sato surname ratio. Estimation Results (1) Growth Rate [Rho] of the Sato Surname It is found that the ratio of the Sato surname x(t) increased from 1.480% in 2013 to 1.530% in 2023, an increase of 0.05 percentage points over more than 10 years. Calculating from the data for the most recent points of 2022 and 2023, the growth rate [Rho] of the Sato surname ratio is (1+[Rho]) = 1.0083.

(2) Future Simulation Assuming that the ratio of the Sato surname to the Japanese population will grow at a rate of 1.0083 each year, starting from 1.530% as of March 2023, and repeating the calculation of x(t+1) = (1+[Rho]) x(t), it was calculated that the ratio will reach 100% in approximately 500 years, in 2531.


And the question is : Is anyone kind enough to write me the Wolfram Language generic code to perform the same ? Ideally, students could get their own country basic data (if possible directly with WolfrmaLanguage functions, if not possible from outside sources) and have the system perform the simulation for them (of course, forcing the Japanese role "select one family name" even if not the case for the specific country)

It is also possible that it will not converge and we will never have a "Rossi only" country not even in year 3500...or maybe yes. But that's also the point of the simulation.

Thanks!!

POSTED BY: Marco Barsotti
4 Replies

Very cool answer! I can add nothing more.

I ran the code, slightly modified such that the upper limit 30 of $n$ is replaced by 1000. I put some sample outputs here. The last one is prophesied future of Japan !? enter image description here enter image description here enter image description here

POSTED BY: Akishi Kato

Just out of curiosity, I wrote a stochastic simulation. I start with a few surnames, give them a random frequency, and build a starting population with those surnames with the given frequencies. The new generation is generated this way: first I make a random permutation of the previous population, then I take half of the remixed population, and finally I redouble this half: this simulates married couples that homogenize their surnames:

Clear[initialFrequencies, startingPopulation, population];
surnames = {"Sato", "Suzuki", "Takahashi",
   "Tanaka", "Watanabe", "Ito"};
initialFrequencies = 
  1000*(Join[{0, 1}, RandomReal[{0, 1}, Length[surnames] - 1]] //
         Sort // Differences) // Sort // Reverse // Round;
startingPopulation = 
  Thread[Inactive[ConstantArray][surnames, initialFrequencies]] //
    Activate // Flatten;
population[0] = startingPopulation;
population[n_] := 
  population[n] = Module[{randomHalfPopulation}, randomHalfPopulation =
     RandomSample[population[n - 1]][[1 ;; -1 ;; 2]]; 
    Join[randomHalfPopulation, randomHalfPopulation]];
ListLinePlot[
 Transpose[
  surnames /. Table[Association @@ (Rule @@@ Tally[population[n]]),
    {n, 0, 30}]],
 PlotLegends -> surnames]

Running the simulation a few times does not show a clear pattern. Sometimes Sato increases, sometimes it decreases. The least frequent surnames are at risk of extinction.

POSTED BY: Gianluca Gorni

hello!

thank you for your kind answer. I suspected that too.

A running wolfram language code to simulate the evolution would be great in fact.

thanks again

POSTED BY: Marco Barsotti

Hi. I am Japanese and my surname is Kato, which is ranked as 10th popular surname.

Perhaps I could help in coding, but the model you cited (geometric series, or, compound interest calculation) is too simplistic to be real. If most people had surname Sato, the percentage of Sato might still increase, but the rate of growth would inevitably slow down. The model totally ignores this fact and percentage will go up over 100%!

Mathematically speaking, such hegemony / supremacy / dominance phenomena are best described by Logistic differential equations. This equation correctly incorporate such slowing-down effect due to saturation.

I don't know why the author of the paper used such oversimplified model. My guess is that his real intention is not in the academic correctness, but to spark a debate on selective separate surnames for married couples, which is a hot political issues in Japan.

POSTED BY: Akishi Kato
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract