Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Time-course variation of statistics embedded in music: Corpus study on implicit learning and knowledge

  • Tatsuya Daikoku

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    daikoku@cbs.mpg.de

    Affiliation Department of Neuropsychology, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany

Abstract

Learning and knowledge of transitional probability in sequences like music, called statistical learning and knowledge, are considered implicit processes that occur without intention to learn and awareness of what one knows. This implicit statistical knowledge can be alternatively expressed via abstract medium such as musical melody, which suggests this knowledge is reflected in melodies written by a composer. This study investigates how statistics in music vary over a composer’s lifetime. Transitional probabilities of highest-pitch sequences in Ludwig van Beethoven’s Piano Sonata were calculated based on different hierarchical Markov models. Each interval pattern was ordered based on the sonata opus number. The transitional probabilities of sequential patterns that are musical universal in music gradually decreased, suggesting that time-course variations of statistics in music reflect time-course variations of a composer’s statistical knowledge. This study sheds new light on novel methodologies that may be able to evaluate the time-course variation of composer’s implicit knowledge using musical scores.

1. Introduction

1.1. Implicit learning and statistical learning

Implicit learning and knowledge is a ubiquitous phenomenon that occurs without the intention to learn and awareness of what we know [1,2]. A number of studies have investigated this form of learning using several paradigms [3,4, 5,6]. In a series of psychological and neurological studies, the learning and knowledge of transitional probability in sequential information such as music and language—called statistical learning and statistical knowledge—have been considered implicit processes that could be performed domain-generally in both music and language regardless of sensory modalities. Thus, the terms implicit learning and statistical learning have been used interchangeably and are regarded as the same phenomenon [6], although some researchers stated statistical knowledge and implicit knowledge may exactly not be same thing. In daily life, using the framework of statistical learning, we benefit immensely from artificial intelligence such as machine learning that gives computers learning ability similar to humans [7,8]. By calculating statistics on what we did in the past, computers can predict what we will do next, and they also give the best information, even when we are unaware of what we statistically want. In several fields of study on machine learning and natural language and music processing [918], the Markov chain has often been used as a model of the artificial grammar of language and music. The Markov chain, which was first reported by Markov [19], is a mathematical system in which the probability of the forthcoming state is statistically defined only by the latest state. The use of the Markov chains embedded in tone sequences allows us to verify the statistical structure of music [1218] and statistical learning and knowledge in humans [2023].

There are a number of evidences suggesting that statistical knowledge is implicit, although some researchers report that knowledge gained from statistical learning may result in both implicit and explicit knowledge [24]. Due to the implicitness of statistical learning in humans, however, we cannot specify the learned statistical knowledge, even when the learning effects are neurophysiologically unveiled [21,25,26] based on the framework of predictive coding [27]. In other words, implicit statistical knowledge can be acquired in our brains even though it does not reach the level of explicit awareness of what we learned. Despite the fact that humans are unaware of the learned statistical knowledge of auditory sequences, they can, however, identify that sequences containing high transitional probabilities sound familiar [20,21]. That is to say, although statistical knowledge itself may be implicit to a certain degree and does not reach the level of explicit awareness, it can be alternatively expressed via an abstract medium such as a musical melody [1218].

1.2. Characteristics of music: Implicit and explicit aspects

Music has numerous domain-specific structures (e.g., isochronic metrical grids, tonal pitch spaces, hierarchical tension, and attraction contours based on the structure of a melody and chord progression), and cannot be fully explained solely by the domain-general statistics of transitional probabilities [28,29]. Thus, the general hypothesis is that music is learned by both statistical learning and music-specific learning [30,31]. In general, musical composers explicitly learn music-specific knowledge such as musical grammar, and intentionally follow these frameworks when composing music. Musical rules and how they should be followed vary over time and among traditions, genres, and composers themselves. Therefore, the characteristics of music can be extracted based on music-specific structures, such as the harmony, tonalities, relative pitches, and musical intervals in the musical scores [3242]. To the best of our knowledge, however, few studies have verified whether the characteristics of music can be extracted based on the statistical structure in the musical scores. Furthermore, no study has investigated how the characteristics of the statistical structure in music varies over a composer’s lifetime.

1.3. Informatics approach of musical implicit learning

According to the previous studies, musical creativities in part depend on statistical learning [43], and the statistical distribution of music could partly represent composer’s statistical knowledge [1218]. That is, in the framework of human statistical learning theory and information theory, a higher-probability sequence may represent one that a composer is more likely to choose compared to lower-probability sequence. Indeed, informatics approach is often used to understand general music acquisition [3342] and the mental representation of implicit knowledge [4446]. Particularly, PARSER [47], Competitive Chunker [48], Information Dynamics of Music (IDyOM) [12,15], and n-gram models [49] underpin the hypothesis that music is acquired based on statistical distribution of music. Rohrmeier and colleague also applied neural network such as simple recurrent network (SRN) to human’s implicit learning model [50]. They showed that strong learning effect for n-gram models and a weaker effect for SRN, partly suggesting that SRN and nGram (i.e., transitional probability based on Markov chain) models may represent the ability of general humans and musical experts, respectively. Thus, the time-course variations of statistical structures in music may represent the time-course variations of a professional composer’s statistical knowledge. Nevertheless, there are few studies that investigate time-course variations of a composer’s statistical knowledge by analyzing musical score. To understand mental representation of statistical knowledge in music, it is important to examine the relationships among statistical structures in music, music composition, and implicit knowledge.

