Re: [R] qqplot for binomial distribution

2017-05-11 Thread Ashim Kapoor
Dear All, when I do : set.seed(123) expected_distribution<-rbinom(1000,100,.05) #Without jitter qqplot(jitter(expected_distribution),count1_vector, xlab="Expected distribution",ylab="Observed values") qqline(count1_vector,distribution = function(probs) { qbinom(probs, size=100, prob=0.05) },co

Re: [R] qqplot for binomial distribution

2017-04-19 Thread Ashim Kapoor
Dear Boris, Many thanks, Ashim On Tue, Apr 18, 2017 at 7:56 PM, Boris Steipe wrote: > As per the help pages, the data samples are expected in the second > argument, "y". > > So try > qqplot(rbinom(n=100, size=100, p=0.05), count1_vector) > > ... and then plot your qqline() > > Alternatively,

Re: [R] qqplot for binomial distribution

2017-04-18 Thread Boris Steipe
As per the help pages, the data samples are expected in the second argument, "y". So try qqplot(rbinom(n=100, size=100, p=0.05), count1_vector) ... and then plot your qqline() Alternatively, try qqline(count1_vector, distribution = function(probs) { qbinom(probs, size=100, prob=0.05)

Re: [R] qqplot for binomial distribution

2017-04-17 Thread Ashim Kapoor
Dear Boris, Thank you for your reply. > dput(count1_vector) c(5, 6, 4, 4, 6, 5, 4, 5, 3, 7, 5, 5, 3, 4, 8, 6, 10, 2, 4, 6, 8, 4, 4, 6, 8, 5, 6, 3, 7, 9, 4, 7, 5, 7, 3, 4, 5, 2, 11, 7, 8, 5, 5, 6, 3, 2, 3, 5, 9, 6, 5, 6, 7, 3, 10, 7, 6, 4, 9, 5, 7, 3, 7, 3, 2, 3, 4, 5, 10, 4, 5, 5, 6, 7, 4, 8, 7,

Re: [R] qqplot for binomial distribution

2017-04-17 Thread Boris Steipe
That's not how qqline() works. The line is drawn with respect to a _reference_distribution_ which is the normal distribution by default. For the binomial distribution, you need to specify the distribution argument. There is an example in the help page that shows you how this is done for qchisq()

Re: [R] qqplot for binomial distribution

2017-04-17 Thread Ashim Kapoor
Dear Boris, Okay and Thanks. Best, Ashim On Mon, Apr 17, 2017 at 6:45 PM, Boris Steipe wrote: > Moreover, setting the seed once, then evaluating two functions means you > are sampling from the same distributions, but you do in fact have different > values. Outliers in the rarefied tails of the

Re: [R] qqplot for binomial distribution

2017-04-17 Thread Boris Steipe
Moreover, setting the seed once, then evaluating two functions means you are sampling from the same distributions, but you do in fact have different values. Outliers in the rarefied tails of the distribution may lie quite considerably off the expected diagonal. Try set.seed(123) qqplot(rbinom(n

Re: [R] qqplot for binomial distribution

2017-04-17 Thread Ashim Kapoor
Dear Spencer, Okay. Many thanks. My next query is how do I use qqline? When I try > qqline(rbinom(n=100,size=100,p=.05)) I don't get the line in the right place. Best Regards, Ashim On Mon, Apr 17, 2017 at 6:31 PM, Spencer Graves < spencer.gra...@effectivedefense.org> wrote: > > > On 2017-04

Re: [R] qqplot for binomial distribution

2017-04-17 Thread Spencer Graves
On 2017-04-17 7:58 AM, Ashim Kapoor wrote: Dear All, set.seed(123) qqplot(rbinom(n=100,size=100,p=.05), rbinom(n=100,size=100,p=.05) ) I expect to see 1 clear line,but I don't. What am I misunderstanding? The distribution is discrete, and points are superimposed. Try the following:

[R] qqplot for binomial distribution

2017-04-17 Thread Ashim Kapoor
Dear All, set.seed(123) qqplot(rbinom(n=100,size=100,p=.05), rbinom(n=100,size=100,p=.05) ) I expect to see 1 clear line,but I don't. What am I misunderstanding? Best Regards, Ashim [[alternative HTML version deleted]] __ R-help@r-project.org