estims_qqplot: QQ-plot of estimator empirical distributions

In EstimDiagnostics: Diagnostic Tools and Unit Tests for Statistical Estimators

Description

Plot QQ-plots of estimators' empirical distributions for different sample sizes.

Usage

estims_qqplot(data, sep = FALSE, ...)

Arguments

data	data frame returned by Estim_diagnost
sep	indicates whether all plots will be stacked together or returned as elements of a list
...	parameters to pass to stat_qq function

Value

ggplot2 object

Examples

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)

EstimDiagnostics documentation built on April 17, 2021, 1:08 a.m.