x <- 5*rnorm(1001)^2+rnorm(1001)^2
x.thin<-sort(x[1+(0:50)*20])
p.invert<-pchisqsum(x.thin,df=c(1,1),a=c(5,1),method="integration")
p.satt<-pchisqsum(x.thin,df=c(1,1),a=c(5,1),method="satt")
plot(p.invert, p.satt,type="l")
abline(0,1,lty=2,col="purple")
Run the code above in your browser using DataCamp Workspace