Package com.nr.ci

Examples of com.nr.ci.Svmpolykernel


      nerror += ((y[i]==1.0) != (linsvm.predict(i) >= 0.0) ? 1 : 0);
    }
    System.out.printf("Errors: %d\n", nerror);

    // Polynomial kernel
    Svmpolykernel polykernel=new Svmpolykernel(data,y,1.0,1.0,2.0);
    Svm polysvm=new Svm(polykernel);
    lambda=10;
    k=0;
    do {
      test=polysvm.relax(lambda,omega);
//      System.out.printf(test);
      k++;
    } while (test > 1.e-3 && k < 100);
    nerror=0;
    for (i=0;i<M;i++) {
      nerror += ((y[i]==1.0) != (polysvm.predict(i) >= 0.0) ? 1 : 0);
    }
    System.out.printf("Errors: %d\n", nerror);

    // Gaussian kernel
    Svmgausskernel gausskernel=new Svmgausskernel(data,y,1.0);
    Svm gausssvm=new Svm(gausskernel);
    lambda=10;
    k=0;
    do {
      test=gausssvm.relax(lambda,omega);
//      System.out.printf(test);
      k++;
    } while (test > 1.e-3 && k < 100);
    nerror=0;
    for (i=0;i<M;i++) {
      nerror += ((y[i]==1.0) != (gausssvm.predict(i) >= 0.0) ? 1 : 0);
    }
    System.out.printf("Errors: %d\n", nerror);

   
    // Need to add tests for harder test case and resolve issue that the two
    // support vectors give an erroneous indication for two of the kernels above

    // Example similar to the book
    Normaldev ndev=new Normaldev(0.0,0.5,17);
    for (j=0;j<4;j++) {   // Four quadrants
      for (i=0;i<M/4;i++) {
        k=(M/4)*j+i;
        if (j == 0) {
          y[k]=1.0;
          data[k][0]=1.0+ndev.dev();
          data[k][1]=1.0+ndev.dev();
        } else if (j == 1) {
          y[k]=-1.0;
          data[k][0]=-1.0+ndev.dev();
          data[k][1]=1.0+ndev.dev();
        } else if (j == 2) {
          y[k]=1.0;
          data[k][0]=-1.0+ndev.dev();
          data[k][1]=-1.0+ndev.dev();
        } else {
          y[k]=-1.0;
          data[k][0]=1.0+ndev.dev();
          data[k][1]=-1.0+ndev.dev();
        }
      }
    }
       
    // Linear kernel
    Svmlinkernel linkernel2=new Svmlinkernel(data,y);
    Svm linsvm2=new Svm(linkernel2);
    System.out.printf("Errors: ");
    for (lambda=0.001;lambda<10000;lambda *= 10) {
      k=0;
      do {
        test=linsvm2.relax(lambda,omega);
//        System.out.printf(test);
        k++;
      } while (test > 1.e-3 && k < 100);
      nerror=0;
      for (i=0;i<M;i++) {
        nerror += ((y[i]==1.0) != (linsvm2.predict(i) >= 0.0) ? 1 : 0);
      }
      System.out.printf("%d ",nerror);
      // Test new data
      nerror=0;
      for (j=0;j<4;j++) {   // Four quadrants
        for (i=0;i<M/4;i++) {
          if (j == 0) {
            yy=1.0;
            x[0]=1.0+ndev.dev();
            x[1]=1.0+ndev.dev();
          } else if (j == 1) {
            yy=-1.0;
            x[0]=-1.0+ndev.dev();
            x[1]=1.0+ndev.dev();
          } else if (j == 2) {
            yy=1.0;
            x[0]=-1.0+ndev.dev();
            x[1]=-1.0+ndev.dev();
          } else {
            yy=-1.0;
            x[0]=1.0+ndev.dev();
            x[1]=-1.0+ndev.dev();
          }
          nerror += ((yy==1.0) != (linsvm2.predict(x) >= 0.0) ? 1 : 0);
        }
      }
      System.out.printf("%d    ",nerror);
    }
    System.out.println();

    // Polynomial kernel
    Svmpolykernel polykernel2 = new Svmpolykernel(data,y,1.0,1.0,4.0);
    Svm polysvm2=new Svm(polykernel2);
    System.out.printf("Errors: ");
    for (lambda=0.001;lambda<10000;lambda *= 10) {
      k=0;
      do {
View Full Code Here

TOP

Related Classes of com.nr.ci.Svmpolykernel

Copyright © 2018 www.massapicom. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.