{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Reliability Scale analysis"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "if(window['d3'] === undefined ||\n",
       "   window['Nyaplot'] === undefined){\n",
       "    var path = {\"d3\":\"https://cdnjs.cloudflare.com/ajax/libs/d3/3.5.5/d3.min\",\"downloadable\":\"https://cdn.rawgit.com/domitry/d3-downloadable/master/d3-downloadable\"};\n",
       "\n",
       "\n",
       "\n",
       "    var shim = {\"d3\":{\"exports\":\"d3\"},\"downloadable\":{\"exports\":\"downloadable\"}};\n",
       "\n",
       "    require.config({paths: path, shim:shim});\n",
       "\n",
       "\n",
       "require(['d3'], function(d3){window['d3']=d3;console.log('finished loading d3');require(['downloadable'], function(downloadable){window['downloadable']=downloadable;console.log('finished loading downloadable');\n",
       "\n",
       "\tvar script = d3.select(\"head\")\n",
       "\t    .append(\"script\")\n",
       "\t    .attr(\"src\", \"https://cdn.rawgit.com/domitry/Nyaplotjs/master/release/nyaplot.js\")\n",
       "\t    .attr(\"async\", true);\n",
       "\n",
       "\tscript[0][0].onload = script[0][0].onreadystatechange = function(){\n",
       "\n",
       "\n",
       "\t    var event = document.createEvent(\"HTMLEvents\");\n",
       "\t    event.initEvent(\"load_nyaplot\",false,false);\n",
       "\t    window.dispatchEvent(event);\n",
       "\t    console.log('Finished loading Nyaplotjs');\n",
       "\n",
       "\t};\n",
       "\n",
       "\n",
       "});});\n",
       "}\n"
      ],
      "text/plain": [
       "\"if(window['d3'] === undefined ||\\n   window['Nyaplot'] === undefined){\\n    var path = {\\\"d3\\\":\\\"https://cdnjs.cloudflare.com/ajax/libs/d3/3.5.5/d3.min\\\",\\\"downloadable\\\":\\\"https://cdn.rawgit.com/domitry/d3-downloadable/master/d3-downloadable\\\"};\\n\\n\\n\\n    var shim = {\\\"d3\\\":{\\\"exports\\\":\\\"d3\\\"},\\\"downloadable\\\":{\\\"exports\\\":\\\"downloadable\\\"}};\\n\\n    require.config({paths: path, shim:shim});\\n\\n\\nrequire(['d3'], function(d3){window['d3']=d3;console.log('finished loading d3');require(['downloadable'], function(downloadable){window['downloadable']=downloadable;console.log('finished loading downloadable');\\n\\n\\tvar script = d3.select(\\\"head\\\")\\n\\t    .append(\\\"script\\\")\\n\\t    .attr(\\\"src\\\", \\\"https://cdn.rawgit.com/domitry/Nyaplotjs/master/release/nyaplot.js\\\")\\n\\t    .attr(\\\"async\\\", true);\\n\\n\\tscript[0][0].onload = script[0][0].onreadystatechange = function(){\\n\\n\\n\\t    var event = document.createEvent(\\\"HTMLEvents\\\");\\n\\t    event.initEvent(\\\"load_nyaplot\\\",false,false);\\n\\t    window.dispatchEvent(event);\\n\\t    console.log('Finished loading Nyaplotjs');\\n\\n\\t};\\n\\n\\n});});\\n}\\n\""
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Analysis 2016-03-26 10:22:36 +0000\n",
      "= Statsample::Reliability\n",
      "  == Reliability Analysis\n",
      "    Summary for Reliability Analysis\n",
      "+------------------------------+-------------------+\n",
      "| Valid Items                  | 20                |\n",
      "| Valid cases                  | 100               |\n",
      "| Sum mean                     | 1.0691            |\n",
      "| Sum sd                       | 19.4551           |\n",
      "| Sum median                   | 4.063620925282448 |\n",
      "+------------------------------+-------------------+\n",
      "| Item mean                    | 0.0535            |\n",
      "| Item sd                      | 0.9728            |\n",
      "+------------------------------+-------------------+\n",
      "| Skewness                     | -0.3883           |\n",
      "| Kurtosis                     | -0.2671           |\n",
      "+------------------------------+-------------------+\n",
      "| Cronbach's alpha             | 0.9979            |\n",
      "| Standarized Cronbach's alpha | 0.9979            |\n",
      "| Mean rpb                     | 0.9788            |\n",
      "| Variances mean               | 0.983765          |\n",
      "| Covariances mean             | 0.944276          |\n",
      "+------------------------------+-------------------+\n",
      "\n",
      "    Items for obtain alpha(0.8) : 0\n",
      "    Items for obtain alpha(0.9) : 0\n",
      "    Items report for Reliability Analysis\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "| item  |  mean   |   sd    | mean if deleted | var if deleted | sd if deleted |  item-total correl. | alpha if deleted |\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "| (v0)  | 0.07570 | 1.02435 | 0.99337         | 340.44186      | 18.45107      | 0.97906             | 0.99780          |\n",
      "| (v1)  | 0.04737 | 1.00360 | 1.02171         | 341.17282      | 18.47086      | 0.97966             | 0.99780          |\n",
      "| (v2)  | 0.09366 | 0.96842 | 0.97541         | 342.58944      | 18.50917      | 0.97556             | 0.99782          |\n",
      "| (v3)  | 0.05185 | 0.97455 | 1.01723         | 342.14884      | 18.49727      | 0.98194             | 0.99779          |\n",
      "| (v4)  | 0.04570 | 0.99565 | 1.02337         | 341.51265      | 18.48006      | 0.97818             | 0.99781          |\n",
      "| (v5)  | 0.06053 | 0.99252 | 1.00855         | 341.53079      | 18.48055      | 0.98092             | 0.99779          |\n",
      "| (v6)  | 0.02392 | 0.98476 | 1.04516         | 341.71873      | 18.48564      | 0.98363             | 0.99777          |\n",
      "| (v7)  | 0.02845 | 0.97928 | 1.04063         | 342.06366      | 18.49496      | 0.97941             | 0.99780          |\n",
      "| (v8)  | 0.06049 | 1.00672 | 1.00859         | 340.94454      | 18.46468      | 0.98291             | 0.99778          |\n",
      "| (v9)  | 0.04812 | 0.96548 | 1.02095         | 342.91228      | 18.51789      | 0.96919             | 0.99786          |\n",
      "| (v10) | 0.03464 | 1.02432 | 1.03444         | 340.43920      | 18.45099      | 0.97916             | 0.99780          |\n",
      "| (v11) | 0.03204 | 0.97632 | 1.03704         | 342.19005      | 18.49838      | 0.97886             | 0.99780          |\n",
      "| (v12) | 0.08685 | 1.01969 | 0.98223         | 340.52967      | 18.45345      | 0.98133             | 0.99779          |\n",
      "| (v13) | 0.05757 | 0.98015 | 1.01150         | 341.99716      | 18.49317      | 0.98042             | 0.99779          |\n",
      "| (v14) | 0.04636 | 0.99898 | 1.02272         | 341.48719      | 18.47937      | 0.97547             | 0.99782          |\n",
      "| (v15) | 0.05206 | 1.01824 | 1.01702         | 340.61883      | 18.45586      | 0.98030             | 0.99780          |\n",
      "| (v16) | 0.08545 | 0.98659 | 0.98362         | 341.94768      | 18.49183      | 0.97511             | 0.99782          |\n",
      "| (v17) | 0.02525 | 0.97607 | 1.04382         | 342.22581      | 18.49935      | 0.97808             | 0.99781          |\n",
      "| (v18) | 0.06765 | 0.95923 | 1.00142         | 342.77599      | 18.51421      | 0.97989             | 0.99780          |\n",
      "| (v19) | 0.04539 | 0.99829 | 1.02368         | 341.43479      | 18.47795      | 0.97767             | 0.99781          |\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "\n",
      "    -50.00 : ***\n",
      "    -40.00 : *****\n",
      "    -30.00 : *******\n",
      "    -20.00 : ************\n",
      "    -10.00 : ******************\n",
      "     0.00 : ********************\n",
      "    10.00 : *********************\n",
      "    20.00 : *********\n",
      "    30.00 : ****\n",
      "    40.00 : *\n",
      "  == Multi Scale analyss\n",
      "    === Reliability analysis of scales\n",
      "      ==== Scale Scale 1\n",
      "        Summary for Scale Scale 1\n",
      "+------------------------------+-------------------+\n",
      "| Valid Items                  | 10                |\n",
      "| Valid cases                  | 100               |\n",
      "| Sum mean                     | 0.4947            |\n",
      "| Sum sd                       | 9.7206            |\n",
      "| Sum median                   | 1.745914154683958 |\n",
      "+------------------------------+-------------------+\n",
      "| Item mean                    | 0.0495            |\n",
      "| Item sd                      | 0.9721            |\n",
      "+------------------------------+-------------------+\n",
      "| Skewness                     | -0.3931           |\n",
      "| Kurtosis                     | -0.2427           |\n",
      "+------------------------------+-------------------+\n",
      "| Cronbach's alpha             | 0.9959            |\n",
      "| Standarized Cronbach's alpha | 0.9960            |\n",
      "| Mean rpb                     | 0.9782            |\n",
      "| Variances mean               | 0.979484          |\n",
      "| Covariances mean             | 0.941068          |\n",
      "+------------------------------+-------------------+\n",
      "\n",
      "        Items for obtain alpha(0.8) : 0\n",
      "        Items for obtain alpha(0.9) : 0\n",
      "        Items report for Scale Scale 1\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "| item  |  mean   |   sd    | mean if deleted | var if deleted | sd if deleted |  item-total correl. | alpha if deleted |\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "| (v1)  | 0.04737 | 1.00360 | 0.44736         | 76.33585       | 8.73704       | 0.97782             | 0.99549          |\n",
      "| (v10) | 0.09366 | 0.96842 | 0.40107         | 76.94542       | 8.77185       | 0.97752             | 0.99550          |\n",
      "| (v2)  | 0.05185 | 0.97455 | 0.44289         | 76.77019       | 8.76186       | 0.98204             | 0.99538          |\n",
      "| (v3)  | 0.04570 | 0.99565 | 0.44903         | 76.48009       | 8.74529       | 0.97732             | 0.99550          |\n",
      "| (v4)  | 0.06053 | 0.99252 | 0.43421         | 76.47154       | 8.74480       | 0.98131             | 0.99540          |\n",
      "| (v5)  | 0.02392 | 0.98476 | 0.47081         | 76.59234       | 8.75170       | 0.98214             | 0.99538          |\n",
      "| (v6)  | 0.02845 | 0.97928 | 0.46629         | 76.74124       | 8.76021       | 0.97863             | 0.99547          |\n",
      "| (v7)  | 0.06049 | 1.00672 | 0.43425         | 76.25645       | 8.73249       | 0.97945             | 0.99545          |\n",
      "| (v8)  | 0.04812 | 0.96548 | 0.44661         | 77.16046       | 8.78410       | 0.96678             | 0.99578          |\n",
      "| (v9)  | 0.03464 | 1.02432 | 0.46010         | 75.96873       | 8.71600       | 0.97855             | 0.99549          |\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "\n",
      "        -30.00 : *\n",
      "        -24.00 : ***\n",
      "        -18.00 : *******\n",
      "        -12.00 : ***********\n",
      "        -6.00 : ***********************\n",
      "         0.00 : **************************\n",
      "         6.00 : ********************\n",
      "        12.00 : ******\n",
      "        18.00 : ***\n",
      "        24.00 : \n",
      "      ==== Scale Scale 2\n",
      "        Summary for Scale Scale 2\n",
      "+------------------------------+--------------------+\n",
      "| Valid Items                  | 9                  |\n",
      "| Valid cases                  | 100                |\n",
      "| Sum mean                     | 0.4986             |\n",
      "| Sum sd                       | 8.7550             |\n",
      "| Sum median                   | 1.5885071766348218 |\n",
      "+------------------------------+--------------------+\n",
      "| Item mean                    | 0.0554             |\n",
      "| Item sd                      | 0.9728             |\n",
      "+------------------------------+--------------------+\n",
      "| Skewness                     | -0.3662            |\n",
      "| Kurtosis                     | -0.3053            |\n",
      "+------------------------------+--------------------+\n",
      "| Cronbach's alpha             | 0.9954             |\n",
      "| Standarized Cronbach's alpha | 0.9954             |\n",
      "| Mean rpb                     | 0.9774             |\n",
      "| Variances mean               | 0.981242           |\n",
      "| Covariances mean             | 0.941926           |\n",
      "+------------------------------+--------------------+\n",
      "\n",
      "        Items for obtain alpha(0.8) : 0\n",
      "        Items for obtain alpha(0.9) : 0\n",
      "        Items report for Scale Scale 2\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "| item  |  mean   |   sd    | mean if deleted | var if deleted | sd if deleted |  item-total correl. | alpha if deleted |\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "| (v11) | 0.03204 | 0.97632 | 0.46660         | 60.84026       | 7.80002       | 0.97543             | 0.99487          |\n",
      "| (v12) | 0.08685 | 1.01969 | 0.41179         | 60.12024       | 7.75372       | 0.97958             | 0.99475          |\n",
      "| (v13) | 0.05757 | 0.98015 | 0.44107         | 60.71265       | 7.79183       | 0.98050             | 0.99470          |\n",
      "| (v14) | 0.04636 | 0.99898 | 0.45228         | 60.48544       | 7.77724       | 0.97605             | 0.99485          |\n",
      "| (v15) | 0.05206 | 1.01824 | 0.44658         | 60.15395       | 7.75590       | 0.97875             | 0.99477          |\n",
      "| (v16) | 0.08545 | 0.98659 | 0.41318         | 60.71287       | 7.79185       | 0.97326             | 0.99494          |\n",
      "| (v17) | 0.02525 | 0.97607 | 0.47338         | 60.81200       | 7.79820       | 0.97779             | 0.99479          |\n",
      "| (v18) | 0.06765 | 0.95923 | 0.43099         | 61.06099       | 7.81415       | 0.97850             | 0.99479          |\n",
      "| (v19) | 0.04539 | 0.99829 | 0.45325         | 60.48152       | 7.77699       | 0.97709             | 0.99481          |\n",
      "+-------+---------+---------+-----------------+----------------+---------------+---------------------+------------------+\n",
      "\n",
      "        -30.00 : \n",
      "        -25.00 : **\n",
      "        -20.00 : ****\n",
      "        -15.00 : ******\n",
      "        -10.00 : *************\n",
      "        -5.00 : ********************\n",
      "         0.00 : **********************\n",
      "         5.00 : ********************\n",
      "        10.00 : **********\n",
      "        15.00 : ***\n",
      "\n"
     ]
    }
   ],
   "source": [
    "require 'statsample'\n",
    "\n",
    "Statsample::Analysis.store(Statsample::Reliability) do \n",
    "  Daru.lazy_update = true\n",
    "  samples=100\n",
    "  a=rnorm(samples)\n",
    "  \n",
    "  ds = Daru::DataFrame.new({})\n",
    "  \n",
    "  20.times do |i|\n",
    "    ds[\"v#{i}\".to_sym]= a + rnorm(samples,0,0.2)\n",
    "  end  \n",
    "  ds.update\n",
    "  \n",
    "  rel=Statsample::Reliability::ScaleAnalysis.new(ds)\n",
    "  summary rel\n",
    "  \n",
    "  ms=Statsample::Reliability::MultiScaleAnalysis.new(:name=>\"Multi Scale analyss\") do |m|\n",
    "    m.scale \"Scale 1\", ds.clone([:v1, :v2, :v3, :v4, :v5, :v6, :v7, :v8, :v9, :v10])\n",
    "    m.scale \"Scale 2\", ds.clone([:v11, :v12, :v13, :v14, :v15, :v16, :v17, :v18, :v19])\n",
    "  end\n",
    "  \n",
    "  summary ms\n",
    "  \n",
    "  Daru.lazy_update = false\n",
    "end\n",
    "\n",
    "Statsample::Analysis.run_batch"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Ruby 2.2.1",
   "language": "ruby",
   "name": "ruby"
  },
  "language_info": {
   "file_extension": ".rb",
   "mimetype": "application/x-ruby",
   "name": "ruby",
   "version": "2.2.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}