package com.imsl.test.example.stat;

import com.imsl.stat.*;
import com.imsl.math.PrintMatrix;
import com.imsl.math.PrintMatrixFormat;
/**
 * <p>Selects a pseudorandom sample from Fisher's iris data.</p>
 * 
 * <p>In this example, the <code>getSamples()</code> method is used to generate
five pseudorandom samples from the set of 150 observations of Fisher's iris
data. The sampled indices are returned via <code>getIndices()</code>.</p>

 * 
 * @see <a href="RandomSamplesEx5.java">Code</a>
 * @see <a href="doc-files/RandomSamplesEx5.out">Output</a>
 */

public class RandomSamplesEx5 {

    public static void main(String[] args) {
        double[][] x = {
            {5.1, 3.5, 1.4, 0.2}, {4.9, 3.0, 1.4, 0.2},
            {4.7, 3.2, 1.3, 0.2}, {4.6, 3.1, 1.5, 0.2},
            {5.0, 3.6, 1.4, 0.2}, {5.4, 3.9, 1.7, 0.4},
            {4.6, 3.4, 1.4, 0.3}, {5.0, 3.4, 1.5, 0.2},
            {4.4, 2.9, 1.4, 0.2}, {4.9, 3.1, 1.5, 0.1},
            {5.4, 3.7, 1.5, 0.2}, {4.8, 3.4, 1.6, 0.2},
            {4.8, 3.0, 1.4, 0.1}, {4.3, 3.0, 1.1, 0.1},
            {5.8, 4.0, 1.2, 0.2}, {5.7, 4.4, 1.5, 0.4},
            {5.4, 3.9, 1.3, 0.4}, {5.1, 3.5, 1.4, 0.3},
            {5.7, 3.8, 1.7, 0.3}, {5.1, 3.8, 1.5, 0.3},
            {5.4, 3.4, 1.7, 0.2}, {5.1, 3.7, 1.5, 0.4},
            {4.6, 3.6, 1.0, 0.2}, {5.1, 3.3, 1.7, 0.5},
            {4.8, 3.4, 1.9, 0.2}, {5.0, 3.0, 1.6, 0.2},
            {5.0, 3.4, 1.6, 0.4}, {5.2, 3.5, 1.5, 0.2},
            {5.2, 3.4, 1.4, 0.2}, {4.7, 3.2, 1.6, 0.2},
            {4.8, 3.1, 1.6, 0.2}, {5.4, 3.4, 1.5, 0.4},
            {5.2, 4.1, 1.5, 0.1}, {5.5, 4.2, 1.4, 0.2},
            {4.9, 3.1, 1.5, 0.2}, {5.0, 3.2, 1.2, 0.2},
            {5.5, 3.5, 1.3, 0.2}, {4.9, 3.6, 1.4, 0.1},
            {4.4, 3.0, 1.3, 0.2}, {5.1, 3.4, 1.5, 0.2},
            {5.0, 3.5, 1.3, 0.3}, {4.5, 2.3, 1.3, 0.3},
            {4.4, 3.2, 1.3, 0.2}, {5.0, 3.5, 1.6, 0.6},
            {5.1, 3.8, 1.9, 0.4}, {4.8, 3.0, 1.4, 0.3},
            {5.1, 3.8, 1.6, 0.2}, {4.6, 3.2, 1.4, 0.2},
            {5.3, 3.7, 1.5, 0.2}, {5.0, 3.3, 1.4, 0.2},
            {7.0, 3.2, 4.7, 1.4}, {6.4, 3.2, 4.5, 1.5},
            {6.9, 3.1, 4.9, 1.5}, {5.5, 2.3, 4.0, 1.3},
            {6.5, 2.8, 4.6, 1.5}, {5.7, 2.8, 4.5, 1.3},
            {6.3, 3.3, 4.7, 1.6}, {4.9, 2.4, 3.3, 1.0},
            {6.6, 2.9, 4.6, 1.3}, {5.2, 2.7, 3.9, 1.4},
            {5.0, 2.0, 3.5, 1.0}, {5.9, 3.0, 4.2, 1.5},
            {6.0, 2.2, 4.0, 1.0}, {6.1, 2.9, 4.7, 1.4},
            {5.6, 2.9, 3.6, 1.3}, {6.7, 3.1, 4.4, 1.4},
            {5.6, 3.0, 4.5, 1.5}, {5.8, 2.7, 4.1, 1.0},
            {6.2, 2.2, 4.5, 1.5}, {5.6, 2.5, 3.9, 1.1},
            {5.9, 3.2, 4.8, 1.8}, {6.1, 2.8, 4.0, 1.3},
            {6.3, 2.5, 4.9, 1.5}, {6.1, 2.8, 4.7, 1.2},
            {6.4, 2.9, 4.3, 1.3}, {6.6, 3.0, 4.4, 1.4},
            {6.8, 2.8, 4.8, 1.4}, {6.7, 3.0, 5.0, 1.7},
            {6.0, 2.9, 4.5, 1.5}, {5.7, 2.6, 3.5, 1.0},
            {5.5, 2.4, 3.8, 1.1}, {5.5, 2.4, 3.7, 1.0},
            {5.8, 2.7, 3.9, 1.2}, {6.0, 2.7, 5.1, 1.6},
            {5.4, 3.0, 4.5, 1.5}, {6.0, 3.4, 4.5, 1.6},
            {6.7, 3.1, 4.7, 1.5}, {6.3, 2.3, 4.4, 1.3},
            {5.6, 3.0, 4.1, 1.3}, {5.5, 2.5, 4.0, 1.3},
            {5.5, 2.6, 4.4, 1.2}, {6.1, 3.0, 4.6, 1.4},
            {5.8, 2.6, 4.0, 1.2}, {5.0, 2.3, 3.3, 1.0},
            {5.6, 2.7, 4.2, 1.3}, {5.7, 3.0, 4.2, 1.2},
            {5.7, 2.9, 4.2, 1.3}, {6.2, 2.9, 4.3, 1.3},
            {5.1, 2.5, 3.0, 1.1}, {5.7, 2.8, 4.1, 1.3},
            {6.3, 3.3, 6.0, 2.5}, {5.8, 2.7, 5.1, 1.9},
            {7.1, 3.0, 5.9, 2.1}, {6.3, 2.9, 5.6, 1.8},
            {6.5, 3.0, 5.8, 2.2}, {7.6, 3.0, 6.6, 2.1},
            {4.9, 2.5, 4.5, 1.7}, {7.3, 2.9, 6.3, 1.8},
            {6.7, 2.5, 5.8, 1.8}, {7.2, 3.6, 6.1, 2.5},
            {6.5, 3.2, 5.1, 2.0}, {6.4, 2.7, 5.3, 1.9},
            {6.8, 3.0, 5.5, 2.1}, {5.7, 2.5, 5.0, 2.0},
            {5.8, 2.8, 5.1, 2.4}, {6.4, 3.2, 5.3, 2.3},
            {6.5, 3.0, 5.5, 1.8}, {7.7, 3.8, 6.7, 2.2},
            {7.7, 2.6, 6.9, 2.3}, {6.0, 2.2, 5.0, 1.5},
            {6.9, 3.2, 5.7, 2.3}, {5.6, 2.8, 4.9, 2.0},
            {7.7, 2.8, 6.7, 2.0}, {6.3, 2.7, 4.9, 1.8},
            {6.7, 3.3, 5.7, 2.1}, {7.2, 3.2, 6.0, 1.8},
            {6.2, 2.8, 4.8, 1.8}, {6.1, 3.0, 4.9, 1.8},
            {6.4, 2.8, 5.6, 2.1}, {7.2, 3.0, 5.8, 1.6},
            {7.4, 2.8, 6.1, 1.9}, {7.9, 3.8, 6.4, 2.0},
            {6.4, 2.8, 5.6, 2.2}, {6.3, 2.8, 5.1, 1.5},
            {6.1, 2.6, 5.6, 1.4}, {7.7, 3.0, 6.1, 2.3},
            {6.3, 3.4, 5.6, 2.4}, {6.4, 3.1, 5.5, 1.8},
            {6.0, 3.0, 4.8, 1.8}, {6.9, 3.1, 5.4, 2.1},
            {6.7, 3.1, 5.6, 2.4}, {6.9, 3.1, 5.1, 2.3},
            {5.8, 2.7, 5.1, 1.9}, {6.8, 3.2, 5.9, 2.3},
            {6.7, 3.3, 5.7, 2.5}, {6.7, 3.0, 5.2, 2.3},
            {6.3, 2.5, 5.0, 1.9}, {6.5, 3.0, 5.2, 2.0},
            {6.2, 3.4, 5.4, 2.3}, {5.9, 3.0, 5.1, 1.8}
        };
        int nSamp = 5;

        RandomSamples rs = new RandomSamples();
        Random r = new Random(123457L);
        r.setMultiplier(16807);
        rs.setRandomObject(r);
        double[][] samples = rs.getSamples(x, nSamp);
        int[] idx = rs.getIndices();

        PrintMatrix pm = new PrintMatrix("Random Samples");
        PrintMatrixFormat mf = new PrintMatrixFormat();
        mf.setNumberFormat(new java.text.DecimalFormat("0.0"));
        mf.setNoColumnLabels();
        pm.print(mf, samples);
        pm = new PrintMatrix("Sampled Indices");
        pm.print(mf, idx);
    }
}