共分散の求め方とは!?共分散の公式

みなさん、こんにちは。今回のテーマはデータ分析の【共分散】とその求め方について解説します。

 

たなか君
たなか君
データ分析って式が複雑わからない


 

データ分析がむずかしそうというイメージを持っているみなさんに向け、共分散をやさしく解説していきます。データ分析は、あまりなじみがないため、はじめは難しいですが、慣れてしまえばかんたんに解くことができます。

 

この記事では、まず共分散とは何かを説明し、共分散の公式と求め方を確認します。共分散は相関係数を求める際にも必要になるので、しっかりと覚えるようにしましょう。

 

この記事を15分で読んでできること・共分散とは何かがわかる

・共分散の求め方がわかる

・自分で実際に共分散を求められる

スポンサーリンク
スポンサーリンク

共分散とは?

共分散とは、2つの変数の偏差の積の平均です。2変数の関係の強さを表す「相関係数」を求めるときに、共分散が必要になります。

 

これだけ聞いても、よくわからないですよね。式を確認してみましょう。

共分散の公式

共分散の公式は下のとおりです。

共分散の公式

$$S_{xy}=\frac{1}{n}\sum_{i=1}^{n}\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)$$

たなか君
たなか君
記号が多すぎてまったく分からない…


 

たしかに公式をみると、見慣れない文字ばかりで、よけいにむずかしく見えるかもしれませんね。それでは、1つずつ確認していきましょう。

 

$S_{xy}$:xとyの共分散

$\sum_{i=1}^{n}$:右の式の$i$に1からnまで代入して、それをすべて足した値
          ※$\sum_{i=1}^{5}k$だと、1+2+3+4+5になります

$x_i$:変数$x$の$i$番目の値

$\bar{x}$:変数$x$の平均

 

なんとなくイメージはつかめてきましたか?次は、実際にデータをつかって、共分散を求める練習をしてみましょう。

共分散の求め方

例題で共分散の求め方を確認しましょう。

例題
x53426
y71435

上のようなデータがあるとします。xとyの共分散を求めましょう。

 

公式は

$$S_{xy}=\frac{1}{n}\sum_{i=1}^{n}\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)$$

でしたね。

 

まずは、xとyそれぞれの平均を求めましょう。

$\bar{x}=\frac{1}{5}\left(5+3+4+2+6\right)=\frac{1}{5}\times 20=4$

$\bar{y}=\frac{1}{5}\left(7+1+4+3+5\right)=\frac{1}{5}\times 20=4$

 

計算に必要になるのが、$x_i-\bar{x}$(xの偏差)と$y_i-\bar{y}$(yの偏差)の積です。次のような表をつくるのがおすすめです。

$x_i-\bar{x}$1-10-22
$y_i-\bar{y}$3-30-11
$\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)$33022

 

$\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)$がわかったので、共分散の公式に代入しましょう。

 

$S_{xy}=\frac{1}{5}\sum_{i=1}^{5}\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)$

$=\frac{1}{5}\left(3+3+0+2+2\right)$

$=\frac{1}{5}\times 10=2$

 

よって、xとyの共分散$S_xy=2$と求められました。

 

共分散の求め方、理解できましたか?さっそく練習問題を解いて、理解を定着させましょう。

問題を解いてみよう

問題

次のデータにおける、xとyの共分散を求めましょう。

x23622
y43215

解答

$\bar{x}=\frac{1}{5}\left(2+3+6+2+2\right)=\frac{1}{5}\times 15=3$

$\bar{y}=\frac{1}{5}\left(4+3+2+1+5\right)=\frac{1}{5}\times 15=3$

$\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)$は以下の表のとおり。

$x_i-\bar{x}$-103-1-1
$y_i-\bar{y}$10-1-22
$\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)$-10-32-2

これを$S_{xy}=\frac{1}{5}\sum_{i=1}^{5}\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)$に代入すると

$S_{xy}=\frac{1}{5}\left(-1+0-3+2-2\right)=\frac{1}{5}\times \left(-4\right)=-0.8$

今回のまとめ

今回は、共分散について解説しました。

 

はじめに述べたとおり、共分散は相関係数を求めるうえで必須の要素です。センター試験でも出題されるので、しっかりと押さえておきましょう。

 

今回もおつかれさまでした。

スポンサーリンク
スポンサーリンク

コメント

