Re: Tes résultats ne sont pas reproductibles

Version plus optimiste que la précédente et sans doute plus juste :

# 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)
vOrig = Inf
aApprox = 0
mApprox = 0
sdApprox = 0
for (a in c(5000:20000)) {
 for (sd in sdTab) {
 for(m in mTab) {
 v = normVariance(italy, m, sd, a)
 if (v < vOrig) {
 sdApprox = sd
 mApprox = m
 aApprox = a
 vOrig = v
 }
 }
 }
}

Les résultats diffèrent un peu de ce que je trouve à la main ... J'analyse et je fais un retours.