Is Latitude Associated with Chronotype?

Daniel Vartanian

University of São Paulo

2025-02-03


Hi there! 👋

This presentation will provide an overview of the thesis objectives, main concepts, methods, and results. Together, we’ll uncover how environment factors can influence the cronotype.

We will cover the following topics:

  1. Introduction
  2. Methods
  3. Results
  4. Discussion
  5. Conclusion

Introduction

Question

This thesis focuses on the expression of sleep and circadian phenotypes (chronotypes), seeking to provide an answer to the following question related to human populations:

Is latitude associated with chronotype?

Chronotype

Entrainment

The Latitude Hypothesis

Leocadio-Miguel et al. 

Sample: \(12,884\) Brazilian participants

Method: HO Questionnaire (Horne & Östberg, 1976)

\(\Delta \ \text{Adjusted} \ \text{R}^{2} = 0.388\%\) (Cohen’s \(f^2 = 0.004143174\))

Methods

MCTQ

The Two Process of Sleep Reg.

Sample

Sample

The dataset used for analysis is made up of \(65,824\) Brazilian individuals aged 18 or older, residing in the UTC-3 timezone, who completed the survey between October 15th and 21st, 2017.

Hypothesis Test

The study tested whether latitude significantly improves model fit when predicting chronotype with nested models.

A restricted model (excluding latitude) with a full model (including latitude) was compared.

To ensure practical significance, a Minimum Effect Size (MES) criterion was applied, following the original Neyman-Pearson framework for hypothesis testing (Neyman & Pearson, 1928a, 1928b; Perezgonzalez, 2015).

The MES was set at Cohen’s threshold for small effects (\(f^2 = 0.02\), equivalent to \(\text{R}^2 = 0.01960784\)) (Cohen, 1988, p. 413). Thus, latitude was considered meaningful only if its inclusion accounted for at least \(1.960784\%\) of the variance in chronotype.

To balance the sample, a weighting procedure was applied to the data. The weights were calculated by cell weighting, using the sex, age group and Brazil’s state as reference.

Hypothesis Test

The hypothesis can be outlined as follows:

Null hypothesis (\(\text{H}_{0}\))
Adding latitude does not meaningfully improve the model’s fit. This implies that the change in adjusted \(\text{R}^{2}\) is negligible, or the F-test is not significant, given a type I error probability (\(\alpha\)) of \(0.05\).
Alternative Hypothesis (\(\text{H}_{a}\))
Adding latitude meaningfully improves the model’s fit. This implies that the change in adjusted \(\text{R}^{2}\) exceeds the expected Minimum Effect Size (MES), and the F-test is significant, given a type I error probability (\(\alpha\)) of \(0.05\).

\[ \begin{cases} \text{H}_{0}: \Delta \ \text{Adjusted} \ \text{R}^{2} \leq \text{MES} \quad \text{or} \quad \text{F-test is not significant} \ (\alpha \geq 0.05) \\ \text{H}_{a}: \Delta \ \text{Adjusted} \ \text{R}^{2} > \text{MES} \quad \text{and} \quad \text{F-test is significant} \ (\alpha < 0.05) \end{cases} \]

Results

Chronotype Distribution

Hypothesis Test

Restricted Model
(Based on Leocadio-Miguel et al. (2017))

  • Predictors (\(4\)):
    Age, Sex, Longitude, and the Average Monthly Global Horizontal Irradiance (GHI) at the time of questionnaire completion.
  • $^2_{} = $
  • \(\text{F}(4, 65818) = 1531.808\)
  • \(p\text{-value} < 1e-05\)

Hypothesis Test

Full Model A
(Based on Leocadio-Miguel et al. (2017))

  • Predictors (\(8\)):
    [restricted predictors] + Average Annual Global Horizontal Irradiance (GHI), Daylight Duration for the Nearest March Equinox, June and December solstices
  • $^2_{} = $
  • \(\text{F}(8, 65814) = 794.12\)
  • \(p\text{-value} < 1e-05\)

Full Model B

  • Predictors (\(5\)):
    [restricted predictors] + Latitude
  • $^2_{} = $
  • \(\text{F}(5, 65817) = 1233.588\)
  • \(p\text{-value} < 1e-05\)

Hypothesis Test

