Login [Register]
Don't have an account? Register now to chat, post, use our tools, and much more.
This is an archived, read-only copy of the United-TI subforum , including posts and topic from May 2003 to April 2012. If you would like to discuss any of the topics in this forum, you can visit Cemetech's TI-BASIC subforum. Some of these topics may also be directly-linked to active Cemetech topics. If you are a Cemetech member with a linked United-TI account, you can link United-TI topics here with your current Cemetech topics.

This forum is locked: you cannot post, reply to, or edit topics. TI-Basic => TI-BASIC
Author Message
thornahawk
μολών λαβέ


Active Member


Joined: 27 Mar 2005
Posts: 569

Posted: 27 Nov 2009 08:46:26 pm    Post subject:

Grubbs's test is a quick method for testing if a certain entry in a small data set is an outlier; i.e. this test answers the question "is this entry very much different from its companions?"

The test is relatively simple to perform compared to an older outlier test, Dixon's Q test, and is now the test recommended by ISO.

Here is a short routine for performing Grubbs's test on a TI-84+SE calculator:


Code:
PROGRAM:GRUBBS
Input "DATA:",L₁
Prompt I
Input "CONFIDENCE",P
abs(L₁(I)mean(L₁))/stdDev(L₁)→Q
dim(L₁)→N
invT(1+.5(P1)/N,N2)→C
(N1)(C/√(C+N2))/√(N)→C
Q≤C


An output of 0 signifies L₁(I) is an outlier at the confidence level P. C is the Grubbs critical value, and Q is the test statistic. A data set to test the program can be found here.

A (slightly longer) version of this program can be done for the TI-83+ and other calculators without the invT( function, which I will post on request.

thornahawk
Back to top
Graphmastur


Advanced Member


Joined: 25 Mar 2009
Posts: 360

Posted: 27 Nov 2009 09:39:35 pm    Post subject:

Oh, so it's like figuring out which point is bad data for a test. Sounds interesting.
Back to top
Galandros


Active Member


Joined: 29 Aug 2008
Posts: 565

Posted: 28 Nov 2009 03:45:51 am    Post subject:

Nice program. Maybe I can use in physics experiments soon...
Back to top
Display posts from previous:   
Register to Join the Conversation
Have your own thoughts to add to this or any other topic? Want to ask a question, offer a suggestion, share your own programs and projects, upload a file to the file archives, get help with calculator and computer programming, or simply chat with like-minded coders and tech and calculator enthusiasts via the site-wide AJAX SAX widget? Registration for a free Cemetech account only takes a minute.

» Go to Registration page
    »
» View previous topic :: View next topic  
Page 1 of 1 » All times are GMT - 5 Hours

 

Advertisement