1.4. Musicological aspects on the statistics of Ludwig van Beethoven’s Piano Sonata

In the transition from eighteenth-century musical classicism to nineteenth-century romanticism in Western music, Ludwig van Beethoven was a German composer and pianist (1770–1882) who was highly influential to other composers. In musicological study, Beethoven’s compositional career is divided into early (around 1802), middle (around 1802 to 1814), and late periods (from about 1814) [51,52]. Compared to his works in the early period, which were strongly influenced by his predecessors in classicism, such as Wolfgang Amadeus Mozart (1756–1792) and Franz Joseph Haydn (1732–1809), his works in the late period, when he was already troubled by deafness and irritability brought on by chronic abdominal pain, were considered to show his personal expression and intellectual depth [51,52]. It is also believed that Beethoven always explored new directions and gradually expanded his scope of music over his lifetime [51,52]. Musicological researchers suggest that Beethoven’s music expressed his vision of life [53,54]. Sullivan [53] claimed that Beethoven’s vision was the product of his character and experience. Thus, his psychological variations on thinking and experience may form the character of his music [55]. The present study hypothesized that the time-course variations of statistical structures in Beethoven’s music over his lifetime might reflect the time-course variations of his statistical knowledge. As such, it would be very interesting to examine if the psychological variations in which Beethoven explored new directions and gradually expanded his scope of music over his lifetime [51,52] were reflected in the statistics of his music. Based on the hypothesis, the present study investigated how statistical structures in music vary over a composer’s lifetime.

1.5. Study purpose

The present study aimed to investigate how statistical structures in music vary over a composer’s lifetime. The transitional probabilities of the sequences containing the highest pitches in Ludwig van Beethoven’s Piano Sonata with all of the movements (Piano Sonata No.1 in F minor, Op.2-1 to No.32 in C minor, Op.111) were calculated based on seven different hierarchical Markov stochastic models (i.e., first- to seven-order Markov chains). Although melody is sometimes not highest pitches e.g. bass melodies), the present study only analyzed the highest pitch because the definition of melody in each title of music is still controversial in musicological study, different melodies could concurrently appear in some titles of music, and melody is often played in highest pitches. To understand how the transitional probability of more general sequences that is consistently used by one composer, the present study just targeted the interval pattern that appear in all 32 sonatas. The time-course variations of the transitional probabilities in each interval pattern that appear in all 32 sonatas were examined. It was hypothesized that there may be two types of time-course variations: transitional probabilities that gradually decrease and those that gradually increase, consistent with the sonata opus number. If so, these findings suggest that the statistical knowledge of music gradually shift over a composer’s lifetime.

2. Methods

The Ludwig van Beethoven’s Piano Sonata and all of its movements (No.1 in F minor, Op.2-1 to No.32 in C minor, Op.111, composed 1795–1822) was used in the present study (Beethoven, Piano Sonata No.1 to No.32, Breitkopf & Härtel, Leiptig). Using a scorewriter (Finale version 25, MI Seven Japan, Inc.), electronic scoring data of the sequences of highest pitch were extracted from the XML files. The highest pitches were chosen based on the following definitions: the highest pitches that can be played at a given point in time, the pitches with slurs can be counted as one, and the grace notes were excluded (S1 Table). Using all the sequences of highest pitches in a movement of a Sonata, the transitional probabilities of the sequences of highest pitches were calculated as a statistic based on Markov chains. The weighted averages of transitional probabilities of all the movements in a Piano Sonata were calculated. The probability of a forthcoming tone was statistically defined by the last tone to seven successive tones, respectively (i.e., first- to seven-order Markov chains).

n-order Markov model is based on the conditional probability of an element ei+1, given the preceding n elements:

(1)

Then, for each type of pitch transition, all pitches were numbered so that the first pitch was 0 in each transition, and an increase or decrease in a semitone was 1 and -1 based on the first pitch, respectively. The representative examples were shown in Fig 1. This revealed interval patterns but not pitch pattern [43,56]. This procedure was employed to eliminate the effects of the change of key on interval patterns. The interpretation of the change of key depends on musicians, and it is difficult to define in an objective manner. Thus, the results in the present study may represent a variation of statistics associated with relative pitch rather than absolute pitch.

thumbnail
Fig 1.

Representative transitions of [0,–2,–3,–5,–7] in Beethoven’s Piano Sonata No.1 in F minor, Op.2-1 (a), and those of [0,–2,–3,–5,–2] in Beethoven’s Piano Sonata No.32 in C minor, Op.111 (b). Based on the fourth-order Markov chain, the forthcoming states with the highest transitional probability defined by the last four states of [0, –2, –3, –5] are -7 in No.1 and -2 in No.32, respectively (see Table 10).

https://doi.org/10.1371/journal.pone.0196493.g001