All coefficients were statistically different from zero (\(p\text{-value} < 0.05\)). Assumption checking and residual diagnostics primarily relied on visual inspection, as formal assumption tests (e.g., Anderson-Darling) are often not recommended for large samples (Shatz, 2024). All validity assumptions were met, and no serious multicollinearity was found among the predictor variables.

An ANOVA for nested models revealed a significant reduction in the residual sum of squares in both tests (A \(\text{F}(4, 65814) = 51.71\), \(p\text{-value} < 1e-05\)) (B \(\text{F}(1, 65817) = 37.325\), \(p\text{-value} < 1e-05\)).


However, similarly to Leocadio-Miguel et al. (2017), when estimating Cohen’s \(f^2\) effect size (\(f^2 = 0.004143174\)), the results were below the MES (i.e., negligible) (A \(f^2 = 0.0031428, 95\% \ \text{CI}[0, 0.012203]\)) (B \(f^2 = 0.0005671, 95\% \ \text{CI}[0, 0.0095426]\)).

Discussion

Chronotype versus Latitude

Chronotype Geo. Distribution

Chronotype Geo. Distribution

The p-value Problem

Large samples and sensitivity

Is a difference of \(0.00001\) valid?

Statistical ritual versus Statistical thinking

Confirmation bias


Comparison of a 95% of confidence level (\(\alpha = 0.05\)) and an n-dependent p-value curve. The parameter \(n_{\alpha}\) represents the minimum sample size to detect statistically significant differences among compared groups. The parameter \(n_{\gamma}\) represents the convergence point of the p-value curve. When the p-value curve expresses practical differences, the area under the red curve (\(A_{p(n)}\)) is smaller than the area under the constant function \(\alpha = 0.05\) (\(A_{\alpha = 0.05}\)) when it is evaluated between \(0\) and \(n_{\gamma}\).

The p-value Problem

SMALL EFFECT SIZE: \(f^2 = .02\). Translated into \(\text{R}^{2}\) or partial \(\text{R}^{2}\) for Case 1, this gives \(.02 / (1 + .02) = .0196\). We thus define a small effect as one that accounts for 2% of the \(\text{Y}\) variance (in contrast with 1% for \(r\)), and translate to an \(\text{R} = \sqrt{0196} = .14\) (compared to .10 for \(r\)). This is a modest enough amount, just barely escaping triviality and (alas!) all too frequently in practice represents the true order of magnitude of the effect being tested. (p. 413)

[…] in many circumstances, all that is intended by “proving” the null hypothesis is that the ES [Effect Size] is not necessarily zero but small enough to be negligible […]. (p. 461)

Nonrealistic Models

Conclusion

Answer

I suggest that it is the aim of science to find satisfactory explanations, of whatever strikes us as being in need of explanation.

This study, using what is arguably one of the largest datasets on chronotype, found no evidence supporting the latitude hypothesis in humans.

Closing remarks

License: MIT License: CC BY 4.0

This thesis was presented to the School of Arts, Sciences and Humanities (EACH) at the University of São Paulo (USP) as a requirement for the degree of Master of Science by the Graduate Program in Complex Systems Modeling (PPGSCX).

I am deeply grateful to my advisor, Prof. Dr. Camilo Rodrigues Neto, for introducing me to complexity science in 2012, guiding my dissertation with patience, and demonstrating remarkable integrity in navigating a challenging supervisory transition.

Closing remarks

License: MIT License: CC BY 4.0

The presentation was created using the Quarto publishing system. All analyses presented are fully reproducible and were conducted using the R programming language.

To explore the code and repository for this thesis, click here. The research compendium is also accessible via The Open Science Framework by clicking here.

Financial support was provided by the Coordination for the Improvement of Higher Education Personnel (CAPES) (Grant number: 88887.703720/2022-00).

References

In accordance with the American Psychological Association (APA) Style, 7th edition.

