Added file tab_5_3.ipynb

Author: AldoCorbelliniUNIPR

cap5/tab_5_3.ipynb

@@ -0,0 +1,102 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "\n",
+    "# <span style=\"color:rgb(213,80,0)\">Theoretical size</span>\n",
+    "\n",
+    "This file creates theoretical parts of Tables 5.2 and 5.3\n",
+    "\n",
+    "## Sample size"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "metadata": {},
+   "source": [
+    "n=[50 100 200 2000];"
+   ],
+   "outputs": []
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## THEORETICAL INDIVIDUAL SIZE"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "metadata": {},
+   "source": [
+    "p=2;\n",
+    "\n",
+    "pro2=1000*(1-betacdf(   n*chi2inv(0.99,1)./(n-p).^2  ,0.5,(n-p-1)/2));\n",
+    "disp(pro2)\n",
+    "\n",
+    "p=10;\n",
+    "\n",
+    "pro10=1000*(1-betacdf(n*chi2inv(0.99,1)./(n-p).^2,0.5,(n-p-1)/2));\n",
+    "disp(pro10)\n",
+    "\n",
+    "varn=\"n=\"+n;\n",
+    "rown=[\"p=2\" \"p=10\"];\n",
+    "TAB52=array2table([pro2; pro10],\"RowNames\",rown,\"VariableNames\",varn);\n",
+    "disp('Theoretical individual size')\n",
+    "disp(TAB52)"
+   ],
+   "outputs": []
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## THEORETICAL SIMULTANEOUS SIZE"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "metadata": {},
+   "source": [
+    "p=2;\n",
+    "pstar=(betacdf(n.*chi2inv(1-0.01./n,1)./((n-p).^2),0.5,(n-p-1)/2));\n",
+    "pstar2=1000*(1-pstar.^n);\n",
+    "% disp(pstar2)\n",
+    "\n",
+    "p=10;\n",
+    "pstar=(betacdf(n.*chi2inv(1-0.01./n,1)./((n-p).^2),0.5,(n-p-1)/2));\n",
+    "pstar10=1000*(1-pstar.^n);\n",
+    "%disp(pstar10)\n",
+    "\n",
+    "TAB53=array2table([pro2; pro10],\"RowNames\",rown,\"VariableNames\",varn);\n",
+    "disp('Theoretical simultaneous size')\n",
+    "disp(TAB53)\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": "25.1.0.2943329"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 4
+}