Hi Peter,
Thank you so much for your explaination. I draw histogram before that, so
maybe that warning messages are meant for that.
From: peter dalgaard
Cc: "r-help@r-project.org"
Sent: Friday, August 2, 2013 3:11 PM
Subject: Re: [R] Kolmogor
On Aug 2, 2013, at 03:24 , Roslina Zakaria wrote:
> Dear r-users,
>
> I am using KS test to test the goodness of fit for my data and the got the
> following output. However, I don't understand about the "warning messages".
> What does it mean by "horizontals" is not a graphical parameter"
>
] Kolmogorov-Smirnov test and the plot of max distance between two
ecdf curves
Sent by:
r-help-boun...@r-project.org
Rui,
Your response nearly answered a similar question of mine except that I
also
have ecdfs of different lengths.
Do you know how I can adjust x <- seq(min(loga, logb), max(l
Hello,
Try length.out = max(length(loga), length(logb))
Note also that all of the previous code and the line above assumes that
we are interested in the max distance, whereas the KS statistic computes
the supremum of the distance. If it's a two sample test then their
values are almost surely
Rui,
Your response nearly answered a similar question of mine except that I also
have ecdfs of different lengths.
Do you know how I can adjust x <- seq(min(loga, logb), max(loga, logb),
length.out=length(loga))
to account for this? It must be in length.out() but I'm unsure how to
proceed.
Just a final correction.
I was wrong, stats::ks.test doesn't use only Marsaglia et al.
It's even clearly written in the help page.
Read the documentation before stating!
Rui Barradas
Em 28-05-2012 11:51, maxbre escreveu:
> thanks for the help: I'll have a look at the papers
> max
>
> Il 28/05/20
thanks for the help: I'll have a look at the papers
max
Il 28/05/2012 12:31, Rui Barradas [via R] ha scritto:
> Hello,
>
> That's a very difficult question. See
>
> Marsaglia, Tsang, Wang (2003)
> http://www.jstatsoft.org/v08/i18/
>
> Simard, L'Ecuyer (2011)
> http://www.jstatsoft.org/v39/i11
>
>
Hello,
That's a very difficult question. See
Marsaglia, Tsang, Wang (2003)
http://www.jstatsoft.org/v08/i18/
Simard, L'Ecuyer (2011)
http://www.jstatsoft.org/v39/i11
R's ks functions are a port of Marsaglia et al. to the .C interface.
Rui Barradas
maxbre wrote
>
> thanks rui
>
> that's what
thanks rui
that's what I was looking for
I have another related question:
- why of the difference between the max distance D calculated with ks.test()
and the max distance D “manually” calculated as in (2)?
I guess it has something to do with the fact that KS is obtained with a
maximisation th
Hello,
Try the following.
(i've changed the color of the first ecdf.)
loga <- log10(a+1) # do this
logb <- log10(b+1) # only once
f.a <- ecdf(loga)
f.b <- ecdf(logb)
# (2) max distance D
x <- seq(min(loga, logb), max(loga, logb), length.out=length(loga))
x0 <- x[which( abs(f.a(x) - f.b(x)) ==
I recently gave a presentation at the 50th Army Operational Research
Symposium at Ft Lee describing an implementation of Conover's exact
calculation method for the KS test applied to discrete distributions. My
implementation was done in Matlab script as opposed to R. Multiple
Monte-Carlo trials wer
801.408.8111
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of Greg Snow
> Sent: Monday, September 26, 2011 11:45 AM
> To: rommel; r-help@r-project.org
> Subject: Re: [R] Kolmogorov-Smirnov test
>
> There
2:30 AM
> To: r-help@r-project.org
> Subject: Re: [R] Kolmogorov-Smirnov test
>
> Dear Dr. Snow,
>
> Thank you for your reply.
>
> 1. Are you doing the 2 sample KS test? Comparing if 2 samples come from
> the same distribution? -Yes, I am doing 2-sample KS test
>
nes you show. You will
need to ask someone who knows the programs you reference to determine what
input they are expecting. R expects the raw data. -Original
Message- From: [hidden email] [mailto: [hidden email] ] On Behalf Of
rommel Sent: Friday, September 23, 2011 7:51 AM To: [hidden e
...@r-project.org] On
Behalf Of rommel
Sent: Friday, September 23, 2011 7:51 AM
To: r-help@r-project.org
Subject: Re: [R] Kolmogorov-Smirnov test
Dear Dr. Snow,
I would like to ask for help on my three questions regarding Kolmogorov
Smirnov test.
1.
'With a sample size over 10,000 you will have power
Dear Dr. Snow,
I would like to ask for help on my three questions regarding Kolmogorov
Smirnov test.
1.
'With a sample size over 10,000 you will have power to detect differences
that are not practically meaningful. '
-Is sample size of 3000 for each sample okay for using Kolmogorov
Smirnov t
Hi,
many thanks for helpful answer.
Best
Marcin M.--
View this message in context:
http://r.789695.n4.nabble.com/Kolmogorov-Smirnov-test-tp3479506p3488364.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing l
l 28, 2011 3:54 PM
> To: r-help@r-project.org
> Subject: Re: [R] Kolmogorov-Smirnov test
>
> Hi,
> thanks for response.
>
> >> The Kolmogorov-Smirnov test is designed for distributions on
> continuous
> >> variable, not discrete like the >> poisson.
This test SnowsPenultimateNormalityTest() is great :)
Best
--
View this message in context:
http://r.789695.n4.nabble.com/Kolmogorov-Smirnov-test-tp3479506p3482401.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org ma
Hi,
thanks for response.
>> The Kolmogorov-Smirnov test is designed for distributions on continuous
>> variable, not discrete like the >> poisson. That is why you are getting
>> some of your warnings.
I read in "Fitting distributions whith R" Vito Ricci page 19 that: "...
Kolmogorov-Smirnov t
A couple of things to consider:
The Kolmogorov-Smirnov test is designed for distributions on continuous
variable, not discrete like the poisson. That is why you are getting some of
your warnings.
With a sample size over 10,000 you will have power to detect differences that
are not practically
gt; ks.test::ks.test(c(0,1), ecdf(c(0,1)))
One-sample Kolmogorov-Smirnov test
data: c(0, 1)
D = 0, p-value = 1
alternative hypothesis: two-sided
Original Message:
Date: Mon, 28 Feb 2011 21:31:26 +1100
From: Glen Barnett
To: tsippel
Cc: r-help@r-project.org
Subject: Re: [R] Kolmogorov-s
It's designed for continuous distributions. See the first sentence here:
http://en.wikipedia.org/wiki/Kolmogorov%E2%80%93Smirnov_test
K-S is conservative on discrete distributions
On Sat, Feb 19, 2011 at 1:52 PM, tsippel wrote:
> Is the kolmogorov-smirnov test valid on both continuous and discr
The KS test was designed for continuous variables. The vcd package has tools
for exploring categorical variables and distributions.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Message-
> From: r-help-boun
.D.
> Statistical Data Center
> Intermountain Healthcare
> greg.s...@imail.org801.408.8111begin_of_the_skype_highlighting 801.408.8111 end_of_the_skype_highlighting
>
>
>
> > -Original Message-
> > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
>
-Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
> project.org] On Behalf Of Kerry
> Sent: Thursday, November 11, 2010 12:02 AM
> To: r-help@r-project.org
> Subject: Re: [R] Kolmogorov Smirnov Test
>
> Thanks for the feedback. My goal is to run a
On 11-Nov-10 04:22:55, Kerry wrote:
> I'm using ks.test (mydata, dnorm) on my data.
I think your problem may lie here! If you look at the documentation
for ks.test, available with the command:
help("ks.test")
or simply:
?ks.test
you will read the following near the beginning:
Usage: ks.test(
Thanks for the feedback. My goal is to run a simple test to show that
the data cannot be rejected as either normally or uniformally
distributed (depening on the variable), which is what a previous K-S
test run using SPSS had shown. The actual distribution I compare to my
sample only matters that it
The way you are running the test the null hypothesis is that the data comes
from a normal distribution with mean=0 and standard deviation = 1. If your
minimum data value is 0, then it seems very unlikely that the mean is 0. So
the test is being strongly influenced by the mean and standard devi
It is not clear what question you are trying to answer. Perhaps if you can
give us an explanation of your overall goal then we can be more helpful.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Message-
>
Hi Roslina,
I believe that you can ignore the warning.
Alternatively, you may add a very small random noise to pairs with ties, i.e.
something like
xobs[which(duplicated(xobs))] <- xobs[which(duplicated(xobs))] +
1.0e-6*sd(xobs)*rnorm(length(which(duplicated(xobs
Regards,
Moshe.
--- On Tu
This is the third homework question you have asked the list to do for you. How
many more should we expect?
The posting guide is pretty clear in that: " Basic statistics and classroom
homework: R-help is not intended for these."
-Original Message-
From: r-help-boun...@r-project.org [m
help.search("kolmogorov")
?ks.test
andydol...@gmail.com
2009/4/29 mathallan
>
> I got a distribution function and a empirical distribution function. How do
> I
> make to Kolmogorov-Smirnov test in R.
>
> Lets call the empirical distribution function >Fn on [0,1]
> and the
33 matches
Mail list logo
