• pmvt(..., df = ) is scalar only