# NOT RUN {
library(ggplot2)
Nmc=500
s<-c(1e3,4e3)
Inference<-function(s){
rrr<-rnorm(n=s)
list(Mn=mean(rrr), Sd=sd(rrr))
}
data <- Estim_diagnost(Nmc, s, Inference)
lisst <- estims_qqplot(data, sep=TRUE)
lisst[2][[1]] + geom_abline(intercept = 1)
pl_joint<-estims_qqplot(data)
pl_joint + geom_abline(slope=1)
pl_joint<-estims_qqplot(data, distribution = stats::qt, dparams = list(df=3, ncp=0.1))
pl_joint + geom_abline(slope=1)
# }
Run the code above in your browser using DataLab