from scipy.stats import kruskal, f_oneway
import pandas as pd
import numpy as np

#if p value less than alpha reject
qiskit_data = np.concatenate((np.full(245, 0.75), np.full(264, 0.5), np.full(243, 0.25), np.full(248, 0)))

cirq_data = np.concatenate((np.full(259, 0.75), np.full(234, 0.5), np.full(249, 0.25), np.full(258, 0)))

qs_data = np.concatenate((np.full(260, 0.75), np.full(234, 0.5), np.full(264, 0.25), np.full(242, 0)))

kruskal(qiskit_data, cirq_data, qs_data)


#if p value less than alpha reject
qiskit_data = np.concatenate((np.full(254, 0.75), np.full(254, 0.5), np.full(230, 0.25), np.full(262, 0)))

cirq_data = np.concatenate((np.full(251, 0.75), np.full(258, 0.5), np.full(269, 0.25), np.full(222, 0)))

qs_data = np.concatenate((np.full(254, 0.75), np.full(232, 0.5), np.full(268, 0.25), np.full(246, 0)))

kruskal(qiskit_data, cirq_data, qs_data)


#if p value less than alpha reject
qiskit_data = np.concatenate((np.full(486, 0.5), np.full(514, 0)))

cirq_data = np.concatenate((np.full(483, 0.5), np.full(517, 0)))

qs_data = np.concatenate((np.full(491, 0.5), np.full(509, 0)))

kruskal(qiskit_data, cirq_data, qs_data)

#if p value less than alpha reject
qiskit_data = np.concatenate((np.full(246, 0.75), np.full(226, 0.5), np.full(237, 0.25), np.full(291, 0)))

cirq_data = np.concatenate((np.full(247, 0.75), np.full(272, 0.5), np.full(230, 0.25), np.full(251, 0)))

qs_data = np.concatenate((np.full(262, 0.75), np.full(243, 0.5), np.full(251, 0.25), np.full(244, 0)))

kruskal(qiskit_data, cirq_data, qs_data)

#if p value less than alpha reject
qiskit_data = np.concatenate((np.full(651, 3+5), np.full(349, 3)))

cirq_data = np.concatenate((np.full(685, 3+5), np.full(315, 3)))

qs_data = np.concatenate((np.full(662, 3+5), np.full(338, 3)))

kruskal(qiskit_data, cirq_data, qs_data)

#if p value less than alpha reject
qiskit_data = np.concatenate((np.full(1000, 3+5), np.full(0, 3)))

cirq_data = np.concatenate((np.full(1000, 3+5), np.full(0, 3)))

qs_data = np.concatenate((np.full(1000, 3+5), np.full(0, 3)))

kruskal(qiskit_data, cirq_data, qs_data)