0

R を使用してヒストグラムに t 密度を重ねるにはどうすればよいですか? これが私の機能です：

simfun <- function(a=56.25102409,b=1.78977412,c=0.08664925,n=18,x1.sd=18.87671,x2.sd=18.87671,e.sd=18.87671) {
   X1 <- rnorm(n, mean=0, sd=x1.sd)
   X2 <- rnorm(n, mean=0, sd=x2.sd) 
   e <-  rnorm(n, mean=0, sd=e.sd)
   Z <- a+b*X1+c*X2+e 
   data.frame(X1,X2,Z)
}

statfun <- function(samples) {
    coef(lm(Z~X1+X2,data=samples))
}

library(plyr)
B=raply(1000,statfun(simfun()))

(hist(B[,2]))
4

2 に答える 2

2

最後の行を次のように変更します。

hist(B[,2], prob=TRUE)

スケーリングを正しくするには、次のようにします

curve( dt(x, df=15), add=TRUE, col='blue' )

と色を任意の値に変更しdfます。

于 2013-09-30T21:09:47.250 に答える