");const o=ie?ie.createHTML(e):e;if(nt===tt)try{t=(new W).parseFromString(o,at)}catch(e){}if(!t||!t.documentElement){t=le.createDocument(nt,"template",null);try{t.documentElement.innerHTML=ot?ae:o}catch(e){}}const i=t.body||t.documentElement;return e&&n&&i.insertBefore(r.createTextNode(n),i.childNodes[0]||null),nt===tt?ue.call(t,Ie?"html":"body")[0]:Ie?t.documentElement:i},bt=function(e){return ce.call(e.ownerDocument||e,e,H.SHOW_ELEMENT|H.SHOW_COMMENT|H.SHOW_TEXT|H.SHOW_PROCESSING_INSTRUCTION|H.SHOW_CDATA_SECTION,null)},St=function(e){return e instanceof B&&("string"!=typeof e.nodeName||"string"!=typeof e.textContent||"function"!=typeof e.removeChild||!(e.attributes instanceof z)||"function"!=typeof e.removeAttribute||"function"!=typeof e.setAttribute||"string"!=typeof e.namespaceURI||"function"!=typeof e.insertBefore||"function"!=typeof e.hasChildNodes)},Rt=function(e){return"function"==typeof b&&e instanceof b},wt=function(e,t,n){pe[e]&&u(pe[e],(e=>{e.call(o,t,n,ut)}))},Ct=function(e){let t=null;if(wt("beforeSanitizeElements",e,null),St(e))return _t(e),!0;const n=st(e.nodeName);if(wt("uponSanitizeElement",e,{tagName:n,allowedTags:Ne}),e.hasChildNodes()&&!Rt(e.firstElementChild)&&_(/<[/\w]/g,e.innerHTML)&&_(/<[/\w]/g,e.textContent))return _t(e),!0;if(e.nodeType===J)return _t(e),!0;if(Me&&e.nodeType===Q&&_(/<[/\w]/g,e.data))return _t(e),!0;if(!Ne[n]||Ce[n]){if(!Ce[n]&&Dt(n)){if(we.tagNameCheck instanceof RegExp&&_(we.tagNameCheck,n))return!1;if(we.tagNameCheck instanceof Function&&we.tagNameCheck(n))return!1}if(Ye&&!qe[n]){const t=re(e)||e.parentNode,n=oe(e)||e.childNodes;if(n&&t){for(let o=n.length-1;o>=0;--o){const r=X(n[o],!0);r.__removalCount=(e.__removalCount||0)+1,t.insertBefore(r,$(e))}}}return _t(e),!0}return e instanceof R&&!Et(e)?(_t(e),!0):"noscript"!==n&&"noembed"!==n&&"noframes"!==n||!_(/<\/no(script|embed|frames)/i,e.innerHTML)?(ke&&e.nodeType===Z&&(t=e.textContent,u([fe,de,he],(e=>{t=g(t,e," ")})),e.textContent!==t&&(p(o.removed,{element:e.cloneNode()}),e.textContent=t)),wt("afterSanitizeElements",e,null),!1):(_t(e),!0)},Lt=function(e,t,n){if(Be&&("id"===t||"name"===t)&&(n in r||n in mt))return!1;if(ve&&!Le[t]&&_(ge,t));else if(De&&_(Te,t));else if(!Se[t]||Le[t]){if(!(Dt(e)&&(we.tagNameCheck instanceof RegExp&&_(we.tagNameCheck,e)||we.tagNameCheck instanceof Function&&we.tagNameCheck(e))&&(we.attributeNameCheck instanceof RegExp&&_(we.attributeNameCheck,t)||we.attributeNameCheck instanceof Function&&we.attributeNameCheck(t))||"is"===t&&we.allowCustomizedBuiltInElements&&(we.tagNameCheck instanceof RegExp&&_(we.tagNameCheck,n)||we.tagNameCheck instanceof Function&&we.tagNameCheck(n))))return!1}else if(Ze[t]);else if(_(Ae,g(n,Ee,"")));else if("src"!==t&&"xlink:href"!==t&&"href"!==t||"script"===e||0!==T(n,"data:")||!Ke[e]){if(Oe&&!_(ye,g(n,Ee,"")));else if(n)return!1}return!0},Dt=function(e){return"annotation-xml"!==e&&h(e,_e)},vt=function(e){wt("beforeSanitizeAttributes",e,null);const{attributes:t}=e;if(!t)return;const n={attrName:"",attrValue:"",keepAttr:!0,allowedAttributes:Se};let r=t.length;for(;r--;){const i=t[r],{name:a,namespaceURI:l,value:c}=i,s=st(a);let p="value"===a?c:y(c);if(n.attrName=s,n.attrValue=p,n.keepAttr=!0,n.forceKeepAttr=void 0,wt("uponSanitizeAttribute",e,n),p=n.attrValue,n.forceKeepAttr)continue;if(At(a,e),!n.keepAttr)continue;if(!xe&&_(/\/>/i,p)){At(a,e);continue}if(Me&&_(/((--!?|])>)|<\/(style|title)/i,p)){At(a,e);continue}ke&&u([fe,de,he],(e=>{p=g(p,e," ")}));const f=st(e.nodeName);if(Lt(f,s,p)){if(!We||"id"!==s&&"name"!==s||(At(a,e),p=Ge+p),ie&&"object"==typeof G&&"function"==typeof G.getAttributeType)if(l);else switch(G.getAttributeType(f,s)){case"TrustedHTML":p=ie.createHTML(p);break;case"TrustedScriptURL":p=ie.createScriptURL(p)}try{l?e.setAttributeNS(l,a,p):e.setAttribute(a,p),St(e)?_t(e):m(o.removed)}catch(e){}}}wt("afterSanitizeAttributes",e,null)},Ot=function e(t){let n=null;const o=bt(t);for(wt("beforeSanitizeShadowDOM",t,null);n=o.nextNode();)wt("uponSanitizeShadowNode",n,null),Ct(n)||(n.content instanceof s&&e(n.content),vt(n));wt("afterSanitizeShadowDOM",t,null)};return o.sanitize=function(e){let t=arguments.length>1&&void 0!==arguments[1]?arguments[1]:{},n=null,r=null,i=null,l=null;if(ot=!e,ot&&(e="\x3c!--\x3e"),"string"!=typeof e&&!Rt(e)){if("function"!=typeof e.toString)throw A("toString is not a function");if("string"!=typeof(e=e.toString()))throw A("dirty is not a string, aborting")}if(!o.isSupported)return e;if(Ue||ft(t),o.removed=[],"string"==typeof e&&(je=!1),je){if(e.nodeName){const t=st(e.nodeName);if(!Ne[t]||Ce[t])throw A("root node is forbidden and cannot be sanitized in-place")}}else if(e instanceof b)n=Nt("\x3c!----\x3e"),r=n.ownerDocument.importNode(e,!0),r.nodeType===V&&"BODY"===r.nodeName||"HTML"===r.nodeName?n=r:n.appendChild(r);else{if(!Fe&&!ke&&!Ie&&-1===e.indexOf("<"))return ie&&ze?ie.createHTML(e):e;if(n=Nt(e),!n)return Fe?null:ze?ae:""}n&&Pe&&_t(n.firstChild);const c=bt(je?e:n);for(;i=c.nextNode();)Ct(i)||(i.content instanceof s&&Ot(i.content),vt(i));if(je)return e;if(Fe){if(He)for(l=se.call(n.ownerDocument);n.firstChild;)l.appendChild(n.firstChild);else l=n;return(Se.shadowroot||Se.shadowrootmode)&&(l=me.call(a,l,!0)),l}let m=Ie?n.outerHTML:n.innerHTML;return Ie&&Ne["!doctype"]&&n.ownerDocument&&n.ownerDocument.doctype&&n.ownerDocument.doctype.name&&_(q,n.ownerDocument.doctype.name)&&(m="\n"+m),ke&&u([fe,de,he],(e=>{m=g(m,e," ")})),ie&&ze?ie.createHTML(m):m},o.setConfig=function(){let e=arguments.length>0&&void 0!==arguments[0]?arguments[0]:{};ft(e),Ue=!0},o.clearConfig=function(){ut=null,Ue=!1},o.isValidAttribute=function(e,t,n){ut||ft({});const o=st(e),r=st(t);return Lt(o,r,n)},o.addHook=function(e,t){"function"==typeof t&&(pe[e]=pe[e]||[],p(pe[e],t))},o.removeHook=function(e){if(pe[e])return m(pe[e])},o.removeHooks=function(e){pe[e]&&(pe[e]=[])},o.removeAllHooks=function(){pe={}},o}();return oe}))
タイトルとURLをコピーしました