少し気になっているのは、

1. randomized時のexplained_variance_ratioの計算式
   (batch sizeが1の時はdaskのPCAと一致した方が良いと思うので)
2. componentsなどのattributesをcompute()していない点
   です。

compute はしたほうがいいんでしょうか。たしかにpartial_fitしまくってグラフが大きくなりすぎないほうがいいような気もするんですが。

ただ、future warningがまだ回避できていません。vstackの部分かと思いましたが別の部分から来ているようです。

239行目のnp.expand_dims の部分じゃないかなと思います。