えろｊｊ - 経済統計練習帳

import numpy as np
import math
import nn
 
number=100
number_eq=5
number_var=2
 
x=[]
for i in range(100):
   r1=np.random.random()
   r2=np.random.random()
   x.append([r1,r2])
 
y=[]
for i in range(100):
    r3=0
    h=0
    if 3*x[i][0]+2*x[i][1]>2.5:h=h+1  
    if h==1:r3=1
    y.append(r3)
 
link=[[3,5,6],[4,5,6],[1,3,4],[2,3,4],[0,1,2]]
 
a=np.zeros([number_eq,number_var+1])
alpha=np.zeros([number_eq,number_var+1])
maxc=0
for de in range(10):
 
 
    for eq1 in range(number_eq): 
        a[eq1][0]=-np.random.random()
        for u in range(number_var):
            a[eq1][u+1]=np.random.random()
 
 
    sub=nn.nn()
    sub.number=number
    sub.number_var=number_var
    sub.number_eq=number_eq
    sub.link=link
    sub.x=x
    sub.y=y
    sub.a=a
    sub.makedata()
    a=sub.a
 
    es=np.zeros([7,number])
    for i in range(number):
        es[5][i]=x[i][0]
        es[6][i]=x[i][1]
 
    for n in range(number_eq):
        for i in range(number):
            f1=a[n][0] 
            for j in range(number_var):
                n1=link[n][j+1]
                f1=f1+a[n][j+1]*es[n1][i]
            n1=link[n][0]
            f2=0
            if f1>0:f2=1
            es[n1][i]=f2
 
    c=0
    for i in range(number):
        if y[i]==es[0][i]:c=c+1
 
    if c>maxc:
        maxc=c
        for eq1 in range(number_eq): 
            for m in range(number_var+1):
                    alpha[eq1][m]=a[eq1][m]
 
    print(c)
 
 
print(alpha)
 
import numpy as np
import math
 
class nn:
 
    def g(self):
 
        g3=0
        for i in range(self.number):
            g3=g3+(self.y[i]-self.e[0][i])*(self.y[i]-self.e[0][i])
 
        return g3
 
    def f(self,ax):
 
        for n in range(self.number_eq):
            for i in range(self.number):
                f1=ax[n][0] 
                for j in range(self.number_var):
                    n1=self.link[n][j+1]
                    f1=f1+ax[n][j+1]*self.e[n1][i]
                n1=self.link[n][0]
                self.e[n1][i]=1/(1+math.exp(-f1)) 
 
    def makedata(self):
 
        self.e=np.zeros([7,self.number])
        df=np.zeros([self.number_eq,self.number_var+1])
        ax=np.zeros([self.number_eq,self.number_var+1])
 
        for i in range(self.number):
            self.e[5][i]=self.x[i][0]
            self.e[6][i]=self.x[i][1]
 
        tr=0
        while tr<100:
            tr=tr+1
 
            self.f(self.a)
            g1=self.g()
 
            for eqs in range(self.number_eq): 
                for mx in range(self.number_var+1):
 
 
 
                    for eq1 in range(self.number_eq): 
                        for m in range(self.number_var+1):
                            ax[eq1][m]=self.a[eq1][m]
 
                    h=0.01
                    ax[eqs][mx]=ax[eqs][mx]+h
                    self.f(ax)
                    g2=self.g()
                    df[eqs][mx]=(g2-g1)/h
 
 
            sig=0
            for eq1 in range(self.number_eq): 
                for m in range(self.number_var+1):
                    sig=sig+df[eq1][m]*df[eq1][m]
 
            for eq1 in range(self.number_eq): 
                for m in range(self.number_var+1):
                    ax[eq1][m]=self.a[eq1][m]-0.01*g1*df[eq1][m]/sig
 
            self.f(ax)
            g2=self.g()
            if g2>g1:tr=1000
            if g2<g1:
                for eq1 in range(self.number_eq): 
                    for m in range(self.number_var+1):
                        self.a[eq1][m]=ax[eq1][m]