similar to: Numerical stability in chisq.test

> On 28 Dec 2017, at 13:08 , Kurt Hornik <Kurt.Hornik at wu.ac.at> wrote: > >>>>>> Jan Motl writes: > >> The chisq.test on line 57 contains following code: >> STATISTIC <- sum(sort((x - E)^2/E, decreasing = TRUE)) > > The preceding 2 lines seem relevant: > > ## Sorting before summing may look strange, but seems to be >

>>>>> Jan Motl writes: > The chisq.test on line 57 contains following code: > STATISTIC <- sum(sort((x - E)^2/E, decreasing = TRUE)) The preceding 2 lines seem relevant: ## Sorting before summing may look strange, but seems to be ## a sensible way to deal with rounding issues (PR#3486): STATISTIC <- sum(sort((x - E) ^ 2 / E,

The chisq.test contains following code: STATISTIC <- sum(sort((x - E)^2/E, decreasing = TRUE)) However, based on book Accuracy and stability of numerical algorithms <http://ftp.demec.ufpr.br/CFD/bibliografia/Higham_2002_Accuracy%20and%20Stability%20of%20Numerical%20Algorithms.pdf> Table 4.1 on page 89, it is better to sort the data in increasing order than in decreasing order, when the

Hi, This report deals with p-values coming from chisq.test using the simulate.p=TRUE option. The issue is numerical accuracy and was brought up in previous bug reports 3486 and 3896. The bug was considered fixed but apparently was only mostly fixed. Just the typical problem of two values that are mathematically equal not ending up numerically equivalent. Consider this series of three 2x2

Hola lista: Una duda rápida (espero). ¿Se puede hacer un "inlcude" o un "input" en markdown? sin emplear Rmarkdown o knitr.... solo puro markdown, de modo que por ejemplo github lo interprete o un visor básico de markdown. la idea es q un file muestre(contenga) todo pero tener la info distribuida en varios files. Gracias. -- Antonio Maurandi López Sec. Apoyo

Hola ¿Quizá por aquí? https://github.com/yihui/knitr/issues/268 El 13 de enero de 2017, 12:26, Rubén Gómez Antolí <lobo en mucharuina.com> escribió: > > > El 13/01/17 a las 11:50, Antonio Maurandi López escribió: > > Hola lista: > > > > Una duda rápida (espero). > > > > ¿Se puede hacer un "inlcude" o un "input" en markdown? >

So far the tutorials I''ve seen are for encrypting text and then decrypting it back. So for example you have a plaintext = "some text" you encrypt it and then you decrypt it depending on how you encrypted it. What I am looking for is to encrypt a whole file and then store it. for example I allow my users to upload .doc/.xls files. I want to encrypt those files so I would not

Full_Name: Tao Shi Version: 1.7.0 OS: Windows XP Professional Submission from: (NULL) (149.142.163.65) > x [,1] [,2] [1,] 149 151 [2,] 1 8 > c2x<-chisq.test(x, simulate.p.value=T, B=100000)$p.value > for(i in (1:20)){c2x<-c(c2x,chisq.test(x, simulate.p.value=T,B=100000)$p.value)} > c2tx<-chisq.test(t(x), simulate.p.value=T, B=100000)$p.value > for(i in

I want to create a general model ''Product'' that I will be able to store several attributes in it. For example the attributes :id :name :description are pretty stantard. So each product will have this attributes. However, I wanted to know if it is possible to create dynamic attributes. In the case of a cell phone device I might want to store the IMEI of the phone, this

Hi all, Hope someone could help me. I am trying to run automatically the conversion of an Rwn file to a tex file. I am using windows 7, and cygwin. I tried to run automatically the Sweave.sh script, in its the most recent version available at R webpage: http://cran.r-project.org/contrib/extra/scripts/Sweave.sh Unfortunately, I got this error message: =========================== Raquel at

Hello all; In my network I have a server with samba4 as AD DC and two domain members as file servers with samba4. One of then work property, but the other not. My samba4 AD DC version is compiled from sources: [root at gtmad ~]# samba -V Version 4.5.5 The samba4 as domain member (files server) are installing from .rpm packages of CentOS7. [root at gtmpve /]# uname --all Linux

.... I have a CentOS 5.n VM running on a Fedora 14 server/host, installed using virt-manager, pretty plain vanilla. I use it to compile binaries to run under RHEL/CentOS 5.n OS. I occasionally notice that when the VM gets paged out by the server, it takes several minutes to get it back in :-/ (see below). On the host, logged in to a shell through a terminal window, this A.M.: [wam at Q6600,

Full_Name: Doug Hershberger Version: 1.2.2 OS: Red Hat-Linux 7.0 Submission from: (NULL) (216.99.65.36) In the R FAQ I find the following entry explaining that R no longer has problems with memory. http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why%20does%20R%20run%20out%20of%20of%20memory%3f However in my installation: R Version 1.2.2 (2001-02-26) Installed from the red hat RPM on your

Hello, I know this is likely to give me a brute force attack hit, but the only thing anyone can accomplish by ssh-ing to my machine is to provide me with a tunnel into your machine. So don't bother. Anyway, my server machine is running this: /usr/bin/ssh -X -R ${port}:localhost:22 -o BatchMode=yes \ -o StrictHostKeyChecking=no ${user}@${my_home_machine} On my local machine: ssh -vvv -X

The Chisquare test is based upon a normal approx of the (essentially) binomial distribution for the data in question. Small EXPECTED (not observed) values (<5) suggest a asymetric distribution and potential errors in inferential conclusions. The alternative is the exact test, which calculates the exact probabilities of the observed distribution, or a more extreme one, given the constraining

Dear all, I have a question about the chisq.test command. As an option one can chose the computation of p-values by Monte-Carlo simulation (simulate.p.value=T). Is there any documentation available how this calculations are done and how this simulation based test behaves in small samples? Thanks Klaus Abberger University of Konstanz, Germany [[alternate HTML version deleted]]

chisq.test with simulate.p.value=TRUE uses the Patefield algorithm, this is not documented, and the original reference is not given, as it ought to be. The reference is: Patefield,W. M. (1981) An efficient method of generating r * c tables with given row and column totals (algorithm AS 159). Applied Statistics 30, 91-97. Kjetil Halvorsen

>>>>> dmurdoch writes: >> Date: Wed, 16 Jul 2003 01:27:25 +0200 (MET DST) >> From: shitao@ucla.edu >>> x >> [,1] [,2] >> [1,] 149 151 >> [2,] 1 8 >>> c2x<-chisq.test(x, simulate.p.value=T, B=100000)$p.value >>> for(i in (1:20)){c2x<-c(c2x,chisq.test(x, >> simulate.p.value=T,B=100000)$p.value)}

On Thu, 11 Dec 2003, Jeffrey Chang wrote: > Hello everybody, > > I'm running R 1.8.1 on both Linux and OS X compiled with gcc 3.2.2 and > 3.3, respectively. The following call seems to freeze the interpreter > on both systems: > > chisq.test(matrix(c(233, 580104, 3776, 5786104), 2, 2), > simulate.p.value=TRUE) > > By freeze, I mean, the function call never

>>>>> "Torsten" == Torsten Hothorn <torsten@hothorn.de> >>>>> on Thu, 11 Dec 2003 18:03:07 +0100 (CET) writes: Torsten> On Thu, 11 Dec 2003, Jeffrey Chang wrote: >> Hello everybody, >> >> I'm running R 1.8.1 on both Linux and OS X compiled with >> gcc 3.2.2 and 3.3, respectively. The following