Each interval pattern was ordered based on the Sonata opus numbers (i.e., Piano Sonata No.1 in F minor, Op.2-1 to No.32 in C minor, Op.111). Using the interval patterns that appear in all 32 sonatas, the time-course variations of the transitional probabilities were analyzed by multiple regression analyses using the stepwise method. Statistical significance levels were set at p = 0.05 for all analyses. The criteria of the variance inflation factor (VIF) and condition index (CI) were set at VIF < 2 and CI < 20 to confirm that there was no multi collinearity [57]. Based on the stepwise method, the regression model that showed greatest numbers of interval patterns were adopted. Next, all of the highest-order transitional probabilities that showed statistical significance were calculated in the first (No.1) and last (No.32) sonatas. Furthermore, the representative phrases that transitioned with the highest probability in the first and last sonatas were decoded as music scores.

3. Results

Twenty-two interval patterns with two tones that appear in all 32 sonatas were detected ([0,–1], [0,–12], [0,–2], [0,–3], [0,–4], [0,–5], [0,–6], [0,–7], [0,–9], [0,0], [0,1], [0,12], [0,2], [0,3], [0,4], [0,5], [0,6], [0,7], [0,8], and [0,9]) (Table 1). Using these interval patterns, a multiple linear regression based on the stepwise method was carried out to predict the sonata opus numbers (i.e., Piano Sonata No.1 in F minor, Op.2-1 to No.32 in C minor, Op.111), based on the transitional probabilities of the 22 interval patterns in the first-order Markov chain. A significant regression equation was found (F(2,29) = 5.96, p = 0.007), with an adjusted R2 of 0.24 (Table 2). The predicted sonata opus number is equal to 18.56–137.64 (transition of [0,1]) + 147.308 (transition of [0,2]). The transitional probabilities of [0,1] and [0,2] gradually decreased and increased consistently with the ascending order of the sonata opus numbers, respectively ([0,1] p = 0.007, [0,2] p = 0.021). These transitional probabilities were significant predictors of the sonata opus numbers.

thumbnail
Table 1. Transitional probabilities calculated using first-order Markov chains for each of the interval patterns.

https://doi.org/10.1371/journal.pone.0196493.t001

thumbnail
Table 2. Multiple regression analyses based on the stepwise method in first-order Markov chain.

https://doi.org/10.1371/journal.pone.0196493.t002

Thirty-seven interval patterns with three tones that appear in all 32 sonatas were detected ([0,–1,–1], [0,–1,–3], [0,–1,0], [0,–1,2], [0,–2,–2], [0,–2,–3], [0,–2,–4], [0,–2,0], [0,–2,2], [0,–2,5], [0,–3,–2], [0,–3,–5], [0,–4,–5], [0,–5,0], [0,0,–1], [0,0,–2], [0,0,0], [0,0,1], [0,0,2], [0,0,3], [0,0,5], [0,1,–1], [0,1,–2], [0,1,–4], [0,1,0], [0,1,1], [0,1,3], [0,2,0], [0,2,3], [0,2,4], [0,3,0], [0,3,1], [0,3,2], [0,3,3], [0,4,2], [0,4,7], and [0,5,9]) (Table 3). Using these interval patterns, a multiple linear regression based on the stepwise method was carried out to predict the numbers of music based on the transitional probabilities of the 37 interval patterns in the second-order Markov chain. A significant regression equation was found (F(4,27) = 10.10, p = 0.0001), with an adjusted R2 of 0.54 (Table 4). The predicted sonata opus number is equal to 32.72–71.15 (transition of [0,–2,–4])– 50.24 (transition of [0,0,–1]) + 36.37 (transition of [0,–4,–5]) + 42.27 (transition of [0,–2,0]). The transitional probabilities of [0,–2,–4] and [0,0,–1] gradually decreased and the transitional probabilities of [0,–4,–5], and [0,–2,0] gradually increased consistently with the ascending order of the sonata opus numbers ([0,–2,–4] p = 0.011, [0,0,–1] p = 0.002, [0,–4,–5] p = 0.008, [0,–2,0] p = 0.020). These transitional probabilities were significant predictors of the sonata opus numbers.

thumbnail
Table 3. Transitional probabilities calculated using second-order Markov chains for each of the interval patterns.

https://doi.org/10.1371/journal.pone.0196493.t003

thumbnail
Table 4. Multiple regression analyses based on the stepwise method in second-order Markov chain.

https://doi.org/10.1371/journal.pone.0196493.t004

Twelve interval patterns with four tones that appear in all 32 sonatas were detected ([0,–1,–3,–5], [0,–1,0,–1], [0,–1,0,2], [0,–2,–3,–2], [0,–2,–3,–5], [0,–2,–4,–5], [0,0,0,0 ], [0,1,3,1 ], [0,1,3,5], [0,2,0,–1], [0,2,3,5], and [0,2,4,5]) (Table 5). Using these interval patterns, a multiple linear regression based on the stepwise method was carried out to predict the sonata opus numbers based on the transitional probabilities of the 12 interval patterns in the third-order Markov chain. A significant regression equation was found (F(2,29) = 9.25, p = 0.001), with an adjusted R2 of 0.35 (Table 6). The predicted sonata opus number is equal to 39.12–28.81 (transition of [0,2,0,–1])– 24.82 (transition of [0,2,4,5]). The transitional probabilities of [0,2,0,–1] and [0,2,4,5] gradually decreased consistently with the ascending order of the sonata opus numbers ([0,2,0,–1] p = 0.003, [0,2,4,5] p = 0.011). These transitional probabilities were significant predictors of the sonata opus numbers.

