相当不错                                                                        
                                       
                                                     
                                                     
                                             
                                                                 
                                                                                 

stdError <- sd(yhat)/sqrt(length(yhat))

> x<- rnorm(50, mean=100, sd=15)

> t.test(x, conf.level=0.99)