solutionsEX/Ex4_2Interactive.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "# <span style=\"color:rgb(213,80,0)\">Hawkins data: brushing residuals from monitoring S residuals</span>\n",
    "\n",
    "This file creates Figures A.7\\-\\-A.9 in an interactive way\n",
    "\n",
    "## Hawkins data loading"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "source": [
    "close all;\n",
    "load('hawkins.txt');\n",
    "y=hawkins(:,9);\n",
    "X=hawkins(:,1:8);\n",
    "\n",
    "n=length(y);\n",
    "prin=0;"
   ],
   "outputs": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Prepare the input for Figures A.7, A.8 and A.9"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "source": [
    "% Perform both the FS and the S residuals monitoring\n",
    "[outLXS]=LXS(y,X,'nsamp',10000);\n",
    "[outFS]=FSReda(y,X,outLXS.bs);\n",
    "\n",
    "% Sregeda with bisquare link brushing\n",
    "outTB=Sregeda(y,X,'plots',1,'rhofunc','bisquare');\n",
    "yl=5;\n",
    "ylim([-yl yl])\n",
    "% The information about the order of entry of the units in the FS is added\n",
    "% to outTB\n",
    "outTB.Un=outFS.Un;"
   ],
   "outputs": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Create Figure A.7\\-A.9\n",
    "\n",
    "Show the plot of minimum deletion residuals\n",
    ""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "source": [
    "mdrplot(outFS,'ylimy',[1 8])\n",
    "\n",
    "standard=struct;\n",
    "standard.ylim=[-yl yl];\n",
    "\n",
    "databrush=struct;\n",
    "databrush.bivarfit='';\n",
    "databrush.selectionmode='Brush'; % Brush selection\n",
    "% databrush.persist='on'; % Enable repeated mouse selections\n",
    "databrush.Label='off'; % Write labels of trajectories while selecting\n",
    "databrush.RemoveLabels='off'; % Do not remove labels after selection\n",
    "\n",
    "disp('Brush the monitoring residual plot')\n",
    "resfwdplot(outTB,'databrush',databrush,'datatooltip','','standard',standard)\n",
    "if prin==1\n",
    "    % print to postscript\n",
    "    print -depsc HDtbexeres6.eps;\n",
    "    print -depsc HDtbexemdr6.eps;\n",
    "\n",
    "    print -depsc HDtbexeres.eps;\n",
    "    print -depsc HDtbexemdr.eps;\n",
    "    print -depsc HDtbexeyXplot.eps;\n",
    "else\n",
    "   % title('Output which depends on your brushing')\n",
    "end\n",
    "\n",
    "%InsideREADME"
   ],
   "outputs": []
}
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "MATLAB (matlabkernel)",
   "language": "matlab",
   "name": "matlab"
  },
  "language_info": {
   "file_extension": ".m",
   "mimetype": "text/matlab",
   "name": "matlab",
   "nbconvert_exporter": "matlab",
   "pygments_lexer": "matlab",
   "version": "24.1.0.2628055"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}