使用R覆蓋t分佈到我的直方圖？

Question

如何使用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])) 

Answer 1

更改最後一行：

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

得到縮放正確，然後做

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

改變df和色彩，以任何你想要的值。