% \documentclass[serif]{beamer} % Serif for Computer Modern math font.\documentclass[serif, handout]{beamer} % Handout to ignore pause statements\hypersetup{colorlinks,linkcolor=,urlcolor=red}\usefonttheme{serif} % Looks like Computer Modern for non-math text -- nice!\setbeamertemplate{navigation symbols}{} % Suppress navigation symbols% \usetheme{Berlin} % Displays sections on top\usetheme{Frankfurt}  % Displays section titles on top: Fairly thin but still swallows some material at bottom of crowded slides%\usetheme{Berkeley}\usepackage[english]{babel}\usepackage{amsmath} % for binom% \usepackage{graphicx} % To include pdf files!% \definecolor{links}{HTML}{2A1B81}% \definecolor{links}{red}\setbeamertemplate{footline}[frame number] \mode<presentation>\title{Using R as a Calculator (on the final exam)\footnote{ This slide show is an open-source document. See last slide for copyright information.}}\subtitle{STA 260 Spring 2020}\date{} % To suppress date\begin{document}\begin{frame}  \titlepage\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}\frametitle{R is a wonderful calculator} %\framesubtitle{} \begin{itemize}    \item You are going to have access to a computer during the final exam anyway.    \item Even a cell phone will do.    \item You may already have it and know how to use it. \pause    \item Free download at \href{https://www.r-project.org} {\small\texttt{https://www.r-project.org}}    \item Run free online at\href{https://rdrr.io/snippets} {\small\texttt{https://rdrr.io/snippets}}\end{itemize}\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{R is a Basic Calculator}\framesubtitle{} \pause{\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> 1+1\end{verbatim}} % End color\begin{verbatim}[1] 2\end{verbatim} \pause{\color{blue}\begin{verbatim}> 1:40\end{verbatim}} % End color\begin{verbatim} [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20[21] 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{R is an advanced calculator}%\framesubtitle{} {\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> n = 50; xbar = 1.56> Gsq = 2*n*(xbar*log(xbar) - (1+xbar)*(log(1+xbar)-log(2))); Gsq\end{verbatim}} % End color\begin{verbatim}[1] 6.174808\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}\frametitle{No More Tables} %\framesubtitle{} We will use R for CDFs and quantiles of all the familiar distributions. \pause        \begin{itemize}            \item Critical values.            \item $p$-values.            \item Posterior probabilities.        \end{itemize} \end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{Normal Distribution}\framesubtitle{Faster than the table} \pause{\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> pnorm(0) # CDF of standard normal\end{verbatim}} % End color\begin{verbatim}[1] 0.5\end{verbatim}  \pause} % End sizeYou can specify $\mu$ and $\sigma$ (not $\sigma^2$). \pause \vspace{2mm}IQ tests are designed to have $\mu=100$ and $\sigma=15$. What's $P(IQ>160)$? \pause{\footnotesize % or scriptsize{\color{blue}\begin{verbatim}1 - pnorm(160,mean=100,sd=15) # Or just pnorm(160,100,15)\end{verbatim}} % End color\begin{verbatim}[1] 3.167124e-05\end{verbatim}  \pause{\color{blue}\begin{verbatim}> options(scipen=999) # Supress scientific notation> 1 - pnorm(160,mean=100,sd=15)\end{verbatim}} % End color\begin{verbatim}[1] 0.00003167124\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{Quantiles}%\framesubtitle{} {\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> qnorm(0.975)\end{verbatim}} % End color\begin{verbatim}[1] 1.959964\end{verbatim} \pause{\color{blue}\begin{verbatim}> # An IQ of ___ is higher than 90% of the population,> qnorm(0.90,100,15) # q for quantile\end{verbatim} } % End color\begin{verbatim}[1] 119.2233\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{$\chi^2$ Distribution}%\framesubtitle{} {\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> n = 50; xbar = 1.56> Gsq = 2*n*(xbar*log(xbar) - (1+xbar)*(log(1+xbar)-log(2))); Gsq\end{verbatim}} % End color\begin{verbatim}[1] 6.174808\end{verbatim}  \pause{\color{blue}\begin{verbatim}> 1-pchisq(Gsq,df=1) # p-value\end{verbatim}} % End color\begin{verbatim}[1] 0.0129582\end{verbatim}  \pause{\color{blue}\begin{verbatim}> qchisq(0.95,1) # Critical value at alpha = 0.05\end{verbatim}} % End color\begin{verbatim}[1] 3.841459\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{$t$ distribution}%\framesubtitle{} {\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> 2*(1 - pt(2.14,df=10) ) # Two-tailed p-value\end{verbatim}} % End color\begin{verbatim}[1] 0.05803497\end{verbatim}  \pause{\color{blue}\begin{verbatim}> qt(0.975,df=10) # Critical value\end{verbatim}} % End color\begin{verbatim}[1] 2.228139\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{$F$ Distribution}%\framesubtitle{} {\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> 1 - pf(3.17,6,114) \end{verbatim}} % End color\begin{verbatim}[1] 0.006504761\end{verbatim} \pause{\color{blue}\begin{verbatim}> qf(0.95,6,114) # Not in the table\end{verbatim}} % End color\begin{verbatim}[1] 2.1791\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{Gamma with $\alpha=$shape, $\lambda=$rate} \pause%\framesubtitle{} {\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> pgamma(1,shape=21,rate=23.35) # P(Lambda < 1 |x)\end{verbatim}} % End color\begin{verbatim}[1] 0.7146466\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{Beta Distribution}%\framesubtitle{} For the coffee taste test, 60/100 consumers chose the new blend of coffee beans, yielding $\alpha^\prime = \alpha + \sum_{i=1}^nx_i = 61$ and $\beta^\prime = \beta + n-\sum_{i=1}^nx_i = 41$. \pause \vspace{5mm}{\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> pbeta(1/2,61,41)\end{verbatim}} % End color\begin{verbatim}[1] 0.02302203\end{verbatim}  } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{Binomial Distribution}%\framesubtitle{} {\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> # 20-question abcd multiple choice, probability of passing> 1 - pbinom(9,20,0.25)\end{verbatim}} % End color\begin{verbatim}[1] 0.01386442\end{verbatim}  \pause{\color{blue}\begin{verbatim}> # Probability of exactly 50 heads in 100 tosses of a fair coin> dbinom(50,100,0.5)\end{verbatim}} % End color\begin{verbatim}[1] 0.07958924\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}[fragile]\frametitle{Poisson Distribution}%\framesubtitle{} If there are really only a population mean of 8 rat hairs in a peanut butter jar, what is the probability of obtaining a sample mean of 9.2 (or more) from a random sample of 30 jars? \pause Distribution of $S = \sum_{i=1}^nX_i$ is Poisson$(30*8)$. \pause{\footnotesize % or scriptsize{\color{blue}\begin{verbatim}> 9.2*30\end{verbatim}} % End color\begin{verbatim}[1] 276\end{verbatim}  \pause{\color{blue}\begin{verbatim}> 1 - ppois(275,240) # P(S geq 276)\end{verbatim}} % End color\begin{verbatim}[1] 0.01226396\end{verbatim} } % End size\end{frame}%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\begin{frame}\frametitle{Copyright Information}This slide show was prepared by  \href{http://www.utstat.toronto.edu/~brunner}{Jerry Brunner},Department of Mathematical and Computational Sciences, University of Toronto Mississauga. It is licensed under a \href{http://creativecommons.org/licenses/by-sa/3.0/deed.en_US}     {Creative Commons Attribution - ShareAlike 3.0 Unported License}. Use any part of it as you like and share the result freely. The \LaTeX~source code is available from the course website:\vspace{5mm}\href{http://www.utstat.toronto.edu/~brunner/oldclass/260s20} {\small\texttt{http://www.utstat.toronto.edu/$^\sim$brunner/oldclass/260s20}}\end{frame}\end{document}\begin{frame}[fragile]\frametitle{Title here}%\framesubtitle{} {\footnotesize % or scriptsize{\color{blue}\begin{verbatim}R input in blue\end{verbatim}} % End color\begin{verbatim}R output in black\end{verbatim} } % End size\end{frame}