Borbély, A. A. (1982). A two process model of sleep regulation. Human Neurobiology, 1(3), 195–204. https://pubmed.ncbi.nlm.nih.gov/7185792
Cohen, J. (1988). Statistical power analysis for the behavioral sciences (2nd ed.). Lawrence Erlbaum Associates.
Cohen, J. (1992). A power primer. Psychological Bulletin, 112(1), 155–159. https://doi.org/10.1037/0033-2909.112.1.155
Comte, A. (1892). Cours de philosophie positive: Tome premier: Les préliminaires généraux et la philosophie mathématique [Course of positive philosophy: Volume one: General preliminaries and mathematical philosophy] (5th ed., Vol. 1). http://archive.org/details/coursdephilosoph01comt
Duffy, J. F., Cain, S. W., Chang, A.-M., Phillips, A. J. K., Münch, M. Y., Gronfier, C., Wyatt, J. K., Dijk, D.-J., Wright, K. P., & Czeisler, C. A. (2011). Sex difference in the near-24-hour intrinsic period of the human circadian timing system. Proceedings of the National Academy of Sciences, 108, 15602–15608. https://doi.org/10.1073/pnas.1010666108
Epstein, J. M. (1999). Agent-based computational models and generative social science. Complexity, 4(5), 41–60. https://doi.org/10.1002/(SICI)1099-0526(199905/06)4:5<41::AID-CPLX9>3.0.CO;2-F
Flanagan, A., Bechtold, D. A., Pot, G. K., & Johnston, J. D. (2021). Chrono-nutrition: From molecular and neuronal mechanisms to human epidemiology and timed feeding patterns. Journal of Neurochemistry, 157(1), 53–72. https://doi.org/10.1111/jnc.15246
Gómez-de-Mariscal, E., Guerrero, V., Sneider, A., Jayatilaka, H., Phillip, J. M., Wirtz, D., & Muñoz-Barrutia, A. (2021). Use of the p-values as a size-dependent function to address practical differences when analyzing large datasets. Scientific Reports, 11(1, 1), 20942. https://doi.org/10.1038/s41598-021-00199-5
Holland, J. H. (2012). Signals and boundaries: Building blocks for complex adaptive systems. MIT Press. https://doi.org/10.7551/mitpress/9412.001.0001
Holland, J. H. (2014). Complexity: A very short introduction (OCLC: ocn888011172). Oxford University Press.
Horne, J. A., & Östberg, O. (1976). A self-assessment questionnaire to determine morningness-eveningness in human circadian rhythms. International Journal of Chronobiology, 4(2), 97–110.
Krakauer, D. (2023, February 17). What is complexity? [Video recording]. Santa Fe Institute. https://www.youtube.com/watch?v=JR93X7xK05o
Kuhlman, S. J., Craig, L. M., & Duffy, J. F. (2018). Introduction to chronobiology. Cold Spring Harbor Perspectives in Biology, 10(9), a033613. https://doi.org/10.1101/cshperspect.a033613
Leocadio-Miguel, M. A., Louzada, F. M., Duarte, L. L., Areas, R. P., Alam, M., Freire, M. V., Fontenele-Araujo, J., Menna-Barreto, L., & Pedrazzoli, M. (2017). Latitudinal cline of chronotype. Scientific Reports, 7(1), 5437. https://doi.org/10.1038/s41598-017-05797-w
Lewin, R. (1993). Complexity: Life at the edge of chaos. Collier Books.
Merrow, M., & Roenneberg, T. (2020). Circadian clocks: How rhythms structure life [Online course]. https://www.coursera.org/learn/circadian-clocks
Neyman, J., & Pearson, E. S. (1928a). On the use and interpretation of certain test criteria for purposes of statistical inference: Part I. Biometrika, 20A(1/2), 175–240. https://doi.org/10.2307/2331945
Neyman, J., & Pearson, E. S. (1928b). On the use and interpretation of certain test criteria for purposes of statistical inference: Part II. Biometrika, 20A(3/4), 263–294. https://doi.org/10.2307/2332112
Nicolis, G., & Nicolis, C. (2012). Foundations of complex systems: Emergence, information and prediction (2nd ed.). World Scientific.
Nobel Prize Outreach AB. (n.d.). Press release. The Nobel Prize. Retrieved September 28, 2023, from https://www.nobelprize.org/prizes/medicine/2017/press-release/
Pereira, E. B., Martins, F. R., Gonçalves, A. R., Costa, R., Lima, F. J. L., Rüther, R., Abreu, S. L., Tiepolo, G. M., Pereira, S. V., & Souza, J. G. (2017). Atlas brasileiro de energia solar [Brazilian atlas of solar energy] (2nd ed.). Instituto Nacional de Pesquisas Espaciais. https://doi.org/10.34024/978851700089
Perezgonzalez, J. D. (2015). Fisher, Neyman-Pearson or NHST? A tutorial for teaching data testing. Frontiers in Psychology, 6. https://doi.org/10.3389/fpsyg.2015.00223
Popper, K. R. (1979). Objective knowledge: An evolutionary approach. Oxford University Press. (Original work published 1972)
Popper, K. R. (2002). Conjectures and refutations: The growth of scientific knowledge. Routledge. (Original work published 1963)
Popper, K. R. (2005). The logic of scientific discovery. Routledge. (Original work published 1934)
Popper, K. R. (2005). Unended quest: An intellectual autobiography. Taylor and Francis. (Original work published 1974)
Rede Globo. (2017, October 15). Metade da população se sente mal no horário de verão, revela pesquisa (TV program Fantástico) [Half of the population feels unwell during daylight saving time, research reveals] [Video recording]. Rede Globo. https://g1.globo.com/fantastico/noticia/2017/10/metade-da-populacao-se-sente-mal-no-horario-de-verao-revela-pesquisa.html
Roenneberg, T., Kuehnle, T., Juda, M., Kantermann, T., Allebrandt, K., Gordijn, M., & Merrow, M. (2007). Epidemiology of the human circadian clock. Sleep Medicine Reviews, 11(6), 429–438. https://doi.org/10.1016/j.smrv.2007.07.005
Roenneberg, T., Wirz-Justice, A., & Merrow, M. (2003). Life between clocks: Daily temporal patterns of human chronotypes (.). Journal of Biological Rhythms, 18(1), 80–90. https://doi.org/10.1177/0748730402239679
Roenneberg, T., Wirz-Justice, A., Skene, D. J., Ancoli-Israel, S., Wright, K. P., Dijk, D.-J., Zee, P., Gorman, M. R., Winnebeck, E. C., & Klerman, E. B. (2019). Why should we abolish daylight saving time? Journal of Biological Rhythms, 34(3), 227–230. https://doi.org/10.1177/0748730419854197
Shatz, I. (2024). Assumption-checking rather than (just) testing: The importance of visualization and effect size in statistical diagnostics. Behavior Research Methods, 56(2), 826–845. https://doi.org/10.3758/s13428-023-02072-x
The Worldwide Experimental Platform. (n.d.). MCTQ. The Worldwide Experimental Platform. Retrieved July 17, 2023, from https://www.thewep.org/documentations/mctq
Wilson, R. (2004, October 2). Emergence isn’t magic. New Scientist, 2467. https://www.newscientist.com/letter/mg18424674-700-emergence-isnt-magic/

