[R] Issues when trying to fit a curve when generating random deviates

Rui Barradas
Mon, 19 Sep 2022 04:44:29 -0000

Hello,

I cannot reproduce the error.

Anyway, if the data is discrete a bar plot might be more appropriate.


barplot(table(x),
         col = "green",
         main = "Bar plot of Binomial Distribution B(20,0.4)")


Or, to have proportions,

barplot(table(x)/length(x), etc)


Hope this helps,

Rui Barradas

?s 05:01 de 19/09/2022, Paul Bernal escreveu:
> Dear friends,
>
> Hope you are doing great.
>
> I first generated random deviates from a binomial distribution with the
> following code (I had to generate from a B(20,0.4) dist):
> #Setting seed
> set.seed(1234567)
>
> #Generating 1000 random deviates from a B(20,0.4) Distribution
> x <- rbinom(1000,20,0.4)
> #Generating histogram for x
> hist(x, col = c("green"),main="Histogram of Binomial
Distribution
> B(20,0.4)",freq=F)
> #Generating the curve for x
> There were 50 or more warnings (use warnings() to see the first 50)
> warnings()
> #the first 5 warnings are shown, but all the warnings are similar
> Warning messages:
> 1: In dbinom(x, 1000, 0.4) : non-integer x = 2.120000
> 2: In dbinom(x, 1000, 0.4) : non-integer x = 2.240000
> 3: In dbinom(x, 1000, 0.4) : non-integer x = 2.360000
> 4: In dbinom(x, 1000, 0.4) : non-integer x = 2.480000
> 5: In dbinom(x, 1000, 0.4) : non-integer x = 2.600000
> and the curve was not generated for the histogram.
>
> Any ideas on what I could be doing wrong?
>
> Best regards,
>
> Paul
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.