• [^] # Re: Tes résultats ne sont pas reproductibles

    Posté par . En réponse au journal Simulation situation en Italie. Évalué à 4.

    Bonjour,

    une proposition de code plus dans l'esprit de R (le langage c'est R, Rstudio n'est qu'un IDE pour R)

    # d: dataset, m: mean, sd: width, a: aplitude
normApplitude <- function(d, m, sd) {
 d[length(d)] / pnorm(length(d):length(d), mean = m, sd = sd)
}
normVariance <- function(d, m, sd, a) {
 var(d - pnorm(1:length(d), mean = m, sd = sd) * a)
}
mTab = c(length(d):(length(d)+40))
sdTab = c(4:10)
print(mTab)
print(sdTab)
dta<-expand.grid(m=mTab,sd=sdTab,a=5000:20000)
dta$v=apply(dta,1,function(x) normVariance(d, x[1], x[2], x[3]))
dta[dta$v==min(dta$v),]