[R] problems hist() and density

jim holtman jholtman at gmail.com
Sun Jul 26 20:27:11 CEST 2009


"sums to one"  I should have said.

On Sun, Jul 26, 2009 at 7:43 AM, Jan
Teichmann<jan.teichmann at googlemail.com> wrote:
> Hello,
>
> I have a problem with the hist() function and showing densities. The
> densities sum to 50 and not to 1! I use R version 2.9.1 (2009-06-26) and
> I load the seqinR library.
>
> My data is the following vector:
> [1] 0.1400000 0.2000000 0.2200000 0.2828283 0.1600000 0.1600000
> 0.3600000
> [8] 0.1600000 0.2200000 0.2600000 0.2000000 0.3000000 0.2200000
> 0.2342342
> [15] 0.1800000 0.2200000 0.1600000 0.2300000 0.2000000 0.2200000
> 0.2400000
> [22] 0.2000000 0.2200000 0.2200000 0.2600000 0.2000000 0.1600000
> 0.2200000
> [29] 0.2342342 0.2000000 0.2200000 0.2000000 0.2000000 0.1400000
> 0.1800000
> [36] 0.2200000 0.1600000 0.1600000 0.1400000 0.2200000 0.2000000
> 0.2871287
> [43] 0.2900000 0.2000000 0.1836735 0.2000000 0.2000000 0.2900000
> 0.2400000
> [50] 0.2200000 0.2800000 0.2000000 0.2745098 0.2200000 0.2300000
> 0.1800000
> [57] 0.2300000 0.1800000 0.2600000 0.2200000 0.2222222 0.2200000
> 0.2600000
> [64] 0.2200000 0.2200000 0.2600000 0.2200000 0.2000000 0.2200000
>
> I use the following command:
> tmp <- hist(data, freq=FALSE, plot=FALSE)
>
> and that's the result:
> $breaks
>  [1] 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 0.36
>
> $counts
>  [1] 10  4 15 19  8  5  2  5  0  0  1
>
> $intensities
>  [1]  7.2463754  2.8985507 10.8695652 13.7681159  5.7971014  3.6231884
>  [7]  1.4492754  3.6231884  0.0000000  0.0000000  0.7246377
>
> $density
>  [1]  7.2463754  2.8985507 10.8695652 13.7681159  5.7971014  3.6231884
>  [7]  1.4492754  3.6231884  0.0000000  0.0000000  0.7246377
>
> $mids
>  [1] 0.15 0.17 0.19 0.21 0.23 0.25 0.27 0.29 0.31 0.33 0.35
>
> $xname
> [1] "data"
>
> $equidist
> [1] TRUE
>
> attr(,"class")
> [1] "histogram"
>
> ______________________________________________
> R-help at r-project.org mailing list
> 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.
>
>



-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem that you are trying to solve?




More information about the R-help mailing list