R help - Aug 2005 - chisq.,test`

Hi
 
 I am trying to do this:
 
chisq.test(c(11, 13, 12, 18, 21, 43, 15, 12, 9, 10, 5, 28, 22, 11, 15,
11, 18, 28, 16, 8, 15, 19, 44, 18, 11, 23, 15, 23, 2, 5, 4, 14, 3, 22,
9, 0, 6, 19, 15, 32, 3, 16, 14, 10, 24, 16, 24, 31, 29, 28, 16, 26, 11,
11, 4, 17, 16, 13, 20, 26, 16, 19, 34, 19, 17, 14, 22, 25, 17, 12, 23,
14, 19, 30, 18, 10, 23, 21, 17, 16, 10, 14, 6, 17, 17, 10, 21, 25, 20,
4, 11, 4, 11, 12, 21, 13, 12, 21, 8, 20, 8, 13, 18, 24, 19, 23, 17, 17,
10, 12, 5, 16, 1, 8, 14, 39, 13, 6, 6, 23, 11, 23, 21, 16, 18, 12, 16,
6, 8, 5, 9, 10, 9, 11, 20, 8, 12, 0, 14, 24, 30, 17, 17, 18, 43, 14, 36,
15, 20, 10, 19, 14, 19, 28, 21, 18, 24, 16, 13, 16, 11, 24, 28, 14, 23,
11, 19, 25, 26, 34, 16, 9, 18, 13, 27, 14, 16, 21, 10, 18, 23, 16, 39,
16, 35, 24, 32, 23, 19, 28, 19, 21, 32, 16, 21, 54, 23, 42, 33, 27, 19,
26, 22, 22, 19, 19, 19, 45, 14, 41, 15, 18, 27, 25, 23, 67, 7, 12, 26,
25, 26, 22, 12, 125, 70, 48, 17, 20, 34, 15, 25, 45, 32, 42, 29, 107,
14, 17, 39, 34, 45, 65, 46, 27, 31, 28, 26, 27, 24, 39, 26, 16, 28, 28,
17, 32, 32, 35, 25, 44, 21, 20, 21, 32, 30, 29, 18, 86, 32, 56, 30, 31,
35, 17, 19, 27, 22, 23, 44, 15, 18, 34, 32, 21, 32, 22, 19, 34, 21, 8,
19, 32, 31, 35, 19, 16, 17, 36, 39, 29, 31, 42, 30, 42, 22, 25, 20, 52,
21, 36, 27, 13, 27, 118, 18, 12, 19, 16, 60, 22, 25, 22, 24, 11, 53, 1,
27, 34, 39, 43, 22, 21, 36, 43, 20, 49, 45, 34, 95, 10, 13, 27, 21, 20,
39, 26, 27, 23, 24, 20, 21, 14, 26, 31, 31, 23, 24, 19, 52, 50, 26, 16,
78, 29, 34, 15, 59, 20, 43, 10, 29, 41, 23, 26, 26, 27, 8, 35, 34, 16,
42, 31, 35, 36, 33, 43, 19, 24, 21, 24, 18, 32, 33, 35, 17, 13, 25, 34,
7, 54, 34, 92, 27, 24, 19, 13, 16, 22, 18, 15, 19, 17, 31, 14, 32),
p=c(0.0016, 0.002752, 0.001728, 0.0016, 0.001792, 0.005953, 0.006081,
0.001792, 0.002432, 0.0016, 0.001984, 0.002624, 0.001984, 0.002432,
0.001728, 0.002304, 0.002176, 0.003008, 0.001792, 0.0016, 0.0016,
0.001664, 0.0032, 0.003328, 0.0016, 0.002112, 0.001472, 0.001728,
0.001728, 0.002432, 0.001856, 0.0016, 0.001536, 0.002624, 0.003072,
0.001664, 0.002624, 0.003264, 0.004417, 0.002816, 0.001536, 0.00288,
0.002048, 0.002368, 0.002624, 0.002048, 0.003008, 0.002432, 0.003008,
0.003584, 0.001856, 0.00256, 0.001728, 0.001856, 0.002048, 0.00288,
0.001984, 0.0016, 0.001728, 0.002432, 0.001472, 0.001792, 0.003328,
0.003136, 0.00256, 0.002496, 0.00192, 0.002048, 0.00192, 0.001792,
0.002432, 0.0016, 0.001728, 0.002624, 0.0016, 0.001728, 0.001728,
0.001664, 0.002304, 0.001472, 0.002176, 0.0016, 0.000512, 0.001792,
0.001856, 0.002432, 0.0016, 0.001856, 0.001856, 0.001664, 0.001856,
0.00192, 0.00128, 0.001664, 0.00256, 0.001728, 0.001856, 0.002176,
0.001472, 0.001728, 0.001664, 0.001408, 0.001728, 0.002688, 0.00224,
0.004289, 0.001728, 0.001408, 0.001984, 0.001344, 0.001664, 0.001024,
0.001984, 0.001536, 0.0016, 0.008066, 0.001152, 0.001472, 0.001536,
0.002048, 0.001792, 0.003969, 0.002432, 0.00192, 0.003008, 0.003008,
0.00224, 0.001856, 0.001344, 0.0016, 0.00224, 0.002176, 0.001664,
0.0016, 0.001664, 0.001664, 0.001664, 0.00128, 0.002048, 0.002432,
0.002944, 0.001728, 0.00192, 0.001984, 0.003328, 0.001536, 0.002624,
0.001984, 0.001728, 0.002176, 0.00192, 0.001856, 0.001984, 0.00224,
0.001856, 0.0016, 0.001792, 0.001344, 0.001728, 0.001856, 0.001536,
0.001792, 0.00256, 0.002112, 0.002816, 0.001792, 0.001728, 0.001984,
0.002048, 0.002752, 0.001984, 0.001728, 0.001536, 0.002432, 0.00224,
0.001664, 0.001664, 0.001856, 0.0016, 0.002176, 0.001792, 0.001472,
0.003648, 0.00192, 0.003264, 0.001792, 0.003136, 0.00224, 0.002112,
0.002368, 0.001664, 0.001792, 0.003328, 0.002368, 0.001792, 0.007553,
0.002944, 0.004033, 0.002752, 0.002944, 0.001856, 0.002304, 0.00224,
0.001728, 0.0016, 0.001344, 0.002368, 0.003392, 0.001856, 0.003584,
0.001408, 0.001984, 0.003136, 0.00192, 0.00288, 0.005633, 0.001408,
0.001664, 0.00288, 0.0032, 0.002752, 0.002304, 0.001664, 0.00973,
0.005633, 0.003392, 0.001664, 0.001792, 0.002432, 0.002048, 0.001984,
0.003392, 0.002368, 0.004929, 0.002176, 0.012739, 0.001984, 0.001664,
0.00352, 0.004161, 0.004993, 0.004865, 0.004417, 0.002368, 0.0032,
0.001792, 0.002112, 0.00224, 0.0016, 0.003328, 0.002304, 0.001344,
0.001856, 0.003072, 0.001664, 0.002432, 0.002496, 0.002688, 0.002176,
0.003008, 0.001664, 0.001856, 0.0016, 0.002624, 0.002176, 0.002368,
0.001472, 0.008962, 0.002496, 0.003904, 0.002368, 0.002496, 0.002624,
0.001472, 0.002048, 0.002048, 0.001536, 0.002304, 0.003072, 0.001536,
0.0016, 0.002368, 0.002304, 0.00192, 0.002304, 0.0016, 0.001536,
0.003008, 0.0016, 0.001984, 0.002176, 0.004993, 0.002368, 0.002304,
0.001344, 0.002368, 0.001728, 0.003136, 0.004033, 0.00224, 0.002112,
0.002944, 0.001984, 0.003392, 0.002048, 0.001984, 0.001728, 0.003776,
0.001856, 0.002752, 0.002624, 0.0016, 0.002816, 0.012163, 0.00352,
0.001344, 0.001728, 0.001856, 0.005505, 0.001664, 0.001728, 0.0016,
0.001664, 0.002112, 0.00384, 0.001792, 0.00192, 0.002368, 0.002752,
0.004097, 0.0016, 0.001536, 0.00256, 0.003264, 0.002112, 0.004353,
0.003264, 0.00224, 0.00845, 0.001408, 0.001536, 0.002112, 0.001472,
0.002368, 0.0032, 0.001728, 0.00224, 0.001472, 0.001728, 0.002112,
0.002176, 0.002112, 0.002304, 0.002816, 0.003072, 0.001728, 0.001856,
0.001216, 0.004353, 0.003264, 0.002496, 0.001536, 0.005633, 0.002112,
0.002624, 0.001792, 0.004161, 0.001792, 0.004353, 0.001344, 0.00256,
0.002752, 0.0016, 0.001856, 0.002112, 0.003456, 0.001408, 0.002432,
0.002496, 0.001472, 0.003392, 0.001984, 0.002496, 0.003136, 0.002304,
0.003264, 0.001344, 0.00224, 0.001536, 0.001984, 0.001728, 0.00256,
0.002944, 0.002688, 0.001856, 0.001216, 0.001856, 0.00224, 0.0016,
0.003969, 0.002176, 0.006529, 0.001856, 0.00256, 0.00192, 0.001536,
0.00192, 0.001664, 0.002112, 0.001792, 0.002304, 0.00224, 0.002112,
0.001664, 0.002304))
 
but I keep on getting ‘syntax error’.  Is this to long, or somehow badly
formed?

-- 
No virus found in this outgoing message.
Checked by AVG Anti-Virus.

22/08/2005
 

	[[alternative HTML version deleted]]

On Fri, 2005-08-26 at 18:15 +1000, Stephen Choularton
wrote:
> Hi
>  
>  I am trying to do this:
I get no syntax error on R 2.1.1-patched. I do get an error though:

Error in chisq.test(c(11, 13, 12, 18, 21, 43, 15, 12, 9, 10, 5, 28,
22,  :
        probabilities must sum to 1.

Which leads us to question the values you provided as argument p (which
I assigned to a vector P first)

sum(P)
[1] 0.999783

So you need to make sure your probabilities add sum to one.

HTH

Ps. It would have been helpful if you'd posted the syntax error you
received and your version of R - type version at the R prompt.

Gav

<snip> 
> but I keep on getting syntax error.  Is this to long, or somehow badly
> formed?
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
-- 
%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
Gavin Simpson                     [T] +44 (0)20 7679 5522
ENSIS Research Fellow             [F] +44 (0)20 7679 7565
ENSIS Ltd. & ECRC                 [E] gavin.simpsonATNOSPAMucl.ac.uk
UCL Department of Geography       [W] http://www.ucl.ac.uk/~ucfagls/cv/
26 Bedford Way                    [W] http://www.ucl.ac.uk/~ucfagls/
London.  WC1H 0AP.
%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%

On 26-Aug-05 Stephen Choularton wrote:
> Hi
>  
>  I am trying to do this:
>  
> chisq.test(c(11, 13, 12, 18, 21, 43, 15, 12, 9, 10, 5, 28, 22, 11, 15,
> [...]
> 7, 54, 34, 92, 27, 24, 19, 13, 16, 22, 18, 15, 19, 17, 31, 14, 32),
> p=c(0.0016, 0.002752, 0.001728, 0.0016, 0.001792, 0.005953, 0.006081,
> [...]
> 0.001664, 0.002304))
>  
> but I keep on getting ?syntax error?.  Is this to long, or somehow
> badly formed?
Having cut&pasted your data as published in your email into R-1.8.0
on Red Hat Linux 9, I encountered no "syntax error" and got the
result:

  data:  c(11, 13, 12, 18, 21, 43, 15, 12, 9, 10, 5, 28, 22, 11, 15,
  [...]

  X-squared = 1094.448, df = 414, p-value = < 2.2e-16

Fee: For each P-value of the form "x.ye-N", N beers.

Hoping this helps ... though your "syntax error" remains mysterious.

Best wishes,
Ted.


--------------------------------------------------------------------
E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk>
Fax-to-email: +44 (0)870 094 0861
Date: 26-Aug-05                                       Time: 09:53:38
------------------------------ XFMail ------------------------------