<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:"Calibri","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-NZ" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoPlainText">Hi Tyler et al.,<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">This is a good discussion topic.<o:p></o:p></p>
<p class="MsoNormal">I do believe the qqplot() has its place; this is why it remains part of the output generated by UniDesc()<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">I think you're being a little harsh on the qqplot. Your comments re the usefulness are correct, but there are instance where the qqplot is an extremely useful tool, such as detecting significant effects in a 2^k factorial design with
no replicates, or a partial fraction of one I plan to implement a VI() summarisation of the qqplot() for this reason alone.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Your comments re the hypothesis tests for normality are true, but they are true for practically all hypothesis tests and not just normality tests.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">The key thing about the normality of a large data set is that it almost certainly means the central limit theorem kicks in and the need for normality is hugely reduced. In such circumstances, the need to test for normality is often redundant
because we "know" what we should be expecting with enough forethought anyway. Our other methods of checking a distribution for the sample and therefore the population can be used instead of qqplot().
<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">In cases where the sample size is small and the population's distribution is unable to be determined, the tests really should not be used as a scape goat for failing to have sufficient data in the first place.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Your comments re skewness and kurtosis are spot on, but the introductory student doesn't know about hypothesis testing when they first encounter the need to determine normality in their exploratory data analysis.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Finally, do be careful when using skewness and kurtosis as there are different measures for both notions. I've been very careful to explicitly state the particular test and measure for each in the UniDesc() output for this reason.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Jonathan<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><span lang="EN-US" style="mso-fareast-language:EN-NZ">-----Original Message-----<br>
From: BlindRUG [mailto:blindrug-bounces@nfbnet.org] On Behalf Of Ensor, Tyler via BlindRUG<br>
Sent: Saturday, 13 June 2015 1:13 a.m.<br>
To: Blind R Users Group<br>
Cc: Ensor, Tyler<br>
Subject: Re: [BlindRUG] Some Questions</span></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">I thought I would add my two cents on assessing normality. Maybe this will help someone. I can't see a Q-Q plot, and used to use hypothesis tests such as the Shapiro-Wilk or Kolmogorov-Smirnov. However, statisticians generally find these
tests unhelpful, as they are too conservative with small sample sizes (i.e., very nonnormal data will still be deemed normal) and too liberal with large samples (i.e., very small deviations from normality will push the p value into the rejection region).<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">The approach I've since adopted is to look at the values of skewness and kurtosis. Skewness and kurtosis both have a standard error associated with them. If the skewness and kurtosis values are both within -2 or +2 standard errors of
0, then you can conclude that data are normal. This isn't a perfect approach, as it is somewhat arbitrary, but inspecting a Q-Q plot also lacks objectivity.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Tyler<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">On 6/12/15, Voica Gavrilut via BlindRUG <blindrug@nfbnet.org> wrote:<o:p></o:p></p>
<p class="MsoPlainText">> Hello,<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">>> 1. Auto Load<o:p></o:p></p>
<p class="MsoPlainText">>> I want BrailleR package to be loaded automatically when starting R.<o:p></o:p></p>
<p class="MsoPlainText">>> How can I do?<o:p></o:p></p>
<p class="MsoPlainText">> For example I use the config file that already exists in the R
<o:p></o:p></p>
<p class="MsoPlainText">> installation - the R_HOME/etc/Rprofile. in this file I added a
<o:p></o:p></p>
<p class="MsoPlainText">> function to change my working directory and to load some of my useful<o:p></o:p></p>
<p class="MsoPlainText">> scripts:<o:p></o:p></p>
<p class="MsoPlainText">> .First <- function() {<o:p></o:p></p>
<p class="MsoPlainText">> setwd("my_path")<o:p></o:p></p>
<p class="MsoPlainText">> source("my_script1.r")<o:p></o:p></p>
<p class="MsoPlainText">> #library("BrailleR")<o:p></o:p></p>
<p class="MsoPlainText">> }<o:p></o:p></p>
<p class="MsoPlainText">> Take care to the dot before the function name and the first letter
<o:p></o:p></p>
<p class="MsoPlainText">> capitalized, are part of the R config requirements.<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> Good luck<o:p></o:p></p>
<p class="MsoPlainText">> --<o:p></o:p></p>
<p class="MsoPlainText">> G.Voica<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> _______________________________________________<o:p></o:p></p>
<p class="MsoPlainText">> BlindRUG mailing list<o:p></o:p></p>
<p class="MsoPlainText">> BlindRUG@nfbnet.org<o:p></o:p></p>
<p class="MsoPlainText">> http://nfbnet.org/mailman/listinfo/blindrug_nfbnet.org<o:p></o:p></p>
<p class="MsoPlainText">> To unsubscribe, change your list options or get your account info for<o:p></o:p></p>
<p class="MsoPlainText">> BlindRUG:<o:p></o:p></p>
<p class="MsoPlainText">> http://nfbnet.org/mailman/options/blindrug_nfbnet.org/tyler.ensor%40mu<o:p></o:p></p>
<p class="MsoPlainText">> n.ca<o:p></o:p></p>
<p class="MsoPlainText">> The list archive can be viewed at:<o:p></o:p></p>
<p class="MsoPlainText">> http://www.nfbnet.org/pipermail/blindrug_nfbnet.org<o:p></o:p></p>
<p class="MsoPlainText">> More information and useful links about using R as a blind person can
<o:p></o:p></p>
<p class="MsoPlainText">> be obtained at:<o:p></o:p></p>
<p class="MsoPlainText">> http://R-Resources.massey.ac.nz<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> Look for help using R commands by reading the accessible e-book "Let's
<o:p></o:p></p>
<p class="MsoPlainText">> Use R Now" compiled by Jonathan Godfrey at:<o:p></o:p></p>
<p class="MsoPlainText">> http://R-Resources.massey.ac.nz/lurn/front.html<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">_______________________________________________<o:p></o:p></p>
<p class="MsoPlainText">BlindRUG mailing list<o:p></o:p></p>
<p class="MsoPlainText">BlindRUG@nfbnet.org<o:p></o:p></p>
<p class="MsoPlainText">http://nfbnet.org/mailman/listinfo/blindrug_nfbnet.org<o:p></o:p></p>
<p class="MsoPlainText">To unsubscribe, change your list options or get your account info for BlindRUG:<o:p></o:p></p>
<p class="MsoPlainText">http://nfbnet.org/mailman/options/blindrug_nfbnet.org/a.j.godfrey%40massey.ac.nz<o:p></o:p></p>
<p class="MsoPlainText">The list archive can be viewed at:<o:p></o:p></p>
<p class="MsoPlainText">http://www.nfbnet.org/pipermail/blindrug_nfbnet.org<o:p></o:p></p>
<p class="MsoPlainText">More information and useful links about using R as a blind person can be obtained at:<o:p></o:p></p>
<p class="MsoPlainText">http://R-Resources.massey.ac.nz<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Look for help using R commands by reading the accessible e-book "Let's Use R Now" compiled by Jonathan Godfrey at:<o:p></o:p></p>
<p class="MsoPlainText">http://R-Resources.massey.ac.nz/lurn/front.html<o:p></o:p></p>
</div>
</body>
</html>