thumbnail
Table 5. Transitional probabilities calculated using third-order Markov chains for each of the interval patterns.

https://doi.org/10.1371/journal.pone.0196493.t005

thumbnail
Table 6. Multiple regression analyses based on the stepwise method in third-order Markov chain.

https://doi.org/10.1371/journal.pone.0196493.t006

Three interval patterns with five tones that appear in all 32 sonatas were detected ([0,–1,–3,–5,–6], [0,–2,–3,–5,–7], and [0,–2,–4,–5,–7]) (Table 7). Using these interval patterns, a multiple linear regression based on the stepwise method was carried out to predict the sonata opus numbers based on the transitional probabilities of the three interval patterns in the fourth-order Markov chain. A significant regression equation was found (F(1,30) = 6.65, p = 0.015), with an adjusted R2 of 0.15 (Table 8). The predicted sonata opus number is equal to 27.13–19.27 (transition of [0,–2,–3,–5,–7]). The transitional probabilities gradually decreased consistently with the ascending order of the sonata opus numbers (p = 0.015). The transitional probabilities were significant predictors of the sonata opus numbers. To understand how the transitional probabilities of [0,–2,–3,–5,–7] were shifted in more detail, the transition matrices based on the fourth-order Markov chain of P(X|0, -2, -3, -5) in the first (No.1) and last (No.32) sonatas are shown in Table 9. Furthermore, the representative phrases that transitioned with the highest probability in the first and last sonata were decoded as music scores (Fig 1).

thumbnail
Table 7. Transitional probabilities calculated using fourth-order Markov chains for each of the interval patterns.

https://doi.org/10.1371/journal.pone.0196493.t007

thumbnail
Table 8. Multiple regression analyses based on the stepwise method in fourth-order Markov chain.

https://doi.org/10.1371/journal.pone.0196493.t008

thumbnail
Table 9. Transition matrices based on the fourth-order Markov chain (P(X|0, -2, -3, -5)).

https://doi.org/10.1371/journal.pone.0196493.t009

One interval pattern with six tones that appeared in all 32 sonatas was detected (0,–2,–4,–5,–7,–9) (Table 10). Using these interval patterns, a multiple linear regression based on the stepwise method was carried out to predict the sonata opus numbers based on the transitional probabilities of the interval patterns in the fifth-order Markov chain. No significant regression equation was detected, however. No sixth- or higher-order interval patterns appeared in all 32 sonatas.

thumbnail
Table 10. Transitional probabilities calculated using first-order Markov chains for each of the interval patterns.

https://doi.org/10.1371/journal.pone.0196493.t010

4. Discussion

4.1. Psychological relationships between variation of statistical structure in music and implicit knowledge

A tone with a higher transitional probability of sequences in a musical score may be one that a composer is more likely to choose compared to other tones with lower transitional probability [2023,5860]. Thus, the transitional probability matrix calculated from music may represent the characteristics of a composer’s statistical knowledge by which a forthcoming tone is implicitly defined. The present study aimed to investigate how the statistical structures in music vary over a composer’s lifetime. I hypothesized that there might be two types of time-course variations: transitional probabilities that gradually decrease and those that gradually increase, consistent with the order of the sonata opus number. If so, these findings suggest that statistical knowledge of music, which are generally considered implicit phenomenon, gradually shifts over a composer’s lifetime.

Consistent with the order of the sonata opus numbers, the transitional probabilities of [0,1], [0,–2,–4], [0,0,–1], [0,2,0,–1], [0,2,4,5], and [0,–2,–3,–5,–7] gradually decreased, and those of [0,2], [0,–4,–5], and [0,–2,0] gradually increased. All of these transitions may be musical universal. The examples of each transition were shown in Table 11. The transition matrices of P(X|0, -2, -3, -5) in the first (No.1) and last sonatas (No.32) are shown in Table 10. The representative transition of [0,–2,–3,–5,–7] in No.1 was decoded as musical scores in Fig 1A. In these familiar sequences, decreasing transitional probabilities were detected more frequently than increasing transitional probabilities. This finding suggests that, in the later period of his life, Beethoven composed statistically novel music in which predictable and familiar transitions did not appear compared to the earlier periods. The sequences with 5 tones based on fourth-order transitional probabilities were the longest in the transitions that showed significance, suggesting that the time-course variation of statistical structure that could implicitly reflect strategy of music composition may depend on until fourth-order transitional probabilities. According to the previous neurophysiological study on statistical learning, fourth-order transitional probabilities in auditory sequences modulate prediction in human’s brain [61]. Furthermore, they suggested that, higher-order transitional probabilities could modulate learner’s prediction when they learn more difficult statistical structure. These findings may suggest that the findings of the present study are in part related to prediction and expectation in the composer, although this study could not detect difference among the orders of transitional probability. On the other hand, the only interval patterns that appear in all 32 sonatas were used to verify the variation of transitional probabilities of the familiar sequences in the 32 sonatas. In total, there are more than several billion kinds of interval patterns in music. Thus, it is difficult to define the statistical significance when we target all of the interval patterns. Future study will need to investigate interval patterns that do not appear in some sonatas.

