- 901 名前：１３２人目の素数さん mailto:sage [2021/04/23(金) 18:15:11.87 ID:Ma96BCT1.net]
- >>882
>871をラジアン表示にしてpdf曲線を重ねて作図。
https://i.imgur.com/UmgkGbh.png
期待値
> mn
[1] 1.741718
> mn*180/pi
[1] 99.79307
pdf <- function(x){
if(pi/3<x & x<pi/2) return(((3*x-pi)*cos(x)-sin(3*x))/sin(x)^2)
if(pi/2<=x & x<pi) return(((pi-x)*cos(x)+sin(x))/sin(x)^2)
else return(0)
}
pdf=Vectorize(pdf)
curve(pdf(x),add=TRUE,lwd=2)
mn=integrate(function(x) x*pdf(x),0,pi)$value
mn
最頻値
mode=optimize(pdf,c(1,3),maximum=TRUE)$max
mode
> mode
[1] 1.447019
> mode*180/pi
[1] 82.90809
|
 |
