formula = 'mood ~ 1+ [negtweet1 ~ badweather_friends]'
mod4 = IV2SLS.from_formula(formula,fb)
model4 = mod4.fit(cov_type='robust')
model4