예제: 비선형 회귀 1

genfit 함수를 사용하여 일반 비선형 적합식의 적합 매개변수를 찾습니다.

1. 데이터 집합을 정의합니다.

2. genfit 함수에 제공할 함수를 만듭니다. 적합식을 입력한 다음 각 매개변수를 고려하여 도함수를 추가합니다. 도함수에는 e 대신 exp 함수를 사용합니다.

매개변수는 지수 연산자를 사용하여 매개변수 벡터의 요소로 지정해야 합니다.

3. 추측값을 입력합니다.

4. genfit 함수를 호출합니다.

<region id="ID0ENA3S" top="883.20000000000016" left="19.200000000000003" height="25.6" width="192.57666666666668" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="52">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <id labels="*" xml:space="preserve">β</id>
      <apply>
        <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">genfit</id>
        <sequence>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">vx</id>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">vy</id>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">guess</id>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">F</id>
        </sequence>
      </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" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

5. 데이터 집합과 일반 비선형 적합식을 도표화합니다.

<region id="ID0ELC3S" top="1132.8000000000002" left="19.200000000000003" height="36.800000000000004" width="145.70333333333335" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="56">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <function>
        <id labels="*" xml:space="preserve">GNfit</id>
        <boundVars>
          <id labels="*" xml:space="preserve">x</id>
        </boundVars>
      </function>
      <apply>
        <indexer />
        <apply>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">F</id>
          <sequence>
            <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">x</id>
            <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">β</id>
          </sequence>
        </apply>
        <real>0</real>
      </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" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

<region id="ID0EZC3S" top="1152.0000000000002" left="0" height="312" width="622" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="58">
          <traceStyle color="#FFFF0000" symbol="circle" line-weight="1" line-style="None">lines</traceStyle>
        </trace>
        <trace resultRef="59" num-of-points="500">
          <traceStyle color="#FF0000FF" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="464.4" height="216" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="1.3">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="69.13" />
          <axisGrid>
            <gridFrequency>14</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 resultRef="64">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">vx</id>
              </math>
              <math resultRef="65">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="66">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">x</id>
              </math>
              <math resultRef="67">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="61">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <endValue resultRef="63">
              <real xmlns="http://schemas.mathsoft.com/math50">1.2</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0.05" end="2.8">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="117.563333333333" />
          <axisGrid>
            <gridFrequency>12</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 resultRef="68">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">vy</id>
              </math>
              <math resultRef="69">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="70">
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id xml:space="preserve">GNfit</id>
                  <id xml:space="preserve">x</id>
                </apply>
              </math>
              <math resultRef="71">
                <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>

6. corr 함수를 사용하여 vy와 GNfit 간의 상관 관계를 확인합니다.

<region id="ID0ECF3S" top="1574.4000000000003" left="38.400000000000006" height="25.6" width="179.84333333333333" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="82">
    <eval xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">corr</id>
        <sequence>
          <apply>
            <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">GNfit</id>
            <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">vx</id>
          </apply>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">vy</id>
        </sequence>
      </apply>
      <unitOverride>
        <placeholder />
      </unitOverride>
    </eval>
    <resultFormat>
      <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

위의 결과에서 알 수 있듯이 적합도가 완벽합니다.

참조

이 예제의 데이터 출처: Lanczos, C., Applied Analysis, Prentice Hall, 1956년, 272-280페이지. 이 문서는 NIST Statistical Reference Dataset Archive에 게시되어 있습니다.