ｍｍｍ＋＋＋ - 経済統計練習帳

import math 
 
def dg(m,a1):
 
    h=0.01
 
    g1=sumg(a1)
 
    a2=[]
    for ap in a1:
        a2.append(ap)
 
    a2[m]=a2[m]+h
 
    g2=sumg(a2)
 
    return (g2-g1)/h
 
 
def sumg(a1):
 
    z=0
    for i in range(len(x)):
        z=z+(y[i]-g(a1,x[i]))*(y[i]-g(a1,x[i]))    
 
    return z
 
def g(a1,x1):
 
    f1=f(a1,x1)
    g1=1/(1+math.exp(-f1))
 
    return g1
 
def f(a1,x1):
    f1=a1[0]
    for i in range(len(x1)):
        f1=f1+a1[i+1]*x1[i]
    return f1
 
y=[0,0,0,1,1]
x=[[5,2],[7,4],[5,6],[1,8],[2,10]]
a=[0.1,0.2,0.3]
 
g1=dg(2,a)
 
print(g1)