similar to: Numerical stability in chisq.test

The chisq.test on line 57 contains following code: STATISTIC <- sum(sort((x - E)^2/E, decreasing = TRUE)) However, based on book "Accuracy and stability of numerical algorithms" available from: 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

>>>>> 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,

> 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 >

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

Hello, I can''t seem to get the network going. I found this message in the log of Domain0: ----------------------------------------------------------------------- Dec 1 15:39:18 zms wait_for_sysfs[4180]: either wait_for_sysfs (udev 039) needs an update to handle the device ''/class/net/br0'' properly (no device symlink) or the sysfs-support of your device''s

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? >

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

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

Hello, This is to announce the first release of the "ant" R package, which has been pushed to CRAN yesterday, and should reach your mirror and your platform soon. The package provides an R-aware version of the famous build tool from the apache project. http://ant.apache.org/ The package ships an R script that can be used to invoke ant with enough plumbing so that it can use R

Hello, This is to announce the first release of the "ant" R package, which has been pushed to CRAN yesterday, and should reach your mirror and your platform soon. The package provides an R-aware version of the famous build tool from the apache project. http://ant.apache.org/ The package ships an R script that can be used to invoke ant with enough plumbing so that it can use R

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

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

> -----Original Message----- > From: Duncan Murdoch [mailto:murdoch at stats.uwo.ca] > Sent: Tuesday, August 04, 2009 5:06 PM > To: Steven McKinney > Cc: 'r-devel at stat.math.ethz.ch'; 'R-bugs at r-project.org' > Subject: Re: [Rd] Wishlist: Navigate to "Index" page of help when no > topic specified (PR#13860) >=20 > On 04/08/2009 7:33 PM,

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

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

>>>>> 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