Bonjour,
j'ai fait un fit avec la fonction nls
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4 setwd("C:/Users/pierre/Desktop/R") a1<-read.table('toshiba.txt',col.names=c('L','DL0')) attach(a1) result_fit <- nls(DL0 ~ Deq*(1-alpha*exp(-4*L/Leq)), data = a1, start = list(Deq = 100, alpha = 0.5, Leq = 40))
ça marche,
mais j'aimerais récupérer les incertitudes sur les valeurs Deq, alpha et Leq qu'il a trouvées.
Vous savez comment faire ?
En gros, j'aimerais comparer la différence entre deux fits et voir si cette différence est significative.
Partager