<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.4 20241031//EN" "JATS-journalpublishing1-4.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="1.4" xml:lang="en">
  <front>
    <journal-meta>
      <journal-id journal-id-type="publisher-id">jpee</journal-id>
      <journal-title-group>
        <journal-title>Journal of Power and Energy Engineering</journal-title>
      </journal-title-group>
      <issn pub-type="epub">2327-5901</issn>
      <issn pub-type="ppub">2327-588X</issn>
      <publisher>
        <publisher-name>Scientific Research Publishing</publisher-name>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.4236/jpee.2026.146002</article-id>
      <article-id pub-id-type="publisher-id">jpee-151769</article-id>
      <article-categories>
        <subj-group>
          <subject>Article</subject>
        </subj-group>
        <subj-group>
          <subject>Engineering</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>Assessment Analysis of Neutral to Earth Voltage (NTEV) in the Distribution System at the Government Building Using Statistical and Analytical Methods Analysis</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <name name-style="western">
            <surname>Ahamad</surname>
            <given-names>Muhammad Edzham bin</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <name name-style="western">
            <surname>Baharuddin</surname>
            <given-names>Zikri Abadi bin</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <name name-style="western">
            <surname>Ibrahim</surname>
            <given-names>Zulkifli bin</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <name name-style="western">
            <surname>Zahirudin</surname>
            <given-names>Mohd Ruzaini bin</given-names>
          </name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <name name-style="western">
            <surname>Rahman</surname>
            <given-names>C Mahbubur</given-names>
          </name>
          <xref ref-type="aff" rid="aff3">3</xref>
        </contrib>
      </contrib-group>
      <aff id="aff1"><label>1</label> Fakulti Teknologi Kejuruteraan Elektrik (FTKE), Universiti Teknikal Malaysia Melaka, Durian Tunggal, Melaka, Malaysia </aff>
      <aff id="aff2"><label>2</label> Bahagian Kejuruteran (BKJ), Agensi Nuklear Malaysia, Kajang, Selangor </aff>
      <aff id="aff3"><label>3</label> C Angstrom Laboratory, Division for Electricity, Department of Engineering Sciences, Uppsala University, Uppsala, Sweeden </aff>
      <author-notes>
        <fn fn-type="conflict" id="fn-conflict">
          <p>The authors declare no conflicts of interest regarding the publication of this paper.</p>
        </fn>
      </author-notes>
      <pub-date pub-type="epub">
        <day>08</day>
        <month>06</month>
        <year>2026</year>
      </pub-date>
      <pub-date pub-type="collection">
        <month>06</month>
        <year>2026</year>
      </pub-date>
      <volume>14</volume>
      <issue>06</issue>
      <fpage>15</fpage>
      <lpage>28</lpage>
      <history>
        <date date-type="received">
          <day>16</day>
          <month>12</month>
          <year>2025</year>
        </date>
        <date date-type="accepted">
          <day>05</day>
          <month>06</month>
          <year>2026</year>
        </date>
        <date date-type="published">
          <day>08</day>
          <month>06</month>
          <year>2026</year>
        </date>
      </history>
      <permissions>
        <copyright-statement>© 2026 by the authors and Scientific Research Publishing Inc.</copyright-statement>
        <copyright-year>2026</copyright-year>
        <license license-type="open-access">
          <license-p> This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license ( <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link> ). </license-p>
        </license>
      </permissions>
      <self-uri content-type="doi" xlink:href="https://doi.org/10.4236/jpee.2026.146002">https://doi.org/10.4236/jpee.2026.146002</self-uri>
      <abstract>
        <p>Elevated neutral-to-earth voltage or simply known as Neutral-To-Earth Voltage (NTEV), and the related phenomenon called stray voltage is typically caused by fundamental frequency currents returning to the source via the neutral conductor and earth. The nature of grounded power systems results in the fact that the neutral conductors are not always at zero potential with respect to the earth. This study also aims to investigate the effect of Neutral to Earth voltage on power distribution. Using the appropriate measurement equipment, the data logger was attached to Main Switch Boards to record data for 4 days for monitoring, recording, and analyzing frequency, harmonics, and voltage fluctuation. In addition, the data logger for this measurement tool have been installed on the main switchboard (MSB1) in Agensi Nuklear Malaysia. This study performed a comparison of peak time (8 am-5 pm) and non-peak time (6 pm-6 am) and also the correlation between each parameter. Based on the result of standard deviation and mean, all parameters deviate very closely to the mean value. The skewness value provide evidence that the measured data was distributed close to the nominal specification. Through correlation analysis, at the end of this study it was found that the correlation between NTEV and THDV of each phase in the first location has a very high correlation.</p>
      </abstract>
      <kwd-group kwd-group-type="author-generated" xml:lang="en">
        <kwd>Neutral-to-Earth Voltage (NTEV)</kwd>
        <kwd>Power Quality Disturbances</kwd>
        <kwd>Total Harmonics Distortion (THD)</kwd>
        <kwd>Pearson Correlation</kwd>
        <kwd>Descriptive Analysis</kwd>
        <kwd>Statistical Analysis</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec1">
      <title>1. Introduction</title>
      <p>Power quality encompasses the key attributes of electrical power, such as voltage, frequency, and waveform integrity, ensuring they meet the necessary standards for optimal device and equipment performance. Power quality disturbances, such as excessive Neutral-to-Earth Voltage (NTEV) and Total Harmonic Distortion (THD), can compromise operational efficiency and lead to equipment failure. Statistical analysis, including Pearson Correlation and Descriptive Analysis, are valuable tools in assessing and quantifying the impact of these disturbances. A detailed understanding of these factors is crucial for maintaining reliable and efficient electrical systems.</p>
      <p>The growing emphasis on power quality has led to notable advancements in monitoring technologies capable of characterizing disturbances and fluctuations in power quality. This paper examines the various types of power quality disturbances and explores methods for effectively analyzing and presenting this data in a meaningful way. </p>
      <sec id="sec1dot1">
        <title>1.1. Overview of NTEV</title>
        <p>Neutral-to-Earth Voltage (NTEV) in distribution systems is a critical concern affecting safety and reliability. Elevated NTEV often results from unbalanced loads, harmonic distortions, and poor grounding, causing voltage differences between the neutral point and earth. This can lead to stray voltage, compromising both human safety and equipment performance. Harmonic currents from nonlinear loads exacerbate the issue by generating triple harmonics that do not cancel out in the neutral conductor. Key factors contributing to elevated NTEV include unbalanced loading, system asymmetries, earth resistivity, and feeder length. By referring to IEE1159 (IEEE Recommended Practice for Monitoring Electric Power Quality), the acceptable thresholds for NTEV or stray voltage is 10 V [<xref ref-type="bibr" rid="B1">1</xref>]. <bold>Table 1</bold> shows the summary table for NTEV level.</p>
        <p><bold>Table 1.</bold> Summary table for NTEV level.</p>
        <table-wrap id="tbl1">
          <label>Table 1</label>
          <table>
            <tbody>
              <tr>
                <td>Condition/Location</td>
                <td>Max. NTEV (V)</td>
                <td>Comments</td>
              </tr>
              <tr>
                <td>Residential Areas</td>
                <td>1.0 V</td>
                <td>Standard for safety in residential system</td>
              </tr>
              <tr>
                <td>Commercial Areas</td>
                <td>1.5 V</td>
                <td>Slightly higher allowance for commercial equipment</td>
              </tr>
              <tr>
                <td>Industrial Areas</td>
                <td>2.0 V</td>
                <td>Higher due to larger electrical systems and loads</td>
              </tr>
              <tr>
                <td>Sensitive Equipment Zones</td>
                <td>0.5 V</td>
                <td>Critical areas such as hospitals or labs</td>
              </tr>
              <tr>
                <td>NTEV with Isolation Transformer</td>
                <td>0.5 V</td>
                <td>Significant reduction when isolation transformers are used</td>
              </tr>
              <tr>
                <td>Acceptable Limits (IEEE 519-2014)</td>
                <td>1.0 V</td>
                <td>Per IEEE guidelines for maintaining safe power quality</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <p>Shafie <italic>et</italic><italic>al.</italic> (2010) examine the use of isolation transformers to mitigate harmonic voltages and Neutral-to-Earth Voltage (NTEV) in distribution systems, focusing on their effectiveness at Universiti Teknologi MARA in Penang. NTEV, or common-mode voltage, represents the potential difference between the neutral conductor and safety ground, which can cause operational issues, equipment failure, and higher losses for linear loads [<xref ref-type="bibr" rid="B2">2</xref>]. The results show that the isolation transformer significantly reduced both NTEV and Total Harmonic Distortion (THD). For example, <xref ref-type="fig" rid="fig1">Figure 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref> show the maximum NTEV in Computer Lab 1 decreased from 1.464 V to 0.977 V, and the THD for phase voltage dropped from 1.927% to 1.501%. These improvements are in line with IEEE 519-2014 standards on harmonic distortion [<xref ref-type="bibr" rid="B3">3</xref>].</p>
        <fig id="fig1">
          <label>Figure 1</label>
          <graphic xlink:href="https://html.scirp.org/file/1771334-rId13.jpeg?20260608020930" />
        </fig>
        <p><bold>Figure 1</bold><bold>.</bold> NTEV profile with Isolation Transformer at Lab 1.</p>
      </sec>
      <sec id="sec1dot2">
        <title>1.2. Objective</title>
        <p>Research on Neutral-to-Earth Voltage (NTEV) in government buildings is of paramount importance, as these structures are subject to stringent safety regulations, possess intricate electrical systems, and have a significant public impact. Conducting such research enables the identification of potential electrical hazards, interference with communication systems, and issues related to legal compliance. A key objective of this research is to evaluate the relationship between Neutral-to-Earth Voltage and Total Harmonic Distortion (THD) in identifying the cause and effect of Neutral-to-Earth Voltage at government buildings. Thereby it aligns the government’s standards with international benchmarks.</p>
        <p>1.2.1. Type of System in This Research</p>
        <p>Recent studies from M.A Shafie highlight the increase in the use of nonlinear loads (like computers, printers, UPS, and electronic lighting) in distribution systems cause harmonic distortion and elevated neutral-to-ground voltages, which can lead to equipment malfunction, higher losses, and reliability issues. This study focuses on Low Voltage system with varies of load that contribute to the NTEV level which is distribution system at government buildings.</p>
        <fig id="fig2">
          <label>Figure 2</label>
          <graphic xlink:href="https://html.scirp.org/file/1771334-rId14.jpeg?20260608020931" />
        </fig>
        <p><bold>Figure 2</bold><bold>.</bold> NTEV profile without Isolation Transformer at Lab 1.</p>
        <p>1.2.2. Parameter</p>
        <p>1) Frequency Stability</p>
        <p>Referring to Malaysia Distribution Code, under clause 5.4.3 system frequency are required to be continuously maintained ±1% from 50 Hz within the range 49.5 Hz to 50.5 Hz this is to make sure the active power output does not decrease by more than 5% [<xref ref-type="bibr" rid="B4">4</xref>].</p>
        <p>2) Harmonic Distortion</p>
        <p>Harmonic distortion is defined as departure of a waveform from sinusoidal shape that is caused by the addition of one or more harmonics to the fundamental waveform. In Malaysian Distribution Code, clause 5.4.6.9 it is stated that Total Harmonic Distortion is to be maintained below 5% from the magnitude of the fundamental wave form [<xref ref-type="bibr" rid="B4">4</xref>]. In principle total harmonic distortion is equal to I harmonic divide by I fundamental [THD = Iн/Iı].</p>
      </sec>
    </sec>
    <sec id="sec2">
      <title>2. Methodology</title>
      <p><xref ref-type="fig" rid="fig3">Figure 3</xref> presented the overall research workflow employed in this study, beginning with site selection and site visits, followed by PQ data logger installation and continuous data measurement and monitoring. The recorded data were then extracted, sorted, processed, and systematically arranged for analysis. Descriptive analysis, correlation analysis, and visual observation were subsequently conducted to evaluate data characteristics and relationships. Finally, the outcomes of these analyses were consolidated and interpreted in the results and discussion stage, marking the completion of the research process.</p>
      <fig id="fig3">
        <label>Figure 3</label>
        <graphic xlink:href="https://html.scirp.org/file/1771334-rId15.jpeg?20260608020932" />
      </fig>
      <p><bold>Figure 3</bold><bold>.</bold> Research methodology flowchart.</p>
      <sec id="sec2dot1">
        <title>2.1. Site Selection &amp; Site Visit</title>
        <p>The site selection and site visit, as well as discussions with the building owner regarding the building’s function and the categories of loads present, are the initial steps of the process.</p>
        <p><bold>Overview</bold><bold>of</bold><bold>THD</bold></p>
        <p>Changes in frequency on each wave whether voltage or current, it can cause the existence of harmonic distortion. This includes neutral voltage and current components [<xref ref-type="bibr" rid="B5">5</xref>]. In a power system, the issue of power quality is very important because it can disrupt the entire system. Among the main components in this determination are harmonics. Therefore, all harmonics is very important for any operation or planning studies.</p>
        <p>According to David from Okanagan College in Kelowna, BC, Canada states that Total Harmonic Distortion (THD) is the ratio for <italic>V</italic><italic><sub>rms</sub></italic> of all the harmonic frequency starting from the 2nd harmonic on, over the <italic>V</italic><italic><sub>rms</sub></italic> fundamental frequency. He also provided the equation to determine the THD value</p>
        <disp-formula id="FD1">
          <mml:math display="inline">
            <mml:mrow>
              <mml:mtext>THD</mml:mtext>
              <mml:mo>=</mml:mo>
              <mml:mfrac>
                <mml:mrow>
                  <mml:msqrt>
                    <mml:mrow>
                      <mml:mstyle displaystyle="true">
                        <mml:msubsup>
                          <mml:mo>∑</mml:mo>
                          <mml:mrow>
                            <mml:mi>n</mml:mi>
                            <mml:mo>=</mml:mo>
                            <mml:mn>2</mml:mn>
                          </mml:mrow>
                          <mml:mi>∞</mml:mi>
                        </mml:msubsup>
                        <mml:mrow>
                          <mml:msubsup>
                            <mml:mi>V</mml:mi>
                            <mml:mrow>
                              <mml:mi>n</mml:mi>
                              <mml:mo>_</mml:mo>
                              <mml:mi>r</mml:mi>
                              <mml:mi>m</mml:mi>
                              <mml:mi>s</mml:mi>
                            </mml:mrow>
                            <mml:mn>2</mml:mn>
                          </mml:msubsup>
                        </mml:mrow>
                      </mml:mstyle>
                    </mml:mrow>
                  </mml:msqrt>
                </mml:mrow>
                <mml:mrow>
                  <mml:msub>
                    <mml:mi>V</mml:mi>
                    <mml:mrow>
                      <mml:mi>f</mml:mi>
                      <mml:mi>u</mml:mi>
                      <mml:mi>n</mml:mi>
                      <mml:mi>d</mml:mi>
                      <mml:mo>_</mml:mo>
                      <mml:mi>r</mml:mi>
                      <mml:mi>m</mml:mi>
                      <mml:mi>s</mml:mi>
                    </mml:mrow>
                  </mml:msub>
                </mml:mrow>
              </mml:mfrac>
            </mml:mrow>
          </mml:math>
        </disp-formula>
        <p>were,</p>
        <p><italic>V</italic><italic><sub>n</sub></italic><italic><sub>_</sub></italic><italic><sub>rms</sub></italic> is the RMS voltage of the nth harmonic;</p>
        <p><italic>V</italic><italic><sub>fund_rms</sub></italic> is the RMS voltage of the fundamental frequency.</p>
        <p><bold>Scope</bold></p>
        <p>Standard and Regulation</p>
        <p>In Malaysia, Total Harmonic Distortion (THD) limits are prescribed through national grid codes and utility-specific technical requirements to ensure power quality and the secure operation of electrical infrastructure. These regulations are primarily documented in the Malaysian Grid Code, the Distribution Code for Peninsular Malaysia, and technical publications by Tenaga Nasional Berhad (TNB), the country’s principal utility provider.</p>
        <p>Voltage THD limits differ based on system voltage levels. According to the Malaysian Grid Code and corroborated by technical analyses, the allowable THD values are as follows:</p>
        <p>High Voltage (HV) Systems (500 kV, 275 kV, 132 kV): Maximum THD is limited to 3.0% [<xref ref-type="bibr" rid="B6">6</xref>].Medium Voltage (MV) Systems (33 kV, 22 kV, 11 kV, 6.6 kV): THD must not exceed 6.5% [<xref ref-type="bibr" rid="B7">7</xref>].Low Voltage (LV) Systems (400 V and below): A THD threshold of 5.0% is applied [<xref ref-type="bibr" rid="B7">7</xref>].</p>
        <p>Maintaining THD within these limits is critical to minimizing harmonic distortion, which can lead to overheating, equipment malfunction, and increased system losses.</p>
        <table-wrap id="tbl2">
          <label>Table 2</label>
          <table>
            <tbody>
              <tr>
                <td>Voltage Level</td>
                <td>Voltage THD Limit</td>
                <td>Current THD Limit(Inverter Output)</td>
              </tr>
              <tr>
                <td>500 kV, 275 kV, 132 kV</td>
                <td>3.0%</td>
                <td>N/A</td>
              </tr>
              <tr>
                <td>33 kV, 22 kV, 11 kV, 6.6 kV</td>
                <td>6.5%</td>
                <td>5.0%</td>
              </tr>
              <tr>
                <td>400 V and below</td>
                <td>5.0%</td>
                <td>5.0%</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
      </sec>
      <sec id="sec2dot2">
        <title>2.2. Type of Analysist</title>
        <p>Two types of analysis are used, which are Descriptive Analysis and Correlation Analysis. </p>
        <p>Descriptive statistical analysis was applied to summarize the fundamental characteristics of the measured variables, including central tendency, variability, and data distribution. This preliminary analysis provides an essential basis for understanding the behavior of Neutral-to-Earth Voltage (NTEV) and Total Harmonic Distortion (THD) prior to inferential evaluation. Pearson correlation analysis was selected because both NTEV and THD are continuous variables, and this method is appropriate for examining the strength and direction of their linear relationship. Consequently, Pearson correlation offers a suitable statistical approach for assessing the association between harmonic distortion and variations in neutral-to-earth voltage.</p>
        <p>The study conducted by Suziyani from the School of Electrical Engineering, College of Engineering, Universiti Teknologi MARA Dungun, Terengganu, demonstrated that the correlation coefficient method is an effective tool for analyzing the relationships between variables [<xref ref-type="bibr" rid="B8">8</xref>].</p>
        <p>Referring to the research by Gabrijela Dimić and colleagues demonstrates the effectiveness of descriptive statistical analysis in educational data mining. By providing detailed insights into data distribution, identifying skewness, and applying histogram discretization, the method improved prediction model accuracy, particularly with the Random Forest classifier, which achieved 91.49% accuracy. The analysis also facilitated the transformation of continuous features into discrete intervals, enhancing algorithm efficiency and model performance. Overall, the study underscores the critical role of descriptive statistics in preprocessing, leading to more accurate and efficient data mining outcomes [<xref ref-type="bibr" rid="B9">9</xref>]. <bold>Table 2</bold>presents the judgement table for the descriptive analysis.</p>
        <p><bold>Table 2</bold><bold>.</bold> Judgement table for descriptive analysis.</p>
        <table-wrap id="tbl3">
          <label>Table 3</label>
          <table>
            <tbody>
              <tr>
                <td>
                  <bold>Measure</bold>
                </td>
                <td>
                  <bold>Description</bold>
                </td>
                <td>
                  <bold>Purpose</bold>
                </td>
                <td>
                  <bold>Interpretation</bold>
                </td>
              </tr>
              <tr>
                <td>
                  <bold>Mean</bold>
                </td>
                <td>The average value of the dataset.</td>
                <td>To identify the central tendency of the data.</td>
                <td>Provides the overall average value of the dataset.</td>
              </tr>
              <tr>
                <td>
                  <bold>Median</bold>
                </td>
                <td>The middle value when the data is arranged in order.</td>
                <td>To determine the central point of the data.</td>
                <td>Useful for skewed distributions or outliers.</td>
              </tr>
              <tr>
                <td>
                  <bold>Mode</bold>
                </td>
                <td>The value that appears most frequently in the dataset.</td>
                <td>To understand the most common value.</td>
                <td>Indicates the most frequent observation.</td>
              </tr>
              <tr>
                <td>
                  <bold>Standard</bold>
                  <bold>Deviation</bold>
                </td>
                <td>The square root of variance, showing data spread in the same units as the data.</td>
                <td>To gauge the dispersion or consistency in the data.</td>
                <td>Helps in understanding the consistency of data.</td>
              </tr>
              <tr>
                <td>
                  <bold>Skewness</bold>
                </td>
                <td>Measures the asymmetry of the data distribution.</td>
                <td>To assess the direction of skew in the data.</td>
                <td>Positive skew: right tail longer; Negative: left tail longer.</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <p><bold>Table 3</bold><bold>.</bold> Pearson correlation coefficient guide table.</p>
        <table-wrap id="tbl4">
          <label>Table 4</label>
          <table>
            <tbody>
              <tr>
                <td>Correlation Coefficient Value</td>
                <td>Interpretation Criteria</td>
              </tr>
              <tr>
                <td>−0.800 to −1.000</td>
                <td>Perfectly Negative</td>
              </tr>
              <tr>
                <td>−0.600 to −0.799</td>
                <td>Strongly Negative</td>
              </tr>
              <tr>
                <td>−0.400 to −0.599</td>
                <td>Moderately Negative</td>
              </tr>
              <tr>
                <td>−0.200 to −0.399</td>
                <td>Weakly Negative</td>
              </tr>
              <tr>
                <td>0 to −0.199</td>
                <td>Very Weakly Negative</td>
              </tr>
              <tr>
                <td>0.000 to 0.199</td>
                <td>Very Weakly Positive</td>
              </tr>
              <tr>
                <td>0.200 to 0.399</td>
                <td>Weakly Positive</td>
              </tr>
              <tr>
                <td>0.400 to 0.599</td>
                <td>Moderately Positive</td>
              </tr>
              <tr>
                <td>0.600 to 0.799</td>
                <td>Strongly Positive</td>
              </tr>
              <tr>
                <td>0.800 to 1.000</td>
                <td>Perfectly Positive</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <p><bold>Table 3</bold> shows a guide to interpreting the Pearson Correlation Coefficient, which measures the strength and direction of the linear relationship between two continuous variables. In practice, the closer the coefficient is to either +1 or −1, the stronger the linear relationship. Values closer to 0 indicate weaker relationships, and the direction of the relationship (positive or negative) is indicated by the sign of the coefficient [<xref ref-type="bibr" rid="B10">10</xref>].</p>
      </sec>
      <sec id="sec2dot3">
        <title>2.3. PQ Data Logger Installation</title>
        <p>In this phase, several measurement locations have been identified for the purpose of acquiring raw data prior to analysis. <xref ref-type="fig" rid="fig4">Figure 4</xref> below illustrates the point diagram, where a Fluke 435 PQ as shows in <xref ref-type="fig" rid="fig5">Figure 5</xref> meter will be installed on the distribution network to collect data according to the desired parameters. Each location will be installed within a period of four days.</p>
        <fig id="fig4">
          <label>Figure 4</label>
          <graphic xlink:href="https://html.scirp.org/file/1771334-rId18.jpeg?20260608020934" />
        </fig>
        <p><bold>Figure 4</bold><bold>.</bold> Point of connection.</p>
        <fig id="fig5">
          <label>Figure 5</label>
          <graphic xlink:href="https://html.scirp.org/file/1771334-rId19.jpeg?20260608020934" />
        </fig>
        <p><bold>Figure 5</bold><bold>.</bold> PQ data logger fluke 435.</p>
      </sec>
      <sec id="sec2dot4">
        <title>2.4. Analytical and Statistical</title>
        <p>This research will utilize statistical methods, such as the Pearson correlation coefficient, to analyze the relationship between frequency, THD, and Voltage variation with NTEV levels. These analyses confirm that as frequency changes, especially with harmonic distortion, there is a significant correlation with fluctuations in NTEV levels. Consequently, Minitab software as shown in <xref ref-type="fig" rid="fig6">Figure 6</xref> was selected as the optimal tool due to its user-friendly interface and extensive statistical capabilities. Within Minitab, two distinct types of analytical methods were employed:</p>
        <p>1) Descriptive Statistical Analysis</p>
        <p>2) Pearson Correlation Analysis</p>
        <fig id="fig6">
          <label>Figure 6</label>
          <graphic xlink:href="https://html.scirp.org/file/1771334-rId20.jpeg?20260608020934" />
        </fig>
        <p><bold>Figure 6</bold><bold>.</bold> Power log software interface</p>
      </sec>
      <sec id="sec2dot5">
        <title>2.5. Visual Observation</title>
        <p>A visual observation for Neutral-to-Earth Voltage (NTEV) graft and Total Harmonic Distortion (THD) graft typically involves comparing two separate graphical representations to assess their respective trends and relationships. </p>
        <p>By comparing these graphs, one can observe if there’s a correlation between NTEV variations and spikes in THD. For instance, an increase in NTEV may coincide with a rise in THD, suggesting that grounding issues could contribute to higher harmonic distortion in the system.</p>
        <p>This type of analysis helps to diagnose power quality problems and take corrective actions to improve the overall stability and performance of the electrical system.</p>
      </sec>
    </sec>
    <sec id="sec3">
      <title>3. Result</title>
      <sec id="sec3dot1">
        <title>3.1. Neutral to Earth Voltage</title>
        <p>Observation</p>
        <p>From the result obtained in <xref ref-type="fig" rid="fig7">Figure 7</xref>, sine wave view, it can be seen that the NTEV pattern is directly proportional to the current value in each phase. High values occur during peak hours while low voltage values occur outside of peak hours.</p>
        <fig id="fig7">
          <label>Figure 7</label>
          <graphic xlink:href="https://html.scirp.org/file/1771334-rId21.jpeg?20260608020936" />
        </fig>
        <p><bold>Figure 7</bold><bold>.</bold> Neutral to earth voltage at MSB1.</p>
        <p><bold>Table 4</bold><bold>.</bold> Minimum and maximum NTEV for MSB1.</p>
        <table-wrap id="tbl5">
          <label>Table 5</label>
          <table>
            <tbody>
              <tr>
                <td>Location</td>
                <td>Neutral Voltage (V)</td>
                <td>Date</td>
                <td>Time</td>
              </tr>
              <tr>
                <td rowspan="4">MSB1</td>
                <td rowspan="3">Min: 0.32</td>
                <td>23/11/2022</td>
                <td rowspan="3">Mostly during off-peak hours from 6.30 pm-7.00 am</td>
              </tr>
              <tr>
                <td>24/11/2022</td>
              </tr>
              <tr>
                <td>25/11/2022</td>
              </tr>
              <tr>
                <td>Max: 1.34</td>
                <td>23/11/2022</td>
                <td>Mostly during peak hours from 7.00 am-6.30 pm</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <p>The Neutral to Earth Voltage’s minimum and maximum values are displayed in <bold>Table 4</bold>. During the recorded measurement period November 22nd to 25th 2022 at the first site, the minimal NTEV value of 0.32 V occurred every day during off-peak hours. The maximum value of 1.34 V happened on November 23, 2022, at 10 am, during peak hour.</p>
      </sec>
      <sec id="sec3dot2">
        <title>3.2. Total Harmonic Distortion (THD)</title>
        <fig id="fig8">
          <label>Figure 8</label>
          <graphic xlink:href="https://html.scirp.org/file/1771334-rId22.jpeg?20260608020936" />
        </fig>
        <p><bold>Figure 8</bold><bold>.</bold> THD for MSB1.</p>
        <p>The percentage THD value at the first location is shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>. The measurement began on November 22, 2022, at 10.40 a.m. and ended on November 25, 2022, at 9.16 a.m. Throughout the measurement record, the THD value remains within the specification range of less than 5%. When compared to the current in each phase, the pattern is nearly same, with high value at peak time and low value at non-peak. This is because VTHD and total current are depended on each other.</p>
        <p><bold>Table 5</bold><bold>.</bold> THD at MSB1.</p>
        <table-wrap id="tbl6">
          <label>Table 6</label>
          <table>
            <tbody>
              <tr>
                <td>Phase</td>
                <td>Min (%)</td>
                <td>Date</td>
                <td>Time</td>
                <td>Max (%)</td>
                <td>Date</td>
                <td>Time</td>
              </tr>
              <tr>
                <td>R-N</td>
                <td>1.38</td>
                <td>24/11/2022</td>
                <td>12.51 am</td>
                <td>2.58</td>
                <td>25/11/2022</td>
                <td>9.14 am</td>
              </tr>
              <tr>
                <td>Y-N</td>
                <td>1.2</td>
                <td>24/11/2022</td>
                <td>1.10 am</td>
                <td>2.95</td>
                <td>24/11/2022</td>
                <td>10.13 am</td>
              </tr>
              <tr>
                <td>B-N</td>
                <td>1.19</td>
                <td>23/11/2022</td>
                <td>7.12 am</td>
                <td>2.82</td>
                <td>24/11/2022</td>
                <td>10.02 am</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <p>According to <bold>Table 5</bold> above, for the Red phase, the lowest value is 1.38% and occurs on November 24, 2022, at 12.51 am, while the highest value occurs at peak time on November 25, 2022, at 9.15 am. On November 24, 2022, at 1.10 am, the yellow phase recorded a lowest value of 1.2%, and at 10.13 am the same day, a maximum value of 2.95%. One of the lowest VTHD values ever recorded was in the blue phase on November 23, 2022, at 7.12 am, while the highest value was at peak time on November 24, 2022, at 10.02 am.</p>
      </sec>
    </sec>
    <sec id="sec4">
      <title>4. Analysis</title>
      <sec id="sec4dot1">
        <title>4.1. Descriptive Statistics</title>
        <p><bold>Table 6</bold><bold>.</bold> Statistic of MSB1.</p>
        <table-wrap id="tbl7">
          <label>Table 7</label>
          <table>
            <tbody>
              <tr>
                <td>
                  <bold>Variable</bold>
                </td>
                <td>
                  <bold>Mean</bold>
                </td>
                <td>
                  <bold>Std.</bold>
                  <bold>Dev.</bold>
                </td>
                <td>
                  <bold>Minimum</bold>
                </td>
                <td>
                  <bold>Median</bold>
                </td>
                <td>
                  <bold>Maximum</bold>
                </td>
                <td>
                  <bold>Skewness</bold>
                </td>
              </tr>
              <tr>
                <td>NTEV</td>
                <td>0.67</td>
                <td>0.31</td>
                <td>0.32</td>
                <td>0.54</td>
                <td>1.34</td>
                <td>0.46</td>
              </tr>
              <tr>
                <td>VTHD R</td>
                <td>1.85</td>
                <td>0.28</td>
                <td>1.38</td>
                <td>1.86</td>
                <td>2.58</td>
                <td>0.08</td>
              </tr>
              <tr>
                <td>VTHD Y</td>
                <td>1.96</td>
                <td>0.50</td>
                <td>1.20</td>
                <td>1.82</td>
                <td>2.95</td>
                <td>0.23</td>
              </tr>
              <tr>
                <td>VTHD B</td>
                <td>1.82</td>
                <td>0.38</td>
                <td>1.19</td>
                <td>1.76</td>
                <td>2.82</td>
                <td>0.17</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <p>NTEV</p>
        <p>Referring to the <bold>Table 6</bold> Statistic of MSB1, NTEV std deviation is 0.31, and mean 0.67. By this two information, the value of standard deviation provides evidence that Neutral to Earth voltage across neutral to ground line deviates very closely to the mean value. While skewness with the value of 0.46 and median of 0.54 provide informative evidence on how the measured data was distributed during measurement period and in this case is skewed to the right. As a result, the data might be interpreted as NTEV voltages tend to be greater than the average value.</p>
        <p>VTHD RYB</p>
        <p>Overall, the VTHD value on each phase is in the range of 1.19% to 2.95%. By this information, the value of standard deviation provides evidence that VTHD each phase deviates very closely to the mean value. While skewness with the value of 0.08, 0.23 and 0.17, provide informative evidence on how the measured data was distributed during measurement period. The data measured in this case is positively skewed. As a result, the data might be interpreted as VTHD each phase tend to be greater than the average value.</p>
      </sec>
      <sec id="sec4dot2">
        <title>4.2. Analysis</title>
        <p><bold>Table 7</bold><bold>.</bold> NTEV and VTHD correlation MSB1.</p>
        <table-wrap id="tbl8">
          <label>Table 8</label>
          <table>
            <tbody>
              <tr>
                <td>
                </td>
                <td>
                  <bold>NTEV</bold>
                </td>
                <td>
                  <bold>VTHD</bold>
                  <bold>Red</bold>
                </td>
                <td>
                  <bold>VTHD</bold>
                  <bold>Yellow</bold>
                </td>
              </tr>
              <tr>
                <td>VTHD Red</td>
                <td>0.845</td>
                <td>
                </td>
                <td>
                </td>
              </tr>
              <tr>
                <td>VTHD Yellow</td>
                <td>0.919</td>
                <td>0.956</td>
                <td>
                </td>
              </tr>
              <tr>
                <td>VTHD Blue</td>
                <td>0.905</td>
                <td>0.952</td>
                <td>0.974</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <p>Interpretation:</p>
        <p><italic>NTEV</italic><italic>and</italic><italic>VTHD</italic><italic>correlat</italic><italic>e</italic><italic>to</italic><italic>each</italic><italic>other</italic>. <italic>Referring</italic><italic>to</italic><bold>Table</bold><bold>7</bold><italic>,</italic><italic>the</italic><italic>analysis</italic><italic>of</italic><italic>Pearson</italic><italic>correlation</italic><italic>coefficient</italic><italic>between</italic><italic>NTEV</italic><italic>to</italic><italic>VTHD</italic><italic>Red</italic>, <italic>VTHD</italic><italic>Yellow</italic><italic>and</italic><italic>VTHD</italic><italic>Blue</italic><italic>are</italic> 0.845, 0.919 <italic>and</italic> 0.905. <italic>According</italic><italic>to</italic><italic>the</italic><italic>Pearson</italic><italic>correlation</italic><italic>guideline</italic>, <italic>a</italic><italic>coefficient</italic><italic>with</italic> 0.7 <italic>to</italic> 0.9 <italic>can</italic><italic>be</italic><italic>interpreted</italic><italic>as</italic><italic>highly</italic><italic>positive</italic><italic>correlation</italic>. <italic>In</italic><italic>this</italic><italic>case</italic><italic>the</italic><italic>NTEV</italic><italic>is</italic><italic>slightly</italic><italic>dependent</italic><italic>on</italic><italic>VTHD</italic><italic>Red</italic><italic>compared</italic><italic>with</italic><italic>VTHD</italic><italic>Yellow</italic><italic>and</italic><italic>Blue</italic>.</p>
        <p>While the size of correlation is more than 0.9, the interpretation is very high correlation. Hence the result can be interpreted as NTEV is dependent to VTHD Red, Yellow and Blue, as VTHD value increase the voltage at Neutral increase. This is due to the effects of harmonics in which it changes the base frequency of the voltage hence changes the voltage magnitude. While within spec value of THD impose limited problem to line stability, severe THD beyond the specification imposes potent risk to the system stability.</p>
      </sec>
    </sec>
    <sec id="sec5">
      <title>5. Discussion</title>
      <p><bold>Table 8</bold><bold>.</bold> Summarize value of VTHD and correlation for all location.</p>
      <table-wrap id="tbl9">
        <label>Table 9</label>
        <table>
          <tbody>
            <tr>
              <td>Location</td>
              <td>Average Max. VTHD (%)</td>
              <td colspan="2">Size of correlation between NTEV &amp; VTHD</td>
            </tr>
            <tr>
              <td>MSB1</td>
              <td>2.78</td>
              <td>0.8 - 0.9</td>
              <td>High positive</td>
            </tr>
          </tbody>
        </table>
      </table-wrap>
      <p>The correlation analysis reveals a strong to very strong positive relationship between NTEV and VTHD across all phases at MSB1(see <bold>Table 8</bold>). This finding indicates that harmonic distortion plays a significant role in influencing NTEV levels in low-voltage distribution systems. From a practical perspective, this suggests that effective harmonic mitigation strategies, such as the installation of passive or active harmonic filters, can contribute to the reduction of NTEV and improve overall power quality in government buildings.</p>
      <p>It should be noted, however, that harmonics are not the sole contributors to elevated NTEV. Other factors, including load imbalance, grounding resistance, and system configuration, may also influence NTEV levels. These aspects were not explicitly analyzed in this study and should be considered in future investigations to provide a more comprehensive understanding of NTEV behavior.</p>
      <p>The study is subject to several limitations, including measurements taken at a single location and over a relatively short monitoring period. Future work should involve multiple measurement points and longer observation durations to capture broader operating conditions.</p>
    </sec>
    <sec id="sec6">
      <title>Acknowledgements</title>
      <p>The corresponding author, on behalf of all authors, declares that there are no conflicts of interest. The authors thank the Faculty of Electrical Engineering Technology (FTKE) and the Centre for Research and Innovation Management (CRIM), Universiti Teknikal Malaysia Melaka (UTeM), for their support and research facilities. Special thanks to Agensi Nuklear Malaysia, Bangi, Selangor, for providing data and the technical assistance.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <title>References</title>
      <ref id="B1">
        <label>1.</label>
        <citation-alternatives>
          <mixed-citation publication-type="web">IEEE (2019) Recommended Practice for Monitoring Electric Power Quality. http://www.ieee.org/web/aboutus/whatis/policies/p9-26.html</mixed-citation>
          <element-citation publication-type="web">
            <year>2019</year>
            <article-title>Recommended Practice for Monitoring Electric Power Quality</article-title>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B2">
        <label>2.</label>
        <citation-alternatives>
          <mixed-citation publication-type="confproc">Shafie, M.A., Singh, H. and Rahman, M.Q.A. (2010). Harmonic and Neutral to Ground Voltage Reduction Using Isolation Transformer. 2010 IEEE International Conference on Power and Energy. Kuala Lumpur, 29 November-1 December 2010, 561-566. https://doi.org/10.1109/PECON.2010.5697645 <pub-id pub-id-type="doi">10.1109/PECON.2010.5697645</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/PECON.2010.5697645">https://doi.org/10.1109/PECON.2010.5697645</ext-link></mixed-citation>
          <element-citation publication-type="confproc">
            <person-group person-group-type="author">
              <string-name>Shafie, M.A.</string-name>
              <string-name>Singh, H.</string-name>
              <string-name>Rahman, M.Q.A.</string-name>
            </person-group>
            <year>2010</year>
            <pub-id pub-id-type="doi">10.1109/PECON.2010.5697645</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B3">
        <label>3.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Mirus International Inc (2021) A Practical and Effective Way of Applying IEEE Std 519-2014 Harmonic Limits.</mixed-citation>
          <element-citation publication-type="other">
            <year>2021</year>
            <article-title>A Practical and Effective Way of Applying IEEE Std 519-2014 Harmonic Limits</article-title>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B4">
        <label>4.</label>
        <citation-alternatives>
          <mixed-citation publication-type="web">Malaysia, S.T. (2017) Distribution Code for Peninsular Malaysia, Sabah &amp; F.T. Labuan (Amendments) 2017. https://www.st.gov.my/</mixed-citation>
          <element-citation publication-type="web">
            <person-group person-group-type="author">
              <string-name>Malaysia, S.T.</string-name>
              <string-name>Malaysia, S</string-name>
            </person-group>
            <year>2017</year>
            <article-title>Distribution Code for Peninsular Malaysia, Sabah &amp; F</article-title>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B5">
        <label>5.</label>
        <citation-alternatives>
          <mixed-citation publication-type="confproc">Variz, A.M., Pereira, J.L.R., Carneiro, S. and Barbosa, P.G. (2009) Harmonic Analysis of the Power Distribution Neutral-To-Earth Voltage (NEV) Test Case Using Four-Wire Three-Phase Harmonic Current Injection Method. 2009 IEEE Power &amp; Energy Society General Meeting. Calgary, 26-30 July 2009, 1-7. https://doi.org/10.1109/PES.2009.5275758 <pub-id pub-id-type="doi">10.1109/PES.2009.5275758</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/PES.2009.5275758">https://doi.org/10.1109/PES.2009.5275758</ext-link></mixed-citation>
          <element-citation publication-type="confproc">
            <person-group person-group-type="author">
              <string-name>Variz, A.M.</string-name>
              <string-name>Pereira, J.L.R.</string-name>
              <string-name>Carneiro, S.</string-name>
              <string-name>Barbosa, P.G.</string-name>
            </person-group>
            <year>2009</year>
            <article-title>Harmonic Analysis of the Power Distribution Neutral-To-Earth Voltage (NEV) Test Case Using Four-Wire Three-Phase Harmonic Current Injection Method</article-title>
            <source>2009 IEEE Power &amp; Energy Society General Meeting. Calgary</source>
            <volume>26</volume>
            <pub-id pub-id-type="doi">10.1109/PES.2009.5275758</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B6">
        <label>6.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Johari, A. and St, J. (2010) Document Control the Malaysian Grid Code.</mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Johari, A.</string-name>
              <string-name>St, J.</string-name>
            </person-group>
            <year>2010</year>
            <article-title>Document Control the Malaysian Grid Code</article-title>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B7">
        <label>7.</label>
        <citation-alternatives>
          <mixed-citation publication-type="web">Al-Shetwi, A.Q. and Sujod, M.Z. (2018) Harmonic Distortion and Voltage Imbalance Study of Photovoltaic Power Plant Connected to the Malaysian Grid Muhamad Zahim Sujod Harmonic Distortion and Voltage Imbalance Study of Photovoltaic Power Plant Connected to the Malaysian Grid. https://www.researchgate.net/publication/323642562</mixed-citation>
          <element-citation publication-type="web">
            <person-group person-group-type="author">
              <string-name>Al-Shetwi, A.Q.</string-name>
              <string-name>Sujod, M.Z.</string-name>
            </person-group>
            <year>2018</year>
            <article-title>Harmonic Distortion and Voltage Imbalance Study of Photovoltaic Power Plant Connected to the Malaysian Grid Muhamad Zahim Sujod Harmonic Distortion and Voltage Imbalance Study of Photovoltaic Power Plant Connected to the Malaysian Grid</article-title>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B8">
        <label>8.</label>
        <citation-alternatives>
          <mixed-citation publication-type="confproc">Rohafauzi, S., Kassim, M., Ja’Afar, H., Rustam, I., Miskon, M.T. and Pakhrudin, N.S.M. (2023) Descriptive Analysis on Data of Stingless Bee IoT Application Monitoring System. 2023 <italic>IEEE</italic> 13 <italic>th International Conference on System Engineering and Technology</italic>( <italic>ICSET</italic>), Shah Alam, 2 October 2023, 159-164. https://doi.org/10.1109/icset59111.2023.10295071 <pub-id pub-id-type="doi">10.1109/icset59111.2023.10295071</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/icset59111.2023.10295071">https://doi.org/10.1109/icset59111.2023.10295071</ext-link></mixed-citation>
          <element-citation publication-type="confproc">
            <person-group person-group-type="author">
              <string-name>Rohafauzi, S.</string-name>
              <string-name>Kassim, M.</string-name>
              <string-name>Afar, H.</string-name>
              <string-name>Rustam, I.</string-name>
              <string-name>Miskon, M.T.</string-name>
              <string-name>Pakhrudin, N.S.M.</string-name>
            </person-group>
            <year>2023</year>
            <article-title>Descriptive Analysis on Data of Stingless Bee IoT Application Monitoring System</article-title>
            <source>2023 IEEE 13th International Conference on System Engineering and Technology (ICSET)</source>
            <volume>2</volume>
            <pub-id pub-id-type="doi">10.1109/icset59111.2023.10295071</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B9">
        <label>9.</label>
        <citation-alternatives>
          <mixed-citation publication-type="confproc">Dimić, G., Rancic, D., Pronić-Rančić, O., <italic>et al</italic>. (2019) Descriptive Statistical Analysis in the Process of Educational Data Mining. 2019 14 <italic>th International Conference on Advanced Technologies</italic>, <italic>Systems and Services in Telecommunications</italic>, Niš Serbia, 23-25 October 2019, 388-391.</mixed-citation>
          <element-citation publication-type="confproc">
            <person-group person-group-type="author">
              <string-name>Rancic, D.</string-name>
              <string-name>Technologies, S</string-name>
              <string-name>Telecommunications, N</string-name>
            </person-group>
            <year>2019</year>
            <article-title>Descriptive Statistical Analysis in the Process of Educational Data Mining</article-title>
            <source>2019 14th International Conference on Advanced Technologies</source>
            <volume>23</volume>
          </element-citation>
        </citation-alternatives>
      </ref>
      <ref id="B10">
        <label>10.</label>
        <citation-alternatives>
          <mixed-citation publication-type="other">Maulana, B.W., Utomo, S.B., Anam, K. and Bilfaqih, I. (2023) Analysis of Movement Patterns Based on Electroencephalograph Signals Using Pearson Correlation Coefficient. 2023 <italic>International Seminar on Intelligent Technology and Its Applications</italic> ( <italic>ISITIA</italic>), Surabaya, 26 July 2023, 593-597. https://doi.org/10.1109/isitia59021.2023.10220971 <pub-id pub-id-type="doi">10.1109/isitia59021.2023.10220971</pub-id><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/isitia59021.2023.10220971">https://doi.org/10.1109/isitia59021.2023.10220971</ext-link></mixed-citation>
          <element-citation publication-type="other">
            <person-group person-group-type="author">
              <string-name>Maulana, B.W.</string-name>
              <string-name>Utomo, S.B.</string-name>
              <string-name>Anam, K.</string-name>
              <string-name>Bilfaqih, I.</string-name>
            </person-group>
            <year>2023</year>
            <article-title>Analysis of Movement Patterns Based on Electroencephalograph Signals Using Pearson Correlation Coefficient</article-title>
            <source>2023 International Seminar on Intelligent Technology and Its Applications (ISITIA)</source>
            <volume>26</volume>
            <pub-id pub-id-type="doi">10.1109/isitia59021.2023.10220971</pub-id>
          </element-citation>
        </citation-alternatives>
      </ref>
    </ref-list>
  </back>
</article>