例: 尖度および歪度

関数 > 統計 > 記述統計 > 例: 尖度および歪度

例: 尖度および歪度

kurt関数とskew関数を使用して、ワイブル分布と正規分布の形状を比較します。

1. 次の分布を定義してプロットします。

<region id="ID0EHHWAB" top="192" left="24" width="600" height="219" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot origin-positioning="false">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="10">
          <traceStyle color="#FF00008B" symbol="none" line-weight="1" line-style="Solid">columns</traceStyle>
        </trace>
        <trace resultRef="11">
          <traceStyle color="#FFF79646" symbol="none" line-weight="1" line-style="Solid">columns</traceStyle>
        </trace>
      </traces>
      <axes>
        <xAxis legend-position="PlotBoundaryBottom" rank="1">
          <axisLine position="bottom" positionticmark="-1" />
          <axisGrid>
            <gridFrequency>11</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <matcol />
                  <id>WH</id>
                  <real>0</real>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <matcol />
                  <id>NH</id>
                  <real>0</real>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis legend-position="PlotBoundaryLeft" rank="1">
          <axisLine position="ticknumberlock" positionticmark="2" />
          <axisGrid>
            <gridFrequency>11</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <matcol />
                  <id>WH</id>
                  <real>1</real>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math>
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <matcol />
                  <id>NH</id>
                  <real>1</real>
                </apply>
              </math>
              <math>
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

2. 分布の尖度を計算して分布のとがり具合を比較します。

ヒストグラムの区間が原因となって正規分布の尖度がゼロからわずかにずれています。ワイブル分布の尖度から、ワイブル分布の方が正規分布よりもわずかにとがっていることがわかります。

次の標本尖度係数によって母集団の尖度を推定できます。

<region id="ID0EYKWAB" top="672" left="48" width="679.516666666667" height="80" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math>
    <define xmlns="http://schemas.mathsoft.com/math50">
      <function>
        <id label="None">kurtosis</id>
        <boundVars>
          <id label="None">v</id>
        </boundVars>
      </function>
      <apply>
        <minus />
        <apply>
          <div />
          <apply>
            <mult />
            <apply>
              <mult />
              <apply>
                <id>rows</id>
                <id>v</id>
              </apply>
              <parens>
                <apply>
                  <plus />
                  <apply>
                    <id>rows</id>
                    <id>v</id>
                  </apply>
                  <real>1</real>
                </apply>
              </parens>
            </apply>
            <apply>
              <summation />
              <lambda>
                <boundVars>
                  <placeholder />
                </boundVars>
                <apply>
                  <pow />
                  <parens>
                    <apply>
                      <minus />
                      <id>v</id>
                      <apply>
                        <id>mean</id>
                        <id>v</id>
                      </apply>
                    </apply>
                  </parens>
                  <real>4</real>
                </apply>
              </lambda>
              <upperBound>
                <placeholder />
              </upperBound>
            </apply>
          </apply>
          <apply>
            <mult />
            <apply>
              <mult />
              <apply>
                <mult />
                <parens>
                  <apply>
                    <minus />
                    <apply>
                      <id>rows</id>
                      <id>v</id>
                    </apply>
                    <real>1</real>
                  </apply>
                </parens>
                <parens>
                  <apply>
                    <minus />
                    <apply>
                      <id>rows</id>
                      <id>v</id>
                    </apply>
                    <real>2</real>
                  </apply>
                </parens>
              </apply>
              <parens>
                <apply>
                  <minus />
                  <apply>
                    <id>rows</id>
                    <id>v</id>
                  </apply>
                  <real>3</real>
                </apply>
              </parens>
            </apply>
            <apply>
              <pow />
              <apply>
                <id>Stdev</id>
                <id>v</id>
              </apply>
              <real>4</real>
            </apply>
          </apply>
        </apply>
        <apply>
          <div />
          <apply>
            <mult />
            <real>3</real>
            <apply>
              <pow />
              <parens>
                <apply>
                  <minus />
                  <apply>
                    <id>rows</id>
                    <id>v</id>
                  </apply>
                  <real>1</real>
                </apply>
              </parens>
              <real>2</real>
            </apply>
          </apply>
          <apply>
            <mult />
            <parens>
              <apply>
                <minus />
                <apply>
                  <id>rows</id>
                  <id>v</id>
                </apply>
                <real>2</real>
              </apply>
            </parens>
            <parens>
              <apply>
                <minus />
                <apply>
                  <id>rows</id>
                  <id>v</id>
                </apply>
                <real>3</real>
              </apply>
            </parens>
          </apply>
        </apply>
      </apply>
    </define>
    <resultFormat>
      <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
      <matrix size="12,12" offset="0,0" show-indices="false" />
    </resultFormat>
  </math>
</region>

3. 分布の歪度を比較します。

ワイブル分布は右に漸減しています。

次の標本歪度係数によって母集団の歪度を推定できます。

<region id="ID0ECQWAB" top="1152" left="24" width="379.883333333333" height="80" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math>
    <define xmlns="http://schemas.mathsoft.com/math50">
      <function>
        <id label="None">skewness</id>
        <boundVars>
          <id label="None">v</id>
        </boundVars>
      </function>
      <apply>
        <div />
        <apply>
          <mult />
          <apply>
            <id>rows</id>
            <id>v</id>
          </apply>
          <apply>
            <summation />
            <lambda>
              <boundVars>
                <placeholder />
              </boundVars>
              <apply>
                <pow />
                <parens>
                  <apply>
                    <minus />
                    <id>v</id>
                    <apply>
                      <id>mean</id>
                      <id>v</id>
                    </apply>
                  </apply>
                </parens>
                <real>3</real>
              </apply>
            </lambda>
            <upperBound>
              <placeholder />
            </upperBound>
          </apply>
        </apply>
        <apply>
          <mult />
          <apply>
            <mult />
            <parens>
              <apply>
                <minus />
                <apply>
                  <id>rows</id>
                  <id>v</id>
                </apply>
                <real>1</real>
              </apply>
            </parens>
            <parens>
              <apply>
                <minus />
                <apply>
                  <id>rows</id>
                  <id>v</id>
                </apply>
                <real>2</real>
              </apply>
            </parens>
          </apply>
          <apply>
            <pow />
            <apply>
              <id>Stdev</id>
              <id>v</id>
            </apply>
            <real>3</real>
          </apply>
        </apply>
      </apply>
    </define>
    <resultFormat>
      <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
      <matrix size="12,12" offset="0,0" show-indices="false" />
    </resultFormat>
  </math>
</region>