在模擬卡方置信區間內繪製卡方（1）的hist的困難

我想通過casella的統計推斷來編碼例子10.3.2。（我在這裏附上了這個例子）。在模擬卡方置信區間內繪製卡方（1）的hist的困難

我有產生相同的情節問題。任何幫助？

模擬數據和對比表：

n=25 
lam<-5 
nsim<-10000 
set.seed(442256)     

poisson<-function(nsim,n,lam){ 
    ratio<-c() 

distributionMean = NULL 
for (i in 1 : nsim) distributionMean = c(distributionMean, mean(rpois(n, lam))) 

d<- 2*n*((lam-distributionMean)-distributionMean*log(lam/distributionMean)) 
     ratio<-c(ratio,d) 
     return(ratio) 
    } 

logLi<-poisson(10000,25,5) 

m<-matrix(0,2, 4) 
m[1,1]=quantile(p1,0.80) 
m[2,1]=qchisq(.80, df=1) 

m[1,2]=quantile(p1,0.90) 
m[2,2]=qchisq(.90, df=1) 

m[1,3]=quantile(p1,0.95) 
m[2,3]=qchisq(.95, df=1) 

m[1,4]=quantile(p1,0.99) 
m[2,4]=qchisq(.99, df=1) 
  row.names(m)<-c("simulated", "Chi-square") 
  colnames(m)<-c("80_perc", "90_perc","95_perc","99_perc")

來源

2017-04-02 sxq2221

它看起來像p1變量被分配之前使用 – tagoma

你想表？嘗試這個。

simulated <- quantile(logLi, c(0.8, 0.9, 0.95, 0.99)) 
chisquare <- qchisq(c(0.8, 0.9, 0.95, 0.99), df = 1) 

rbind(simulated, chisquare)

該圖的代碼如下。

a <- hist(logLi, freq=FALSE, xlim = c(0,4), 
      breaks = seq(0, ceiling(max(logLi)), by = 0.1)) 
lines(a$mids, dchisq(a$mids, df = 1))

來源

2017-04-02 12:15:04

非常感謝您清理我的表和幫助，與圖！有用***：） – sxq2221

在模擬卡方置信區間內繪製卡方（1）的hist的困難

回答

相關問題