Thank you!

Appendices

(AP) Mimosa Pudica

(AP) Entrainment

(AP) Emergence

Stable macroscopic patterns arising from local interaction of agents (Epstein, 1999).

When the aggregate exhibits properties not attained by summation (Holland, 2014).

(AP) Emergence

(AP) Emergence

An aggregate behavior emerges from the interactions of the parts (CAS) (Holland, 2012).

(AP) Earth’s Rotation Axis

(AP) Corr. Matrix

(AP) Chronotype Comparison

Brazil

Europe

(AP) Age versus Chronotype

Brazil

Europe

(AP) Chronotype Geo. Dist.

(AP) Chronotype Geo. Dist.

(AP) Cohen’s Benchmarks

(AP) Directions for Fut. Research

WorldClim 2.1 in NetLogo (NetLogo/Scala-Java).

(AP) Popper’s hypothetico-deductive method

The activity can be represented by a general schema of problem- solving by the method of imaginative conjectures and criticism, or, as I have often called it, by the method of conjecture and refutation. The schema (in its simplest form) is this:

flowchart LR
  A(P1) --> B(TT)
  B --> C(EE)
  C --> D(P2)

flowchart LR
  A(P) --> B(TT)
  B --> C(EE)
  C --> A

(AP) Popper’s hypothetico-deductive method

