import numpy as np
import scipy.stats as stats
np.random.seed(2019)
rvs1 = stats.norm.rvs(loc=5, scale=10, size=(500,4))
rvs2 = stats.norm.rvs(loc=5, scale=10, size=(300,4))
t_statistic, p_val = stats.ttest_ind(rvs1, rvs2)
print(t_statistic, p_val)
t_statistic, p_val = stats.ttest_ind(rvs1[:, 0], rvs2[:, 0])
print(t_statistic, p_val)
aW1wb3J0IG51bXB5IGFzIG5wCmltcG9ydCBzY2lweS5zdGF0cyBhcyBzdGF0cwpucC5yYW5kb20uc2VlZCgyMDE5KQpydnMxID0gc3RhdHMubm9ybS5ydnMobG9jPTUsIHNjYWxlPTEwLCBzaXplPSg1MDAsNCkpCnJ2czIgPSBzdGF0cy5ub3JtLnJ2cyhsb2M9NSwgc2NhbGU9MTAsIHNpemU9KDMwMCw0KSkKdF9zdGF0aXN0aWMsIHBfdmFsID0gc3RhdHMudHRlc3RfaW5kKHJ2czEsIHJ2czIpCnByaW50KHRfc3RhdGlzdGljLCBwX3ZhbCkKdF9zdGF0aXN0aWMsIHBfdmFsID0gc3RhdHMudHRlc3RfaW5kKHJ2czFbOiwgMF0sIHJ2czJbOiwgMF0pCnByaW50KHRfc3RhdGlzdGljLCBwX3ZhbCkK