thumbnail
Table 11. The examples of transitions in which transitional probabilities were gradually increased and decreased.

https://doi.org/10.1371/journal.pone.0196493.t011

To understand how the highest-order (i.e., fourth-order) transitional probabilities of [0,–2,–3,–5,–7] shifted in more detail, the transition matrices based on the fourth-order Markov chain of P(X|0, -2, -3, -5) in the first (No.1) and last (No.32) sonatas were investigated (Table 10). In the transition matrices of Sonata No.1 (Table 10), the distribution of transitional probabilities was biased: the transitional probability of [0,–2,–3,–5,–7] was markedly higher than that of any other transitional probabilities. In contrast, in the transition matrices of Sonata No.32 (Table 10), the transitional probabilities were equally distributed compared to the transition matrices of Sonata No.1: there was no obvious highest transitional probability. The representative phrases with the highest transitional probability in Sonatas No.1 and No.32 were decoded as musical scores in Fig 1A and 1B. These findings are consistent with the hypothesis that, in the last sonata, Beethoven tried several composition methods in which he avoided familiar interval patterns. The previous studies indicated that statistical learning contributes to many types of mental representation in music: the comprehension and production of music [62], intuitive decision-making [63,64] and creativity involved in musical composition [15]. The several studies also suggest that, compared to language [65], musical representation including tonality is mainly formed by a tacit knowledge [6672]. Furthermore, the neurophysiological studies indicates that the distribution of pitch frequency that mainly forms musical melody, rather than other acoustic features such as formant frequencies that play an important role in forming speech sounds instead of melody, facilitate auditory statistical learning in human’s brain [22]. Thus, it is widely accepted that statistical knowledge is tied to musical expression such as composition, playing, and creativity [1217], even if learners cannot realize the acquired knowledge due to implicitness of statistical learning [21]. Given these previous findings, the results of the present study may suggest that the time-course variations of statistical structures in music represent the time-course variations of a composer’s statistical knowledge.

The previous studies demonstrated that the transitional probabilities of relative-pitch patterns could be learned [21,22,73]. Based on these findings, this study analyzed pitch-interval patterns but not absolute-pitch pattern like the methodologies of the previous studies [43,56] in order to eliminate the effects of the change of key on interval patterns. In other words, the present study verifies statistical distribution of relative- but not absolute-pitch structures.

4.2. Musicological relationships between Beethoven’s lifetime and the statistical variations of music

The numbers of Beethoven’s Piano Sonata indicate the order in which the music was published, not when it was composed. Although most of the numbers are consistent with the order in which the music was composed, Piano Sonata No. 19 in G minor and Piano Sonata No. 20 in G major, which were published around 1805, are considered to have been composed a decade earlier, between 1795 and 1798 [74,75]. This means that Piano Sonata No.19 and No.20 may have been written around the time that Beethoven composed Piano Sonata No.3 and No.4. The specific time of composition remains controversial. Therefore, in the present study, the time courses in which the music was published are used, as they are clearer than the time courses in which the music was composed. The present study revealed the time-course variations of statistical structures in Beethoven’s music over his lifetime. In other words, the statistically novel music was composed during Beethoven’s significant life events in late periods when he was already troubled by deafness and irritability brought on by chronic abdominal pain and were considered to show his personal expression and intellectual depth [51,52]. According to the previous studies, statistical learning could modulate strategies of musical composition, and the musical training and experience is associated with cognitive model of probabilistic structure in music involved in statistical learning [43,1218]. The neurophysiological studies also demonstrated that musical training modulates the abilities of statistical learning [76,77,78]. These previous findings may suggest that musical experience allows a composer to change the strategies of musical composition based on the acquired statistical knowledge over his/her lifetime. The results of time-course variation of statistical structure in the present study may imply that, in later compositional periods, because Beethoven experienced a lot of composition strategies and explored new directions, he might compose statistically novel music in which predictable transitions did not appear compared to that of earlier periods. According to the previous studies, corpus analyses can detect the historical characteristics of music and in part distinguish them based on the era (e.g., [37,40,79]), suggesting that variation of statistical structure can be detected over long periods of time. The present corpus study also detected the variation of characteristics within a composer in his/her lifetime. It is of note that the present study did not directly demonstrate that the implicit statistical knowledge of music varied, as only the statistics of musical scores were analyzed. This suggests that there are other possible explanations for the results in the present study. For instance, it might be part of his plan to compose the sonatas from familiar to increasingly unfamiliar based on the statistical structure of music. It cannot be excluded the possibility that the findings in the present study do not necessarily reflect that Beethoven’s statistical knowledge changed. The human’s ability to generate random sequences of numbers [80] has been associated with creativity in human [81]. In addition, to understand variation of transitional probability of more general phrase that is consistently used by Beethoven, the present study only analyzed the sequences that appeared in all pieces of music. One of the reasons is that the phrases that do not appear in all sonatas and showed transitional probability of 0% in some pieces of music is difficult to define composer’s knowledge. Compared with these phrases, the consistently appeared phrases may depend on his knowledge. However, there may be another specific variation in phrases that showed transitional probability of 0% in some pieces of music. Future study should investigate the effects of statistical learning and knowledge of all types of phrase on music compositions. The present study shed new light on novel methodologies that may be able to evaluate the time-course variation of composer’s implicit knowledge using musical scores.

