Ich versuche, ein Gammamischungsmodell (zwei Gammaverteilungen) an ein Alter-Fruchtbarkeits-Profil anzupassen. Ich habe einen Datensatz, der die altersspezifischen Fertilitätsraten und das Alter enthält, und ich möchte zwei Gammas anpassen, um die entsprechenden Parameter zu finden (am Ende verwende ich Fruchtbarkeitsprofile aus verschiedenen Jahren und versuche zu sehen, wie sich die Parameter über die Zeit entwickeln). Ich habe bisher versucht, Mixtools Library (GammamixEM) zu verwenden, aber ohne Erfolg. Für Hilfe wäre ich sehr dankbar.

Ale

a<- structure(list(EDAD = structure(1:45, .Label = c("11", "12", "13", "14", "15",   
"16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", 
"30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", 
"44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "Total" ),
class = "factor"), value = c(0, 0, 0, 0, 0.002761668, 0.006712018, 0.010820244, 
0.017867778, 0.029533765, 0.034055242, 0.036665669, 0.043840421, 0.042949584, 
0.042344612, 0.050588917, 0.050187588, 0.054114728, 0.057258792, 0.059280324, 
0.062566731, 0.062369629, 0.062154767, 0.063734337, 0.058236776, 0.052623842, 
0.046330921, 0.040639027, 0.033707865, 0.02531141, 0.017651534, 0.010953808, 
0.007463863, 0.003224766, 0.002190101, 0.001117443, 0.000465116, 0.000363901, 
0.00012647, 0.000267326, 0.000280308, 0, 0, 0, 0, 0)), .Names = c("EDAD", "value"), 
class = "data.frame", row.names = 79596:79640)