Here \(\text{P}_1\), is the problem from which we start, \(\text{TT}\) (the ‘tentative theory’) is the imaginative conjectural solution which we first reach, for example our first tentative interpretation. \(\text{EE}\) (‘error- elimination’) consists of a severe critical examination of our conjecture, our tentative interpretation: it consists, for example, of the critical use of documentary evidence and, if we have at this early stage more than one conjecture at our disposal, it will also consist of a critical discussion and comparative evaluation of the competing conjectures. \(\text{P}_2\) is the problem situation as it emerges from our first critical attempt to solve our problems. It leads up to our second attempt (and so on).

flowchart LR
  A(P1) --> B(TT)
  B --> C(EE)
  C --> D(P2)

(AP) Popper’s hypothetico-deductive method

The history of ideas teaches us very clearly that ideas emerge in logical or, if the term is preferred, in dialectical contexts. My various schemata such as

flowchart LR
  A(P1) --> B(TT)
  B --> C(EE)
  C --> D(P2)

may indeed be looked upon as improvements and rationalizations of the Hegelian dialectical schema: they are rationalizations because they operate entirely within the classical logical organon of rational criticism, which is based upon the so-called law of contradiction; that is to say, upon the demand that contradictions, whenever we discover them, must be eliminated. Critical error-elimination on the scientific level proceeds by way of a conscious search for contradictions.

(AP) Popper against positivism (or the problem of induction)

Notre activité intellectuelle est suffisamment excitée par le pur espoir de découvrir les lois des phénomènes, par le simple désir de confirmer ou d’infirmer une théorie. […] la philosophie positive est le véritable état définitif de l’intelligence humaine […] (Comte, 1892, pp. 9–10).


Why post-positivist?

  • Anti-inductivist (deduction, satisfactory answers)
  • Anti-verificationist (falsifiability as demarcation criterion)
  • Anti-historicist (rebuttable presumptions)

(AP) Popper against positivism (or the problem of induction)

But I shall certainly admit a system as empirical or scientific only if it is capable of being tested by experience. These considerations suggest that not the verifiability but the falsifiability of a system is to be taken as a criterion of demarcation. In other words: I shall not require of a scientific system that it shall be capable of being singled out, once and for all, in a positive sense; but I shall require that its logical form shall be such that it can be singled out, by means of empirical tests, in a negative sense: it must be possible for an empirical scientific system to be refuted by experience.

(AP) Popper against positivism (or the problem of induction)

Everybody knows nowadays that logical positivism is dead. But nobody seems to suspect that there may be a question to be asked here—the question “Who is responsible?” or, rather, the question “Who has done it?”. (Passmore’s excellent historical article does not raise this question.) I fear that I must admit responsibility.

(AP) The 7 conclusions of Popper on science

  1. It is easy to obtain confirmations, or verifications, for nearly every theory—if we look for confirmations.
  2. Confirmations should count only if they are the result of risky predictions; that is to say, if, unenlightened by the theory in question, we should have expected an event which was incompatible with the theory—an event which would have refuted the theory.

(AP) The 7 conclusions of Popper on science

  1. Every ‘good’ scientific theory is a prohibition: it forbids certain things to happen. The more a theory forbids, the better it is.
  2. A theory which is not refutable by any conceivable event is nonscientific. Irrefutability is not a virtue of a theory (as people often think) but a vice.

(AP) The 7 conclusions of Popper on science

  1. Every genuine test of a theory is an attempt to falsify it, or to refute it. Testability is falsifiability; but there are degrees of testability: some theories are more testable, more exposed to refutation, than others; they take, as it were, greater risks.
  2. Confirming evidence should not count except when it is the result of a genuine test of the theory; and this means that it can be presented as a serious but unsuccessful attempt to falsify the theory. (I now speak in such cases of ‘corroborating evidence’).

(AP) The 7 conclusions of Popper on science

  1. Some genuinely testable theories, when found to be false, are still upheld by their admirers—for example by introducing ad hoc some auxiliary assumption, or by re-interpreting the theory ad hoc in such a way that it escapes refutation. Such a procedure is always possible, but it rescues the theory from refutation only at the price of destroying, or at least lowering, its scientific status. (I later described such a rescuing operation as a ‘conventionalist twist’ or a ‘conventionalist stratagem’).

(AP) The 7 conclusions of Popper on science

One can sum up all this by saying that the criterion of the scientific status of a theory is its falsifiability, or refutability, or testability.