import numpy as np
from scipy import stats

def calculate_probability(data_a, data_b):
    # Calculate relative likelihood
    prob_a = stats.norm(data_a.mean(), data_a.std())
    prob_b = stats.norm(data_b.mean(), data_b.std())
    
    return prob_a.pdf(0.5) / prob_b.pdf(0.5)

result = calculate_probability(dataset_1, dataset_2)
print(f"Relative Probability: {result:.4f}")

# Output:
Relative Probability