18/12/04 18:16:27.40 nfS0Ph7I.net
直接二次モーメントをRに計算させても同じ。
d2w8 <- function(x){
c(mean=mean(x),var=var(x),w8=mean(x^2)/mean(x)/2) # w8=E[X2]/2E[X]
}
> d2w8(diff(tt))
mean var w8
6.666667 9.250000 3.950000
18/12/04 18:16:27.40 nfS0Ph7I.net
直接二次モーメントをRに計算させても同じ。
d2w8 <- function(x){
c(mean=mean(x),var=var(x),w8=mean(x^2)/mean(x)/2) # w8=E[X2]/2E[X]
}
> d2w8(diff(tt))
mean var w8
6.666667 9.250000 3.950000