For文を使ったグラフの書き方

前提

for文を使って統計量G1を計算したのですが、サンプルサイズn=10とn=500の時のグラフだけ表記させたいのですが、自分で書いたコードだとfor文を使った意味がないと思います
正規分布に従うことを示すために正規分布のグラフも追加しています

実現したいこと

for文を使ってn=10とn=500のグラフだけを作成したいです.

該当のソースコード

python
1from scipy.stats import norm, skew
2import numpy as np
3import matplotlib.pyplot as plt
4
5n = [10, 20, 30, 50, 100, 300, 500, 1000] #samplesize
6N =10000 
7alfa = 0.05
8for i in range(len(n)):
9    X1 = norm.rvs(loc = 0, scale =1, size=(n[i],N))
10    G1 = np.sqrt(n[i] / 6) * skew(X1)
11n1=10
12X1 = norm.rvs(loc = 0, scale =1, size=(n1,N))    
13G1 = np.sqrt(n1 / 6) * skew(X1)
14# histogram
15plt.hist(G1, bins=100, edgecolor = 'black',density=True)
16x = np.linspace(-3, 3, 100)
17norm_pdf = norm.pdf(x)
18plt.plot(x, norm_pdf, lw = 3, color='r')
19plt.show()
20n2=500
21X1 = norm.rvs(loc = 0, scale =1, size=(n2,N))
22G1 = np.sqrt(n2 / 6) * skew(X1)
23plt.hist(G1, bins=100,edgecolor = 'black', density=True)
24x = np.linspace(-4, 4, 100)
25norm_pdf = norm.pdf(x)
26plt.plot(x, norm_pdf, lw = 3, color='r')
27plt.show()
28

行動規範の内容に同意します

回答1件

ベストアンサー

一例。

python
1from scipy.stats import norm, skew
2import numpy as np
3import matplotlib.pyplot as plt
4
5N = 10000
6n = [(10, 3), (500, 4)]
7for i, j in n:
8    X1 = norm.rvs(loc=0, scale=1, size=(i, N))
9    G1 = np.sqrt(i/6) * skew(X1)
10    # histogram
11    plt.hist(G1, bins=100, edgecolor = 'black',density=True)
12    x = np.linspace(-j, j, 100)
13    norm_pdf = norm.pdf(x)
14    plt.plot(x, norm_pdf, lw = 3, color='r')
15    plt.show()