[R] Issues when trying to fit a curve when generating random deviates
Rui Barradas
ru|pb@rr@d@@ @end|ng |rom @@po@pt
Mon Sep 19 06:44:29 CEST 2022
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 using 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.
More information about the R-help
mailing list