使用R将t分布覆盖到我的直方图上? [英] Overlay a t-distribution to my histogram using R?

问题描述

如何使用R在我的直方图中叠加t密度?这是我的功能:

How can I overlay a t-density to my histogram using R? Here's my function:

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]))

推荐答案

将最后一行更改为:

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

要获得正确的缩放比例，然后做

to get the scaling correct, then do

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

将df和颜色更改为所需的任何值.

changing the df and color to whatever values you want.

这篇关于使用R将t分布覆盖到我的直方图上?的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！