Supporting information

S1 Table. The coded sequence in each piece of music.

https://doi.org/10.1371/journal.pone.0196493.s001

(XLSX)

Acknowledgments

This work was supported by Grant-in-Aid for Nakayama Foundation for Human Science. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

  1. 1. Reber AS. Implicit learning of artificial grammars. J. Verbal Learn. Verbal Behav. 1967; 6: 855–863.
  2. 2. Reber AS. Implicit Learning and Tacit Knowledge: An Essay on the Cognitive Unconscious, Oxford University Press, 1993.
  3. 3. Cleeremans A., et al. Implicit learning: News from the front. Trends Cogn. Sci. 1998; 2: 406–416. pmid:21227256
  4. 4. Shanks DR. Implicit learning. In Handbook of Cognition (Lamberts K and Goldstone R, eds), pp. 202–220, Sage Publications, 2005.
  5. 5. Saffran JR, Aslin RN, Newport EL. Statistical learning by 8-month-old infants. Science. 1996; 274: 1926–1928. pmid:8943209
  6. 6. Perruchet P, Pacton S. Implicit learning and statistical learning: one phenomenon, two approaches. Trends Cogn Sci. 2006; 10: 233–238. pmid:16616590
  7. 7. Blackburn P, Bos J. “Representation and Inference for Natural Language: A First Course,” in Computational Semantics (Studies in Computational Linguistics), Center for the Study of Language and Information, 2005.
  8. 8. Clark A, Fox F. The Handbook of Computational Linguistics and Natural Language Processing,Wiley-Blackwell, 2012.
  9. 9. Singla P, Domingos P. Entity resolution with markov logic. In ICDM ‘06: Proceedings of the sixth international conference on data mining (pp. 572–582). Washington, DC, USA: IEEE Computer Society, 2006.
  10. 10. Richardson M, Domingos P. Markov logic networks. Machine Learning. 2006; 62: 107–136.
  11. 11. Poon H, Domingos P. Joint inference in information extraction. In Proceedings of the twenty-second national conference on artificial intelligence (pp. 913–918). Vancouver, Canada: AAAI Press, 2007.
  12. 12. Pearce M. The Construction and Evaluation of Statistical Models of Melodic Structure in Music Perception and Composition. PhD thesis, School of Informatics, City University, London, 2005.
  13. 13. Pearce M, Wiggins GA. Expectation in melody: The influence of context and learning. Music Perception. 2006; 23: 377–405.
  14. 14. Pearce M, Müllensiefen D, Wiggins G. The role of expectation and probabilistic learning in auditory boundary perception: A model comparison. Perception. 2010; 39: 1365–1389. pmid:21180358
  15. 15. Pearce M, Wiggins G. Auditory expectation: The information dynamics of music perception and cognition. Topics in Cognitive Science. 2012; 4: 625–652. pmid:22847872
  16. 16. Omigie D, Pearce M, Stewart L. Tracking of pitch probabilities in congenital amusia. Neuropsychologia. 2012; 50: 1483–1493. pmid:22414591
  17. 17. Omigie D, Pearce M, Williamson V, Stewart L. Electrophysiological correlates of melodic processing in congenital amusia. Neuropsychologia. 2013; 51: 1749–1762. pmid:23707539
  18. 18. Hansen NC, Pearce M. Predictive Uncertainty in Auditory Sequence Processing. Frontiers in Psychology. 2014; 5: 1052. pmid:25295018
  19. 19. Markov AA. Extension of the limit theorems of probability theory to a sum of variables connected in a chain. Markov chains (Vol. 1). John Wiley and Sons (Reprinted in Appendix B of: R. Howard D. Dynamic Probabilistic Systems), 1971.
  20. 20. Furl N1, Kumar S, Alter K, Durrant S, Shawe-Taylor J, Griffiths TD. Neural prediction of higher-order auditory sequence statistics. NeuroImage. 2011; 54: 2267–2277. pmid:20970510
  21. 21. Daikoku T, Yatomi Y, Yumoto M. Implicit and explicit statistical learning of tone sequences across spectral shifts. Neuropsychologia. 2014; 63: 194–204. pmid:25192632
  22. 22. Daikoku T, Yatomi Y, Yumoto M. Statistical learning of music- and language-like sequences and tolerance for spectral shifts. Neurobiol Learn Mem. 2015; 118: 8–19. pmid:25451311
  23. 23. Daikoku T, Yatomi Y, Yumoto M. Pitch-class distribution modulates the statistical learning of atonal chord sequences. Brain and Cognition. 2016; 108: 1–10. pmid:27429093
  24. 24. Batterink LJ, Reber PJ, Neville HJ, Paller KA. Implicit and explicit contributions to statistical learning. J Mem Lang. 2015; 83: 62–78. pmid:26034344
  25. 25. Abla D, Katahira K, Okanoya K. On-line Assessment of Statistical Learning by Event-related Potentials, Journal of Cognitive Neuroscience. 2008; 20: 952–964. pmid:18211232
  26. 26. Koelsch S., Busch T., Jentschke S., Rohrmeier M. Under the hood of statistical learning: a statistical MMN reflects the magnitude of transitional probabilities in auditory sequences. Sci. Rep. 2016; 6: 19741. pmid:26830652
  27. 27. Friston K. A theory of cortical responses. Phil. Trans. R. Soc. 2005; B 360: 815–836.
  28. 28. Hauser MD, Chomsky N, Fitch WT. The faculty of language: What is it, who has it, and how did it evolve? Science. 2002; 298: 1569–1579. pmid:12446899
  29. 29. Jackendoff R, Lerdahl F. The capacity for music: What is it, and what’s special about it? Cognition. 2006; 100: 33–72. pmid:16384553
  30. 30. Jonaitis EM, Saffran JR. Learning harmony: The role of serial statistics. Cognitive Science. 2009; 33: 951–968. pmid:21585492
  31. 31. Tillmann B, Bharucha JJ, Bigand E. Implicit learning of tonality: A selforganizing approach. Psychological Review. 2000; 107: 885–913. pmid:11089410
  32. 32. Pollastri E, Simoncelli G. Classification of melodies by composer with hidden Markov models. Web Delivering of Music. Proceedings. First International Conference, 88–95, 2001.
  33. 33. Cilibrasi R, Vitanyi P, Wolf R. Algorithmic Clustering of Music Based on String Compression. Computer Music Journal. 2004; 28: 49–67.
  34. 34. Hardoon DR, Saunders C, Shawe-Taylor J. Using Fisher Kernels and Hidden Markov Models for the Identification of Famous Composers from their Sheet Music. PASCAL 2005.
  35. 35. Backer E, van Kranenburg P. “On musical stylometry––a pattern recognition approach”, in Pattern Recognition Let- ters. 2005; 26: 299–309.
  36. 36. Kranenburg P. Composer attribution by quantifying compositional strategies. Proc ISMIR, 2006.
  37. 37. Gjerdingen RO. Music Perception: An Interdisciplinary Journal. 2012; 31: 192–204.
  38. 38. Ito JP. Koch’s Metrical Theory and Mozart’s Music: A Corpus Study. Music Perception: An Interdisciplinary Journal. 2012; 31: 205–222.
  39. 39. Prince JB, Schmuckler MA. The Tonal-Metric Hierarchy: A Corpus Analysis. Music Perception: An Interdisciplinary Journal. 2012; 31: 254–270.
  40. 40. Albrecht J, Huron D. A Statistical Approach to Tracing the Historical Development of Major and Minor Pitch Distributions, 1400–1750. Music Perception: An Interdisciplinary Journal. 2012; 31: 223–243.
  41. 41. Albrecht J, Shanahan D. The Use of Large Corpora to Train a New Type of Key-Finding Algorithm: An Improved Treatment of the Minor Mode. Music Perception: An Interdisciplinary Journal. 2013; 31: 59–67.
  42. 42. London J. Building a Representative Corpus of Classical Music. Music Perception: An Interdisciplinary Journal. 2013; 31: 68–90.
  43. 43. Norgaard M. How jazz musicians improvise: The central role of auditory and motor pa erns. Music Perception. 2014; 31: 271–287.
  44. 44. Rohrmeier M, Rebuschat P. Implicit learning and acquisition of music. Top Cogn Sci. 4:525–53 (2012). pmid:23060126
  45. 45. Dubnov S. Musical information dynamics as models of auditory anticipation. In Wang W.(Ed.), Machine audition: Principles, algorithms and systems (pp. 371–397). Hershey, PA: IGI Global Publication (2010).
  46. 46. Wang W. (Ed.) Machine audition: Principles, algorithms and systems. Hershey, PA: IGI Global Publication (2010).
  47. 47. Perruchet P, Vinter A. PARSER: A model of word segmentation. Journal of Memory and Language, 39, 246–263 (1998).
  48. 48. Servan-Schreiber D, Anderson JR. Learning artificial grammars with competitive chunking. Journal of Experimental Psychology: Learning, Memory, and Cognition, 16, 592–608 (1990).
  49. 49. Pearce MT, Wiggins GA. Improved methods for statistical modelling of monophonic music. Journal of New Music Research 33, 367–385 (2004).
  50. 50. Rohrmeier MA, Cross I. Modelling unsupervised online-learning of artificial grammars: linking implicit and statistical learning. Conscious Cogn. 2014; 27:155–67. pmid:24905545
  51. 51. Adorno-Wiesengrund TW. Beethoven: The Philosophy of Music; Fragments and Texts, 1993.
  52. 52. Dahlhaus C. Ludwig van Beethoven: Approaches to His Music. Translated by Mary Whittall. Oxford: Clarendon Press; New York: Oxford University Press, 1991.
  53. 53. Sullivan JWN. Beethoven: His Spiritual Development. New York: A.A. Knopf, 1927.
  54. 54. Boucourechliev A. Beethoven (in French). Seuil, 1963 (ISBN: 978–2020002431).
  55. 55. Johnson D, Tyson A, Wnter R (1985). The Beethoven Sketchbooks, Oxford: Clarendon.
  56. 56. Lynch MP, Eilers RE, Oller KD, Urbano RC, Wilson P. Influences of acculturation and musical sophistication on perception of musical interval patterns. J Exp Psychol Hum Percept Perform. 1991; 17: 967–75. pmid:1837307
  57. 57. Cohen H, Cohen O, West SG. Aiken LS. Applied multiple regression/correlation analysis for the behavioral scineces. 3rd eds. LEA, 2003.
  58. 58. Yumoto M, Daikoku T. IV Auditory system. 5 basic function. In Tobimatsu S& Kakigi R (Eds.). Clinical Applications of Magnetoencephalography. Springer (ISBN 4431557288) (2016).
  59. 59. Daikoku T, Yatomi Y, Yumoto M. Statistical learning of an auditory sequence and reorganization of acquired knowledge: A time course of word segmentation and ordering. Neuropsychologia. 2017; 95: 1–10. pmid:27939187
  60. 60. Daikoku T, Yumoto M. Single, but not dual, attention facilitates statistical learning of two concurrent auditory sequences. Sci Rep. 2017;7: 10108. pmid:28860466
  61. 61. Daikoku T, Okano T, Yumoto M. Relative difficulty of auditory statistical learning based on tone transition diversity modulates chunk length in the learning strategy. Biomagnetic Sendai, T-16 (2017).
  62. 62. Rohrmeier M, Rebuschat P. Implicit learning and acquisition of music. Top Cogn Sci. 4:525–53 (2012). pmid:23060126
  63. 63. Berry DC, Dienes Z. Implicit learning: Theoretical and empirical issues. Hove, England: Lawrence Erlbaum (1993).
  64. 64. Perkovic S, Orquin JL. Implicit Statistical Learning in Real-World Environments Leads to Ecologically Rational Decision Making. Psychol Sci. 1:956797617733831 (2017).
  65. 65. Chomsky N. Syntactic structures. The Hague, The Netherlands: Mouton (1957).
  66. 66. Bigand E, Poulin-Charronnat B. Are we ‘‘experienced listeners”? A review of the musical capacities that do not depend on formal musical training. Cognition, 100, 100–130 (2006). pmid:16412412
  67. 67. Ettlinger M, Margulis EH, Wong PCM. Implicit memory in music and language. Frontiers in Psychology, 211(2), (2011). pmid:21927608
  68. 68. Deliège I, Mèlen M, Stammers D, Cross I. Musical schemata in real time listening to a piece of music. Music Perception, 14(2), 117–160 (1996).
  69. 69. Deliège I. Prototype effects in music listening: An empirical approach to the notion of imprint. Music Perception 18, 371–407 (2001).
  70. 70. Huron D. Two challenges in cognitive musicology. Topics in Cognitive Science 4, 678–684 (2012). pmid:23060127
  71. 71. Koelsch S, Gunter T, Friederici A. Brain indices of music processing: ‘‘Nonmusicians” are musical. Journal of Cognitive Neuroscience, 12(3), 520–541 (2000). pmid:10931776
  72. 72. Koelsch S. Toward a neural basis of music perception—a review and updated model. Front Psychology, 2, 110. (2011). pmid:21713060
  73. 73. Saffran JR, Reeck K, Niebuhr A, Wilson D. Changing the tune: the structure of the input affects infants' use of absolute and relative pitch. Dev Sci, 8:1–7 (2005). pmid:15647061
  74. 74. Nottebohm G. Beethoveniana, Leipzig: Peters, 1872 (r: New York: Johnson Reprint 1970).
  75. 75. Nottebohm G. Zweite Beethoveniana, Leipzig: Peters, 1887 (r: New York: Johnson Reprint 1970)
  76. 76. Paraskevopoulos E, Kuchenbuch A., Herholz SC, & Pantev C. Statistical learning effects in musicians and non-musicians: An MEG study. Neuropsychologia, 50, 341–349 (2012). pmid:22197571
  77. 77. Paraskevopoulos E, Chalas N, Bamidis P. Functional connectivity of the cortical network supporting statistical learning in musicians and non-musicians: an MEG study. Sci Rep. 2017, 24;7:16268.
  78. 78. Mandikal Vasuki PR, Sharma M, Ibrahim R, Arciuli J. Statistical learning and auditory processing in children with music training: An ERP study. Clin Neurophysiol. 2017 Jul;128(7):1270–1281. pmid:28545016
  79. 79. White CM. Changing Styles, Changing Corpora, Changing Tonal Models. Music Perception: An Interdisciplinary Journal. 2012; 31: 244–253.
  80. 80. Wagenaar WA. Generation of random sequences by human subjects: A critical survey of the literature. Psychological Bulletin. 1972; 77: 65–72.
  81. 81. Bains W. Random number generation and creativity. Med Hypotheses. 2008;70:186–90. pmid:17920778