ANOVA¶
Analysis of Variance models
Examples¶
In [1]: import statsmodels.api as sm
In [2]: from statsmodels.formula.api import ols
In [3]: moore = sm.datasets.get_rdataset("Moore", "car",
...: cache=True) # load data
...:
---------------------------------------------------------------------------
URLError Traceback (most recent call last)
<ipython-input-3-924d52d36060> in <module>()
1 moore = sm.datasets.get_rdataset("Moore", "car",
----> 2 cache=True) # load data
/build/statsmodels-ungkPp/statsmodels-0.6.1/debian/python-statsmodels/usr/lib/python2.7/dist-packages/statsmodels/datasets/utils.pyc in get_rdataset(dataname, package, cache)
284 "master/doc/"+package+"/rst/")
285 cache = _get_cache(cache)
--> 286 data, from_cache = _get_data(data_base_url, dataname, cache)
287 data = read_csv(data, index_col=0)
288 data = _maybe_reset_index(data)
/build/statsmodels-ungkPp/statsmodels-0.6.1/debian/python-statsmodels/usr/lib/python2.7/dist-packages/statsmodels/datasets/utils.pyc in _get_data(base_url, dataname, cache, extension)
215 url = base_url + (dataname + ".%s") % extension
216 try:
--> 217 data, from_cache = _urlopen_cached(url, cache)
218 except HTTPError as err:
219 if '404' in str(err):
/build/statsmodels-ungkPp/statsmodels-0.6.1/debian/python-statsmodels/usr/lib/python2.7/dist-packages/statsmodels/datasets/utils.pyc in _urlopen_cached(url, cache)
206 # not using the cache or didn't find it in cache
207 if not from_cache:
--> 208 data = urlopen(url).read()
209 if cache is not None: # then put it in the cache
210 _cache_it(data, cache_path)
/usr/lib/python2.7/urllib2.pyc in urlopen(url, data, timeout, cafile, capath, cadefault, context)
152 else:
153 opener = _opener
--> 154 return opener.open(url, data, timeout)
155
156 def install_opener(opener):
/usr/lib/python2.7/urllib2.pyc in open(self, fullurl, data, timeout)
427 req = meth(req)
428
--> 429 response = self._open(req, data)
430
431 # post-process response
/usr/lib/python2.7/urllib2.pyc in _open(self, req, data)
445 protocol = req.get_type()
446 result = self._call_chain(self.handle_open, protocol, protocol +
--> 447 '_open', req)
448 if result:
449 return result
/usr/lib/python2.7/urllib2.pyc in _call_chain(self, chain, kind, meth_name, *args)
405 func = getattr(handler, meth_name)
406
--> 407 result = func(*args)
408 if result is not None:
409 return result
/usr/lib/python2.7/urllib2.pyc in https_open(self, req)
1239 def https_open(self, req):
1240 return self.do_open(httplib.HTTPSConnection, req,
-> 1241 context=self._context)
1242
1243 https_request = AbstractHTTPHandler.do_request_
/usr/lib/python2.7/urllib2.pyc in do_open(self, http_class, req, **http_conn_args)
1196 except socket.error, err: # XXX what error?
1197 h.close()
-> 1198 raise URLError(err)
1199 else:
1200 try:
URLError: <urlopen error [Errno -2] Name or service not known>
In [4]: data = moore.data
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-4-07b15baaa84d> in <module>()
----> 1 data = moore.data
NameError: name 'moore' is not defined
In [5]: data = data.rename(columns={"partner.status" :
...: "partner_status"}) # make name pythonic
...:
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-5-aaf036e8c484> in <module>()
----> 1 data = data.rename(columns={"partner.status" :
2 "partner_status"}) # make name pythonic
3
NameError: name 'data' is not defined
In [6]: moore_lm = ols('conformity ~ C(fcategory, Sum)*C(partner_status, Sum)',
...: data=data).fit()
...:
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-6-0a291eef3047> in <module>()
1 moore_lm = ols('conformity ~ C(fcategory, Sum)*C(partner_status, Sum)',
----> 2 data=data).fit()
NameError: name 'data' is not defined
In [7]: table = sm.stats.anova_lm(moore_lm, typ=2) # Type 2 ANOVA DataFrame
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-7-4e32df33effd> in <module>()
----> 1 table = sm.stats.anova_lm(moore_lm, typ=2) # Type 2 ANOVA DataFrame
NameError: name 'moore_lm' is not defined
In [8]: print table
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-8-97e179e9235e> in <module>()
----> 1 print table
NameError: name 'table' is not defined
A more detailed example can be found here: