scipyのstatsで確率密度関数が使用できないこととseabornをimportすると解消する理由

Question

前提・実現したいこと

プログラミング初心者です。
基本的な内容だと思いますが、Pythonにおける確率密度関数の使用方法について2点質問させていただきます。
(JupyterNotebookのPython3系を用いて実施)

以下のようにベータ分布の出力をscipyのstatsを用いて試しましたが、下記エラーが発生してしまいました。
ただ、下記「試したこと」のようにseabornをimportすると何故か実行できました。
質問させて頂きたいことは下記の2点です。
1.scipyをimportしているにもかかわらず、なぜstatsは使用できないのか？
2.seabornをimportするとなぜ解決できるのか？

よろしくお願いいたします。

発生している問題・エラーメッセージ

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-2-0dd6a0d1e596> in <module>
      3 
      4 x = np.linspace(0,1,100)
----> 5 y = sp.stats.beta.pdf(x,2,2)

AttributeError: module 'scipy' has no attribute 'stats'

該当のソースコード

import numpy as np
import scipy as sp

x = np.linspace(0,1,100)
y = sp.stats.beta.pdf(x,2,2)

試したこと

import numpy as np
import scipy as sp
import seaborn as sns

x = np.linspace(0,1,100)
y = sp.stats.beta.pdf(x,2,2)
の場合は正しく実行されることを確認しています。

Answer 1

1.scipy.statsについて
以下のようにするとうまく動くはずです。

import scipy.stats as stats

x = np.linspace(0,1,100)
y = stats.beta.pdf(x,2,2)