以前の astropy バージョンでは、次のような方法で不確実性の伝播を処理することができました。
from astropy.nddata import NDData, StdDevUncertainty
x = NDData( 16.0, uncertainty=StdDevUncertainty( 4.0 ))
y = NDData( 361.0, uncertainty=StdDevUncertainty( 19.0 ))
print x.add(y)
NDData への変更により、この機能が削除されたようです。「AttributeError: 'NDData' object has no attribute 'add'」というメッセージが表示され、ドキュメントに役立つアドバイスが見つかりません。エラー伝播は現在どのように処理されていますか?