out.R

In artitzco/prueba:

function(fun){
fit<-list()
prob<-list()
fit[[1]]<-function(x) 0.2+0.652900429496863*fun(x, shape=0.929303051939395, scale=141.461421999742)
prob[[1]]<-function(x) fun(x, shape=0.929303051939395, scale=141.461421999742)
fit[[2]]<-function(x) 0.2+0.533671248768375*fun(x, shape=0.847708511937514, scale=106.368640252579)
prob[[2]]<-function(x) fun(x, shape=0.847708511937514, scale=106.368640252579)
fit[[3]]<-function(x) 0.2+0.799876259716282*fun(x, shape=0.984630372719014, scale=158.000796916026)
prob[[3]]<-function(x) fun(x, shape=0.984630372719014, scale=158.000796916026)
fit[[4]]<-function(x) 0.2+0.799943836105268*fun(x, shape=0.792643710111436, scale=208.438833722765)
prob[[4]]<-function(x) fun(x, shape=0.792643710111436, scale=208.438833722765)
return(list(fit=fit, prob=prob))
}